Insight derived from molecular dynamics simulation into substrate-induced changes in protein motions of proteinase K

Computational Insights into SARS-CoV-2 Main Protease Mutations and Nirmatrelvir Efficacy: The Effects of P132H and P132H-A173V

Nirmatrelvir, a pivotal component of the oral antiviral Paxlovid for COVID-19, targets the SARS-CoV-2 main protease (Mpro) as a covalent inhibitor. Here, we employed combined computational methods to explore how the prevalent Omicron variant mutation P132H, alone and in combination with A173V (P132H-A173V), affects nirmatrelvir's efficacy. Our findings suggest that P132H enhances the noncovalent binding affinity of Mpro for nirmatrelvir, whereas P132H-A173V diminishes it. Although both mutants catalyze the rate-limiting step more efficiently than the wild-type (WT) Mpro, P132H slows the overall rate of covalent bond formation, whereas P132H-A173V accelerates it. Comprehensive analysis of noncovalent and covalent contributions to the overall binding free energy of the covalent complex suggests that P132H likely enhances Mpro sensitivity to nirmatrelvir, while P132H-A173V may confer resistance. Per-residue decompositions of the binding and activation free energies pinpoint key residues that significantly affect the binding affinity and reaction rates, revealing how the mutations modulate these effects. The mutation-induced conformational perturbations alter drug-protein local contact intensities and the electrostatic preorganization of the protein, affecting noncovalent binding affinity and the stability of key reaction states, respectively. Our findings inform the mechanisms of nirmatrelvir resistance and sensitivity, facilitating improved drug design and the detection of resistant strains.

A sub-100 nm thickness flat jet for extreme ultraviolet to soft X-ray absorption spectroscopy

Experimental characterization of the structural, electronic and dynamic properties of dilute systems in aqueous solvents, such as nanoparticles, molecules and proteins, are nowadays an open challenge. X-ray absorption spectroscopy (XAS) is probably one of the most established approaches to this aim as it is element-specific. However, typical dilute systems of interest are often composed of light elements that require extreme-ultraviolet to soft X-ray photons. In this spectral regime, water and other solvents are rather opaque, thus demanding radical reduction of the solvent volume and removal of the liquid to minimize background absorption. Here, we present an experimental endstation designed to operate a liquid flat jet of sub-micrometre thickness in a vacuum environment compatible with extreme ultraviolet/soft XAS measurements in transmission geometry. The apparatus developed can be easily connected to synchrotron and free-electron-laser user-facility beamlines dedicated to XAS experiments. The conditions for stable generation and control of the liquid flat jet are analyzed and discussed. Preliminary soft XAS measurements on some test solutions are shown.

Interaction between the antioxidant compound safranal and α-chymotrypsin in spectroscopic fields and molecular modeling approaches

Among various herbal plants, saffron has been the subject of study in various medical and food fields. Among the compounds of saffron, safranal is one of them. Safranal is a monoterpene aldehyde. The precursor of safranal is called picrocrocin, whose hydrolysis leads to the production of safranal. picrocrocin has two sugar components and aglycone. sugar component was separated during the drying process of saffron and safranal is produced. Saffron is the cause of the saffron aroma. Previous studies have shown that safranal offers many benefits such as antioxidants, blood pressure regulation and anti-tumor qualities. On the other hand, α-Chy is an enzyme secreted by the pancreas into the intestine and then acts as an efficient protease. In this study, various methods, such as molecular dynamics (MD) simulation and molecular binding, and different spectroscopic techniques, as well as protein stability techniques, were used to investigate the possible interactions between safranal and α-Chy. UV spectroscopic studies were showing that the existence of safranal decreased α-Chy absorption intensity. safranal caused the intrinsic fluorescence of α-Chy to be quenched too. According to the Stern-Volmer equation, the interaction between safranal and α-Chy was of the static type. In thermodynamic calculations, the interaction between safranal and α-Chy was stabilized by hydrophobic forces. And it was found that this interaction continued spontaneously. These results were, thus, consistent with the Docking data simulation (with the negative ΔG° number and positive changes in enthalpy and entropy). The thermal stability of α-Chy was also measured, showing that its melting point was shifted to a higher threshold as a result of the interaction. also, MD simulation indicated that α-Chy became more stable in the presence of safranal. In this paper, all the results of the laboratory techniques were confirmed by molecular dynamic simulations, so the correctness of the results was confirmed. From this research, we hope to carefully observe the possible changes in the behavior and structure of the enzyme in the presence of safranal.Communicated by Ramaswamy H. Sarma.

Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network

Psychrophilic proteases have attracted enormous attention in past decades, due to their high catalytic activity at low temperatures in a wide range of industrial processes, especially in the detergent and leather industries. Among them, H5 is an alkaline protease mutant, which featuring psychrophilic-like behavior, but the reasons that H5 with higher activity at low temperatures are still poorly understood. Herein, the molecular dynamics (MD) simulations combined with residue interaction network (RIN) were utilized to investigate the mechanisms of the cold-adaption of mutant H5. The results demonstrated that two loops involved in the substrate binding G100-S104 and S125-S129 in H5 had higher mobility, and the distance enlargement between the two loops modulated the substrate's accessibility compared with wild type counterpart. Besides, H5 enhanced conformational flexibility by weakening salt bridges and increasing interaction with the solvent. In particular, the absence of Lys251-Asp197-Arg247 salt bridge network may contribute to the structural mobility. Based on the free energy landscape and molecular mechanics Poisson-Boltzmann surface area of the wild type and H5, it was elucidated that H5 possessed a large population of interconvertible conformations, resulting in the weaker substrate binding free energy. The calculated RIN topology parameters such as the average degree, average cluster coefficient, and average path length further verified that the mutant H5 attenuated residue-to-residue interactions. The investigation of the mechanisms by which how the residue mutation affects the stability and activity of enzymes provides a theoretical basis for the development of cold-adapted protease.

