Flexible Residues Research Articles

Preliminary Structural Data Revealed That the SARS-CoV-2 B.1.617 Variant's RBD Binds to ACE2 Receptor Stronger Than the Wild Type to Enhance the Infectivity.

The evolution of new SARS‐CoV‐2 variants around the globe has made the COVID‐19 pandemic more worrisome, further pressuring the health care system and immunity. Novel variations that are unique to the receptor‐binding motif (RBM) of the receptor‐binding domain (RBD) spike glycoprotein, i. e. L452R‐E484Q, may play a different role in the B.1.617 (also known as G/452R.V3) variant's pathogenicity and better survival compared to the wild type. Therefore, a thorough analysis is needed to understand the impact of these mutations on binding with host receptor (RBD) and to guide new therapeutics development. In this study, we used structural and biomolecular simulation techniques to explore the impact of specific mutations (L452R‐E484Q) in the B.1.617 variant on the binding of RBD to the host receptor ACE2. Our analysis revealed that the B.1.617 variant possesses different dynamic behaviours by altering dynamic‐stability, residual flexibility and structural compactness. Moreover, the new variant had altered the bonding network and structural‐dynamics properties significantly. MM/GBSA technique was used, which further established the binding differences between the wild type and B.1.617 variant. In conclusion, this study provides a strong impetus to develop novel drugs against the new SARS‐CoV‐2 variants.

Effects of Ribosomal Protein S10 Flexible Loop Mutations on Tetracycline and Tigecycline Susceptibility of Escherichia coli.

Tigecycline is a tetracycline derivative that is being used as an antibiotic of last resort. Both tigecycline and tetracycline bind to the small (30S) ribosomal subunit and inhibit translation. Target mutations leading to resistance to these antibiotics have been identified both in the 16S ribosomal RNA and in ribosomal proteins S3 and S10 (encoded by the rpsJ gene). Several different mutations in the S10 flexible loop tip residue valine 57 (V57) have been observed in tigecycline-resistant Escherichia coli isolates. However, the role of these mutations in E. coli has not yet been characterized in a defined genetic background. In this study, we chromosomally integrated 10 different rpsJ mutations into E. coli, resulting in different exchanges or a deletion of S10 V57, and investigated the effects of the mutations on growth and tigecycline/tetracycline resistance. While one exchange, V57K, decreased the minimal inhibitory concentration (MIC) (Etest) to tetracycline to 0.75 μg/ml (compared to 2 μg/ml in the parent strain) and hence resulted in hypersensitivity to tetracycline, most exchanges, including the ones reported previously in resistant isolates (V57L, V57D, and V57I) resulted in slightly increased MICs to tigecycline and tetracycline. The strongest increase was observed for the V57L mutant, with a MIC (Etest) to tigecycline of 0.5 μg/ml (compared to 0.125 μg/ml in the parent strain) and a MIC to tetracycline of 4.0 μg/ml. Nevertheless, none of these exchanges increased the MIC to the extent observed in previously described clinical tigecycline-resistant isolates. We conclude that, next to S10 mutations, additional mutations are necessary in order to reach high-level tigecycline resistance in E. coli. In addition, our data reveal that mutants carrying S10 V57 exchanges or deletion display growth defects and, in most cases, also thermosensitivity. The defects are particularly strong in the V57 deletion mutant, which is additionally cold-sensitive. We hypothesize that the S10 loop tip residue is critical for the correct functioning of S10. Both the S10 flexible loop and tigecycline are in contact with helix h31 of the 16S rRNA. We speculate that exchanges or deletion of V57 alter the positioning of h31, thereby influencing both tigecycline binding and S10 function.

Dual-Knockout of Mutant Isocitrate Dehydrogenase 1 and 2 Subtypes Towards Glioma Therapy: Structural Mechanistic Insights on the Role of Vorasidenib.

