Thermal Stability of Globins: Implications of Flexibility and Heme Coordination Studied by Molecular Dynamics Simulations

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

Interactions between carbon nanotubes and external structures of SARS-CoV-2 using molecular docking and molecular dynamics

Molecular modeling techniques are used to describe the process of interaction between nanotubes and the main structures of the Covid-19 virus: the envelope protein, the main protease, and the Spike glycoprotein. Molecular docking studies show that the ligands have interaction characteristics capable of adsorbing the structures. Molecular dynamics simulations provide information on the mean squared deviation of atomic positions ​​between 0.5 and 3.0 Å. The Gibbs free energy model and solvent accessible surface area approaches are used. Through the results obtained through molecular dynamics simulations, it is noted that the zig-zag nanotube prefers to interact with E-pro, M-pro, and S-gly, respectively. Molecular couplings and free energy showed that the S-gly active site residues strongly interact with zigzag, chiral, and armchair nanotubes, in this order. The interactions demonstrated in this manuscript may predict some promising candidates for virus antagonists, which may be confirmed through experimental approaches.

Architectural digest: Thermodynamic stability and domain structure of a consensus monomeric globin

Artificial proteins representing the consensus of a set of homologous sequences have attracted attention for their increased thermodynamic stability and conserved activity. Here, we applied the consensus approach to a b-type heme-binding protein to inspect the contribution of a dissociable cofactor to enhanced stability and the chemical consequences of creating a generic heme environment. We targeted the group 1 truncated hemoglobin (TrHb1) subfamily of proteins for their small size (∼120 residues) and ease of characterization. The primary structure, derived from a curated set of ∼300 representative sequences, yielded a highly soluble consensus globin (cGlbN) enriched in acidic residues. Optical and NMR spectroscopies revealed high-affinity heme binding in the expected site and in two orientations. At neutral pH, proximal and distal iron coordination was achieved with a pair of histidine residues, as observed in some natural TrHb1s, and with labile ligation on the distal side. As opposed to studied TrHb1s, which undergo additional folding upon heme binding, cGlbN displayed the same extent of secondary structure whether the heme was associated with the protein or not. Denaturation required guanidine hydrochloride and showed that apo- and holoprotein unfolded in two transitions-the first (occurring with a midpoint of ∼2 M) was shifted to higher denaturant concentration in the holoprotein (∼3.7 M) and reflected stabilization due to heme binding, while the second transition (∼6.2 M) was common to both forms. Thus, the consensus sequence stabilized the protein but exposed the existence of two separately cooperative subdomains within the globin architecture, masked as one single domain in TrHb1s with typical stabilities. The results suggested ways in which specific chemical or thermodynamic features may be controlled in artificial heme proteins.

Computational Exploration of Anti-cancer Potential of Flavonoids against Cyclin-Dependent Kinase 8: An In Silico Molecular Docking and Dynamic Approach

Over the centuries, cancer has been considered one of the significant health threats. It holds the position in the list of deadliest diseases over the globe. In women, breast cancer is the most common among many cancers and is the second most common cancer all over the world, while lung cancer is the first. Cyclin-dependent kinase 8 (CDK8) has been identified as a critical oncogenic driver that is found in breast cancer and associated with tumor progression. Flavonoids were virtually screened against CDK8 using molecular docking, drug-likeness, ADMET prediction, and a molecular dynamics (MD) simulation approach to determine the potential flavonoid structure against CDK8. The results indicated that ZINC000005854718 showed the highest negative binding affinity of -10.7 kcal/mol with the targeted protein and passed all the drug-likeness parameters. Performed molecular dynamics simulation showed that docked complex systems have good conformational stability over 100 ns in different temperatures (298, 300, 305, 310, and 320 K). The comparison between calculated binding free energy via MM/PB(GB)SA methods and binding affinity calculated via molecular docking suggested tight binding of ZINC000005854718 with targeted protein. The results concluded that ZINC000005854718 has drug-like properties with tight and stable binding with the targeted protein.

Oxidative Implications of Substituting a Conserved Cysteine Residue in Sugar Beet Phytoglobin BvPgb 1.2

Phytoglobins (Pgbs) are plant-originating heme proteins of the globin superfamily with varying degrees of hexacoordination. Pgbs have a conserved cysteine residue, the role of which is poorly understood. In this paper, we investigated the functional and structural role of cysteine in BvPgb1.2, a Class 1 Pgb from sugar beet (Beta vulgaris), by constructing an alanine-substituted mutant (Cys86Ala). The substitution had little impact on structure, dimerization, and heme loss as determined by X-ray crystallography, size-exclusion chromatography, and an apomyoglobin-based heme-loss assay, respectively. The substitution significantly affected other important biochemical properties. The autoxidation rate increased 16.7- and 14.4-fold for the mutant versus the native protein at 25 °C and 37 °C, respectively. Thermal stability similarly increased for the mutant by ~2.5 °C as measured by nano-differential scanning fluorimetry. Monitoring peroxidase activity over 7 days showed a 60% activity decrease in the native protein, from 33.7 to 20.2 U/mg protein. When comparing the two proteins, the mutant displayed a remarkable enzymatic stability as activity remained relatively constant throughout, albeit at a lower level, ~12 U/mg protein. This suggests that cysteine plays an important role in BvPgb1.2 function and stability, despite having seemingly little effect on its tertiary and quaternary structure.