A VEGFB-Based Peptidomimetic Inhibits VEGFR2-Mediated PI3K/Akt/mTOR and PLCγ/ERK Signaling and Elicits Apoptotic, Antiangiogenic, and Antitumor Activities

Vascular endothelial growth factor receptor 2 (VEGFR2) mediates VEGFA signaling mainly through the PI3K/AKT/mTOR and PLCγ/ERK1/2 pathways. Here we unveil a peptidomimetic (VGB3) based on the interaction between VEGFB and VEGFR1 that unexpectedly binds and neutralizes VEGFR2. Investigation of the cyclic and linear structures of VGB3 (named C-VGB3 and L-VGB3, respectively) using receptor binding and cell proliferation assays, molecular docking, and evaluation of antiangiogenic and antitumor activities in the 4T1 mouse mammary carcinoma tumor (MCT) model showed that loop formation is essential for peptide functionality. C-VGB3 inhibited proliferation and tubulogenesis of human umbilical vein endothelial cells (HUVECs), accounting for the abrogation of VEGFR2, p-VEGFR2 and, subsequently, PI3K/AKT/mTOR and PLCγ/ERK1/2 pathways. In 4T1 MCT cells, C-VGB3 inhibited cell proliferation, VEGFR2 expression and phosphorylation, the PI3K/AKT/mTOR pathway, FAK/Paxillin, and the epithelial-to-mesenchymal transition cascade. The apoptotic effects of C-VGB3 on HUVE and 4T1 MCT cells were inferred from annexin-PI and TUNEL staining and activation of P53, caspase-3, caspase-7, and PARP1, which mechanistically occurred through the intrinsic pathway mediated by Bcl2 family members, cytochrome c, Apaf-1 and caspase-9, and extrinsic pathway via death receptors and caspase-8. These data indicate that binding regions shared by VEGF family members may be important in developing novel pan-VEGFR inhibitors that are highly relevant in the pathogenesis of angiogenesis-related diseases.

Enhancement of thermal stability of proteinase K by biocompatible cholinium-based ionic liquids

Proteinase K (PK) is a proteolytic enzyme that has been widely used in nucleic acid purification, leather production, environmental protection, and other industrial applications. However, this biocatalyst cannot tolerate high temperatures which has severely restricted its wider application. As reported in previous studies, cholinium-based ionic liquids (ILs) have gained tremendous attention serving as a promising media to stabilize and preserve proteins, DNA, and other biomolecules due to their environmentally benign nature and biocompatibility. In this work, we chose 13 different kinds of cholinium-based ILs to examine their effects on the thermal stability and enzymatic activity of PK. We found that biocompatible cholinium-based ions with appropriately chosen anions can greatly improve the thermal stability of PK, whose melting temperature (T_m) is increased from ∼74.4 °C to 87.7 °C. However, the enzymatic activity is slightly reduced in the presence of ILs. Further comparison of our results with other literature findings suggests that kosmotropic anions of cholinium-based ILs are crucial to maintain the thermal stability of proteins. However, to achieve the best performance, the choice of IL anions is protein specific.

Structural localization of pathogenic mutations in the central nucleotide-binding domain (NBD) of nucleotide-binding oligomerization domain-2 (NOD2) protein and their inference in inflammatory disorders

The human NBD domain which is centrally located in the NOD2 protein displays an essential role in oligomerization and initiates the immune response via CARD-RIPK2 interaction. The mutations associated with the NBD domain have been largely implicated in inflammatory disorders such as Blau syndrome and sarcoidosis. This study aims to determine the structural and phenotypic effect of a lethal mutation that occurs in the NBD domain which has an axiomatic impact on protein dysfunction. Initially, the most deleterious missense mutations were screened through various in silico analysis. Out of 33 variants, I-Mutant 3.0, SIFT, PolyPhen 2, Align GVGD, PHD SNP and SNP&GO have statistically identified 5 variants (R42W, D90E, E91K, G189D & W198L) as less stable, deleterious and damaging. Our predicted models have paved the way to understand the various structural properties such as physiochemical, secondary structural arrangements and stabilizing residues in folding associated with the native and mutant NBD domain especially of the functionally important regions. From the aforementioned results, R42W and G189D were found to be the more predominant among the mutants. Precisely, through molecular simulation, we have strongly justified the significant conformational disruption of R42W and G189D through the stabilization factors, folding and essential dynamics. Conclusively, these regions (α₃41-44, α₁₃185-191 and β₆133-143β₇) seem to adopt such structures that are not conducive to wild-type-like functionality. Our prediction and validation of lethal mutations based on structural stability may be useful for conducting experimental studies in detail to uncover the protein deregulation leading to inflammatory disorders.

Recombinant Production and Characterization of an Extracellular Subtilisin-Like Serine Protease from Acinetobacter baumannii of Fermented Food Origin

