Temperature-induced unfolding of epidermal growth factor (EGF): insight from molecular dynamics simulation

Solar-driven strongly coupled plasmonic Au nanoarrays on mesoporous silica nanodisks enable selective fungal and bacterial inactivation in well water

Well water is an important water source in isolated rural areas but easily suffers from microbial contamination. Herein, we anchored periodic Au nanoarrays on mesoporous silica nanodisks (Au-MSN) to fabricate a solar-driven nano-stove for well water disinfection. The solar/Au-MSN process completely inactivated 3.98, 6.55, 7.11 log10 cfu/mL, and 3.37 log10 pfu/mL of Aspergillus niger spores, Escherichia coli, chlorine-resistant Spingopyxis sp. BM1-1, and bacteriophage MS2 within 5 min, respectively. Moreover, the complete inactivation of various microorganisms (even at a viable but nonculturable state) was achieved in the flow-through reactor under natural solar light in real well water matrixes. Thorough characterizations and theoretical simulations verified that the densely anchoring strategy of Au-MSN's nanoarray worked on broadband absorption via the photon confinement effect, and trace amounts of Au can induce strong electromagnetic fields and collective localized heating. The resulting surge of 1O2 and heat synergically destroyed membranes, dysfunction cellular self-defense and metabolic system, induced intracellular oxidative stress, and ultimately inactivated microorganisms. Additionally, the 1O2-dominated oxidation and cell adhesion facilitated the selective disinfection in real well water matrixes. This study provides a cost-effective and practical solution for efficient well water disinfection, which assists isolated rural areas in getting safe drinking water.

An efficient method to predict protein thermostability in alanine mutation

The relationship between protein sequence and its thermodynamic stability is a critical aspect of computational protein design. In this work, we present a new theoretical method to calculate the free energy change (ΔΔG) resulting from a single-point amino acid mutation to alanine in a protein sequence. The method is derived based on physical interactions and is very efficient in estimating the free energy changes caused by a series of alanine mutations from just a single molecular dynamics (MD) trajectory. Numerical calculations are carried out on a total of 547 alanine mutations in 19 diverse proteins whose experimental results are available. The comparison between the experimental ΔΔG_exp and the calculated values shows a generally good correlation with a correlation coefficient of 0.67. Both the advantages and limitations of this method are discussed. This method provides an efficient and valuable tool for protein design and engineering.

Refolding of bioactive human epidermal growth factor from E. coli BL21(DE3) inclusion bodies & evaluations on its in vitro & in vivo bioactivity

Human epidermal growth factor (hEGF) is a mitogenic protein widely used in pharmaceutical and cosmetic industries, thus recombinant DNA technology has been applied to meet the high demand for hEGF. The overexpression of recombinant protein in E. coli often leads to the formation of inclusion bodies (IBs). Mild solubilisation preserves the native secondary protein structure in IBs, thereby the high recovery of active protein from IBs. The redox system also plays a pivotal role in the formation of disulphide bonds during refolding of disulphide bond-containing protein. This study aimed to recover hEGF from bacterial IBs through freeze-thawing solubilisation and glutathione-based oxidative refolding. CBD-Ssp DnaB-hEGF fusion protein was expressed as IBs in E. coli, washed with Triton X-100 and urea to remove most protein contaminants, then the solubilised fusion protein was obtained by freeze-thawing with the addition of 2 M urea. The solubilised protein was subsequently refolded by intein cleavage via a glutathione-based redox system. The refolded hEGF demonstrated heat-resistant properties, interacted with specific antibodies on ELISA, stimulated keratinocyte proliferation and possessed significant in vivo wound healing properties on the 8th day, confirming that hEGF was correctly folded. In summary, the protocol described is suitable for the recovery of refolded hEGF from bacterial IBs by mild solubilisation and oxidative refolding.

Computational Investigation of Protein Photoinactivation by Molecular Hyperthermia