Effects of the Temperature and Salt Concentration on the Structural Characteristics of the Protein (PDB Code 1BBL)

The effect of the temperature and salt solution on the structural characteristics of the protein 1BBL was investigated by molecular dynamics simulations. The paper presents simulation results regarding the non-bonded energy and the structural stability of the protein immersed in salt solutions with different concentrations and temperatures. Our work demonstrates that the electrostatic potential energy and van der Waals energy of the system show the opposite changes with the influence of the external environment. Since the electrostatic potential energy changes more obviously, it is dominated in the non-bonding interactions. The structural parameters, such as the root mean square deviation and the radius of gyration, increased initially and decreased afterward with the increase of the salt concentration. The protein presented the loose structure with a relative low stability when it was immersed in a monovalent solution with a salt concentration of 0.8 mol/L. The salt concentration corresponding to the maximum value of structural parameters in the monovalent salt solution was double that in the divalent salt solution. It was also concluded that the protein presented a compact and stable structure when immersed in salt solutions with a high concentration of 2.3 mol/L. The analysis of the root mean square deviation and root mean square fluctuation of the protein sample also exhibited that the structural stability and chain flexibility are strongly guided by the effect of the temperature. These conclusions help us to understand the structural characteristics of the protein immersed in the salt solutions with different concentrations and temperatures.

An impedimetric approach for COVID-19 detection

In this study, an electrochemical approach for the determination of coronavirus disease (COVID-19) was developed. The biosensor system relied on the spike protein (S-protein) based infection mechanism of the virus and included separate interactions of receptors like angiotensin-converting enzyme 2 (ACE2) and CD147. After the optimization of experimental parameters, the analytical characteristics of both receptors ACE2 and CD147 were investigated. For ACE2 receptor, the linear detection ranges of the S-protein were found in the range of 700 ng mL^-1 to 1500 ng mL^-1 and from 1500 ng mL^-1 to 7000 ng mL^-1 with a limit of detection (LOD) value of 299.30 ng mL^-1. Meanwhile, for CD147 receptor the linear range was in the range of 500 ng mL^-1 to 5000 ng mL^-1 with a LOD value of 38.99 ng mL^-1. After the examination of analytical characteristics, the developed electrochemical approach was applied for severe acute respiratory syndrome coronavirus 2 samples and the obtained results were validated with real time polymerase chain reaction method.

Temperature effect on the SARS-CoV-2: A molecular dynamics study of the spike homotrimeric glycoprotein

