Protein folding: Molecular dynamics simulations and in vitro studies for probing mechanism of urea- and guanidinium chloride-induced unfolding of horse cytochrome-c

Investigating pH-induced conformational switch in PIM-1: An integrated multi spectroscopic and MD simulation study

PIM-1 is a Ser/Thr kinase, which has been extensively studied as a potential target for cancer therapy due to its significant roles in various cancers, including prostate and breast cancers. Given its importance in cancer, researchers are investigating the structure of PIM-1 for pharmacological inhibition to discover therapeutic intervention. This study examines structural and conformational changes in PIM-1 across different pH using various spectroscopic and computational techniques. Spectroscopic results indicate that PIM-1 maintains its secondary and tertiary structure within the pH range of 7.0-9.0. However, protein aggregation occurs in the acidic pH range of 5.0-6.0. Additionally, kinase assays suggested that PIM-1 activity is optimal within the pH range of 7.0-9.0. Subsequently, we performed a 100 ns all-atom molecular dynamics (MD) simulation to see the effect of pH on PIM-1 structural stability at the molecular level. MD simulation analysis revealed that PIM-1 retains its native conformation in alkaline conditions, with some residual fluctuations in acidic conditions as well. A strong correlation was observed between our MD simulation, spectroscopic, and enzymatic activity studies. Understanding the pH-dependent structural changes of PIM-1 can provide insights into its role in disease conditions and cellular homeostasis, particularly regarding protein function under varying pH conditions.

Photothermally Heated Asymmetric Thin Nanopores Suggest the Influence of Temperature on the Intermediate Conformational State of Cytochrome c in an Electric Field

Nanopore sensing is a label-free single-molecule technique that enables the study of the dynamical structural properties of proteins. Here, we detect the translocation of cytochrome c (Cyt c) through an asymmetric thin nanopore with photothermal heating to evaluate the influence of temperature on Cyt c conformation during its translocation in an electric field. Before Cyt c translocates through an asymmetric thin SiNx nanopore, ∼1 ms trapping events occur due to electric field-induced denaturation. These trapping events were corroborated by a control analysis with a transmission electron microscopy-drilled pore and denaturant buffer. Cyt c translocation events exhibited markedly greater broad current blockade when the pores were photothermally heated. Collectively, our molecular dynamics simulation predicted that an increased temperature facilitates denaturation of the α-helical structure of Cyt c, resulting in greater blockade current during Cyt c trapping. Our photothermal heating method can be used to study the influence of temperature on protein conformation at the single-molecule level in a label-free manner.

Evidence for paraquat-pepsin interaction: In vitro and silico study

The use of the herbicide paraquat (PQ) has raised concerns about potential environmental consequences due to its toxicity and persistence in the environment. Considering the affinity of dangerous compounds to biological molecules, it is necessary to know their binding properties. This article focuses on the behavior of the pepsin enzyme following its contact with paraquat poison, and the interaction between paraquat and pepsin has been investigated in laboratory conditions and simulated physiological conditions using multispectral techniques. Fluorescence experiments showed that PQ uses a static method to quench pepsin's intrinsic fluorescence. By causing structural damage to pepsin, PQ may be detrimental as it alters its conformational function based on FT-IR spectroscopy. The coupling reaction is a spontaneous process caused by hydrogen bonding and van der Waals forces according to the analysis of the thermodynamic parameters of each system at three different temperatures. The molecular structure of pepsin changes when it binds to PQ. Also, the results showed that PQ is a pepsin inhibitor that changes the function of the enzyme.

Suitability of Adenosine Derivatives in Improving the Activity and Stability of Cytochrome c under Stress: Insights into the Effect of Phosphate Groups

It is well known that adenosine and its phosphate derivatives play a crucial role in biological phenomena such as apoptosis and cell signaling and act as the energy currency of the cell. Although their interactions with various proteins and enzymes have been described, the focus of this work is to demonstrate the effect of the phosphate group on the activity and stability of the native heme metalloprotein cytochrome c (Cyt c), which is important from both biological and industrial aspects. In situ and in silico characterizations are used to correlate the relationship between the binding affinity of adenosine and its phosphate groups with unfolding behavior, corresponding peroxidase activities, and stability factors. Interaction of adenosine (ADN), adenosine monophosphate (AMP), adenosine 5'-diphosphate (ADP), and adenosine 5'-triphosphate (ATP) with Cyt c increases peroxidase-like activity by up to 1.8-6.5-fold compared to native Cyt c. This activity is significantly maintained even after multiple stress conditions such as oxidative stress and the presence of a chaotropic agent such as guanidine hydrochloride (GuHCl). With binding affinities on the order of ADN < AMP < ADP < ATP, adenosine derivatives were found to stabilize Cyt c by varying the secondary structural features of the protein. Thus, in addition to being a fundamental study, the current work also proposes a way of stabilizing protein systems to be used for real-time biocatalytic applications.