Recently, Vorasidenib (AG-881) has been reported as a therapeutic alternative that exerts potent dual inhibitory activity against mIDH1/2 towards the treatment of low-grade glioma. However, structural and dynamic events associated with its dual inhibition mechanism remain unclear. As such, we employ integrative computer-assisted atomistic techniques to provide thorough structural and dynamic insights. Our analysis proved that the dual-targeting ability of AG-881 is mediated by Val255/Val294 within the binding pockets of both mIDH1 and mIDH2 which are shown to elicit a strong intermolecular interaction, thus favoring binding affinity. The structural orientations of AG-881 within the respective hydrophobic pockets allowed favorable interactions with binding site residues which accounted for its high binding free energy of -28.69 kcal/mol and -19.89 kcal/mol towards mIDH1 and mIDH2, respectively. Interestingly, upon binding, AG-881 was found to trigger systemic alterations of mIDH1 and mIDH2 characterized by restricted residue flexibility and a reduction in exposure of residues to the solvent surface area. As a result of these structural alterations, crucial interactions of the mutant enzymes were inhibited, a phenomenon that results in a suppression of the production of oncogenic stimulator 2-HG. Findings therefore provide thorough structural and dynamic insights associated with the dual inhibitory activity of AG-881 towards glioma therapy.

Structure and Thermal Stability of wtRop and RM6 Proteins through All-Atom Molecular Dynamics Simulations and Experiments

In the current work we study, via molecular simulations and experiments, the folding and stability of proteins from the tertiary motif of 4-α-helical bundles, a recurrent motif consisting of four amphipathic α-helices packed in a parallel or antiparallel fashion. The focus is on the role of the loop region in the structure and the properties of the wild-type Rop (wtRop) and RM6 proteins, exploring the key factors which can affect them, through all-atom molecular dynamics (MD) simulations and supporting by experimental findings. A detailed investigation of structural and conformational properties of wtRop and its RM6 loopless mutation is presented, which display different physical characteristics even in their native states. Then, the thermal stability of both proteins is explored showing RM6 as more thermostable than wtRop through all studied measures. Deviations from native structures are detected mostly in tails and loop regions and most flexible residues are indicated. Decrease of hydrogen bonds with the increase of temperature is observed, as well as reduction of hydrophobic contacts in both proteins. Experimental data from circular dichroism spectroscopy (CD), are also presented, highlighting the effect of temperature on the structural integrity of wtRop and RM6. The central goal of this study is to explore on the atomic level how a protein mutation can cause major changes in its physical properties, like its structural stability.

Expanding the SiMPl Plasmid Toolbox for Use with Spectinomycin/Streptomycin.

We recently developed the SiMPl plasmid toolbox, which is constituted by pairs of plasmids, generically indicated as pSiMPlx_N and pSiMPlx_C, which can be stably maintained in Escherichia coli with a single antibiotic x. The method exploits the split intein gp41–1 to reconstitute the enzyme conferring resistance toward the antibiotic x, whereby each enzyme fragment is expressed from one of the plasmids in the pair. pSiMPl plasmids are currently available for use with ampicillin, kanamycin, chloramphenicol, hygromycin, and puromycin. Here, we introduce another pair for use with spectinomycin/streptomycin, broadening the application spectrum of the SiMPl toolbox. To find functional splice sites in aminoglycoside adenylyltransferase, we apply a streamlined strategy looking exclusively at the flexibility of native cysteine and serine residues, which we first validated splitting the enzymes conferring resistance toward ampicillin, kanamycin, chloramphenicol, and hygromycin. This strategy could be used in the future to split other enzymes conferring resistance toward antibiotics.

Structural recognition of the MYC promoter G-quadruplex by a quinoline derivative: insights into molecular targeting of parallel G-quadruplexes.

