Probing immobilization mechanism of alpha-chymotrypsin onto carbon nanotube in organic media by molecular dynamics simulation

Molecular Insight into the Effect of HIV-TAT Protein on Amyloid-β Peptides

Increased deposition of amyloid-β (Aβ) plaques in the brain is a frequent pathological feature observed in human immunodeficiency virus (HIV)-positive patients. Emerging evidence indicates that HIV regulatory proteins, particularly the transactivator of transcription (TAT) protein, could interact with Aβ peptide, accelerating the formation of Aβ plaques in the brain and potentially contributing to the onset of Alzheimer's disease in individuals with HIV infection. Nevertheless, the molecular mechanisms underlying these processes remain unclear. In the present study, we have used long all-atom molecular dynamics simulations to probe the direct interactions between the TAT protein and Aβ peptide at the molecular level. Sampling over 28.0 μs, our simulations show that TAT protein induces a shift in the Aβ monomer ensemble toward elongated conformations, exposing aggregation-prone regions on the surface and thereby inducing subsequent aggregation. TAT protein also appears to enhance the stability of preformed Aβ fibrils, while increasing the β-sheet content within these fibrils. Our atomistically detailed simulations qualitatively agree with previous in vitro and in vivo studies. Importantly, our simulations identify key interactions between Aβ and the TAT protein that drive the Aβ aggregation process and stabilize the preformed Aβ aggregates, which are particularly challenging to obtain through current experimental techniques.

The stability and dynamics of the Aβ40/Aβ42 interlaced mixed fibrils

The accumulation of fibrillar amyloid-β (Aβ) aggregates in the brain, predominantly comprising 40- and 42-residue amyloid-β (Aβ40 and Aβ42), is a major pathological hallmark of Alzheimer's disease (AD). Aβ40 and Aβ42 naturally coexist in the brain under normal physiological conditions, and their interplay is generally considered to be a critical factor in the progression of AD. In addition to forming homogeneous oligomers and fibrils, Aβ40 and Aβ42 are also reported to co-assemble into hetero-oligomers and interlaced mixed fibrils, as evidenced by solid-state nuclear magnetic resonance spectroscopy (NMR), high molecular weight mass spectrometry and cross-seeding experiments. However, the exact molecular mechanisms underlying these processes remain unclear. In this study, we have used a recently resolved structurally uniform 1:1 mixture of Aβ40/Aβ42 interlaced mixed fibril as a prototype to gain insights into the molecular-level interactions between Aβ40 and Aβ42. We employed fully atomistic molecular dynamics simulation and compared the results with a homogeneous U-shaped Aβ40 fibrillar model. Our simulations using two different force fields provide conclusive evidence that the Aβ40/Aβ42 interlaced mixed fibril is energetically more favorable than the homogeneous Aβ40 fibrillar model. Furthermore, we also show that the increase in stability observed in the mixed model stems primarily from the packing interfaces and the stacking interfaces between C-termini. Our simulation results provide valuable mechanistic insights that are not readily accessible in experiment and could have significant implications for both the pathogenesis of AD and the development of current therapeutic strategies.Communicated by Ramaswamy H. Sarma.

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.

Assessment of Enzyme Functionality at Metal-Organic Framework Interfaces Developed through Molecular Simulations

The catalytic efficiency and unrivaled selectivity with which enzymes convert substrates to products have been tapped for widespread chemical transformations within biomedical technology, biofuel production, gas sensing, and the upgrading of commodity chemicals, just to name a few. However, the feasibility of enzymes implementation is challenged by the lack of reusability and loss of native catalytic activity due to the irreversible biocatalyst denaturation at high temperatures and in the presence of industrial solvents. Enzyme immobilization, a prerequisite for enzyme reusability, offers controllable strategies for increased functional viability of the biocatalyst in a synthetic environment. Herein we used molecular dynamics (MD) simulations and probed the noncovalent interactions between model enzymes of technological interest, i.e., carbonic anhydrase (CA) and myeloperoxidase (MPO), with selected metal-organic frameworks (MOFs; MIL-160 and ZIF-8) of proven industrial implementation. We found that the CA and MPO can bind to MIL-160 at optimal binding energies of 201 and 501 kJ mol^-1, respectively, that are strongly influenced by the increased incidence of hydrogen bonding between enzymes and the frameworks. The free energy of binding of enzymes to ZIF-8, on the other hand, was found to be less strongly influenced by hydrogen bonding networks relative to the occurrence of hydrophobic-hydrophobic interactions that yielded 106 kJ mol^-1 for CA and 201 kJ mol^-1 for MPO.

