Effect of external stresses on protein conformation: a computer modelling study

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

Progress in Simulation Studies of Insulin Structure and Function

Insulin is a peptide hormone known for chiefly regulating glucose level in blood among several other metabolic processes. Insulin remains the most effective drug for treating diabetes mellitus. Insulin is synthesized in the pancreatic β-cells where it exists in a compact hexameric architecture although its biologically active form is monomeric. Insulin exhibits a sequence of conformational variations during the transition from the hexamer state to its biologically-active monomer state. The structural transitions and the mechanism of action of insulin have been investigated using several experimental and computational methods. This review primarily highlights the contributions of molecular dynamics (MD) simulations in elucidating the atomic-level details of conformational dynamics in insulin, where the structure of the hormone has been probed as a monomer, dimer, and hexamer. The effect of solvent, pH, temperature, and pressure have been probed at the microscopic scale. Given the focus of this review on the structure of the hormone, simulation studies involving interactions between the hormone and its receptor are only briefly highlighted, and studies on other related peptides (e.g., insulin-like growth factors) are not discussed. However, the review highlights conformational dynamics underlying the activities of reported insulin analogs and mimetics. The future prospects for computational methods in developing promising synthetic insulin analogs are also briefly highlighted.

Effect of thermal and microwave processing on secondary structure of bovine β-lactoglobulin: A molecular modeling study

Milk allergy is known to cause severe allergic reactions in hypersensitive patients, especially in infants and children. β-Lactoglobulin is one of the major allergens in bovine milk. The influence of thermal and microwave processing on the structural deviations of β-lactoglobulin protein have been studied using molecular modeling techniques. The structural deviations are studied using root mean square deviations, radius of gyration, dipole moment, and solvent accessible surface area. STRIDE analysis showed significant changes in the β-lactoglobulin, especially when oscillating electric fields were applied along with heat. Root mean square fluctuations (RMSF) has been assessed for known epitopes in the β-lactoglobulin molecule. This showed that when the protein is exposed to certain thermal stress, it compacts by burying hydrophobic residues in the core. However, few allergic epitope residues also exhibit increased RMSF leading to higher reactive sites on the surface of the protein molecule. PRACTICAL APPLICATIONS: This study showed that molecular modeling can be used to gain valuable insights regarding the structural changes during processing. In the future, with more computational capacity, it can be used to make comparison between results obtained from simulations and real-time experiments. The current techniques used in food industries such as Nuclear Magnetic Resonance Imaging, Fourier Transformation Infrared Spectroscopy, X-ray diffraction can analyze pre- and post-processing effects. Hence, it become necessary to understand the changes that takes place during the processing techniques. Molecular dynamic simulation could be a useful technique in analyzing the changes occurring during the processing.

The Enigma of Amyloid Forming Proteins: Insights From Molecular Simulations

Molecular level insight into the interplay between protein sequence, structure, and conformational dynamics is crucial for the comprehensive understanding of protein folding, misfolding, and aggregation phenomena that are pertinent to the formation of amyloid fibrils implicated in several degenerative diseases. Computational modelling provides insight into protein behaviour at spatial and temporal resolution still largely outside the reach of experiments. Herein we present an account of our theoretical modelling research conducted in collaboration with several experimental groups where we explored the effects of local environment on the structure and aggregation propensity of several types of amyloidogenic peptides and proteins, including apolipoprotein C-II, insulin, amylin, and amyloid-β using a variety of computational approaches.

Deconstruction of the human connexin 26 hemichannel due to an applied electric field; A molecular dynamics simulation study

Connexins are a 21-member membrane protein family constituting channels evolved in direct communication between adjacent cells by passaging cytoplasmic molecules and ions. Hexametrical assembly of connexin proteins in plasma membrane forms a wide aqueous pore known as connexin hemichannel. These hemichannels mediate cytoplasm and extracellular milieu communication both in many external tissues and in the central nervous system. In this study, a series of molecular dynamics simulations has been performed to investigate the effect of applied static and alternating electric fields on the stability and conformation of human connexin26 hemichannel. The root mean square deviations of C-alpha atoms, the dipole moment distribution, the number of inter-protein hydrogen bonds and the number of water-protein hydrogen bonds were used to assess connexin26 hemichannel stability. In the static field case, our results show that although the lowest field used in this study (0.1V/nm) does not lead to the hemichannel deconstruction, stronger fields (>0.1V/nm), however, disrupt the protein structure during the simulations time period. Furthermore, in the alternating case, compared to static field case, field effects on the connexin26 hemichannel conformation are reduced and consequently the protein maintains its native structure for longer times. Specifically, for the highest frequency used in this study (50GHz), the hemichannel keeps its structure even under the effect of the strongest field (0.4V/nm). According to our results, the protein secondary structure is preserved in the characteristic times determined for the protein deconstruction. Consequently, we suggest that the protein deconstruction is due to the tertiary and quaternary structure loss.