DNA G-Quadruplexes (G4s) formed in oncogene promoters regulate transcription. The oncogene MYC promoter G4 (MycG4) is the most prevalent G4 in human cancers. However, the most studied MycG4 sequence bears a mutated 3′-residue crucial for ligand recognition. Here, we report a new drug-like small molecule PEQ without a large aromatic moiety that specifically binds MycG4. We determined the NMR solution structures of the wild-type MycG4 and its 2:1 PEQ complex, as well as the structure of the 2:1 PEQ complex of the widely used mutant MycG4. Comparison of the two complex structures demonstrates specific molecular recognition of MycG4 and shows the clear effect of the critical 3′-mutation on the drug binding interface. We performed a systematic analysis of the four available complex structures involving the same mutant MycG4, which can be considered a model system for parallel G4s, and revealed for the first time that the flexible flanking residues are recruited in a conserved and sequence-specific way, as well as unused potential for selective ligand-G4 hydrogen-bond interactions. Our results provide the true molecular basis for MycG4-targeting drugs and new critical insights into future rational design of drugs targeting MycG4 and parallel G4s that are prevalent in promoter and RNA G4s.

Characterization of SdGA, a cold-adapted glucoamylase from Saccharophagus degradans

We investigated the structural and functional properties of SdGA, a glucoamylase (GA) from Saccharophagus degradans, a marine bacterium which degrades different complex polysaccharides at high rate. SdGA is composed mainly by a N-terminal GH15_N domain linked to a C-terminal catalytic domain (CD) found in the GH15 family of glycosylhydrolases with an overall structure similar to other bacterial GAs. The protein was expressed in Escherichia coli cells, purified and its biochemical properties were investigated. Although SdGA has a maximum activity at 39 °C and pH 6.0, it also shows high activity in a wide range, from low to mild temperatures, like cold-adapted enzymes. Furthermore, SdGA has a higher content of flexible residues and a larger CD due to various amino acid insertions compared to other thermostable GAs. We propose that this novel SdGA, is a cold-adapted enzyme that might be suitable for use in different industrial processes that require enzymes which act at low or medium temperatures.

An iterative scheme of flexibility‐based component mode synthesis with higher‐order residual modal compensation

AbstractAn Iterative Flexibility‐based Component Mode Synthesis (IF‐CMS) is presented for the model reduction of partitioned structural dynamic systems via localized Lagrange multipliers. A distinct IF‐CMS feature is the inclusion of hierarchical residual flexibility at each iteration, resulting in considerably improved accuracy compared with the classical Craig–Bampton (CB) CMS method. The present IF‐CMS method does not increase the number of substructural modes during the hierarchical iteration process. In particular, to alleviate ill‐conditioning during the iteration steps, the proposed IF‐CMS method adopts the so‐called (qd, , ub)‐formulation by condensing Lagrange multipliers from the original F‐CMS method. The performance of the proposed method is evaluated through numerical examples, which illustrates improved accuracy over existing CMS methods.

Role of Intrinsic Disorder in Alternatively Spliced Sortilin Variants in 3T3L1 Adipocytes

Type 2 diabetes mellitus (T2DM) is a metabolic disorder defined by systemic insulin resistance affecting more than 400 million people worldwide. Adipose tissue is an important regulator of glucose metabolism as well as an endocrine organ secreting hormones and cytokines. Defective signaling in insulin resistant adipocytes leads to impaired trafficking of main glucose transporter 4 (Glut4) thus leading to reduced glucose uptake and hyperglycemia. Notably, trafficking protein sortilin has an important role in regulating Glut4 movement in adipocytes thereby affecting glucose levels. We demonstrate that in mouse 3T3L1 adipocytes, alternative splicing generates a sortilin splice variant with an included exon 17b (Sort17b) expressed along with SortWT. Insulin resistance induced by high serum increases expression levels of Sort17b. Further, our results demonstrate that increased Sort17b decreases Glut4 translocation and glucose uptake. Using bioinformatic analysis we show inclusion of exon 17b introduces a novel intrinsic disorder region in a crucial ligand binding regulation site. Since intrinsic disorder is a measure of lack of protein folding potentially affecting ligand binding and protein interactions, we performed molecular modeling of both sortilin variants (Sort17b and SortWT) and analyzed flexibility of protein backbone (RMSD) and individual residues (RMSF). Our results show increased flexibility in Sort17b within the amino acids of exon 17b. Using co-immunoprecipitation assays we demonstrate that Sort17b binds robustly to Glut4 and the binding efficiency is increased in insulin resistant adipocytes. Overall, the results demonstrate that the novel splice variant of sortilin in mouse 3T3L1 adipocytes contains an intrinsically disordered region that affects glucose metabolism and may be a potential therapeutic target in T2DM.

