Electrostatic interactions of enzymes in non-aqueous conditions: insights from molecular dynamics simulations.

Electrostatic interactions of enzymes and their effects on enzyme activity and stability are poorly understood in non-aqueous conditions. Here, we investigate the contribution of the electrostatic interactions on the stability and activity of enzymes in the non-aqueous environment using molecular dynamics simulations. Lipase was selected as active and lysozyme as inactive model enzymes in non-aqueous media. Hexane was used as a common non-aqueous solvent model. In agreement with the previous experiments, simulations show that lysozyme has more structural instabilities than lipase in hexane. The number of hydrogen bonds and salt bridges of both enzymes is dramatically increased in hexane. In contrast to the other opinions, we show that the increase of the electrostatic interactions in non-aqueous media is not so favorable for enzymatic function and stability. In this condition, the newly formed hydrogen bonds and salt bridges can partially denature the local structure of the enzymes. For lysozyme, the changes in electrostatic interactions occur in all domains including the active site cleft, which leads to enzyme inactivation and destabilization. Interestingly, most of the changes in electrostatic interactions of lipase occur far from the active site regions. Therefore, the active site entrance regions remain functional in hexane. The results of this study reveal how the changes in electrostatic interactions can affect enzyme stability and activity in non-aqueous conditions. Moreover, we show for the first time how some enzymes, such as lipase, remain active in a non-aqueous environment.Communicated by Ramaswamy H. Sarma.

Supercritical CO2/Deep Eutectic Solvent Biphasic System as a New Green and Sustainable Solvent System for Different Applications: Insights from Molecular Dynamics Simulations.

Deep eutectic solvents (DESs) are one of the most interesting research subjects in green chemistry nowadays. Due to their low toxicity, simple synthesis, and lower prices, they have gradually taken the place of other green solvents such as ionic liquids (ILs) in sustainable processes. However, problems such as high viscosity and high polarity limit the applications of DESs in areas such as extraction, catalysis, and biocatalysis. In this work, we introduce and evaluate the potential application of scCO2/DES for the first time. Molecular dynamics simulations were used to examine the phase behavior, polarity, molecular mobilities, and microstructure of this system. Results show that CO2 molecules can significantly diffuse to the DES phase, while DES components do not appear in the scCO2 phase. The diffused CO2 molecules significantly enhanced the molecular mobility of the DES components. The presence of CO2 molecules changes the DES polarity so that hexane can be solubilized and dispersed in the DES phase. Radial distribution functions show that the solubilized CO2 molecules have negligible effects on the microstructure of DES. It was shown that chloride and urea are the main interaction sites of CO2 in DES. The results of this study show that scCO2/DES as a new class of green and versatile solvents can open a new promising window for research in sustainable chemistry and engineering.

Large interfacial relocation in RBD-ACE2 complex may explain fast-spreading property of Omicron.

The Omicron variant of SARS-CoV-2 emerged in South African in late 2021. This variant has a large number of mutations, and regarded as fastest-spreading Covid variant. The spike RBD region of SARS-CoV-2 and its interaction with human ACE2 play fundamental role in viral infection and transmission. To explore the reason of fast-spreading properties of Omicron variant, we have modeled the interactions of Omicron RBD and human ACE2 using docking and molecular dynamics simulations. Results show that RBD-ACE2 binding site may drastically relocate with an enlarged interface. The predicted interface has large negative binding energies and shows stable conformation in molecular dynamics simulations. It was found that the interfacial area in Omicron RBD-ACE2 complex is increased up to 40% in comparison to wild-type Sars-Cov-2. Moreover, the number of hydrogen bonds significantly increased up to 80%. The key interacting residues become also very different in Omicron variant. The new binding interface can significantly accommodate R403, as a key RBD residue, near ACE2 surface which leads to two new strong salt bridges. The exploration of the new binding interface can help to understand the reasons of high transmission rate of Omicron.

The protein-stabilizing effects of TMAO in aqueous and non-aqueous conditions.

