Molecular simulation study of feruloyl esterase adsorption on charged surfaces: effects of surface charge density and ionic strength

Simulation Study of Polyethylene Terephthalate Hydrolase Adsorption on Self-Assembled Monolayers

Polyethylene terephthalate (PET) hydrolase, discovered in Ideonella sakaiensis (IsPETase), is a promising agent for the biodegradation of PET under mild reaction conditions, yet the thermal stability is poor. The efficient immobilization and orientation of IsPETase on different solid substrates are essential for its application. In this work, the combined parallel tempering Monte Carlo simulation with the all-atom molecular dynamics simulation approach was adopted to reveal the adsorption mechanism, orientation, and conformational changes of IsPETase adsorbed on charged self-assembled monolayers (SAMs), including COOH-SAM and NH2-SAM with different surface charge densities (SCDs). The results show that the protein adsorption orientation was determined not only by attraction interactions but also by repulsion interactions. IsPETase is adsorbed on the COOH-SAM surface with an "end-on" orientation, which favors the exposure of the catalyzed triplet to the solution. In addition, the entrance to the catalytic active center is larger on the COOH-SAM surface with a low SCD. This work reveals the controlled orientation and conformational information on IsPETase on charged surfaces at the atomistic level. This study would certainly promote our understanding of the mechanism of IsPETase adsorption and provide theoretical support for the design of substrates for IsPETase immobilization.

Enzyme immobilization studied through molecular dynamic simulations

In recent years, simulations have been used to great advantage to understand the structural and dynamic aspects of distinct enzyme immobilization strategies, as experimental techniques have limitations in establishing their impact at the molecular level. In this review, we discuss how molecular dynamic simulations have been employed to characterize the surface phenomenon in the enzyme immobilization procedure, in an attempt to decipher its impact on the enzyme features, such as activity and stability. In particular, computational studies on the immobilization of enzymes using i) nanoparticles, ii) self-assembled monolayers, iii) graphene and carbon nanotubes, and iv) other surfaces are covered. Importantly, this thorough literature survey reveals that, while simulations have been primarily performed to rationalize the molecular aspects of the immobilization event, their use to predict adequate protocols that can control its impact on the enzyme properties is, up to date, mostly missing.

Molecular Insights of Cellobiose Dehydrogenase Adsorption on Self-Assembled Monolayers

Cellobiose dehydrogenase (CDH) is capable of direct electron transfer (DET) on electrodes and is a promising redox enzyme for bioelectrochemical applications. Its unique two-domain structure makes the function of CDH adsorbed on the surface of the electrode deeply affected by the external environment, such as ion species, strength, pH, and surface charge density. To date, however, the exact mechanism of how the external environment tailors the structure and dynamics of CDH adsorbed on the electrode surface still remains poorly understood. Here, multiscale simulations were performed to look for insight into the effect of Na⁺ and Ca²⁺ ions on the activation of CDH on oppositely charged self-assembled monolayer (NH₂-SAM and COOH-SAM) surfaces with different surface charge densities (SCDs). Both Na⁺ and Ca²⁺ can promote CDH conformation switch from the open state to the closed state, while the promotion effect of Ca²⁺ is stronger than that of Na⁺ at the same conditions. However, the high ionic strength (IS) of Ca²⁺ renders the cytochrome (CYT) domain of CDH away from the NH₂-SAM with low SCD. In contrast, whatever the IS, the NH₂-SAM surface with high SCD can not only enhance the CYT-surface interaction but also achieve a closed-state conformation due to a similar role of Ca²⁺. Overall, this study gains molecular-level insights into the role of ion species and surface charge in modulating the structure and conformation of CDH on the SAM surface, thereby tailoring its activity.

Developing G value as an indicator for assessing the molecular status of immobilized antibody