To precisely control protein activity in a living system is a challenging yet long-pursued objective in biomedical sciences. Recently, we have developed a new approach named molecular hyperthermia (MH) to photoinactivate protein activity of interest without genetic modification. MH utilizes nanosecond laser pulse to create nanoscale heating around plasmonic nanoparticles to inactivate adjacent protein in live cells. Here we use a numerical model to study important parameters and conditions for MH to efficiently inactivate proteins in nanoscale. To quantify the protein inactivation process, the impact zone is defined as the range where proteins are inactivated by the nanoparticle localized heating. Factors that reduce the MH impact zone include the laser pulse duration, temperature-dependent thermal conductivity (versus constant properties), and nonspherical nanoparticle geometry. In contrast, the impact zone is insensitive to temperature-dependent material density and specific heat, as well as thermal interface resistance based on reported data in the literature. The low thermal conductivity of cytoplasm increases the impact zone. Different proteins with various Arrhenius kinetic parameters have significantly different impact zones. This study provides guidelines to design the protein inactivation process by MH.

Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis

Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.

Non-Arrhenius Reaction-Diffusion Kinetics for Protein Inactivation over a Large Temperature Range

Understanding protein folding and unfolding has been a long-standing fundamental question and has important applications in manipulating protein activity in biological systems. Experimental investigations of protein unfolding have been predominately conducted by small temperature perturbations (e.g., temperature jump), while molecular simulations are limited to small time scales (microseconds) and high temperatures to observe unfolding. Thus, it remains unclear how fast a protein unfolds irreversibly and loses function (i.e., inactivation) across a large temperature range. In this work, using nanosecond pulsed heating of individual plasmonic nanoparticles to create precise localized heating, we examine the protein inactivation kinetics at extremely high temperatures. Connecting this with protein inactivation measurements at low temperatures, we observe that the kinetics of protein unfolding is less sensitive to temperature change at the higher temperatures, which significantly departs from the Arrhenius behavior extrapolated from low temperatures. To account for this effect, we propose a reaction-diffusion model that modifies the temperature-dependence of protein inactivation by introducing a diffusion limit. Analysis of the reaction-diffusion model provides general guidelines in the behavior of protein inactivation (reaction-limited, transition, diffusion-limited) across a large temperature range from physiological temperature to extremely high temperatures. We further demonstrate that the reaction-diffusion model is particularly useful for designing optimal operating conditions for protein photoinactivation. The experimentally validated reaction-diffusion kinetics of protein unfolding is an important step toward understanding protein-inactivation kinetics over a large temperature range. It has important applications including molecular hyperthermia and calls for future studies to examine this model for other protein molecules.

Molecular Hyperthermia: Spatiotemporal Protein Unfolding and Inactivation by Nanosecond Plasmonic Heating

Spatiotemporal control of protein structure and activity in biological systems has important and broad implications in biomedical sciences as evidenced by recent advances in optogenetic approaches. Here, this study demonstrates that nanosecond pulsed laser heating of gold nanoparticles (GNP) leads to an ultrahigh and ultrashort temperature increase, coined as "molecular hyperthermia", which causes selective unfolding and inactivation of proteins adjacent to the GNP. Protein inactivation is highly dependent on both laser pulse energy and GNP size, and has a well-defined impact zone in the nanometer scale. It is anticipated that the fine control over protein structure and function enabled by this discovery will be highly enabling within a number of arenas, from probing the biophysics of protein folding/unfolding to the nanoscopic manipulation of biological systems via an optical trigger, to developing novel therapeutics for disease treatment without genetic modification.

Computational modeling of acrylodan-labeled cAMP dependent protein kinase catalytic subunit unfolding

Structure of the cAMP-dependent protein kinase catalytic subunit, where the asparagine residue 326 was replaced with acrylodan-cystein conjugate to implement this fluorescence reporter group into the enzyme, was modeled by molecular dynamics (MD) method and the positioning of the dye molecule in protein structure was characterized at temperatures 300K, 500K and 700K. It was found that the acrylodan moiety, which fluorescence is very sensitive to solvating properties of its microenvironment, was located on the surface of the native protein at 300K that enabled its partial solvation with water. At high temperatures the protein structure significantly changed, as the secondary and tertiary structure elements were unfolded and these changes were sensitively reflected in positioning of the dye molecule. At 700K complete unfolding of the protein occurred and the reporter group was entirely expelled into water. However, at 500K an intermediate of the protein unfolding process was formed, where the fluorescence reporter group was directed towards the protein interior and buried in the core of the formed molten globule state. This different positioning of the reporter group was in agreement with the two different shifts of emission spectrum of the covalently bound acrylodan, observed in the unfolding process of the protein.