Acinetobacter baumannii is a ubiquitous bacteria that is increasingly becoming a formidable nosocomial pathogen. Due to its clinical relevance, studies on the bacteria's secretory molecules especially extracellular proteases are of interest primarily in relation to the enzyme's role in virulence. Besides, favorable properties that extracellular proteases possess may be exploited for commercial use thus there is a need to investigate extracellular proteases from Acinetobacter baumannii to gain insights into their catalytic properties. In this study, an extracellular subtilisin-like serine protease from Acinetobacter baumannii designated as SPSFQ that was isolated from fermented food was recombinantly expressed and characterized. The mature catalytically active form of SPSFQ shared a high percentage sequence identity of 99% to extracellular proteases from clinical isolates of Acinetobacter baumannii and Klebsiella pneumoniae as well as a moderately high percentage identity to other bacterial proteases with known keratinolytic and collagenolytic activity. The homology model of mature SPSFQ revealed its structure is composed of 10 β-strands, 8 α-helices, and connecting loops resembling a typical architecture of subtilisin-like α/β motif. SPSFQ is catalytically active at an optimum temperature of 40 °C and pH 9. Its activity is stimulated in the presence of Ca²⁺ and severely inhibited in the presence of PMSF. SPSFQ also displayed the ability to degrade several tissue-associated protein substrates such as keratin, collagen, and fibrin. Accordingly, our study shed light on the catalytic properties of a previously uncharacterized extracellular serine protease from Acinetobacter baumannii that warrants further investigations into its potential role as a virulence factor in pathogenicity and commercial applications.

Exploring the Cold-Adaptation Mechanism of Serine Hydroxymethyltransferase by Comparative Molecular Dynamics Simulations

Cold-adapted enzymes feature a lower thermostability and higher catalytic activity compared to their warm-active homologues, which are considered as a consequence of increased flexibility of their molecular structures. The complexity of the (thermo)stability-flexibility-activity relationship makes it difficult to define the strategies and formulate a general theory for enzyme cold adaptation. Here, the psychrophilic serine hydroxymethyltransferase (pSHMT) from Psychromonas ingrahamii and its mesophilic counterpart, mSHMT from Escherichia coli, were subjected to μs-scale multiple-replica molecular dynamics (MD) simulations to explore the cold-adaptation mechanism of the dimeric SHMT. The comparative analyses of MD trajectories reveal that pSHMT exhibits larger structural fluctuations and inter-monomer positional movements, a higher global flexibility, and considerably enhanced local flexibility involving the surface loops and active sites. The largest-amplitude motion mode of pSHMT describes the trends of inter-monomer dissociation and enlargement of the active-site cavity, whereas that of mSHMT characterizes the opposite trends. Based on the comparison of the calculated structural parameters and constructed free energy landscapes (FELs) between the two enzymes, we discuss in-depth the physicochemical principles underlying the stability-flexibility-activity relationships and conclude that (i) pSHMT adopts the global-flexibility mechanism to adapt to the cold environment and, (ii) optimizing the protein-solvent interactions and loosening the inter-monomer association are the main strategies for pSHMT to enhance its flexibility.

Structural Studies on an Anti-Angiogenic Peptide Using Molecular Modeling

BACKGROUND

Development of VEGF antagonists, which inhibit its interaction with the receptors, is a widely used strategy for the inhibition of angiogenesis and tumor growth.

OBJECTIVES

In the present study, a VEGFR-1 antagonistic peptide was designed and its potential for binding to VEGFR-1 and VEGFR-2 was evaluated by theoretical studies.

MATERIALS AND METHODS

Based on the X-ray structure of VEGF-B/VEGFR-1 D2 (PDB ID: 2XAC), an antagonistic peptide (known as VGB1) was designed, and its model structure was constructed using homology modeling in the MODELLER, version 9.16. The validity of the modeled structures was estimated employing several web tools. Finally, one model was chosen and molecular dynamics (MD) simulation was applied using the GROMACS package, version 5.1.4, to allow conformational relaxation of the structure. Next, docking process of the peptide with VEGFR-1 and VEGFR-2 was performed by HADDOCK web server and the docking structures were optimized by MD simulation for 20 ns. The far-UV circular dichroism (CD) spectrum of VGB1 was recorded to evaluate the overall structure of the peptide.

RESULTS

The far-UV CD spectrum indicated that VGB1 contains α helix structure. The results from docking studies suggested that Van der Waals and nonpolar interactions play the most important role in the peptide binding to VEGFR-1. In addition, our results implicated the relevance of both Van der Waals and electrostatic interactions in the formation of complex between VGB1 and VEGFR-2. In addition to the common binding residues in the corresponding region of VEGF-A and VEGF-B, additional binding residues also were predicted for the interaction of VGB1 with VEGFR-1 and VEGFR-2.

CONCLUSIONS

The results of MD and molecular docking simulations predicted that VGB1 recognizes both VEGFR-1 and VEGFR-2, which may lead to the prevention of the downstream signaling triggered by these receptors.

Molecular dynamics simulations reveal distinct differences in conformational dynamics and thermodynamics between the unliganded and CD4-bound states of HIV-1 gp120

The entry of human immunodeficiency virus type I (HIV-1) into host cells is initiated by binding to the cell-surface receptor CD4, which induces a conformational transition of the envelope (Env) glycoprotein gp120 from the closed, unliganded state to the open, CD4-bound state. Despite many available structures in these two states, detailed aspects on the dynamics and thermodynamics of gp120 remain elusive. Here, we performed microsecond-scale (μs-scale) multiple-replica molecular dynamics (MD) simulations to explore the differences in the conformational dynamics, protein motions, and thermodynamics between the unliganded and CD4-bound/complexed forms of gp120. Comparative analyses of MD trajectories reveal that CD4 binding promotes the structural deviations/changes and conformational flexibility, loosens the structural packing, and complicates the molecular motions of gp120. Comparison of the constructed free energy landscapes (FELs) reveals that the CD4-complexed gp120 has more conformational substates, larger conformational entropy, and lower thermostability than the unliganded form. Therefore, the unliganded conformation represents a structurally and energetically stable "ground state" for the full-length gp120. The observed great increase in the mobility of V1/V2 and V3 along with their more versatile movement directions in the CD4-bound gp120 compared to the unliganded form suggests that their orientations with respect to each other and to the structural core determine the differences in the conformational dynamics and thermodynamics between the two gp120 forms. The results presented here provide a basis by which to better understand the functional and immunological properties of gp120 and, furthermore, to deploy appropriate strategies for the development of anti-HIV-1 drugs or vaccines.