Antibody-functionalized nanoparticles (Ab-NPs) are widely used in bioassays due to their excellent affinity, specificity toward antigen, and ease of operation. However, the uncontrollable molecular status of antibodies on NPs severely limits their applications. This work aims at developing a simple method to evaluate the antigen-binding activity of Ab-NPs using two parameters, i.e., antibody adsorption amount and antigen-binding strength. Herein, we proposed a mathematical expression, G, to quantitively describe the amount and strength of Ab-NPs. G value could be used to assess the antigen-binding performance of NPs influenced by surface and solution factors. Seven types of polymers with different surface properties, including four positively and three negatively charged polymer brushes, were grown from silica NPs via surface-initiated atom transfer radical polymerization (SI-ATRP). A pair of antigen and antibody, human chorionic gonadotropin (hCG) and anti-hCG, were selected to screen the antibody immobilization property of polymer brushes. Among them, the G values of 2 polymer-NPs with opposite charges reached maximum, resulting in low detection limits for hCG, where pDMAEA-NP and pMMA-NP represent Poly[N,N-(dimethylamino)ethyl acrylate]-NP and poly(methyl methacrylate)-NP, respectively. The G value of Ab-NPs makes it feasible to estimate the molecular status of the adsorbed antibodies on surfaces, thus showing great potential for in vitro biosensing and bioseparation.

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.

Simulated revelation of the adsorption behaviours of acetylcholinesterase on charged self-assembled monolayers

An acetylcholinesterase (AChE)-based electrochemical biosensor, as a promising alternative to detect organophosphates (OPs) and carbamate pesticides, has gained considerable attention in recent years, due to the advantages of simplicity, rapidity, reliability and low cost. The bio-activity of AChE immobilized on the surface and the direct electron transfer (DET) rate between an enzyme and an electrode directly determined the analytical performances of the AChE-based biosensor, and experimental studies have shown that the charged surfaces have a strong impact on the detectability of the AChE-based biosensor. Therefore, it is very important to reveal the behaviour of AChE in bulk solution and on charged surfaces at the molecular level. In this work, the adsorption orientation and conformation of AChE from Torpedo californica (TcAChE) on oppositely charged self-assembled monolayers (SAMs), COOH-SAM and NH₂-SAM with different surface charge densities, were investigated by parallel tempering Monte Carlo (PTMC) and all-atom molecular dynamics simulations (AAMD). Simulation results show that TcAChE could spontaneously and stably adsorb on two oppositely charged surfaces by the synergy of an electric dipole and charged residue patch, and opposite orientations were observed. The active-site gorge of TcAChE is oriented toward the surface with the "end-on" orientation and the active sites are close to the surface when it is adsorbed on the positively charged surface and the tunnel cost for the substrate is lower than that on the negatively charged surface and in bulk solution, while for TcAChE adsorbed on the negatively charged surface, the active site of TcAChE is far away from the surface and the active-site gorge is oriented toward the solution with a "back-on" orientation. It suggests that the positively charged surface could provide a better microenvironment for the efficient bio-catalytic reaction and quick DET between TcAChE and the electrode surface. Moreover, the RMSD, RMSF, dipole moment, gyration radius, eccentricity and superimposed structures show that only a slight conformational change occurred on the relatively flexible structure of TcAChE during simulations, and the native conformation is well preserved after adsorption. This work helps us better comprehend the adsorption mechanism of TcAChE on charged surfaces and might provide some guidelines for the development of new TcAChE-based amperometric biosensors for the detection of organophosphorus pesticides.

Feruloyl esterases: Biocatalysts to overcome biomass recalcitrance and for the production of bioactive compounds

Ferulic acid and its hydroxycinnamate derivatives represent one of the most abundant forms of low molecular weight phenolic compounds in plant biomass. Feruloyl esterases are part of a microorganism's plant cell wall-degrading enzymatic arsenal responsible for cleaving insoluble wall-bound hydroxycinnamates and soluble cytosolic conjugates. Stimulated by industrial requirements, accelerating scientific discoveries and knowledge transfer, continuous improvement efforts have been made to identify, create and repurposed biocatalysts dedicated to plant biomass conversion and biosynthesis of high-added value molecules. Here we review the basic knowledge and recent advances in biotechnological characteristics and the gene content encoding for feruloyl esterases. Information about several enzymes is systematically organized according to their function, biochemical properties, substrate specificity, and biotechnological applications. This review contributes to further structural, functional, and biotechnological R&D both for obtaining hydroxycinnamates from agricultural by-products as well as for lignocellulose biomass treatments aiming for production of bioethanol and other derivatives of industrial interest.

