Preferential Adsorption of Prominent Amino Acids in the Urease Enzyme of Sporosarcina pasteurii on Arid Soil Components: A Periodic DFT Study

The urease enzyme is commonly used in microbially induced carbonate precipitation (MICP) and enzyme-induced carbonate precipitation (EICP) to heal and strengthen soil. Improving our understanding of the adsorption of the urease enzyme with various soil surfaces can lead to advancements in the MICP and EICP engineering methods as well as other areas of soil science. In this work, we use density functional theory (DFT) to investigate the urease enzyme's binding ability with four common arid soil components: quartz, corundum, albite, and hematite. As the urease enzyme cannot directly be simulated with DFT due to its size, the amino acids comprising at least 5% of the urease enzyme were simulated instead. An adsorption model incorporating the Gibbs free energy was used to determine the existence of amino acid-mineral binding modes. It was found that the nine simulated amino acids bind preferentially to the different soil components. Alanine favors corundum, glycine and threonine favor hematite, and aspartic acid favors albite. It was found that, under the standard environmental conditions considered here, amino acid binding to quartz is unfavorable. In the polymeric form where the side chains would dominate the binding interactions, hematite favors aspartic acid through its R-OH group and corundum favors glutamic acid through its R-Ket group. Overall, our model predicts that the urease enzyme produced by Sporosarcina pasteurii can bind to various oxides found in arid soil through its alanine, glycine, aspartic/glutamic acid, or threonine residues.

A review of clay based photocatalysts: Role of phyllosilicate mineral in interfacial assembly, microstructure control and performance regulation

Over the past decades, inspired by the outstanding properties of clay minerals such as abundance, low-cost, environmental benignity, high stability, and regularly arranged silica-alumina framework, researchers put much efforts on the interface assembly and surface modification of natural minerals with bare photocatalysts, i.e. TiO₂, g-C₃N₄, ZnO, MoS₂, etc. The clay-based hybrid photocatalysts have resulted in a rich database for their tailor-designed microstructures, characterizations, and environmental-related applications. Therefore, in this study, we took a brief introduction of three representative minerals, i.e. kaolinite, montmorillonite and rectorite, and discussed their basic merits in photocatalysis applications. After that, we summarized the recent advances in construction of stable visible-light driven photocatalysts based on these minerals. The structure-activity relationships between the properties of clay types, pore structure, distribution/dispersion and light absorption, carrier separation efficiency as well as redox performance were illustrated in detail. Such representative information would provide theoretical basis and scientific support for the application of clay based photocatalysts. Finally, we pointed out the major challenges and future directions at the end of this review. Undoubtedly, control and preparation of novel photocatalysts based on clays will continue to witness many breakthroughs in the arena of solar-driven technologies.

Buffer Influence on the Amino Acid Silica Interaction

Protein-surface interactions are exploited in various processes in life sciences and biotechnology. Many of such processes are performed in presence of a buffer system, which is generally believed to have an influence on the protein-surface interaction but is rarely investigated systematically. Combining experimental and theoretical methodologies, we herein demonstrate the strong influence of the buffer type on protein-surface interactions. Using state of the art chromatographic experiments, we measure the interaction between individual amino acids and silica, as a reference to understand protein-surface interactions. Among all the 20 proteinogenic amino acids studied, we found that arginine (R) and lysine (K) bind most strongly to silica, a finding validated by free energy calculations. We further measured the binding of R and K at different pH in presence of two different buffers, MOPS (3-(N-morpholino)propanesulfonic acid) and TRIS (tris(hydroxymethyl)aminomethane), and find dramatically different behavior. In presence of TRIS, the binding affinity of R/K increases with pH, whereas we observe an opposite trend for MOPS. These results can be understood using a multiscale modelling framework combining molecular dynamics simulation and Langmuir adsorption model. The modelling approach helps to optimize buffer conditions in various fields like biosensors, drug delivery or bio separation engineering prior to the experiment.

Role of Hydrogen-Bonding and OH-π Interactions in the Adhesion of Epoxy Resin on Hydrophilic Surfaces