Electromagnetic-field effects on structure and dynamics of amyloidogenic peptides

Electromagnetic fields (EMFs) are ever-present, and so is the need to better understand their influence on human health and biological matter in general. The interaction between a molecular system and external EMF can alter the structure, and dynamical behaviour, and, hence, biological function of proteins with uncertain health consequences. This urges a detailed investigation of EMF-induced effects on basic protein biophysics. Here, we used all-atom non-equilibrium molecular dynamics simulations to understand and quantify the response mechanisms of the amyloidogenic apoC-II(60-70) peptides to non-ionising radiation by modelling their behaviour under external electromagnetic and electric fields of different strengths. Our simulations show high strength fields (>0.04 V/nm) cause structural changes in apoC-II(60-70) due to the peptide dipole alignment along the applied field direction, which disrupts the inherent β-hairpin conformation known to be the intermediate state for fibril formation. The intermediate field-strength range (0.04-0.004 V/nm) causes a significant acceleration in peptide dynamics, which leads to the increased population of structures with fibril-inhibiting characteristics, such as the separated N- and C-termini and colocation of the aromatic residues at the same peptide face. In contrast, lower field strengths (<0.004 V/nm) promote the formation of the amyloid-prone hairpin structures relative to the ambient conditions. These findings suggest that intermediate-strength electromagnetic fields could be considered for designing alternative treatments of amyloid diseases, while the very high and low field strengths could be employed for engineering well-ordered fibrillar aggregates for non-medicinal applications.

Theoretical and computational studies of peptides and receptors of the insulin family

Synergistic interactions among peptides and receptors of the insulin family are required for glucose homeostasis, normal cellular growth and development, proliferation, differentiation and other metabolic processes. The peptides of the insulin family are disulfide-linked single or dual-chain proteins, while receptors are ligand-activated transmembrane glycoproteins of the receptor tyrosine kinase (RTK) superfamily. Binding of ligands to the extracellular domains of receptors is known to initiate signaling via activation of intracellular kinase domains. While the structure of insulin has been known since 1969, recent decades have seen remarkable progress on the structural biology of apo and liganded receptor fragments. Here, we review how this useful structural information (on ligands and receptors) has enabled large-scale atomically-resolved simulations to elucidate the conformational dynamics of these biomolecules. Particularly, applications of molecular dynamics (MD) and Monte Carlo (MC) simulation methods are discussed in various contexts, including studies of isolated ligands, apo-receptors, ligand/receptor complexes and intracellular kinase domains. The review concludes with a brief overview and future outlook for modeling and computational studies in this family of proteins.

Electric field effects on alanine tripeptide in sodium halide solutions

The electric field effects on conformational properties of trialanine in different halide solutions were explored with long-scale molecular dynamics simulations. NaF, NaCl, NaBr and NaI solutions of low (0.2 M) and high (2 M) concentrations were exposed to a constant electric field of 1000 V/m. Generally, the electric field does not disturb trialanine's structure. Large structural changes appear only in the case of the supersaturated 2.0 M NaF solution containing NaF crystals. Although the electric field affects in a complex way, all the ions-water-peptide interactions, it predominantly affects the electroselectivity effect, which describes specific interactions such as the ion-pair formation.

A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations

The anti-apoptotic B-cell lymphoma 2 (Bcl-2) protein interacts with several proteins that regulate the apoptotic properties of cells. In this research, we conduct several all-atom molecular dynamics (MD) simulations under high-temperature unfolding conditions, from 400 to 800 K, for 25 ns. These simulations were performed using a model of an engineered Bcl-2 human protein (Bcl-2-Δ22Σ3), which lacks 22 C-terminal residues of the transmembrane domain. The aim of this study is to gain insight into the structural behavior of Bcl-2-Δ22Σ3 by mapping the conformational movements involved in Bcl-2 stability and its biological function. To build a Bcl-2-Δ22Σ3 three-dimensional model, the protein core was built by homology modeling and the flexible loop domain (FLD, residues 33-91) by ab initio methods. Further, the entire protein model was refined by MD simulations. Afterwards, the production MD simulations showed that the FLD at 400 and 500 K has several conformations reaching into the protein core, whereas at 600 K some of the alpha-helices were lost. At 800 K, the Bcl-2 core is destabilized suggesting a possible mechanism for protein unfolding, where the alpha helices 1 and 6 were the most stable, and a reduction in the number of hydrogen bonds initially occurs. In conclusion, the structural changes and the internal protein interactions suggest that the core and the FLD are crucial components of Bcl-2 in its function of regulate ng access to the recognition sites of kinases and caspases.