Visualization of Protein-Drug Interactions for the Analysis of Drug Resistance in Lung Cancer.

Non-small cell lung cancer (NSCLC) caused by mutation of the epidermal growth factor receptor (EGFR) is a major cause of death worldwide. Tyrosine kinase inhibitors (TKIs) of EGFR have been developed and show promising results at the initial stage of therapy. However, in most cases, their efficacy becomes limited due to the emergence of secondary mutations causing drug resistance after about a year. In this work, we investigated the mechanism of drug resistance due to these mutations. We performed molecular dynamics (MD) simulations of EGFR-drug interactions to obtain Euclidean distance and binding free energy values to analyse drug resistance and visualize drug-protein interactions. A PCA-based method is proposed to find normal, rigid, flexible, and critical residues. We have established a systematic method for the visualization of protein-drug interactions, which provides an effective framework for the analysis of drug resistance in lung cancer at the atomic level.

Molecular modeling and molecular dynamics simulation studies on thyroid hormone receptor from Rattus norvegicus: role of conserved water molecules.

Thyroid hormone receptor (THR) belongs to the nuclear receptor (NR) superfamily that is activated by binding of appropriate ligand molecules (thyroid hormones). These receptors directly bind to specific DNA sequences for gene expression, which is essential for metabolism, homeostasis, and the development of organisms, making it an important drug target. Extensive MD-simulation studies of triiodothyronine (T3) docked modeled rnTHRβ1 structures have indicated the presence of twelve conserved water molecules at the DNA-DBD (DNA binding domain) interface. The W1-W5 water centers have been involved in the recognition between the A-chain of DBD to C-chain of DNA, W6 and W7 mediated the interaction between A-chain of DBD and D-chain of DNA, W8 and W9 recognized the B-chain of DBD and C-chain of DNA, and W9-W12 centers conjugated the residues of B-chain of DBD to D-chain of DNA through hydrogen bonds. The conformation flexibility of Phe272 and Met313 residues in the absence of T3 at the LBD (ligand-binding domain) region have been observed and reported.

Flexible Auto-Weighted Local-Coordinate Concept Factorization: A Robust Framework for Unsupervised Clustering

Concept Factorization (CF) and its variants may produce inaccurate representation and clustering results due to the sensitivity to noise, hard constraint on the reconstruction error, and pre-obtained approximate similarities. To improve the representation ability, a novel unsupervised Robust Flexible Auto-weighted Local-coordinate Concept Factorization (RFA-LCF) framework is proposed for clustering high-dimensional data. Specifically, RFA-LCF integrates the robust flexible CF by clean data space recovery, robust sparse local-coordinate coding, and adaptive weighting into a unified model. RFA-LCF improves the representations by enhancing the robustness of CF to noise and errors, providing a flexible constraint on the reconstruction error and optimizing the locality jointly. For robust learning, RFA-LCF clearly learns a sparse projection to recover the underlying clean data space, and then the flexible CF is performed in the projected feature space. RFA-LCF also uses a L2,1-norm based flexible residue to encode the mismatch between the recovered data and its reconstruction, and uses the robust sparse local-coordinate coding to represent data using a few nearby basis concepts. For auto-weighting, RFA-LCF jointly preserves the manifold structures in the basis concept space and new coordinate space in an adaptive manner by minimizing the reconstruction errors on clean data, anchor points and coordinates. By updating the local-coordinate preserving data, basis concepts and new coordinates alternately, the representation abilities can be potentially improved. Extensive results on public databases show that RFA-LCF delivers the state-of-the-art clustering results compared with other related methods.

The Role of Rigid Residues in Modulating TEM-1 β-Lactamase Function and Thermostability.