Epoxy resin adhesives are widely used for joining metal alloys in various industrial fields. To elucidate the adhesion mechanism microscopically, we investigated the interfacial interactions of epoxy resin with hydroxylated silica (0 0 1) and γ-alumina (0 0 1) surfaces using periodic density functional theory calculations as well as density of states (DOS) and crystal orbital Hamilton population (COHP) analyses. To better understand the interfacial interactions, we employed and analyzed water and benzene molecules as hydrophilic and hydrophobic adsorbates, respectively. Structural features and calculated adhesion energies reveal that these small adsorbates have a higher affinity for the γ-alumina surface than that for the silica surface, while a fragmentary model for the epoxy resin exhibits a strong interaction with the silica surface. This discrepancy suggests that the structural features of the hydroxylated silica surface dictate its affinity to a specific species. Partial DOS and COHP curves provide evidence for the presence of OH-π interactions between the OH groups on the surfaces and the benzene rings of the epoxy resin fragments. The orbital interaction energies of the H-bonding and OH-π interactions evaluated from the integrated COHP indicate that the OH-π interaction is a nonnegligible origin of the adhesion interaction, even when polymers with hydrophobic benzene rings are adsorbed on hydroxylated surfaces.

Silicate-mediated interstellar water formation: A theoretical study

Water is one of the most abundant molecules in the form of solid ice phase in the different regions of the interstellar medium (ISM). This large abundance cannot be properly explained by using only traditional low temperature gas-phase reactions. Thus, surface chemical reactions are believed to be major synthetic channels for the formation of interstellar water ice. Among the different proposals, hydrogenation of atomic O (i.e., 2H + O → H2O) is a chemically "simple" and plausible reaction toward water formation occurring on the surfaces of interstellar grains. Here, novel theoretical results concerning the formation of water adopting this mechanism on the crystalline (010) Mg2SiO4 surface (a unequivocally identified interstellar silicate) are presented. The investigated reaction aims to simulate the formation of the first water ice layer covering the silicate core of dust grains. Adsorption of the atomic O as a first step of the reaction has been computed, results indicating that a peroxo (O 2 2 - ) group is formed. The following steps involve the adsorption, diffusion and reaction of two successive H atoms with the adsorbed O atom. Results indicate that H diffusion on the surface has barriers of 4-6 kcal mol-1, while actual formation of OH and H2O present energy barriers of 22-23 kcal mol-1. Kinetic study results show that tunneling is crucial for the occurrence of the reactions and that formation of OH and H2O are the bottlenecks of the overall process. Several astrophysical implications derived from the theoretical results are provided as concluding remarks.

When the Surface Matters: Prebiotic Peptide-Bond Formation on the TiO2 (101) Anatase Surface through Periodic DFT-D2 Simulations

The mechanism of the peptide-bond formation between two glycine (Gly) molecules has been investigated by means of PBE-D2* and PBE0-D2* periodic simulations on the TiO₂ (101) anatase surface. This is a process of great relevance both in fundamental prebiotic chemistry, as the reaction univocally belongs to one of the different organizational events that ultimately led to the emergence of life on Earth, as well as from an industrial perspective, since formation of amides is a key reaction for pharmaceutical companies. The efficiency of the surface catalytic sites is demonstrated by comparing the reactions in the gas phase and on the surface. At variance with the uncatalyzed gas-phase reaction, which involves a concerted nucleophilic attack and dehydration step, on the surface these two steps occur along a stepwise mechanism. The presence of surface Lewis and Brönsted sites exerts some catalytic effect by lowering the free energy barrier for the peptide-bond formation by about 6 kcal mol^-1 compared to the gas-phase reaction. Moreover, the co-presence of molecules acting as proton-transfer assistants (i.e., H₂ O and Gly) provide a more significant kinetic energy barrier decrease. The reaction on the surface is also favorable from a thermodynamic standpoint, involving very large and negative reaction energies. This is due to the fact that the anatase surface also acts as a dehydration agent during the condensation reaction, since the outermost coordinatively unsaturated Ti atoms strongly anchor the released water molecules. Our theoretical results provide a comprehensive atomistic interpretation of the experimental results of Martra et al. (Angew. Chem. Int. Ed. 2014, 53, 4671), in which polyglycine formation was obtained by successive feedings of Gly vapor on TiO₂ surfaces in dry conditions and are, therefore, relevant in a prebiotic context envisaging dry and wet cycles occurring, at mineral surfaces, in a small pool.

The Adhesion Mechanism of Marine Mussel Foot Protein: Adsorption of L-Dopa onα- and β-Cristobalite Silica Using Density Functional Theory