Rapid spread of SARS-CoV-2 virus have boosted the need of knowledge about inactivation mechanisms to minimize the impact of COVID-19 pandemic. Recent studies have shown that SARS-CoV-2 virus can be disabled by heating, the exposure time for total inactivation depending on the reached temperature (e.g. more than 45 min at 329 K or less than 5 min at 373 K. In spite of recent crystallographic structures, little is known about the molecular changes induced by the temperature. Here, we unravel the molecular basis of the effect of the temperature over the SARS-CoV-2 spike glycoprotein, which is a homotrimer with three identical monomers, by executing atomistic molecular dynamics (MD) simulations at 298, 310, 324, 338, 358 and 373 K. Furthermore, both the closed down and open up conformational states, which affect the accessibility of receptor binding domain, have been considered. Our results suggest that the spike homotrimer undergoes drastic changes in the topology of the hydrogen bonding interactions and important changes on the secondary structure of the receptor binding domain (RBD), while electrostatic interactions (i.e. salt bridges) are mainly preserved. The proposed inactivation mechanism has important implications for engineering new approaches to fight the SARS-CoV-2 coronavirus, as for example, cleaving or reorganizing the hydrogen bonds through chaotropic agents or nanoparticles with local surface resonant plasmon effect.

Investigation of the Effect of Temperature on the Structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations

Statistical and epidemiological data imply temperature sensitivity of the SARS-CoV-2 coronavirus. However, the molecular level understanding of the virus structure at different temperature is still not clear. Spike protein is the outermost structural protein of the SARS-CoV-2 virus which interacts with the Angiotensin Converting Enzyme 2 (ACE2), a human receptor, and enters the respiratory system. In this study, we performed an all atom molecular dynamics simulation to study the effect of temperature on the structure of the Spike protein. After 200 ns of simulation at different temperatures, we came across some interesting phenomena exhibited by the protein. We found that the solvent exposed domain of Spike protein, namely S1, is more mobile than the transmembrane domain, S2. Structural studies implied the presence of several charged residues on the surface of N-terminal Domain of S1 which are optimally oriented at 10-30°C. Bioinformatics analyses indicated that it is capable of binding to other human receptors and should not be disregarded. Additionally, we found that receptor binding motif (RBM), present on the receptor binding domain (RBD) of S1, begins to close around temperature of 40°C and attains a completely closed conformation at 50°C. We also found that the presence of glycan moieties did not influence the observed protein dynamics. Nevertheless, the closed conformation disables its ability to bind to ACE2, due to the burying of its receptor binding residues. Our results clearly show that there are active and inactive states of the protein at different temperatures. This would not only prove beneficial for understanding the fundamental nature of the virus, but would be also useful in the development of vaccines and therapeutics.

Non-Hydrolyzable Plastics - An Interdisciplinary Look at Plastic Bio-Oxidation

Enzymatic plastic conversion has emerged recently as a potential adjunct and alternative to conventional plastic waste management technology. Publicity over progress in the enzymatic degradation of polyesters largely neglects that the majority of commercial plastics, including polyethylene, polypropylene, polystyrene and polyvinyl chloride, are still not biodegradable. Details about the mechanisms used by enzymes and an understanding of macromolecular factors influencing these have proved to be vital in developing biodegradation methods for polyesters. To expand the application of enzymatic degradation to other more recalcitrant plastics, extensive knowledge gaps need to be addressed. By drawing on interdisciplinary knowledge, we suggest that physicochemical influences also have a crucial impact on reactions in less well-studied types of plastic, and these need to be investigated in detail.

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

Investigation of the Effect of Temperature on the Structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations

Statistical and epidemiological data imply temperature sensitivity of the SARS-CoV-2 coronavirus. However, the molecular level understanding of the virus structure at different temperature is still not clear. Spike protein is the outermost structural protein of the SARS-CoV-2 virus which interacts with the Angiotensin Converting Enzyme 2 (ACE2), a human receptor, and enters the respiratory system. In this study, we performed an all atom molecular dynamics simulation to study the effect of temperature on the structure of the Spike protein. After 200 ns of simulation at different temperatures, we came across some interesting phenomena exhibited by the protein. We found that the solvent exposed domain of Spike protein, namely S1, is more mobile than the transmembrane domain, S2. Structural studies implied the presence of several charged residues on the surface of N-terminal Domain of S1 which are optimally oriented at 10-30°C. Bioinformatics analyses indicated that it is capable of binding to other human receptors and should not be disregarded. Additionally, we found that receptor binding motif (RBM), present on the receptor binding domain (RBD) of S1, begins to close around temperature of 40°C and attains a completely closed conformation at 50°C. We also found that the presence of glycan moieties did not influence the observed protein dynamics. Nevertheless, the closed conformation disables its ability to bind to ACE2, due to the burying of its receptor binding residues. Our results clearly show that there are active and inactive states of the protein at different temperatures. This would not only prove beneficial for understanding the fundamental nature of the virus, but would be also useful in the development of vaccines and therapeutics.

Structural and biochemical insights into an engineered high-redox potential laccase overproduced in Aspergillus

Fungal laccases have great potential as biocatalysts oxidizing a variety of aromatic compounds using oxygen as co-substrate. Here, the crystal structure of 7D5 laccase (PDB 6H5Y), developed in Saccharomyces cerevisiae and overproduced in Aspergillus oryzae, is compared with that of the wild type produced by basidiomycete PM1 (Coriolopsis sp.), PDB 5ANH. SAXS showed both enzymes form monomers in solution, 7D5 laccase with a more oblate geometric structure due to heavier and more heterogeneous glycosylation. The enzyme presents superior catalytic constants towards all tested substrates, with no significant change in optimal pH or redox potential. It shows noticeable high catalytic efficiency with ABTS and dimethyl-4-phenylenediamine, 7 and 32 times better than the wild type, respectively. Computational simulations demonstrated a more favorable binding and electron transfer from the substrate to the T1 copper due to the introduced mutations. PM1 laccase is exceptionally stable to thermal inactivation (t_1/2 70 °C = 1.2 h). Yet, both enzymes display outstanding structural robustness at high temperature. They keep folded during 2 h at 100 °C though, thereafter, 7D5 laccase unfolds faster. Rigidification of certain loops due to the mutations added on the protein surface would diminish the capability to absorb temperature fluctuations leading to earlier protein unfolding.