Integrated spectroscopic and MD simulation approach to decipher the effect of pH on the structure function of Staphylococcus aureus thymidine kinase

Staphylococcus aureus is a major human pathogen responsible for a variety of clinical infections, becoming increasingly resistant to antibiotics. To address this challenge, there is a need to identify new cellular targets and innovative approaches to expand treatment options. One such target is thymidine kinase (TK), a crucial enzyme in the pyrimidine salvage pathway, which plays a key role in the phosphorylation of thymidine, an essential component in DNA synthesis and repair. In this study, we have successfully cloned, expressed, and purified the TK protein. A comprehensive investigation into how different pH levels affect the structure and functional activity of TK, using a combination of spectroscopy, classical molecular dynamics simulations, and enzyme activity assays was conducted. Our study revealed that variation in pH disrupts secondary and tertiary structures of TK with noticeable aggregate formation at pH 5.0. Enzyme activity studies demonstrated that TK exhibited its maximum kinase activity within the physiological pH range. These findings strongly suggest a connection between structural changes and enzymatic activity, which was further supported by the agreement between the spectroscopic features we measured and the results of our MD simulations. Our study provides a deeper insight into the structural features of TK, which could potentially be harnessed for the development of therapeutic strategies aimed at combatting infectious diseases. Conformational dynamics plays an essential role in the design and development of effective inhibitors. Considering the effects of pH on the conformational dynamics of TK, our findings may be implicated in the development of potent and selective inhibitors.Communicated by Ramaswamy H. Sarma.

Evaluating the Cysteine-Rich and Catalytic Subdomains of Human Tyrosinase and OCA1-Related Mutants Using 1 μs Molecular Dynamics Simulation

The inherited disorder oculocutaneous albinism type 1 (OCA1) is caused by mutations in the TYR gene encoding tyrosinase (Tyr), an enzyme essential to producing pigments throughout the human body. The intramelanosomal domain of Tyr consists of the cysteine-rich and tyrosinase catalytic subdomains, which are essential for enzymatic activity. In protein unfolding, the roles of these subdomains are not well established. Here, we performed six molecular dynamics simulations at room temperature for Tyr and OCA1-related mutant variants P406L and R402Q intramelanosomal domains. The proteins were simulated for 1 μs in water and urea to induce unfolding. In urea, we observed increases in surface area, decreases in intramolecular hydrogen bonding, and decreases in hydrophobic interactions, suggesting a 'molten globule' state for each protein. Between all conditions, the cysteine-rich subdomain remains stable, whereas the catalytic subdomain shows increased flexibility. This flexibility is intensified by the P406L mutation, while R402Q increases the catalytic domain's rigidity. The cysteine-rich subdomain is rigid, preventing the protein from unfolding, whereas the flexibility of the catalytic subdomain accommodates mutational changes that could inhibit activity. These findings match the conclusions from our experimental work suggesting the function alteration by the P406L mutation, and the potential role of R402Q as a polymorphism.

Apo-metallothionein-3 cooperatively forms tightly compact structures under physiological conditions

Metallothioneins (MTs) are essential mammalian metal chaperones. MT isoform 1 (MT1) is expressed in the kidneys and isoform 3 (MT3) is expressed in nervous tissue. For MTs, the solution-based NMR structure was determined for metal-bound MT1 and MT2, and only one X-ray diffraction structure on a crystallized mixed metal-bound MT2 has been reported. The structure of solution-based metalated MT3 is partially known using NMR methods; however, little is known about the fluxional de novo apo-MT3 because the structure cannot be determined by traditional methods. Here, we used cysteine modification coupled with electrospray ionization mass spectrometry, denaturing reactions with guanidinium chloride, stopped-flow methods measuring cysteine modification and metalation, and ion mobility mass spectrometry to reveal that apo-MT3 adopts a compact structure under physiological conditions and an extended structure under denaturing conditions, with no intermediates. Compared with apo-MT1, we found that this compact apo-MT3 binds to a cysteine modifier more cooperatively at equilibrium and 0.5 times the rate, providing quantitative evidence that many of the 20 cysteines of apo-MT3 are less accessible than those of apo-MT1. In addition, this compact apo-MT3 can be identified as a distinct population using ion mobility mass spectrometry. Furthermore, proposed structural models can be calculated using molecular dynamics methods. Collectively, these findings provide support for MT3 acting as a noninducible regulator of the nervous system compared with MT1 as an inducible scavenger of trace metals and toxic metals in the kidneys.