The relationship between protein motions (i.e., dynamics) and enzymatic function has begun to be explored in β-lactamases as a way to advance our understanding of these proteins. In a recent study, we analyzed the dynamic profiles of TEM-1 (a ubiquitous class A β-lactamase) and several ancestrally reconstructed homologues. A chief finding of this work was that rigid residues that were allosterically coupled to the active site appeared to have profound effects on enzyme function, even when separated from the active site by many angstroms. In the present work, our aim was to further explore the implications of protein dynamics on β-lactamase function by altering the dynamic profile of TEM-1 using computational protein design methods. The Rosetta software suite was used to mutate amino acids surrounding either rigid residues that are highly coupled to the active site or to flexible residues with no apparent communication with the active site. Experimental characterization of ten designed proteins indicated that alteration of residues surrounding rigid, highly coupled residues, substantially affected both enzymatic activity and stability; in contrast, native-like activities and stabilities were maintained when flexible, uncoupled residues, were targeted. Our results provide additional insight into the structure-function relationship present in the TEM family of β-lactamases. Furthermore, the integration of computational protein design methods with analyses of protein dynamics represents a general approach that could be used to extend our understanding of the relationship between dynamics and function in other enzyme classes.

Deciphering the Mechanism of Inhibition of SERCA1a by Sarcolipin Using Molecular Simulations.

SERCA1a is an ATPase calcium pump that transports Ca2+ from the cytoplasm to the sarco/endoplasmic reticulum lumen. Sarcolipin (SLN), a transmembrane peptide, regulates the activity of SERCA1a by decreasing its Ca2+ transport rate, but its mechanism of action is still not well-understood. To decipher this mechanism, we have performed normal mode analysis in the all-atom model, with the SERCA1a-SLN complex, or the isolated SERCA1a, embedded in an explicit membrane. The comparison of the results allowed us to provide an explanation at the atomic level for the action of SLN that is in good agreement with experimental observations. In our analyses, the presence of SLN locally perturbs the TM6 transmembrane helix and as a consequence modifies the position of D800, one of the key metal-chelating residues. Additionally, it reduces the flexibility of the gating residues, V304, and E309 in TM4, at the entrance of the Ca2+ binding sites, which would decrease the affinity for Ca2+. Unexpectedly, SLN has also an effect on the ATP binding site more than 35 Å away, due to the straightening of TM5, a long helix considered as the spine of the protein. The straightening of TM5 modifies the structure of the P-N linker that sits above it, and which comprises the 351DKTG354 conserved motif, resulting in an increase of the distance between ATP and the phosphorylation site. As a consequence, the turn-over rate could be affected. All this gives SERCA1a the propensity to go toward a Ca2+ low-affinity E2-like state in the presence of SLN and toward a Ca2+ high-affinity E1-like state in the absence of SLN. In addition to a general mechanism of inhibition of SERCA1a regulatory peptides, this study also provides an insight into the conformational transition between the E2 and E1 states.

Molecular dynamics study of in silico mutations in the auto-inhibitory loop of human endothelial nitric oxide synthase FMN sub-domain.