A cyclic peptide reproducing the α1 helix of VEGF-B binds to VEGFR-1 and VEGFR-2 and inhibits angiogenesis and tumor growth

Vascular endothelial growth factors (VEGFs) and their receptors (VEGFRs) are pivotal regulators of angiogenesis. The VEGF–VEGFR system is therefore an important target of anti-angiogenesis therapy. Based on the X-ray structure of VEGF-B/VEGFR-1 D2, we designed a cyclic peptide (known as VGB1) reproducing the α1 helix and its adjacent region to interfere with signaling through VEGFR-1. Unexpectedly, VGB1 bound VEGFR-2 in addition to VEGFR-1, leading to inhibition of VEGF-stimulated proliferation of human umbilical vein endothelial cells and 4T1 murine mammary carcinoma cells, which express VGEFR-1 and VEGFR-2, and U87 glioblastoma cells that mostly express VEGFR-2. VGB1 inhibited different aspects of angiogenesis, including proliferation, migration and tube formation of endothelial cells stimulated by VEGF-A through suppression of extracellular signal-regulated kinase 1/2 and AKT (Protein Kinase B) phosphorylation. In a murine 4T1 mammary carcinoma model, VGB1 caused regression of tumors without causing weight loss in association with impaired cell proliferation (decreased Ki67 expression) and angiogenesis (decreased CD31 and CD34 expression), and apoptosis induction (increased TUNEL staining and p53 expression, and decreased Bcl-2 expression). According to far-UV circular dichroism (CD) and molecular dynamic simulation data, VGB1 can adopt a helical structure. These results, for the first time, demonstrate that α1 helix region of VEGF-B recognizes both VEGFR-1 and VEGFR-2.

Catalytic activity, structure and stability of proteinase K in the presence of biosynthesized CuO nanoparticles

Here, CuO nanoparticles were synthesized using Sambucus nigra (elderberry) fruit extract. Further, the binding of proteinase K, as a model enzyme with green synthesized nanoparticles was investigated. The results demonstrated that the structural changes in enzyme were induced by the binding of nanoparticles. These changes were accompanied by the decrease in the Michaelis-Menten constant at 298 K. This means that the enzyme affinity for the substrate was increased. Thermodynamic parameters of protein stability and protein-ligand binding were estimated from the spectroscopic measurements at 298-333 K. Depending on the temperature, CuO nanoparticles showed a dual effect on the thermodynamic stability and binding affinity of enzyme. Nanoparticles increase the stability of the native state of enzyme at room temperature. On the other hand, nanoparticles stabilize the unfolded state of enzyme at 310-333 K. An overall favorable Gibbs energy change was observed for the binding process at 298-333 K. The enzyme-nanoparticle binding is enthalpically driven at room temperature. It was concluded that hydrogen bonding plays a key role in the interaction of enzyme with nanoparticles at 298-310 K. At higher temperatures, the protein-ligand binding is entropically driven. This means that hydrophobic association plays a major role in the proteinase K-CuO binding at 310-333 K.

Theoretical Study on Zearalenol Compounds Binding with Wild Type Zearalenone Hydrolase and V153H Mutant

Zearalenone hydrolase (ZHD) is the only reported α/β-hydrolase that can detoxify zearalenone (ZEN). ZHD has demonstrated its potential as a treatment for ZEN contamination that will not result in damage to cereal crops. Recent researches have shown that the V153H mutant ZHD increased the specific activity against α-ZOL, but decreased its specific activity to β-ZOL. To understand whyV153H mutation showed catalytic specificity for α-ZOL, four molecular dynamics simulations combining with protein network analysis for wild type ZHD α-ZOL, ZHD β-ZOL, V153H α-ZOL, and V153H β-ZOL complexes were performed using Gromacs software. Our theoretical results indicated that the V153H mutant could cause a conformational switch at the cap domain (residues Gly161–Thr190) to affect the relative position catalytic residue (H242). Protein network analysis illustrated that the V153H mutation enhanced the communication with the whole protein and residues with high betweenness in the four complexes, which were primarily assembled in the cap domain and residues Met241 to Tyr245 regions. In addition, the existence of α-ZOL binding to V153H mutation enlarged the distance from the O_AE atom in α-ZOL to the NE2 atom in His242, which prompted the side chain of H242 to the position with catalytic activity, thereby increasing the activity of V153H on the α-ZOL. Furthermore, α-ZOL could easily form a right attack angle and attack distance in the ZHD and α-ZOL complex to guarantee catalytic reaction. The alanine scanning results indicated that modifications of the residues in the cap domain produced significant changes in the binding affinity for α-ZOL and β-ZOL. Our results may provide useful theoretical evidence for the mechanism underlying the catalytic specificity of ZHD.

Structural insights on bioremediation of polycyclic aromatic hydrocarbons using microalgae: a modelling-based computational study