Nanostructured materials as a host matrix to develop robust peroxidases-based nanobiocatalytic systems

Nanostructured materials constitute an interesting and novel class of support matrices for the immobilization of peroxidase enzymes. Owing to the high surface area, robust mechanical stability, outstanding optical, thermal, and electrical properties, nanomaterials have been rightly perceived as immobilization matrices for enzyme immobilization with applications in diverse areas such as nano-biocatalysis, biosensing, drug delivery, antimicrobial activities, solar cells, and environmental protection. Many nano-scale materials have been employed as support matrices for the immobilization of different classes of enzymes. Nanobiocatalysts, enzymes immobilized on nano-size materials, are more stable, catalytically robust, and could be reused and recycled in multiple reaction cycles. In this review, we illustrate the unique structural/functional features and potentialities of nanomaterials-immobilized peroxidase enzymes in different biotechnological applications. After a comprehensive introduction to the immobilized enzymes and nanocarriers, the first section reviewed carbonaceous nanomaterials (carbon nanotube, graphene, and its derivatives) as a host matrix to constitute robust peroxidases-based nanobiocatalytic systems. The second half covers metallic nanomaterials (metals, and metal oxides) and some other novel materials as host carriers for peroxidases immobilization. The next section vetted the potential biotechnological applications of the resulted nanomaterials-immobilized robust peroxidases-based nanobiocatalytic systems. Concluding remarks, trends, and future recommendations for nanomaterial immobilized enzymes are also given.

Structure-function-dynamics of α-chymotrypsin based conjugates as a function of polymer charge

The field of protein-polymer conjugates has suffered from a lack of predictive tools and design guidelines to synthesize highly active and stable conjugates. In order to develop this type of information, structure-function-dynamics relationships must be understood. These relationships depend strongly on protein-polymer interactions and how these influence protein dynamics and conformations. Probing nanoscale interactions is experimentally difficult, but computational tools, such as molecular dynamics simulations, can easily obtain atomic resolution. Atomistic molecular dynamics simulations were used to study α-chymotrypsin (CT) densely conjugated with either zwitterionic, positively charged, or negatively charged polymers. Charged polymers interacted with the protein surface to varying degrees and in different regions of the polymer, depending on their flexibilities. Specific interactions of the negatively charged polymer with CT caused structural deformations in CT's substrate binding pocket and active site while no deformations were observed for zwitterionic and positively charged polymers. Attachment of polymers displaced water molecules from CT's surface into the polymer phase and polymer hydration correlated with the Hofmeister series.

Mechanics and dynamics of lysozyme immobilisation inside nanotubes

Lysozyme is an enzyme often used as an antibacterial agent in food industries and biochemical and pharmaceutical laboratories. Immobilisation of lysozyme by encapsulating in a nanotube has received much interest as it can enhance stability of the enzyme in ambient condition. Experimentally, various types of nanotubes have been proposed as a host for lysozyme. Here, we mathematically model the immobilisation process and the interaction between lysozyme and various types of nanotubes in order to compare the effectiveness of different nanotube materials. In this paper, we consider boron nitride, carbon, silicon, silicon carbide and titania nanotubes. For each type of nanotubes, we determine the critical radius that will maximise the interaction between the lysozyme molecule and the nanotube. Our results suggest that titania nanotube stands out as the most promising candidate for lysozyme storage and delivery. The model presented here can be extended to further investigate the interaction between different types of nanotube materials and protein structures for the development of effective molecular storage.