Structural flexibility of the peptide linker connecting two domains is essential for the functioning of multi-domain complex. Nitric oxide synthase (NOS) isoforms contain the oxygenase and the reductase domains connected by calmodulin binding linker (CBL) region. Additionally, the endothelial NOS (eNOS) isoform contain an auto-inhibitory loop (AI) in the FMN reductase sub-domain which represses the inter-domain electron transfer process. Binding of Ca2+-Calmodulin complex on the CBL region relieves the AI loop repression and facilitates electron transfer from FMN in the reductase domain to the heme in the oxygenase domain. Few experimental studies have reported that in vitro mutation of Serine-615 (S615D) and Serine-633 (S633D) in the FMN reductase sub-domain to aspartic acid increased NO production and increased Ca2+ sensitivity. To understand the role of AI loop in eNOS repression and activation in serine mutants (S615D and S633D), we modelled the FMN reductase sub-domain of human eNOS protein with and without the CBL region. Molecular dynamics simulations performed indicated that the mutant protein AI loop structure was stabilized by salt bridge formed between D615 and R602. It was also found that mutation increased the flexibility of C-terminal residues of eNOS CBL region. The hinge-like movement of the AI loop allowed rotation of the FMN sub-domain clockwise which may favour electron-transfer in the mutant protein. This study provides insight on mutation (S615D and S633D) induced changes in AI loop and increased flexibility of CBL region which may lead to the protein activation and may also facilitate Calmodulin binding at physiological Ca2+ concentration. Graphical Abstract Mutation of amino acid residues contribute to structural changes at molecular level leading to alteration in protein dynamics and its function. Serine-615 and Serine-633 in the auto-inhibitory loop of human eNOS reductase model was mutated to aspartic acid in silico and molecular dynamics simulations of the protein showed that steric hindrance due to mutation altered the auto-inhibitory loop rearrangement and the FMN sub-domain movement favouring electron transfer.

Structure determination and analyses of the GAPDH from the parasite Schistosoma mansoni, the first one from a platyhelminth

The enzyme Glyceraldehyde-3-Phosphate Dehydrogenase from Schistosoma mansoni (SmGAPDH) is characterized as a therapeutical target for schistosomiasis. In this context, we report here the experimental structure, structural analyses and comparisons of SmGAPDH, the first one from a Platyhelminth. The enzyme was expressed, purified and assayed for crystallization, what allowed the obtainment of crystals of sufficient quality to collect X-ray diffraction data up to 2.51 Å resolution. SmGAPDH is the only GAPDH to present the sequence NNR (its residues 114–116) which leads to (especially R116) a hydrogen bond network that possibly reflects on the flexibility of residues to interact with the adenine part of NAD+, speculated to be important for differential drug design.

Oncogenic mutations on Rac1 affect global intrinsic dynamics underlying GTP and PAK1 binding

Rac1 is a small member of the Rho GTPase family. One of the most important downstream effectors of Rac1 is a serine/threonine kinase, p21-activated kinase 1 (PAK1). Mutational activation of PAK1 by Rac1 has oncogenic signaling effects. Here, although we focus on Rac1-PAK1 interaction by atomic-force-microscopy-based single-molecule force spectroscopy experiments, we explore the effect of active mutations on the intrinsic dynamics and binding interactions of Rac1 by Gaussian network model analysis and molecular dynamics simulations. We observe that Rac1 oncogenic mutations are at the hinges of three global modes of motion, suggesting the mechanical changes as potential markers of oncogenicity. Indeed, the dissociation of wild-type Rac1-PAK1 complex shows two distinct unbinding dynamic states that are reduced to one with constitutively active Q61L and oncogenic Y72C mutant Rac1, as revealed by single-molecule force spectroscopy experiments. Q61L and Y72C mutations change the mechanics of the Rac1-PAK1 complex by increasing the elasticity of the protein and slowing down the transition to the unbound state. On the other hand, Rac1’s intrinsic dynamics reveal more flexible GTP and PAK1-binding residues on switches I and II with Q61L, Y72C, oncogenic P29S and Q61R, and negative T17N mutations. The cooperativity in the fluctuations of GTP-binding sites around the p-loop and switch I decreases in all mutants, mostly in Q61L, whereas some PAK1-binding residues display enhanced coupling with GTP-binding sites in Q61L and Y72C and within each other in P29S. The predicted binding free energies of the modeled Rac1-PAK1 complexes show that the change in the dynamic behavior likely means a more favorable PAK1 interaction. Overall, these findings suggest that the active mutations affect intrinsic functional dynamic events and alter the mechanics underlying the binding of Rac1 to GTP and upstream and downstream partners including PAK1.

Computational driven molecular dynamics simulation of keratinocyte growth factor behavior at different pH conditions