Effect of External Electric Field Stress on Gliadin Protein Conformation

A molecular dynamic (MD) modeling approach was applied to evaluate the effect of external electric field on gliadin protein structure and surface properties. Static electric field strengths of 0.001 V/nm and 0.002 V/nm induced conformational changes in the protein but had no significant effect on its surface properties. The study of hydrogen bond evolution during the course of simulation revealed that the root mean square deviation, radius of gyration and secondary structure formation, all depend significantly on the number hydrogen bonds formed. This study demonstrated that it is necessary to gain insight into protein dynamics under external electric field stress, in order to develop the novel food processing techniques that can be potentially used to reduce or eradicate food allergens.

Soybean hydrophobic protein response to external electric field: a molecular modeling approach

The molecular dynamic (MD) modeling approach was applied to evaluate the effect of an external electric field on soybean hydrophobic protein and surface properties. Nominal electric field strengths of 0.002 V/nm and 0.004 V/nm had no major effect on the structure and surface properties of the protein isolate but the higher electric field strength of 3 V/nm significantly affected the protein conformation and solvent accessible surface area. The response of protein isolate to various external field stresses demonstrated that it is necessary to gain insight into protein dynamics under electromagnetic fields in order to be able to develop the techniques utilizing them for food processing and other biological applications.

Insulin dimer dissociation and unfolding revealed by amide I two-dimensional infrared spectroscopy

A structurally sensitive probe of the monomer/dimer equilibrium of insulin was developed using 2DIR spectroscopy and interpreted using calculated spectra.

Exploring the folding free energy landscape of insulin using bias exchange metadynamics

The bias exchange metadynamics (BE-META) technique was applied to investigate the folding mechanism of insulin, one of the most studied and biologically important proteins. The BE-META simulations were performed starting from an extended conformation of chain B of insulin, using only eight replicas and seven reaction coordinates. The folded state, together with the intermediate states along the folding pathway were identified and their free energy was determined. Three main basins were found separated from one another by a large free energy barrier. The characteristic native fold of chain B was observed in one basin, while the other two most populated basins contained "molten-globule" conformations stabilized by electrostatic and hydrophobic interactions, respectively. Transitions between the three basins occur on the microsecond time scale. The implications and relevance of this finding to the folding mechanisms of insulin were investigated.

Electric field effects on insulin chain-B conformation

The response of proteins to different forms of stress continues to be a topic of major interest, especially with the proliferation of electromagnetic devices conjectured to have detrimental effects on human health. In this paper, we have performed molecular dynamics simulations on insulin chain-B under the influence of both static and oscillating electric fields, ranging from 10(7) to 10(9) V/m. We have found that both variants have an effect on the normal behavior of the protein, with oscillating fields being more disruptive to the structure as compared to static fields of similar effective strength. The application of a static field had a stabilizing effect on the secondary structure, restricting the inherent flexibility that is crucial for insulin's biological activity.

Effect of frequency on insulin response to electric field stress

There are many unanswered questions regarding the precise way in which proteins respond to external stress. Since the function of proteins is critically linked to their three-dimensional structures, exposure to any form of stress which may induce changes in conformation can potentially initiate severe cellular dysfunction. This is particularly relevant with regard to the increasing presence of electromagnetic devices in today's environment and the possible effects on human health. Previously, we investigated the effect of electric field of various strengths on insulin chain-B under static and oscillating conditions. This paper expands on our previous work by subjecting the peptide to an oscillating electric field of different frequencies. We observed a frequency-dependent effect where the application of lower-frequency oscillating fields resulted in static-field-like behavior of the peptide, whereby the intrinsic flexibility of the protein is constrained, thus potentially restricting access to the protein's active state.

Protein flexibility: multiple molecular dynamics simulations of insulin chain B

Multiple molecular dynamics simulations totaling more than 100 ns were performed on chain B of insulin in explicit solvent at 300 K and 400 K. Despite some individual variations, a comparison of the protein dynamics of each simulation showed similar trends and most structures were consistent with NMR experimental values, even at the elevated temperature. The importance of packing interactions in determining the conformational transitions of the protein was observed, sometimes resulting in conformations induced by localized hydrophobic interactions. The high temperature simulation generated a more diverse range of structures with similar elements of secondary structure and populated conformations to the simulations at room temperature. A broad sampling of the conformational space of insulin chain B illustrated a wide range of conformational states with many transitions at room temperature in addition to the conformational states observed experimentally. The T-state conformation associated with insulin activity was consistently present and a possible mechanism of behavior was suggested.