Bilirubin Oxidase Adsorption onto Charged Self-Assembled Monolayers: Insights from Multiscale Simulations

The efficient immobilization and orientation of bilirubin oxidase (BOx) on different solid substrates are essential for its application in biotechnology. The T1 copper site within BOx is responsible for the electron transfer. In order to obtain quick direct electron transfer (DET), it is important to keep the distance between the T1 copper site and electrode surface small and to maintain the natural structure of BOx at the same time. In this work, the combined parallel tempering Monte Carlo simulation with the all-atom molecular dynamics simulation approach was adopted to reveal the adsorption mechanism, orientation, and conformational changes of BOx from Myrothecium verrucaria (MvBOx) adsorbed on charged self-assembled monolayers (SAMs), including COOH-SAM and NH₂-SAM with different surface charge densities (±0.05 and ±0.19 C·m^-2). The results show that MvBOx adsorbs on negatively charged surfaces with a "back-on" orientation, whereas on positively charged surfaces, MvBOx binds with a "lying-on" orientation. The locations of the T1 copper site are closer to negatively charged surfaces. Furthermore, for negatively charged surfaces, the T1 copper site prefers to orient closer to the surface with lower surface charge density. Therefore, the negatively charged surface with low surface charge density is more suitable for the DET of MvBOx on electrodes. Besides, the structural changes primarily take place on the relatively flexible turns, coils, and α-helix. The native structure of MvBOx is well preserved when it adsorbs on both charged surfaces. This work sheds light on the controlling orientation and conformational information on MvBOx on charged surfaces at the atomistic level. This understanding would certainly promote our understanding of the mechanism of MvBOx immobilization and provide theoretical support for BOx-based bioelectrode design.

Charged Surface Regulates the Molecular Interactions of Electrostatically Repulsive Peptides by Inducing Oriented Alignment

Regulation of molecular orientation of charged dipeptides and involved interactions by electrostatic repulsion from like-charged surfaces were studied using all-atom molecular dynamics simulations. It was found that a charged surface can induce oriented alignment of like-charged peptides, and the oriented alignment leads to enhanced electrostatic repulsion between the peptide molecules. The findings are consistent with previous experimental results about the inhibition of charged protein aggregation using like-charged ion-exchange resin. Furthermore, the simulations provided molecular insights into this process, and demonstrated the distinct regulation effect of like-charged surfaces on the molecular interactions between peptides that possess an electric dipole structure. Both the charged surface and the electric dipole structure of peptides were confirmed to be crucial for the regulation. The research is expected to facilitate the rational design of surfaces or devices to regulate the behavior of amphoteric molecules such as proteins for both in vivo and in vitro applications, which would contribute to the regulation of protein-protein interactions and its application in life science and biotechnology.

Understanding the curvature effect of silica nanoparticles on lysozyme adsorption orientation and conformation: a mesoscopic coarse-grained simulation study

In nanobiotechnology applications, curvature of nanoparticles has a significant effect on protein activities. In this work, lysozyme adsorption on different-sized silica nanoparticles (SNPs) was simulated at the microsecond timescale by using mesoscopic coarse-grained molecular dynamics simulations. It is found that, with the increase of nanoparticle size, which indicates a decrease of surface curvature, adsorbed lysozyme shows a narrower orientation distribution and a greater conformation change, as the electrostatic attraction dominates lysozyme adsorption, and this trend is more pronounced on larger SNPs. Interestingly, the effect induced by different SNP surface curvatures is not related to the direct contact area between lysozyme and SNPs, but to the interfacial hydration layer above the silica surface, since a smaller curvature can lead to a stronger interfacial hydration and make the distribution of interfacial water molecules more ordered. Besides, at higher ionic strength, lysozyme conformation is less affected by strongly negatively charged SNPs, especially for larger nanoparticles. This work might shed some light on how to prepare protein coronas with higher bioactivities in nanobiotechnology.

Hamiltonian replica exchange simulations of glucose oxidase adsorption on charged surfaces