Model architectures for bacterial membranes

The complex composition of bacterial membranes has a significant impact on the understanding of pathogen function and their development towards antibiotic resistance. In addition to the inherent complexity and biosafety risks of studying biological pathogen membranes, the continual rise of antibiotic resistance and its significant economical and clinical consequences has motivated the development of numerous in vitro model membrane systems with tuneable compositions, geometries, and sizes. Approaches discussed in this review include liposomes, solid-supported bilayers, and computational simulations which have been used to explore various processes including drug-membrane interactions, lipid-protein interactions, host-pathogen interactions, and structure-induced bacterial pathogenesis. The advantages, limitations, and applicable analytical tools of all architectures are summarised with a perspective for future research efforts in architectural improvement and elucidation of resistance development strategies and membrane-targeting antibiotic mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-021-00913-7.

Insights to Human γD-Crystallin Unfolding by NMR Spectroscopy and Molecular Dynamics Simulations

Human γD-crystallin (HGDC) is an abundant lens protein residing in the nucleus of the human lens. Aggregation of this and other structural proteins within the lens leads to the development of cataract. Much has been explored on the stability and aggregation of HGDC and where detailed investigation at the atomic resolution was needed, the X-ray structure was used as an initial starting conformer for molecular modeling. In this study, we implemented NMR-solution HGDC structures as starting conformers for molecular dynamics simulations to provide the missing pieces of the puzzle on the very early stages of HGDC unfolding leading up to the domain swap theories proposed by past studies. The high-resolution details of the conformational dynamics also revealed additional insights to possible early intervention for cataractogenesis.

Effects of natural mutations (L94I and L94V) on the stability and mechanism of folding of horse cytochrome c: A combined in vitro and molecular dynamics simulations approach

Known crystal structures of 10 cytochromes (cyts) c from different sources led to the conclusion that natural mutations in these proteins does not affect their 3D structure, hence evolution preserved structure for function. A sequence alignment of horse cyt c with all other 284 cyts c led to two important conclusions: (i) Leu at position 94 is conserved in all 30 mammalian known sequences, and (ii) there are 14 other species which have either Val or Ile at 94th position. We asked a question: Is the avoidance of substitution by Val or Ile at position 94 in the mammalian cyts c by design or by chance? To answer this question, we introduced natural substitutes of Leu94 by Val and Ile in horse cyt c using site-directed mutagenesis. Here, from our in vitro and molecular dynamic simulation studies on L94V and L94I mutants, we concluded that (i) although the natural mutations destabilize the wild type cyt c, it does not significantly affect the mechanism of folding of the protein, (ii) urea-induced denaturation of WT cyt c and its mutants is a two-state process, and (iii) denaturation of WT cyt c and its mutants by guanidinium chloride is not a two-state process.

Impact of amino acid substitution in the kinase domain of Bruton tyrosine kinase and its association with X-linked agammaglobulinemia

X-linked agammaglobulinemia (XLA) is a rare disease that affects the immune system, characterized by a serial development of bacterial infection from the onset of infantile age. Bruton tyrosine kinase (BTK) is a non-receptor cytoplasmic kinase that plays a crucial role in the B-lymphocyte maturation. The altered expression, mutation and/or structural variations of BTK are responsible for causing XLA. Here, we have performed extensive sequence and structure analyses of BTK to find deleterious variations and their pathogenic association with XLA. First, we screened the pathogenic variations in the BTK from a pool of publicly available resources, and their pathogenicity/tolerance and stability predictions were carried out. Finally, two pathogenic variations (E589G and M630K) were studied in detail and subjected to all-atom molecular dynamics simulation for 200 ns. Intramolecular hydrogen bonds (H-bonds), secondary structure, and principal component analysis revealed significant conformational changes in variants that support the structural basis of BTK dysfunction in XLA. The free energy landscape analysis revealed the presence of multiple energy minima, suggests that E589G brings a large destabilization and consequently unfolding behavior compared to M630K. Overall, our study suggests that amino acid substitutions, E589G, and M630K, significantly alter the structural conformation and stability of BTK.