T7 RNA Polymerase Discriminates Correct and Incorrect Nucleoside Triphosphates by Free Energy

RNA polymerase (RNAP) is the primary machine responsible for transcription. Its ability to distinguish between correct (cognate) and incorrect (noncognate) nucleoside triphosphates (NTPs) is important for fidelity control in transcription. In this work, we investigated the substrate selection mechanism of T7 RNAP from the perspective of energetics. The dissociation free energies were determined for matched and unmatched base pairs in the preinsertion complex using the umbrella sampling method. A clear hydrogen-bond-rupture peak is observed in the potential of mean force curve for a matched base pair, whereas no such peaks are present in the position of mean force profiles for unmatched ones. The free-energy barrier could prevent correct substrates from being separated from the active site. Therefore, when NTPs diffuse into the active site, correct ones will stay for chemistry once they establish effective base pairing contacts with the template nucleotide, whereas incorrect ones will be withdrawn from the active site and rejected back to solution. This result provides an important energy evidence for the substrate selection mechanism of RNAP. Then we elucidated energetics and molecular details for correct NTP binding to the active site of the insertion complex. Our observations reveal that strong interactions act on the triphosphate of NTP to constrain its movement, whereas relatively weak interactions serve to position the base in the correct conformation. Triple interactions, hydrophobic contacts from residues M635 and Y639, base stacking from the 3' RNA terminal nucleotide, and base pairing from the template nucleotide act together to position the NTP base in a catalytically competent conformation. At last, we observed that incorrect NTPs cannot be as well-stabilized as the correct one in the active site when they are misincorporated in the insertion site. It is expected that our work can be helpful for comprehensively understanding details of this basic step in genetic transcription.

Interactions between Endohedral Metallofullerenes and Proteins: The Gd@C60-Lysozyme Model

Endohedral metallofullerenes (EMFs) have great potential as radioisotope carriers for nuclear medicine and as contrast agents for X-ray and magnetic resonance imaging. EMFs have still important restrictions for their use due to low solubility in physiological environments, low biocompatibility, nonspecific cellular uptake, and a strong dependence of their peculiar properties on physiological parameters, such as pH and salt content. Conjugation of the EMFs with proteins can overcome many of these limitations. Here we investigated the thermodynamics of binding of a model EMF (Gd@C60) with a protein (lysozyme) that is known to act as a host for the empty fullerene. As a rule, even if the shape of an EMF is exactly the same as that of the related fullerene, the interactions with a protein are significantly different. The estimated interaction energy (ΔG binding) between Gd@C60 and lysozyme is -18.7 kcal mol-1, suggesting the possibility of using proteins as supramolecular carriers for EMFs. π-π stacking, hydrophobic interactions, surfactant-like interactions, and electrostatic interactions govern the formation of the hybrid between Gd@C60 and lysozyme. The comparison of the energy contributions to the binding between C60 or Gd@C60 and lysozyme suggests that, although shape complementarity remains the driving force of the binding, the presence of electron transfer from the gadolinium atom to the carbon cage induces a charge distribution on the fullerene cage that strongly affects its interaction with the protein.

Phosphorylation promotes binding affinity of Rap-Raf complex by allosteric modulation of switch loop dynamics