Research on bioremediation of polycyclic aromatic hydrocarbons (PAHs) has established that several remediating microbial species are capable of degrading only low molecular weight (LMW)-PAHs, whereas high molecular weight (HMW)-PAHs are hardly degradable. In the present study, the efficiency of degradation of both LMW and HMW-PAHs by cytochrome P450 monooxygenase (CYP) of microalgae was studied. CYP have a key role in the detoxification of xenobiotics. So far, the structure of CYP in microalgae is not predicted; the protein structure was constructed by molecular modelling in the current study using the available template of microalgal CYP. Modelled microalgae 3D structures were docked against 38 different PAH compounds, and the information regarding the interaction between protein and PAHs viz. binding sites along with mode of interactions was investigated. We report that CYP from the microalgae Haematococcus pluvialis and Parachlorella kessleri was found to possess broad oxidising capability towards both LMW and HMW-PAHs. P. kessleri showed a least value with extra precision glide score of - 10.23 and glide energy of - 23.48 kcal/mol. PAHs bind to CYP active sites at Lys ⁶⁹, Trp ⁹⁶, Gln ³⁹⁷ and Arg ³⁹⁸ through intermolecular hydrogen bonding. Also, study revealed that PAHs interacted with CYP active sites through intermolecular hydrogen bonding, hydrophobic bonding, π-π interactions and van der waals interactions. Thus, structural elucidation study confirms that microalgae Parachlorella kessleri have the capacity to remediate HMW more efficiently than other microorganisms. Our results provide a framework in understanding the structure and the possible binding sites of CYP protein for degradation of PAH and that could be a screening tool in identifying the phycoremediation potential.

Insights into the role of electrostatics in temperature adaptation: a comparative study of psychrophilic, mesophilic, and thermophilic subtilisin-like serine proteases

To investigate the role of electrostatics in different temperature adaptations, we performed a comparative study on subtilisin-like serine proteases from psychrophilic Vibrio sp. PA-44 (VPR), mesophilic Engyodontium album (Tritirachium album) (PRK), and thermophilic Thermus aquaticus (AQN) using multiple-replica molecular dynamics (MD) simulations combined with continuum electrostatics calculations. The results reveal that although salt bridges are not a crucial factor in determining the overall thermostability of these three proteases, they on average provide the greatest, moderate, and least electrostatic stabilization to AQN, PRK, and VPR, respectively, at the respective organism growth temperatures. Most salt bridges in AQN are effectively stabilizing and thus contribute to maintaining the overall structural stability at 343 K, while nearly half of the salt bridges in VPR interconvert between being stabilizing and being destabilizing, likely aiding in enhancing the local conformational flexibility at 283 K. The individual salt bridges, salt-bridge networks, and calcium ions contribute differentially to local stability and flexibility of these three enzyme structures, depending on their spatial distributions and electrostatic strengths. The shared negatively charged surface potential at the active center of the three enzymes may provide the active-center flexibility necessary for nucleophilic attack and proton transfer. The differences in distributions of the electro-negative, electro-positive, and electro-neutral potentials, particularly over the back surfaces of the three proteases, may modulate/affect not only protein solubility and thermostability but also structural stability and flexibility/rigidity. These results demonstrate that electrostatics contributes to both heat and cold adaptation of subtilisin-like serine proteases through fine-tuning, either globally or locally, the structural stability and conformational flexibility/rigidity, thus providing a foundation for further engineering and mutagenesis studies.

Differently charged surface patches contribute to temperature adaptation of subtilisin-like serine proteases through affecting/modulating the protein solubility and thermostability and the structural flexibility/rigidity/stability.

Depuration of Cadmium from Blue Mussel (Mytilus edulis) by Protein Hydrolysate-Fe2+ Complex: The Role of Metallothionein

Metallothionein (MT) plays a major role in the detoxification of toxic metal ions in mussel. This study served to investigate the effects of prepared protein hydrolysate-Fe²⁺ (PH-Fe²⁺ ) on depuration of cadmium salt (Cd²⁺ ) from blue mussel (Mytilus edulis). The results indicated that Cd²⁺ concentrations in control ranged from 46.1 to 47.1 μg/g during 15 d of depuration. While, 40 mg/L PH-Fe²⁺ feed mussels exhibited obvious changes in Cd²⁺ concentration variables, which decreased by 22.8 μg/g after 15 d of depuration, making them significantly lower than the values of the control. Our assumption was that Cd²⁺ could be dissociated effectively from the complex of MT-Cd²⁺ in mussels affected by the incorporation of PH-Fe²⁺ during the feeding period. Further results of homology modeling and molecular dynamics (MD) confirmed that the combined power between MT and Cd²⁺ weakened significantly by PH-Fe²⁺ . This condition affected the charge density and/or the loop flexibility of MT and decreased the interaction energy within MT-Cd²⁺ complex and resulted in the release of Cd²⁺ from the complex, thereby exhibiting excretion detoxification. Finally, by comparing the experimental results to MD results, significant positive correlations were observed between PH-Fe²⁺ and the depuration of Cd²⁺ from MT-Cd²⁺ complex. Overall, the findings of this study may help better understand the depuration mechanisms of Cd²⁺ combined with MT, and the PH-Fe²⁺ can be recommended as a depuration agent to decrease Cd²⁺ concentrations in mussels.

PRACTICAL APPLICATION

Metallothionein (MT) is a low-molecular-weight protein with high metal-ion affinity. If the intracellular concentrations of metals are too high or if toxic metals are present within the cell, then the synthesis of MTs is induced and generated. In our previous work, it was found that the prepared hydrolysate-metal element complex showed obvious depuration activity of heavy metals (Cd²⁺ ) from blue mussel (Mytilus edulis). This study provided further the depuration mechanisms of Cd²⁺ from mussel (M. edulis), in particular to the role of MT and its chelate during the depuration process.

Evaluation of peptide designing strategy against subunit reassociation in mucin 1: A steered molecular dynamics approach