We present a new water-dependent molecular mechanism for the widely-used protein stabilizing osmolyte, trimethylamine N-oxide (TMAO), whose mode of action has remained controversial. Classical interpretations, such as osmolyte exclusion from the vicinity of protein, cannot adequately explain the behavior of this osmolyte and were challenged by recent data showing the direct interactions of TMAO with proteins, mainly via hydrophobic binding. Solvent effect theories also fail to propose a straightforward mechanism. To explore the role of water and the hydrophobic association, we disabled osmolyte-protein hydrophobic interactions by replacing water with hexane and using lipase enzyme as an anhydrous-stable protein. Biocatalysis experiments showed that under this non-aqueous condition, TMAO does not act as a stabilizer, but strongly deactivates the enzyme. Molecular dynamics (MD) simulations reveal that TMAO accumulates near the enzyme and makes many hydrogen bonds with it, like denaturing osmolytes. Some TMAO molecules even reach the active site and interact strongly with the catalystic traid. In aqueous solvent, the enzyme functions well: the extent of TMAO interactions is reduced and can be divided into both polar and non-polar terms. Structural analysis shows that in water, some TMAO molecules bind to the enzyme surface like a surfactant. We show that these interactions limit water-protein hydrogen bonds and unfavorable water-hydrophobic surface contacts. Moreover, a more hydrophobic environment is formed in the solvation layer, which reduces water dynamics and subsequently, rigidifies the backbone in aqueous solution. We show that osmolyte amphiphilicity and protein surface heterogeneity can address the weaknesses of exclusion and solvent effect theories about the TMAO mechanism.

The effects of mutation and modification on the structure and stability of human lysozyme: A molecular link between carbamylation and atherosclerosis.

Amino acid mutations in some proteins such as lysozyme lead to genetically disorder variants and adverse pathogenic consequences. Recently, amino acid modifications were known as a risk factor in many related diseases such as uremia and atherosclerosis, showing the importance of these surface-structure changes. Although the structural consequences of the hereditary proteins have been examined extensively, such effects for the protein modifications are known to a lesser extent. One drawback in the examination of protein modifications is hardness in experimental detection of modifications by techniques such as NMR and crystallography. Molecular modeling and simulation can help to understand such phenomena at the molecular levels. It is more rational that the effects of both mutation and modification can be compared in a single protein model. Here, molecular dynamics simulation is used to compare the effects of a disease-related carbamylation modification and an amyloidogenic mutation (D67H) in human lysozyme as a model protein. The results show that the carbamylation adversely effects on the tertiary structure, leading to the similar unfolding pathway to the hereditary amyloidogenic form. The carbamylation leads to the instability of the overall protein conformation, especially on the ß-domain, which is a characteristic of hereditary amyloidosis in human lysozymes. The aggregation behaviors of both modified and mutant lysozyme were examined by molecular docking calculations. The results showed that the partially unfolded lysozyme might form tight protein aggregates upon carbamylation similar to the amyloidogenic variant. Both single and all-residues carbamylations impose serious conformational changes to the tertiary structure of lysozyme. It was obtained that carbamylation of lysozyme strongly effects on the stability of N-terminal ß-sheet, which can produce a highly unstable conformation. The results of this study not only show the adverse structural consequences of a disease-associated post-translational modification, but it also may be very helpful to understand the molecular basis for many carbamylation-related diseases such as atherosclerosis in ESRD patients. The results show that non-native post-translational modifications may be as structurally important as hereditary mutations.

Interaction of camel Lactoferrin derived peptides with DNA: a molecular dynamics study.