The effects of phosphorylation of a serine residue on the structural and dynamic properties of Ras-like protein, Rap, and its interactions with effector protein Ras binding domain (RBD) of Raf kinase, in the presence of GTP, are investigated via molecular dynamics simulations. The simulations show that phosphorylation significantly effects the dynamics of functional loops of Rap which participate in the stability of the complex with effector proteins. The effects of phosphorylation on Rap are significant and detailed conformational analysis suggest that the Rap protein, when phosphorylated and with GTP ligand, samples different conformational space as compared to non-phosphorylated protein. In addition, phosphorylation of SER11 opens up a new cavity in the Rap protein which can be further explored for possible drug interactions. Residue network analysis shows that the phosphorylation of Rap results in a community spanning both Rap and RBD and strongly suggests transmission of allosteric effects of local alterations in Rap to distal regions of RBD, potentially affecting the downstream signalling. Binding free energy calculations suggest that phosphorylation of SER11 residue increases the binding between Rap and Raf corroborating the network analysis results. The increased binding of the Rap-Raf complex can have cascading effects along the signalling pathways where availability of Raf can influence the oncogenic effects of Ras proteins. These simulations underscore the importance of post translational modifications like phosphorylation on the functional dynamics in proteins and can be an alternative to drug-targeting, especially in notoriously undruggable oncoproteins belonging to Ras-like GTPase family.

Shining Light on Molecular Mechanism for Odor-selectivity of CNT-immobilized Olfactory Receptor

Olfactory receptor (OR)-based bioelectronic nose is a new type of bio-affinity sensor applied for detecting numerous odorant molecules. In order to elucidate the effect of the adsorption of nanomaterial carriers on the receptor structure and its selectivity to odors, we used a systematic computation-scheme to study two OR models immobilized onto carbon nanotube. Our result indicates that there is a multistep OR-adsorption process driven by hydrophobic interaction. Many allosteric communication pathways exist between the absorbed residues and the pocket ones, leading to a significant shrinkage of the pocket. Consequently, the size-selectivity of the receptor to the odors is changed to some extent. But, the odor size and its hydrophobicity, rather than specific functional groups of the odor, still play a determinant role in binding OR, at least for the 132 odors under study. Regardless of the limitation for the odor size in initial recognition, the different-size odors could induce significant changes in the pocket conformation so that it could better match the pocket space, indicating the importance of the ligand-fit binding. Due to the CNT-induced shrinkage of the pocket, the CNT immobilization could increase the binding affinity through enhancing van der Waals interaction, in particular for the large odors.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Mechanism of NTP Binding to the Active Site of T7 RNA Polymerase Revealed by Free-Energy Simulation

In genetic transcription, molecular dynamic details and energetics of NTP binding to the active site of RNA polymerase (RNAP) are poorly understood. In this article, we investigated the NTP binding process in T7 RNAP using all-atom MD simulation combined with the umbrella sampling technique. Based on our simulations, a two-step mechanism was proposed to explain NTP binding: first, substrate NTP in aqueous solution, which carries a magnesium ion, diffuses through a secondary channel of RNAP to attain a pore region, where it undergoes conformational changes to give a correct orientation; next, the NTP establishes initial basepairing contacts with the template nucleoside (TN). Our free-energy calculations suggest that both steps are spontaneous. This mechanism can easily explain the problem of NTP binding with different orientations. Moreover, it is found that the nascent NTP:TN basepair is fragile and easily broken by thermal disturbance. Therefore, we speculate that the fingers domain will be triggered to close, so as to create a steady environment for the next chemical step. The observations from the work provide valuable information for comprehensively understanding the mechanism of the basic step in genetic transcription.

Influence of Hyperglycemic Conditions on Self-Association of the Alzheimer's Amyloid β (Aβ1-42) Peptide

Clinical studies have identified a correlation between type-2 diabetes mellitus and cognitive decrements en route to the onset of Alzheimer's disease (AD). Recent studies have established that post-translational modifications of the amyloid β (Aβ) peptide occur under hyperglycemic conditions; particularly, the process of glycation exacerbates its neurotoxicity and accelerates AD progression. In view of the assertion that macromolecular crowding has an altering effect on protein self-assembly, it is crucial to characterize the effects of hyperglycemic conditions via crowding on Aβ self-assembly. Toward this purpose, fully atomistic molecular dynamics simulations were performed to study the effects of glucose crowding on Aβ dimerization, which is the smallest known neurotoxic species. The dimers formed in the glucose-crowded environment were found to have weaker associations as compared to that of those formed in water. Binding free energy calculations show that the reduced binding strength of the dimers can be mainly attributed to the overall weakening of the dispersion interactions correlated with substantial loss of interpeptide contacts in the hydrophobic patches of the Aβ units. Analysis to discern the differential solvation pattern in the glucose-crowded and pure water systems revealed that glucose molecules cluster around the protein, at a distance of 5-7 Å, which traps the water molecules in close association with the protein surface. This preferential exclusion of glucose molecules and resulting hydration of the Aβ peptides has a screening effect on the hydrophobic interactions, which in turn diminishes the binding strength of the resulting dimers. Our results imply that physical effects attributed to crowded hyperglycemic environments are incapable of solely promoting Aβ self-assembly, indicating that further mechanistic studies are required to provide insights into the self-assembly of post-translationally modified Aβ peptides, known to possess aggravated toxicity, under these conditions.