Subunit reassociation in mucin 1, a breast cancer tumor marker, is reported as one of the critical factors for its cytoplasmic activation. Inhibition of its heterodimeric association would therefore result in loss of its function and alter disease progression. The present study aimed at evaluating peptide inhibitor designing strategies that may serve as antagonist against this receptor-ligand alliance. Several peptides and their derivatives were designed based on native residues, subunit interface, hydrogen bonding and secondary structure. Docking studies with the peptides were carried on the receptor subunit and their binding affinities were evaluated using steered molecular dynamics simulation and umbrella sampling. Our results showed that among all the different classes of peptides evaluated, the receptor based peptide showed the highest binding affinity. This result was concurrent with the experimental observation that the receptor-ligand alliance in mucin 1 is highly specific. Our results also show that peptide ligand against this subunit association is only stabilized through native residue inter-protein interaction irrespective of the peptide structure, peptide length and number of hydrogen bonds. Consistency in binding affinity, pull force and free energy barrier was observed with only the receptor derived peptides which resulted in favorable interprotein interactions at the interface. Several observations were made and discussed which will eventually lead to designing efficient peptide inhibitors against mucin 1 heterodimeric subunit reassociation.

Insights into Cryoprotective Roles of Carrageenan Oligosaccharides in Peeled Whiteleg Shrimp (Litopenaeus vannamei) during Frozen Storage

The cryoprotective effects of carrageenan oligosaccharides on peeled whiteleg shrimp were investigated and compared with sodium pyrophosphate treatment during frozen storage, primarily the interaction mechanisms between oligosaccharides and shrimp myosin. Data revealed significant profitable effects on water-holding capacity and textural variables in oligosaccharide-treated shrimp compared to the control. Chemical analyses showed that these saccharides maintained a higher myofibrillar protein content and Ca²⁺-ATPase activity in frozen shrimp. Additionally, the hematoxylin and eosin staining results indicated that the saccharides significantly slowed the damage to muscle tissue structures. The assumption was that water replacement hypothesis played a leading role in cryoprotection of frozen shrimp. Furthermore, the homology modeling and molecular dynamics simulations confirmed that the saccharides substituted water molecules around the shrimp myosin surface by forming hydrogen bonds with polar residues of amino acids, thereby stabilizing the structures in the absence of water, leading to an increase in protein stability during frozen storage.

Substrate-induced changes in dynamics and molecular motions of cuticle-degrading serine protease PL646: a molecular dynamics study

Cuticle-degrading serine proteases secreted by nematophagous fungi can degrade the nematode cuticle during the infection processes.

Molecular Dynamics and Free Energy Simulations of Phenylacetate and CO2 Release from AMDase and Its G74C/C188S Mutant: A Possible Rationale for the Reduced Activity of the Latter

Arylmalonate decarboxylase (AMDase) catalyzes the decarboxylation of α-aryl-α-methyl malonates to produce optically pure α-arylpropionates of industrial and medicinal importance. Herein, atomistic molecular dynamics simulations have been carried out to delineate the mechanism of the release of product molecules phenylacetate (PAC) and carbon dioxide (CO₂), from the wild-type (WT) and its G74C/C188S mutant enzymes. Both of the product molecules follow a crystallographically characterized solvent-accessible channel to come out of the protein interior. A higher free energy barrier for the release of PAC from G74C/C188S compared to that in the WT is consistent with the experimentally observed compromised efficiency of the mutant. The release of CO₂ precedes that of PAC; free energy barriers for CO₂ and PAC release in the WT enzyme are calculated to be ∼1-2 and ∼23 kcal/mol, respectively. Postdecarboxylation, CO₂ moves toward a hydrophobic pocket formed by Pro 14, Leu 38, Leu 40, Leu 77, and the side chain of Tyr 48 which serves as its temporary "reservoir". CO₂ releases following a channel mainly decorated by apolar residues, unlike in the case of oxalate decarboxylase where polar residues mediate its transport.

Casting the critical regions in nucleotide binding oligomerization domain 2 protein: a signature mediated structural dynamics approach

Nucleotide binding oligomerization domain 2 (NOD2), a protein involved in the first line defence mechanism has a pivotal role in innate immunity. Impaired function of this protein is implicated in disorders such as Blau syndrome and Crohn's disease. Since an altered function is linked to protein's structure, we framed a systematic strategy to interpret the structure-function relationship of the protein. Initiated with mutation-based pattern prediction and identified a distant ortholog (DO) of NOD2 from which the intra-residue interaction network was elucidated. The network was used to identify hotspots that serve as critical points to maintain the stable architecture of the protein. Structural comparison of NOD2 domains with a DO revealed the minimal number of intra-protein interactions required by the protein to maintain the structural fold. In addition, the conventional molecular dynamics simulation emphasized the conformational transitions at hot spot residues between native NOD2 domains and its respective mutants (G116R, R42W and R54A) structures. The analysis of intra-protein interactions globally and the displacement of residues locally around the mutational site revealed loss of several critical bonds and residues vital for the protein's function. Conclusively we report, about 10 residues in leucine-rich repeat, 13 residues in NOD and 6 residues in CARD domain are required by the NOD2 to maintain its function. This protocol will help the researchers to achieve for more prospective studies to attest druggable site utility in discovering novel drug candidates.

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.

Influence of the R823W mutation on the interaction of the ANKS6-ANKS3: insights from molecular dynamics simulation and free energy analysis