Study on effects of co-solvents on the structure of DhaA by molecular dynamics simulation

With the increasing application of enzymes in various research fields, the choices of co-solvents in enzymatic preparations which directly related to the catalytic activity have been attracted attention. Thus, researching on the stabilization or destabilization behaviors of enzymes in different solvents is extremely essential. In this study, the structural changes of DhaA in two typical aprotic co-solvents (acetonitrile and tetrahydrofuran) were firstly investigated by molecular dynamics (MD) simulation. The simulation results revealed the strong van der Waals force between co-solvents and DhaA which could induce the structural change of enzyme. Interestingly, the differences of molecular size and the electrostatic force with enzyme of two co-solvents led to quite different influences on DhaA. As for acetonitrile, solvent molecules could penetrate into the catalytic site of DhaA which promoted by the electrostatic interaction. On the contrary, tetrahydrofuran molecules were mainly distributed around the catalytic site due to the relative weak electrostatic interaction and steric resistance effect. It can be concluded that different co-solvent can affect the key domains, substrate pathway and catalytic pocket of DhaA.Communicated by Ramaswamy H. Sarma.

Effects of the deglycosylation on the structure and activity of chloroperoxidase: Molecular dynamics simulation approach

Chloroperoxidase (CPO) is a versatile fungal heme-thiolate protein that catalyzes a variety of one electron and two-electron oxidations. Chloroperoxidase is a versatile fungal heme-thiolate protein that catalyzes a variety of oxidations. CPO enzyme contains thirteen sugars, including five N-acetyl D-glucosamines (NAG) and eight mannoses (MAN), which are attached to the protein via the glycosidic bonds. Removal of the sugars from CPO leads to increase the hydrophobicity of the enzyme, as well as the reduction of the alkylation reactions. However, due to the lack of the proper force field for the sugars, they are ignored in the theoretical studies. The present study aims to assess the effects of the sugar segments on the structure and activity of CPO through the simulation of the halo structure and the structures without the sugar segment. Despite the difficulty of the process and being time-consuming, the suitable force field is introduced successfully for the sugars. According to molecular dynamics simulation (MD), seven channels and fifteen cavities are identified in the CPO structure. Two of the channels provide the substrate access to the active site. The MD simulation results reveal that the removal of NAG decreases the number of the cavities from fifteen to eleven. Besides, the removal of NAG is associated with removing the channel providing the substrate access. The number of the cavities decreases from fifteen to fourteen through the removal of MAN; however, channel providing the substrate access to the active site is partly preserved. The MD simulation results indicate that the structures without the sugar units are more compact in comparison with the halo structures. The removal of the sugar segments induces the significant changes in the flexibility of the residues that affect the catalytic activity of the enzyme. As a result, the enzyme activities, such as the oxidation, alkylation, halogenation, and epoxidation cannot occur when the sugar segments of the enzyme are removed.

Effect of pH on the structure and function of pyruvate dehydrogenase kinase 3: Combined spectroscopic and MD simulation studies

Pyruvate dehydrogenase kinase-3 (PDK3) plays important role in the glucose metabolism and is associated with cancer progression, and thus being considered as an attractive target for cancer therapy. In this study, we employed spectroscopic techniques to study the structural and conformational changes in the PDK3 at varying pH conditions ranging from pH 2.0 to 12.0. UV/Vis, fluorescence and circular dichroism spectroscopic measurements revealed that PDK3 maintains its native-like structure (both secondary and tertiary) in the alkaline conditions (pH 7.0-12.0). However, a significant loss in the structure was observed under acidic conditions (pH 2.0-6.0). The propensity of aggregate formation at pH 4.0 was estimated by thioflavin T fluorescence measurements. To further complement structural data, kinase activity assay was performed, and maximum activity of PDK3 was observed at pH 7.0-8.0 range; whereas, its activity was lost under acidic pH. To further see conformational changes at atomistic level we have performed all-atom molecular dynamics at different pH conditions for 150 ns. A well defined correlation was observed between experimental and computational studies. This work highlights the significance of structural dependence of pH for wide implications in protein-protein interaction, biological function and drug design procedures.