Marine mussels strongly adhere to various surfaces and endure their attachment under a variety of conditions. In order to understand the basic mechanism involved, we study the adsorption of L-dopa molecule on hydrophilic geminal and terminal isolated silanols of silica (001) surface. High content of modified amino acid L-dopa is found in the glue-like material secreted by the mussels through which it sticks to various surfaces under water. To understand the adsorption behavior, we have made use of periodic Density Functional Theory (DFT) study. The L-dopa molecule adheres to silica surfaces terminated with geminal and terminal silanols via its catechol part. In both cases, the adhesion is achieved through the formation of 4 H-bonds. A binding energy of 29.48 and 31.67 kcal/mol has been estimated, after the inclusion of dispersion energy, for geminal and terminal silanols of silica, respectively. These results suggest a relatively stronger adhesion of dopa molecule for surface with terminal isolated silanols.

Does hydrohalic acid HX (X = F, Cl) form true N-protonated twisted amide salts? Effects of anions on the ion-pair interactions and on the amide moiety in N-protonated tricyclic twisted amide salts

The [BF₄]⁻ and [RSO₃]⁻ anions interact with N-protonated amide cations through N–H⋯F and N–H⋯O strong hydrogen bonds and hydrohalic acids form very weak N⋯H–X hydrogen bonds.

Investigating the adsorption mechanism of glycine in comparison with catechol on cristobalite surface using density functional theory for bio-adhesive materials

Amino acid proteins exist in Mussel's adhesive (mefp's) of which glycine has a significant amount. A density functional theory simulation study was performed in a belief that all the proteins in mefp's are responsible for the versatile adhesion.

Theoretical insights into structure, bonding, reactivity and importance of ion-pair interactions in Kirby's tetrafluoroboric acid salts of twisted amides

The hydrolysis of amide 1 is more exothermic and is more favorable than amides 2 and 3 with bridgehead methyl.

Computational investigation of hole mobilities in organic semiconductors: comparison of single crystal structures and surface adsorbed clusters

We report hole mobilities obtained computationally based on both single crystal geometries and those obtained from crystal fragments optimised on a model surface. Such computational estimates can differ considerably from experimentally measured thin film mobilities. One source of this discrepancy is due to a difference in the morphology of the thin film compared with that of the crystal. Here, predictions of thin film hole mobilities based on optimised structures are given. A model surface is used to provide an inert geometric platform for the formation of an organic monolayer. The model is tested on pentacene and TIPS-pentacene for which experimental information of the surface morphology exists. The model has also been applied to four previously uninvestigated structures. Two of the compounds studied had fairly low predicted mobilities in their single crystal structures, which were vastly improved post-optimisation. This is in accord with experiment.

Zwitterionic versus canonical amino acids over the various defects in zeolites: a two-layer ONIOM calculation

Defects are often considered as the active sites for chemical reactions. Here a variety of defects in zeolites are used to stabilize zwitterionic glycine that is not self-stable in gas phase; in addition, effects of acidic strengths and zeolite channels on zwitterionic stabilization are demonstrated. Glycine zwitterions can be stabilized by all these defects and energetically prefer to canonical structures over Al and Ga Lewis acidic sites rather than Ti Lewis acidic site, silanol and titanol hydroxyls. For titanol (Ti-OH), glycine interacts with framework Ti and hydroxyl sites competitively, and the former with Lewis acidity predominates. The transformations from canonical to zwitterionic glycine are obviously more facile over Al and Ga Lewis acidic sites than over Ti Lewis acidic site, titanol and silanol hydroxyls. Charge transfers that generally increase with adsorption energies are found to largely decide the zwitterionic stabilization effects. Zeolite channels play a significant role during the stabilization process. In absence of zeolite channels, canonical structures predominate for all defects; glycine zwitterions remain stable over Al and Ga Lewis acidic sites and only with synergy of H-bonding interactions can exist over Ti Lewis acidic site, while automatically transform to canonical structures over silanol and titanol hydroxyls.

Does Dispersion Dominate over H-Bonds in Drug-Surface Interactions? The Case of Silica-Based Materials As Excipients and Drug-Delivery Agents

Amorphous silica is widely employed in pharmaceutical formulations both as a tableting, anticaking agent and as a drug delivery system, whereas MCM-41 mesoporous silica has been recently proposed as an efficient support for the controlled release of drugs. Notwithstanding the relevance of this topic, the atomistic details about the specific interactions between the surfaces of the above materials and drugs and the energetic of adsorption are almost unknown. In this work, we resort to a computational ab initio approach, based on periodic Density Functional Theory (DFT), to study the adsorption behavior of two popular drugs (aspirin and ibuprofen) on two models of an amorphous silica surface characterized by different hydrophilic/hydrophobic properties due to different SiOH surface groups' density. Particular effort was devoted to understand the role of dispersive (vdW) interactions in the adsorption mechanism and their interplay with H-bond interactions. On the hydrophilic silica surface, the H-bond pattern of the Si-OH groups rearranges to comply with the formation of new H-bond interactions triggered by the adsorbed drug. The interaction energy of ibuprofen with the hydrophilic model of the silica surface is computed to be very close to the sublimation energy of the ibuprofen molecular crystal, accounting for the experimental evidence of ibuprofen crystal amorphization induced by the contact with the mesoporous silica material. For both surface models, dispersion interactions play a crucial role in dictating the features of the drug/silica system, and they become dominant for the hydrophobic surface. It was proved that a competition may exist between directional H-bonds and nonspecific dispersion driven interactions, with important structural and energetic consequences for the adsorption. The results of this work emphasize the inadequacy of plain DFT methods to model adsorption processes involving inorganic surfaces and drugs of moderate size, due to the missing term accounting for London dispersion interactions.

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.