Electrostatics-mediated α-chymotrypsin inhibition by functionalized single-walled carbon nanotubes

Electrostatics-mediated α-chymotrypsin inhibition by functionalized single-walled carbon nanotubes.

Molecular mechanism of carbon nanotube to activate Subtilisin Carlsberg in polar and non-polar organic media

In the work, we mainly used molecular dynamics (MD) simulation and protein structure network (PSN) to study subtilisin Carlsberg (SC) immobilized onto carbon nanotube (CNT) in water, acetonitrile and heptane solvents, in order to explore activation mechanism of enzymes in non-aqueous media. The result indicates that the affinity of SC with CNT follows the decreasing order of water > acetonitrile > heptane. The overall structure of SC and the catalytic triad display strong robustness to the change of environments, responsible for the activity retaining. However, the distances between two β-strands of substrate-binding pocket are significantly expanded by the immobilization in the increasing order of water < acetonitrile < heptane, contributing to the highest substrate-binding energy in heptane media. PSN analysis further reveals that the immobilization enhances structural communication paths to the substrate-binding pocket, leading to its larger change than the free-enzymes. Interestingly, the increase in the number of the pathways upon immobilization is not dependent on the absorbed extent but the desorbed one, indicating significant role of shifting process of experimental operations in influencing the functional region. In addition, some conserved and important hot-residues in the paths are identified, providing molecular information for functional modification.

α-Chymotrypsin superactivity in quaternary ammonium salt solution: kinetic and computational studies

Ammonium salts determine an increase of the hydrophobicity of the α-chymotrypsin catalytic site and therefore an improvement of its activity.

Unraveling origins of the heterogeneous curvature dependence of polypeptide interactions with carbon nanostructures

Underlying causes of the differential polypeptide interactions on carbon nanosurfaces of varying curvatures emerge from a synchronized computational study.

Cloning, Expression, Mutagenesis Library Construction of Glycerol Dehydratase, and Binding Mode Simulation of Its Reactivase with Ligands

The production of 1, 3-propanediol (1, 3-PD) and 3-hydroxypropionaldehyde (3-HPA) by enzyme reaction has been a hot field, and glycerol dehydratase (GDHt) is the key and rate-limiting enzyme involved in their biosynthesis. The gldABC gene encoding GDHt was cloned from Klebsiella pneumoniae, and the activity of the corresponding proteins expressed extracellularly and intracellularly was 6.8 and 3.2 U/mg, respectively, about six and three times higher than that of the wild strain. The change of amino acids for the β subunit can adjust the length of the Co-N bond and affect the homolysis rate of the Co-C bond to change GDHt activity. The expression plasmid, pET-32a-gldAC (containing no gldB which encodes the β subunit of GDHt), was constructed to build the mutagenesis library to improve the GDHt activity. The binding models of glycerol dehydratase reactivation factor (GDHtR) with ATP, CTP, or GTP were simulated by semi-flexible docking, respectively, and there was almost no difference between them. This research provided the basis for studying the quantitative structure-activity relationships between GDHtR and its ligands, as well as searching inexpensive ligands to replace ATP. These results and methods are of great use in economical and highly efficient production of 3-HPA and 1, 3-PD by the enzyme method.