BACKGROUND: Lactoferrampin (LFampin), Lactoferricin (LFcin), and LFchimera are three well-known antimicrobial peptides derived from Lactoferrin and proposed as alternatives for antibiotics. Although the intracellular activity of these peptides has been previously demonstrated, their mode of action is not yet fully understood. Here, we performed a molecular dynamics simulation study to understand the molecular interactions between camel Lactoferrin derived peptides, including CLFampin, CLFcin, and CLFchimera, and DNA as an important intracellular target. RESULTS: Our results indicate that all three peptides bind to DNA, albeit with different propensities, with CLFchimera showing the highest binding affinity. The secondary structures of the peptides, modeled on Lactoferrin, did not undergo significant changes during simulation, supporting their functional relevance. Main residues involved in the peptide-DNA interaction were identified based on binding free energy estimates calculated over 200 ns, which, as expected, confirmed strong electrostatic interactions between DNA phosphate groups and positively charged peptide side chains. Interaction between the different concentrations of CLFchimera and DNA revealed that after binding of four copies of CLFchimera to DNA, hydrogen bonds between the two strands of DNA start to break from one of the termini. CONCLUSIONS: Importantly, our results revealed that there is no DNA-sequence preference for peptide binding, in line with a broad antimicrobial activity. Moreover, the results showed that the strength of the interaction between DNA and CLFchimera is concentration dependent. The insight provided by these results can be used for the rational redesign of natural antimicrobial peptides targeting the bacterial DNA.

Molecular mechanism of enzyme tolerance against organic solvents: Insights from molecular dynamics simulation.

Biocatalysis in presence of organic solvents has numerous industrially attractive advantages in comparison to traditional aqueous solvents. In some cases, the presence of organic molecules such as methanol in the processes such as enzymatic production of biodiesel is inevitable. However, enzyme inactivation and/or instability in organic solvents limits such biotechnological processes. Although it was found that some enzymes are more and others are less tolerant against organic solvents, the structural basis of such differences is relatively unknown. In this work, using molecular dynamics simulations, we have investigated the structural behavior of enzymes with completely different structural architecture including lipase, laccase and lysozyme in the presence of methanol as polar and hexane as non-polar organic solvents. In agreement with the previous experimental observations, simulations showed that lipase is more tolerant against both polar and non-polar organic solvents. It is found that lipase has high stability in pure hexane even higher than that obtained in the aqueous solvent. In contrast, laccase shows better stability in the aqueous conditions. To obtain general mechanism of enzyme inactivation in the presence of methanol and hexane, we have treated lysozyme as model enzyme in the different percentages of these solvents in long MD simulations. It is found that lysozyme is completely denatured at high concentration- of methanol, but it remains native at low concentration of this solvent. Interestingly, the concentration-dependence structural behavior of enzyme was completely different in the presence of hexane. It was obtained that low concentrations of hexane may impose more instability on the enzyme conformation than higher percentages. Results also showed that presence of water is determining factor in the enzyme stability at high concentrations of hexane. Pure hexane may also lead to the surface denaturation of the enzymes. Both methanol and hexane denaturation mechanisms were initiated by diffusion of organic solvent in hydrophobic core. However, enzyme denaturation in hexane was continued by a collapse of hydrophobic core and entering hexane molecules to the core, but in methanol it was completed by decomposition of the secondary structures. In both cases it was found that beta structures are more prone to destabilize than helix structures. This may be a reason for obtained results about lower stability of laccase with ß-barrel architecture than lipase with multiple helixes at it surface. In total, by our extensive structural data, it was found that the forces which stabilize tertiary structure have pivotal role in enzyme tolerance against both polar and non-polar organic solvents.

NGR (Asn-Gly-Arg)-targeted delivery of coagulase to tumor vasculature arrests cancer cell growth.