Hamiltonian replica exchange Monte Carlo simulations efficiently identify the lowest-energy orientations of proteins on charged surfaces at variable ionic strengths.

Albumin (BSA) adsorption onto graphite stepped surfaces

Nanomaterials are good candidates for the design of novel components with biomedical applications. For example, nano-patterned substrates may be used to immobilize protein molecules in order to integrate them in biosensing units. Here, we perform long MD simulations (up to 200 ns) using an explicit solvent and physiological ion concentrations to characterize the adsorption of bovine serum albumin (BSA) onto a nano-patterned graphite substrate. We have studied the effect of the orientation and step size on the protein adsorption and final conformation. Our results show that the protein is stable, with small changes in the protein secondary structure that are confined to the contact area and reveal the influence of nano-structuring on the spontaneous adsorption, protein-surface binding energies, and protein mobility. Although van der Waals (vdW) interactions play a dominant role, our simulations reveal the important role played by the hydrophobic lipid-binding sites of the BSA molecule in the adsorption process. The complex structure of these sites, that incorporate residues with different hydrophobic character, and their flexibility are crucial to understand the influence of the ion concentration and protein orientation in the different steps of the adsorption process. Our study provides useful information for the molecular engineering of components that require the immobilization of biomolecules and the preservation of their biological activity.

Hydrolysis-controlled protein adsorption and antifouling behaviors of mixed charged self-assembled monolayer: A molecular simulation study

Understanding the mechanism of the antimicrobial and antifouling properties of mixed charged materials is of great significance. The interactions between human gamma fibrinogen (γFg) and mixed carboxylic methyl ether-terminated (COOCH3-) and trimethylamino-terminated (N(CH3)3(+)-) SAMs and the influence of hydrolysis were studied by molecular simulations. After hydrolysis, the mixed SAMs exhibit behaviors from antimicrobial to antifouling, since the COOCH3-thiols were translated into carboxylic acid (COO(-)-) terminated thiols, which carried a net charge of -1 e. Simulation results showed that the main differences between COOCH3-/N(CH3)3(+)-SAM and COO(-)-/N(CH3)3(+)-SAM are the charged property and the hydration layer above the surface. γFg could stably adsorb on the positively-charged COOCH3-/N(CH3)3(+)-SAM. The adsorption behavior is mainly induced by the strong electrostatic attraction. There is a single hydration layer bound to the surface, which is related to the N(CH3)3(+) groups. The van der Waals repulsion between γFg and the single hydration layer are not strong enough to compensate the strong electrostatic attraction. After hydrolysis, the positively-charged SAM was transferred to a neutral mixed charged surface, the electrostatic attraction between γFg and the surface disappears. Meanwhile, the SAM surface is covered by double hydration layers, which is induced by the N(CH3)3(+) and COO(-) groups; water molecules around COO(-) groups are obviously denser than that around N(CH3)3(+) groups. With the combined contribution from double hydration layers and the vanishment of electrostatic attraction, γFg is forced to desorb from the surface. After hydrolysis, the internal structure of mixed SAM appears more ordered due to the electrostatic interactions between charged groups on the top of SAMs.

STATEMENT OF SIGNIFICANCE

The antimicrobial and antifouling materials are of great importance in many biological applications. The strong hydration property of surfaces and the interactions between proteins and surfaces play a key role in resisting protein adsorption. The mixed SAMs, constructed from a 1:1 combination of COOCH3- and N(CH3)3(+)-terminated thiols, can induce protein adsorption mainly through the electrostatic interaction. When the COOCH3-terminated thiols were hydrolyzed to negatively charged COO(-)-terminated thiols, the mixed-charged SAMs switched from antimicrobial to antifouling. Due to the strong hydration property of the mixed charged SAMs, the adsorbed γFg moved away from the surface. Understanding the interactions between protein and mixed-charged SAMs in the atomistic level is important for the practical design and development of new antimicrobial and antifouling materials.

Molecular simulations of cytochrome c adsorption on positively charged surfaces: the influence of anion type and concentration

The influence of anion type and concentration on the adsorption of cytochrome c onto the positively charged NH₂-SAM surface.