Probing the Folding-Unfolding Transition of a Thermophilic Protein, MTH1880

The folding mechanism of typical proteins has been studied widely, while our understanding of the origin of the high stability of thermophilic proteins is still elusive. Of particular interest is how an atypical thermophilic protein with a novel fold maintains its structure and stability under extreme conditions. Folding-unfolding transitions of MTH1880, a thermophilic protein from Methanobacterium thermoautotrophicum, induced by heat, urea, and GdnHCl, were investigated using spectroscopic techniques including circular dichorism, fluorescence, NMR combined with molecular dynamics (MD) simulations. Our results suggest that MTH1880 undergoes a two-state N to D transition and it is extremely stable against temperature and denaturants. The reversibility of refolding was confirmed by spectroscopic methods and size exclusion chromatography. We found that the hyper-stability of the thermophilic MTH1880 protein originates from an extensive network of both electrostatic and hydrophobic interactions coordinated by the central β-sheet. Spectroscopic measurements, in combination with computational simulations, have helped to clarify the thermodynamic and structural basis for hyper-stability of the novel thermophilic protein MTH1880.

Correlated parameter fit of arrhenius model for thermal denaturation of proteins and cells

Thermal denaturation of proteins is critical to cell injury, food science and other biomaterial processing. For example protein denaturation correlates strongly with cell death by heating, and is increasingly of interest in focal thermal therapies of cancer and other diseases at temperatures which often exceed 50 °C. The Arrhenius model is a simple yet widely used model for both protein denaturation and cell injury. To establish the utility of the Arrhenius model for protein denaturation at 50 °C and above its sensitivities to the kinetic parameters (activation energy E a and frequency factor A) were carefully examined. We propose a simplified correlated parameter fit to the Arrhenius model by treating E a, as an independent fitting parameter and allowing A to follow dependently. The utility of the correlated parameter fit is demonstrated on thermal denaturation of proteins and cells from the literature as a validation, and new experimental measurements in our lab using FTIR spectroscopy to demonstrate broad applicability of this method. Finally, we demonstrate that the end-temperature within which the denaturation is measured is important and changes the kinetics. Specifically, higher E a and A parameters were found at low end-temperature (50 °C) and reduce as end-temperatures increase to 70 °C. This trend is consistent with Arrhenius parameters for cell injury in the literature that are significantly higher for clonogenics (45-50 °C) vs. membrane dye assays (60-70 °C). Future opportunities to monitor cell injury by spectroscopic measurement of protein denaturation are discussed.

Theoretical volume profiles as a tool for probing transition states: folding kinetics

The mechanism by which conformational changes, particularly folding and unfolding, occur in proteins and other biopolymers has been widely discussed in the literature. Molecular dynamics (MD) simulations of protein folding present a formidable challenge since these conformational changes occur on a time scale much longer than what can be afforded at the current level of computational technology. Transition state (TS) theory offers a more economic description of kinetic properties of a reaction system by relating them to the properties of the TS, or for flexible systems, the TS ensemble (TSE). The application of TS theory to protein folding is limited by ambiguity in the definition of the TSE for this process. We propose to identify the TSE for conformational changes in flexible systems by comparison of its experimentally determined volumetric property, known as the volume of activation, to the structure-specific volume profile of the process calculated using MD. We illustrate this approach by its successful application to unfolding of a model chain system.

A modified rat model of exercise-induced renal injury and the protective effects of losartan and yishen huanji decoction