The sterile alpha motif (SAM) domain of the protein ANKS6, a protein-protein interaction domain, is responsible for autosomal dominant polycystic kidney disease. Although the disease is the result of the R823W point mutation in the SAM domain of the protein ANKS6, the molecular details are still unclear. We applied molecular dynamics simulations, the principal component analysis, and the molecular mechanics Poisson-Boltzmann surface area binding free energy calculation to explore the structural and dynamic effects of the R823W point mutation on the complex ANKS6-ANKS3 (PDB ID: 4NL9) in comparison to the wild proteins. The energetic analysis presents that the wild type has a more stable structure than the mutant. The R823W point mutation not only disrupts the structure of the ANKS6 SAM domain but also negatively affects the interaction of the ANKS6-ANKS3. These results further clarify the previous experiments to understand the ANKS6-ANKS3 interaction comprehensively. In summary, this study would provide useful suggestions to understand the interaction of these proteins and their fatal action on mediating kidney function.

Effect of the Solvent Temperatures on Dynamics of Serine Protease Proteinase K

To obtain detailed information about the effect of the solvent temperatures on protein dynamics, multiple long molecular dynamics (MD) simulations of serine protease proteinase K with the solute and solvent coupled to different temperatures (either 300 or 180 K) have been performed. Comparative analyses demonstrate that the internal flexibility and mobility of proteinase K are strongly dependent on the solvent temperatures but weakly on the protein temperatures. The constructed free energy landscapes (FELs) at the high solvent temperatures exhibit a more rugged surface, broader spanning range, and higher minimum free energy level than do those at the low solvent temperatures. Comparison between the dynamic hydrogen bond (HB) numbers reveals that the high solvent temperatures intensify the competitive HB interactions between water molecules and protein surface atoms, and this in turn exacerbates the competitive HB interactions between protein internal atoms, thus enhancing the conformational flexibility and facilitating the collective motions of the protein. A refined FEL model was proposed to explain the role of the solvent mobility in facilitating the cascade amplification of microscopic motions of atoms and atomic groups into the global collective motions of the protein.

The effect of spermidine on the structure, kinetics and stability of proteinase K: spectroscopic and computational approaches

Polyamines (such as spermidine) are low molecular weight compounds which can be used as cosolvents in biological and industrial applications.

Elucidation, functional clustering and structural characterization of βTrCP1 substrates through a molecular dynamics study

Structural knowledge of substrate recognition by SKP1–CUL1–βTrCP1 complex for targeted cancer therapeutic strategy.

Engineering trypsin for inhibitor resistance