Induction of selective thrombosis and infarction in tumor-feeding vessels represents an attractive strategy to combat cancer. Here we took advantage of the unique coagulation properties of staphylocoagulase and genetically engineered it to generate a new fusion protein with novel anti-cancer properties. This novel bi-functional protein consists of truncated coagulase (tCoa) and an NGR (GNGRAHA) motif that recognizes CD13 and αvß3 integrin receptors, targeting it to tumor endothelial cells. Herein, we report that tCoa coupled by its C-terminus to an NGR sequence retained its normal binding activity with prothrombin and avß3 integrins, as confirmed in silico and in vitro. Moreover, in vivo biodistribution studies demonstrated selective accumulation of FITC-labeled tCoa-NGR fusion proteins at the site of subcutaneously implanted PC3 tumor xenografts in nude mice. Notably, systemic administration of tCoa-NGR to mice bearing 4T1 mouse mammary xenografts or PC3 human prostate tumors resulted in a significant reduction in tumor growth. These anti-tumor effects were accompanied by massive thrombotic occlusion of small and large tumor vessels, tumor infarction and tumor cell death. From these findings, we propose tCoa-NGR mediated tumor infarction as a novel and promising anti-cancer strategy targeting both CD13 and integrin αvß3 positive tumor neovasculature.

The influence of two imidazolium-based ionic liquids on the structure and activity of glucose oxidase: Experimental and theoretical studies.

The search for ionic liquids (ILs) with biochemical and biomedical applications has recently gained great attention. IL containing solvents can change the structure, stability and function of proteins. The study of protein conformation in ILs is important to understand enzymatic activity. In this work, conformational stability and activity of the enzyme in two imidazolium-based ILs (1-butyl 3-methyl-imidozolium and 1-hexyl 3-methyl-imidozoliumbromides) were investigated. We treated glucose oxidase as dimer-active enzyme in different IL concentration and seen that GOx activity was inhibited in the presence of ILs. Our experimental data showed that inhibition of activity and reduction of enzyme tertiary structure are more for hexyl than butyl derivative. These experimental results are in agreement with foregoing observations. To find a possible mechanism, a series of molecular dynamics simulation of the enzyme were performed at different IL concentration. The structure parameters obtained from MD simulation showed that conformational changes at the active site and FAD-binding site support the hypothesis of enzyme inhibition at the presence of ILs. Root mean square deviation and fluctuation calculations indicated that the enzyme has stable conformation at higher IL concentration, in agreement with experimental observation. But hexyl derivative has a much stronger stabilization effect on the protein structure. In summary, the present study could improve our understanding of the molecular mechanism about the ionic liquid effects on the structure and activity of proteins.

How a multimeric macromolecule is affected by divalent salts? Experimental and simulation study.

Salts exist in any cell and living organism in contact with biological macromolecules. How these salts affect biomolecules such as enzyme is important from both basic sciences and practical technologies. It was observed that divalent salts can change structure and function of protein at higher concentrations. Here, we investigated the effect of divalent salt on the behavior of a multimeric enzyme. We treated glucose oxidase as dimer-active enzyme in different CaCl2 concentration and seen that the enzyme become inactive at high concentration of salt. These experimental results are in agreement with recently published researches. To find a possible mechanism, a series of molecular dynamics simulation of the enzyme were performed at different salt concentration. According to the MD simulation, the conformational changes at the active site and FAD-binding site support the hypothesis of enzyme inactivation at high CaCl2 concentration. MD simulations also showed that enzyme has an unstable conformation at higher salt concentration which is in agreement with our experimental data. Detailed structural properties of the enzyme have been analyzed under different conditions. To the best of our knowledge, this is the first study that bears detailed structural mechanism about the salt effects on multimeric macromolecules.

RGD delivery of truncated coagulase to tumor vasculature affords local thrombotic activity to induce infarction of tumors in mice.