Abstract Chronic wounds with a long healing period are considered a worldwide problem. As an important player in the wound healing process, keratinocyte growth factor (KGF) is one of the most effective methods of wound treatment. To prevent the wounds from becoming chronic, effective treatments should be applied at the first steps of wound generation. However, the low stability of KGF at acidic pH conditions at the first stages of the wound prevents its application in wound treatment. Lack of knowledge of the KGF behavior at different pHs has made it difficult to design the pH stable variants of this protein. In this study, we tried to simulate the KGF behavior at different pH conditions using molecular dynamics. To do this, the protonation and deprotonation states of KGF at different pH(s) were analyzed with PROPKA webserver. Then, the PDB files generated at different pHs were subjected to the molecular dynamic simulation at temperatures 300, 400, and 500 K using GROMACS software. The conformational stability, residue flexibility, protein compactness, and hydrogen-bond analysis were performed at the atomic level. The results showed that by decreasing the total charge of the protein from 10 at neutral pH to 7 at alkaline pH (pH: 9.0) at 500K, the KGF stability was significantly increased. Furthermore, analysis of the solvent-accessible surface area (SASA) showed that the KGF SASA is decreased with reduced protein charge. Analysis of the intra-protein and protein-solvent hydrogen bonds (H-bonds) showed that not only the number of H-bonds at alkaline pHs are higher compared to the neutral and acidic pHs, but also, they did not significantly decrease with increasing temperature at alkaline pH(s). Our results showed that the presence of more positively charged residues and their repulsion makes the KGF an unstable protein. Therefore, it seems that targeted mutations to reduce protein charge can lead to increased stability of KGF.

Probing the physico-chemical, antioxidant and anticancer influence of β-lactoglobulin on dietary flavonoid daidzein

The molecular nature, antioxidant, and anticancer actions of the dietary flavonoid daidzein binding to two types of β-lactoglobulin (βLG) are investigated in this work by using different kinds of spectroscopic techniques, electrochemical analysis, molecular docking, and molecular dynamics simulation. The results showed that hydrogen bonds and van der Waals interactions played a major role in the development of the stable complex. It may be deduced from molecular docking studies and competitive assays that daidzein binds to the outer surface of βLG with appropriate binding energies. The results of molecular dynamics showed the creation of a relatively stable complex in which daidzein binding reduced the flexibility of key residues in the binding site. In terms of the masking effect of βLG, the best antioxidant activity was found at low concentrations of the βLG-daidzein complex. After 48 h of exposure to a 50 μM mixture of βLG and daidzein, the complexes showed increased anticancer activity with maximal mortalities of 38.62 and 47.08% for HeLa and MCF-7 cells, respectively. The improvement of daidzein's solubility and anticancer activity as a result of its complexation with βLG, as well as its light and heat stability, could be a beneficial approach in the food and pharmaceutical industries.

Molecular Dynamics Simulation to Uncover the Mechanisms of Protein Instability During Freezing.

Freezing is a common process applied in the pharmaceutical industry to store and transport biotherapeutics. Herewith, multi-scale molecular dynamics simulations of Lactate dehydrogenase (LDH) protein in phosphate buffer with/without ice formation performed to uncover the still poorly understood mechanisms and molecular details of protein destabilization upon freezing. Both fast and slow ice growing conditions were simulated at 243K from one or two-side of the simulation box, respectively. The rate of ice formation at all-atom simulations was crucial to LDH stability, as faster freezing rates resulted in enhanced structural stability maintained by a higher number of intramolecular hydrogen bonds, less flexible protein's residues, lower solvent accessibility and greater structural compactness. Further, protein aggregation investigated by coarse-grained simulations was verified to be initiated by extended protein structures and retained by electrostatic interactions of the salt bridges between charged residues and hydrogen bonds between polar residues of the protein. Lastly, the study of free energy of dissociation through steered molecular dynamics simulation revealed LDH was destabilized by the solvation of the hydrophobic core and the loss of hydrophobic interactions. For the first time, experimentally validated molecular simulations revealed the detailed mechanisms of LDH destabilization upon ice formation and cryoconcentration of solutes.