The development of effective protease therapeutics requires that the proteases be more resistant to naturally occurring inhibitors while maintaining catalytic activity. A key step in developing inhibitor resistance is the identification of key residues in protease-inhibitor interaction. Given that majority of the protease therapeutics currently in use are trypsin-fold, trypsin itself serves as an ideal model for studying protease-inhibitor interaction. To test the importance of several trypsin-inhibitor interactions on the prime-side binding interface, we created four trypsin single variants Y39A, Y39F, K60A, and K60V and report biochemical sensitivity against bovine pancreatic trypsin inhibitor (BPTI) and M84R ecotin. All variants retained catalytic activity against small, commercially available peptide substrates [kcat /KM  = (1.2 ± 0.3) × 10(7) M(-1 ) s(-1) . Compared with wild-type, the K60A and K60V variants showed increased sensitivity to BPTI but less sensitivity to ecotin. The Y39A variant was less sensitive to BPTI and ecotin while the Y39F variant was more sensitive to both. The relative binding free energies between BPTI complexes with WT, Y39F, and Y39A were calculated based on 3.5 µs combined explicit solvent molecular dynamics simulations. The BPTI:Y39F complex resulted in the lowest binding energy, while BPTI:Y39A resulted in the highest. Simulations of Y39F revealed increased conformational rearrangement of F39, which allowed formation of a new hydrogen bond between BPTI R17 and H40 of the variant. All together, these data suggest that positions 39 and 60 are key for inhibitor binding to trypsin, and likely more trypsin-fold proteases.

Investigate the binding of catechins to trypsin using docking and molecular dynamics simulation

To explore the inhibitory mechanism of catechins for digestive enzymes, we investigated the binding mode of catechins to a typical digestive enzyme-trypsin and analyzed the structure-activity relationship of catechins, using an integration of molecular docking, molecular dynamics simulation and binding free energy calculation. We found that catechins with different structures bound to a conservative pocket S1 of trypsin, which is comprised of residues 189–195, 214–220 and 225–228. In the trypsin-catechin complexes, Asp189 by forming strong hydrogen bonding, and Gln192, Trp215 and Gly216 through hydrophobic interactions, all significantly contribute to the binding of catechins. The number and the position of hydroxyl and aromatic groups, the structure of stereoisomers, and the orientation of catechins in the binding pocket S1 of trypsin all affect the binding affinity. The binding affinity is in the order of Epigallocatechin gallate (EGCG) > Epicatechin gallate (ECG) > Epicatechin (EC) > Epigallocatechin (EGC), and 2R-3R EGCG shows the strongest binding affinity out of other stereoisomers. Meanwhile, the synergic conformational changes of residues and catechins were also analyzed. These findings will be helpful in understanding the knowledge of interactions between catechins and trypsin and referable for the design of novel polyphenol based functional food and nutriceutical formulas.

Physicochemical bases for protein folding, dynamics, and protein-ligand binding

Proteins are essential parts of living organisms and participate in virtually every process within cells. As the genomic sequences for increasing number of organisms are completed, research into how proteins can perform such a variety of functions has become much more intensive because the value of the genomic sequences relies on the accuracy of understanding the encoded gene products. Although the static three-dimensional structures of many proteins are known, the functions of proteins are ultimately governed by their dynamic characteristics, including the folding process, conformational fluctuations, molecular motions, and protein-ligand interactions. In this review, the physicochemical principles underlying these dynamic processes are discussed in depth based on the free energy landscape (FEL) theory. Questions of why and how proteins fold into their native conformational states, why proteins are inherently dynamic, and how their dynamic personalities govern protein functions are answered. This paper will contribute to the understanding of structure-function relationship of proteins in the post-genome era of life science research.

Construction of an intermediate-resolution lattice model and re-examination of the helix-coil transition: a dynamic Monte Carlo simulation

In protein modeling, spatial resolution and computational efficiency are always incompatible. As a compromise, an intermediate-resolution lattice model has been constructed in the present work. Each residue is decomposed into four basic units, i.e. the α-carbon group, the carboxyl group, the imino group, and the side-chain group, and each basic coarse-grained unit is represented by a minimum cubic box with eight lattice sites. The spacing of the lattice is about 0.56 Å, holding the highest spatial resolution for the present lattice protein models. As the first report of this new model, the helix-coil transition of a polyalanine chain was examined via dynamic Monte Carlo simulation. The period of formed α-helix was about 3.68 residues, close to that of a natural α-helix. The resultant backbone motion was found to be in the realistic regions of the conformational space in the Ramachandran plot. Helix propagation constant and nucleation constant were further determined through the dynamic hydrogen bonding process and torsional angle variation, and the results were used to make comparison between classical Zimm-Bragg theory and Lifson-Roig theory based on the Qian-Schellman relationship. The simulation results confirmed that our lattice model can reproduce the helix-coil transition of polypeptide and construct a moderately fine α-helix conformation without significantly weakening the priority in efficiency for a lattice model.

Effects of Temperature Control Algorithms on Transport Properties and Kinetics in Molecular Dynamics Simulations

Temperature control algorithms in molecular dynamics (MD) simulations are necessary to study isothermal systems. However, these thermostatting algorithms alter the velocities of the particles and thus modify the dynamics of the system with respect to the microcanonical ensemble, which could potentially lead to thermostat-dependent dynamical artifacts. In this study, we investigate how six well-established thermostat algorithms applied with different coupling strengths and to different degrees of freedom affect the dynamics of various molecular systems. We consider dynamic processes occurring on different times scales by measuring translational and rotational self-diffusion as well as the shear viscosity of water, diffusion of a small molecule solvated in water, and diffusion and the dynamic structure factor of a polymer chain in water. All of these properties are significantly dampened by thermostat algorithms which randomize particle velocities, such as the Andersen thermostat and Langevin dynamics, when strong coupling is used. For the solvated small molecule and polymer, these dampening effects are reduced somewhat if the thermostats are applied to the solvent alone, such that the solute's temperature is maintained only through thermal contact with solvent particles. Algorithms which operate by scaling the velocities, such as the Berendsen thermostat, the stochastic velocity rescaling approach of Bussi and co-workers, and the Nosé-Hoover thermostat, yield transport properties that are statistically indistinguishable from those of the microcanonical ensemble, provided they are applied globally, i.e. coupled to the system's kinetic energy. When coupled to local kinetic energies, a velocity scaling thermostat can have dampening effects comparable to a velocity randomizing method, as we observe when a massive Nose-Hoover coupling scheme is used to simulate water. Correct dynamical properties, at least those studied in this paper, are obtained with the Berendsen thermostat applied globally, despite the fact that it yields the wrong kinetic energy distribution.

Protein dynamics and motions in relation to their functions: several case studies and the underlying mechanisms

Proteins are dynamic entities in cellular solution with functions governed essentially by their dynamic personalities. We review several dynamics studies on serine protease proteinase K and HIV-1 gp120 envelope glycoprotein to demonstrate the importance of investigating the dynamic behaviors and molecular motions for a complete understanding of their structure–function relationships. Using computer simulations and essential dynamic (ED) analysis approaches, the dynamics data obtained revealed that: (i) proteinase K has highly flexible substrate-binding site, thus supporting the induced-fit or conformational selection mechanism of substrate binding; (ii) Ca²⁺ removal from proteinase K increases the global conformational flexibility, decreases the local flexibility of substrate-binding region, and does not influence the thermal motion of catalytic triad, thus explaining the experimentally determined decreased thermal stability, reduced substrate affinity, and almost unchanged catalytic activity upon Ca²⁺ removal; (iii) substrate binding affects the large concerted motions of proteinase K, and the resulting dynamic pocket can be connected to substrate binding, orientation, and product release; (iv) amino acid mutations 375 S/W and 423 I/P of HIV-1 gp120 have distinct effects on molecular motions of gp120, facilitating 375 S/W mutant to assume the CD4-bound conformation, while 423 I/P mutant to prefer for CD4-unliganded state. The mechanisms underlying protein dynamics and protein–ligand binding, including the concept of the free energy landscape (FEL) of the protein–solvent system, how the ruggedness and variability of FEL determine protein's dynamics, and how the three ligand-binding models, the lock-and-key, induced-fit, and conformational selection are rationalized based on the FEL theory are discussed in depth.

ITP Adjuster 1.0: A New Utility Program to Adjust Charges in the Topology Files Generated by the PRODRG Server

The suitable computation of accurate atomic charges for the GROMACS topology *.itp files of small molecules, generated in the PRODRG server, has been a tricky task nowadays because it does not calculate atomic charges using an ab initio method. Usually additional steps of structure optimization and charges calculation, followed by a tedious manual replacement of atomic charges in the *.itp file, are needed. In order to assist this task, we report here the ITP Adjuster 1.0, a utility program developed to perform the replacement of the PRODRG charges in the *.itp files of small molecules by ab initio charges.