We modified a rat model of exercise-induced renal injury by forcing 4-week-old male Sprague-Dawley (SD) rats (the MOD group) running on a treadmill for 8 weeks under conditions of high temperature, high humidity, bearing weight with some additional stimulations. Compared with the control (CON) group, the traditional running group (TRA), the losartan potassium intervention running group (LOS) and the traditional Chinese medicine "Yishen Huanji Decoction" intervention running group (CHI), the urinary, biochemistry indicators, the concentrations of angiotensin II (Ang II) were significantly higher in the MOD group than in the TRA group. After 3--4 weeks and 8-week training program, the 24-h urine protein and NAG levels in the LOS group and CHI group were lower than in the MOD group respectively. The BUN and SCr levels in the CHI group were lower than in the MOD group and higher than in the LOS group. AngII concentrations in the LOS group were higher than the MOD group. The modified rat renal injury model induced greater lesions than the traditional model. High temperatures, humidity and weight bearing were critical factors to induce Ang II activation, which can aggravate renal injury. Losartan potassium and the "Yishen Huanji Decoction" can protect against renal injury.

Computational Molecular Engineering as an Emerging Technology in Process Engineering

The present level of development of molecular force field methods is assessed from the point of view of simulation-based engineering, outlining the immediate perspective for further development and highlighting the newly emerging discipline of “Computational Molecular Engineering (CME)” which makes basic research in soft matter physics fruitful for industrial applications. Within the coming decade, major breakthroughs can be reached if a research focus is placed on processes at interfaces, combining aspects where an increase in the accessible length and time scales due to massively parallel high-performance computing will lead to particularly significant improvements.

Predicting folding free energy changes upon single point mutations

MOTIVATION

The folding free energy is an important characteristic of proteins stability and is directly related to protein's wild-type function. The changes of protein's stability due to naturally occurring mutations, missense mutations, are typically causing diseases. Single point mutations made in vitro are frequently used to assess the contribution of given amino acid to the stability of the protein. In both cases, it is desirable to predict the change of the folding free energy upon single point mutations in order to either provide insights of the molecular mechanism of the change or to design new experimental studies.

RESULTS

We report an approach that predicts the free energy change upon single point mutation by utilizing the 3D structure of the wild-type protein. It is based on variation of the molecular mechanics Generalized Born (MMGB) method, scaled with optimized parameters (sMMGB) and utilizing specific model of unfolded state. The corresponding mutations are built in silico and the predictions are tested against large dataset of 1109 mutations with experimentally measured changes of the folding free energy. Benchmarking resulted in root mean square deviation = 1.78 kcal/mol and slope of the linear regression fit between the experimental data and the calculations was 1.04. The sMMGB is compared with other leading methods of predicting folding free energy changes upon single mutations and results discussed with respect to various parameters.

AVAILABILITY

All the pdb files we used in this article can be downloaded from http://compbio.clemson.edu/downloadDir/mentaldisorders/sMMGB_pdb.rar.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Insights into the structural stability of Bax from molecular dynamics simulations at high temperatures

Bax is a member of the Bcl-2 protein family that participates in mitochondrion-mediated apoptosis. In the early stages of the apoptotic pathway, this protein migrates from the cytosol to the outer mitochondrial membrane, where it is inserted and usually oligomerizes, making cytochrome c-compatible pores. Although several cellular and structural studies have been reported, a description of the stability of Bax at the molecular level remains elusive. This article reports molecular dynamics simulations of monomeric Bax at 300, 400, and 500 K, focusing on the most relevant structural changes and relating them to biological experimental results. Bax gradually loses its α-helices when it is submitted to high temperatures, yet it maintains its globular conformation. The resistance of Bax to adopt an extended conformation could be due to several interactions that were found to be responsible for maintaining the structural stability of this protein. Among these interactions, we found salt bridges, hydrophobic interactions, and hydrogen bonds. Remarkably, salt bridges were the most relevant to prevent the elongation of the structure. In addition, the analysis of our results suggests which conformational movements are implicated in the activation/oligomerization of Bax. This atomistic description might have important implications for understanding the functionality and stability of Bax in vitro as well as within the cellular environment.