Induction of thrombosis in tumor vasculature represents an appealing strategy for combating cancer. Herein, we combined unique intrinsic coagulation properties of staphylocoagulase with new acquired functional potentials introduced by genetic engineering, to generate a novel bi-functional fusion protein consisting of truncated coagulase (tCoa) bearing an RGD motif on its C-terminus for cancer therapy. We demonstrated that free coagulase failed to elicit any significant thrombotic activity. Conversely, RGD delivery of coagulase retained coagulase activity and afforded favorable interaction of fusion proteins with prothrombin and αvß3 endothelial cell receptors, as verified by in silico, in vitro, and in vivo experiments. Although free coagulase elicited robust coagulase activity in vitro, only targeted coagulase (tCoa-RGD) was capable of producing extensive thrombosis, and subsequent infarction and massive necrosis of CT26 mouse colon, 4T1 mouse mammary and SKOV3 human ovarian tumors in mice. Additionally, systemic injections of lower doses of tCoa-RGD produced striking tumor growth inhibition of CT26, 4T1 and SKOV3 solid tumors in animals. Altogether, the nontoxic nature, unique shortcut mechanism, minimal effective dose, wide therapeutic window, efficient induction of thrombosis, local effects and susceptibility of human blood to coagulase suggest tCoa-RGD fusion proteins as a novel and promising anticancer therapy for human trials.

Effects of Natural Osmolytes on the Protein Structure in Supercritical CO2: Molecular Level Evidence.

Protein instability in supercritical CO2 limits the application of this green solvent in enzyme-catalyzed reactions. CO2 molecules act as a protein denaturant at high pressure under supercritical conditions. Here, for the first time, we show that natural osmolytes could stabilize protein conformation in supercritical CO2. Molecular dynamics simulation is used to monitor the effects of adding different natural osmolytes on the conformation and dynamics of chymotrypsin inhibitor 2 (CI2) in supercritical CO2. Simulations showed that CI2 is denatured at 200 bar in supercritical CO2, which is in agreement with experimental observations. Interestingly, the protein conformation remains native after addition of ∼1 M amino acid- and sugar-based osmolyte models. These molecules stabilize protein through the formation of supramolecular self-assemblies resulting from macromolecule-osmolyte hydrogen bonds. Nevertheless, trimethylamine N-oxide, which is known as a potent osmolyte for protein stabilization in aqueous solutions, amplifies protein denaturation in supercritical CO2. On the basis of our structural analysis, we introduce a new mechanism for the osmolyte effect in supercritical CO2, an "inclusion mechanism". To the best of our knowledge, this is the first study that introduces the application of natural osmolytes in a supercritical fluid and describes mechanistic insights into osmolyte action in nonaqueous media.

How a protein can remain stable in a solvent with high content of urea: insights from molecular dynamics simulation of Candida antarctica lipase B in urea : choline chloride deep eutectic solvent.

Deep eutectic solvents (DESs) are utilized as green and inexpensive alternatives to classical ionic liquids. It has been known that some of DESs can be used as solvent in the enzymatic reactions to obtain very green chemical processes. DESs are quite poorly understood at the molecular level. Moreover, we do not know much about the enzyme microstructure in such systems. For example, how some hydrolase can remain active and stable in a deep eutectic solvent including 9 M of urea? In this study, the molecular dynamics of DESs as a liquid was simulated at the molecular level. Urea : choline chloride as a well-known eutectic mixture was chosen as a model DES. The behavior of the lipase as a biocatalyst was studied in this system. For comparison, the enzyme structure was also simulated in 8M urea. The thermal stability of the enzyme was also evaluated in DESs, water, and 8M urea. The enzyme showed very good conformational stability in the urea : choline chloride mixture with about 66% urea (9 M) even at high temperatures. The results are in good agreement with recent experimental observations. In contrast, complete enzyme denaturation occurred in 8M urea with only 12% urea in water. It was found that urea molecules denature the enzyme by interrupting the intra-chain hydrogen bonds in a "direct denaturation mechanism". However, in a urea : choline chloride deep eutectic solvent, as a result of hydrogen bonding with choline and chloride ions, urea molecules have a low diffusion coefficient and cannot reach the protein domains. Interestingly, urea, choline, and chloride ions form hydrogen bonds with the surface residues of the enzyme which, instead of lipase denaturation, leads to greater enzyme stability. To the best of our knowledge, this is the first study in which the microstructural properties of a macromolecule are examined in a deep eutectic solvent.