Dispersion Interactions with Density-Functional Theory: Benchmarking Semiempirical and Interatomic Pairwise Corrected Density Functionals

We present a comparative assessment of the accuracy of two different approaches for evaluating dispersion interactions: interatomic pairwise corrections and semiempirical meta-generalized-gradient-approximation (meta-GGA)-based functionals. This is achieved by employing conventional (semi)local and (screened-)hybrid functionals, as well as semiempirical hybrid and nonhybrid meta-GGA functionals of the M06 family, with and without interatomic pairwise Tkatchenko-Scheffler corrections. All of those are tested against the benchmark S22 set of weakly bound systems, a representative larger molecular complex (dimer of NiPc molecules), and a representative dispersively bound solid (hexagonal boron nitride). For the S22 database, we also compare our results with those obtained from the pairwise correction of Grimme (DFT-D3) and nonlocal Langreth-Lundqvist functionals (vdW-DF1 and vdW-DF2). We find that the semiempirical kinetic-energy-density dependence introduced in the M06 functionals mimics some of the nonlocal correlation needed to describe dispersion. However, long-range contributions are still missing. Pair-wise interatomic corrections, applied to conventional semilocal or hybrid functionals, or to M06 functionals, provide for a satisfactory level of accuracy irrespectively of the underlying functional. Specifically, screened-hybrid functionals such as the Heyd-Scuseria-Ernzerhof (HSE) approach reduce self-interaction errors in systems possessing both localized and delocalized orbitals and can be applied to both finite and extended systems. Therefore, they serve as a useful underlying functional for dispersion corrections.

Do H-bond features of silica surfaces affect the H2O and NH3 adsorption? Insights from periodic B3LYP calculations

The adsorption of a single H(2)O and NH(3) molecule on different fully hydroxylated α-quartz, cristobalite, and tridymite surfaces has been studied at the B3LYP level of theory, within a periodic approach using basis sets of polarized triple-ζ quality and accounting for basis set superposition error (BSSE). Fully hydroxylated crystalline silica exhibits SiOH as terminal groups whose distribution and H-bond features depend on both the considered silica polymorph and the crystallographic plane, which gives rise to isolated, H-bond interacting SiOH pairs or infinitely connected H-bond chains. A key point of the present study is to understand how the H-bond features of a dry crystalline silica surface influence its adsorption properties. Results reveal that the silica-adsorbate (H(2)O and NH(3)) interaction energy anticorrelates with the density of SiOH groups at the surface. This counterintuitive observation arises from the fact that pre-existing H-bonds of the dry surface need to be broken to establish new H-bonds between the surface and the adsorbate, which manifests in a sizable energy cost due to surface deformation. A simple method is also proposed to estimate the strength of the pre-existing H-bonds at the dry surfaces, which is shown to anticorrelate with the adsorbate interaction energy, in agreement with the above trends.

Local MP2 with Density Fitting for Periodic Systems: A Parallel Implementation

A parallel implementation is presented for the evaluation of local second-order Møller-Plesset perturbation theory (LMP2) energies in periodic, nonconducting crystalline systems with a density-fitting approximation of two-electron repulsion integrals. Peculiarities of the periodic case with respect to parallel LMP2 implementations in molecular codes, such as the use of translational and point symmetry, impose different strategies in order to achieve good parallel performance. The implementation is benchmarked on a few systems, representing a choice of the most interesting solid state quantum-chemistry problems where the MP2 approach can be decisive. Good parallel efficiency of the algorithms is demonstrated for up to 54 processors. Test systems include a metal organic framework (MOF-5) 3D crystalline structure with a triple-ζ-quality basis set: this is the largest calculation performed so far with 106 atoms, 532 correlated electrons, and 2884 atomic orbitals per unit cell.