Biophysical Insight into the SARS-CoV2 Spike-ACE2 Interaction and Its Modulation by Hepcidin through a Multifaceted Computational Approach.

At the center of the SARS-CoV2 infection, the spike protein and its interaction with the human receptor ACE2 play a central role in the molecular machinery of SARS-CoV2 infection of human cells. Vaccine therapies are a valuable barrier to the worst effects of the virus and to its diffusion, but the need of purposed drugs is emerging as a core target of the fight against COVID19. In this respect, the repurposing of drugs has already led to discovery of drugs thought to reduce the effects of the cytokine storm, but still a drug targeting the spike protein, in the infection stage, is missing. In this work, we present a multifaceted computational approach strongly grounded on a biophysical modeling of biological systems, so to disclose the interaction of the SARS-CoV2 spike protein with ACE2 with a special focus to an allosteric regulation of the spike-ACE2 interaction. Our approach includes the following methodologies: Protein Contact Networks and Network Clustering, Targeted Molecular Dynamics, Elastic Network Modeling, Perturbation Response Scanning, and a computational analysis of energy flow and SEPAS as a protein-softness and monomer-based affinity predictor. We applied this approach to free (closed and open) states of spike protein and spike-ACE2 complexes. Eventually, we analyzed the interactions of free and bound forms of spike with hepcidin (HPC), the major hormone in iron regulation, recently addressed as a central player in the COVID19 pathogenesis, with a special emphasis to the most severe outcomes. Our results demonstrate that, compared with closed and open states, the spike protein in the ACE2-bound state shows higher allosteric potential. The correspondence between hinge sites and the Allosteric Modulation Region (AMR) in the S-ACE complex suggests a molecular basis for hepcidin involvement in COVID19 pathogenesis. We verify the importance of AMR in different states of spike and then study its interactions with HPC and the consequence of the HPC-AMR interaction on spike dynamics and its affinity for ACE2. We propose two complementary mechanisms for HPC effects on spike of SARS-CoV-2; (a) HPC acts as a competitive inhibitor when spike is in a preinfection state (open and with no ACE2), (b) the HPC-AMR interaction pushes the spike structure into the safer closed state. These findings need clear molecular in vivo verification beside clinical observations.

Effect of Lithium Drug on Binding Affinities of Glycogen Synthase Kinase-3 ß to Its Network Partners: A New Computational Approach.

Finding new methods to study the effect of small molecules on protein interaction networks provides us with invaluable tools in the fields of pharmacodynamics and drug design. Lithium is an antimanic drug that has been used for the treatment of bipolar disorder for more than 60 years. Here, we utilized a new approach to study the effect of lithium as a drug on the protein interaction network of GSK-3ß as a hub protein and computed the affinities of GSK-3ß to its partners in the presence of lithium or sodium ions. For this purpose, ensembles of GSK-3ß protein structures were created in the presence of either lithium or sodium ions using adaptive tempering molecular dynamics simulations. The protein binding patches of GSK-3ß for its partners were determined, and finally, the affinity of each binding patch to the related partner was computed for structures of ensembles using a monomer-based approach. Besides, by comparing structural dynamics of GSK-3ß during MD simulations in the presence of LiCl and NaCl, we suggested a new mechanism for the inhibitory effect of lithium on GSK-3ß.

Studying the Effects of ACE2 Mutations on the Stability, Dynamics, and Dissociation Process of SARS-CoV-2 S1/hACE2 Complexes.

A highly infectious coronavirus, SARS-CoV-2, has spread in many countries. This virus recognizes its receptor, angiotensin-converting enzyme 2 (ACE2), using the receptor binding domain of its spike protein subunit S1. Many missense mutations are reported in various human populations for the ACE2 gene. In the current study, we predict the affinity of many ACE2 variants for binding to S1 protein using different computational approaches. The dissociation process of S1 from some variants of ACE2 is studied in the current work by molecular dynamics approaches. We study the relation between structural dynamics of ACE2 in closed and open states and its affinity for S1 protein of SARS-CoV-2.

The Discovery of a Putative Allosteric Site in the SARS-CoV-2 Spike Protein Using an Integrated Structural/Dynamic Approach.

SARS-CoV-2 has caused the largest pandemic of the twenty-first century (COVID-19), threatening the life and economy of all countries in the world. The identification of novel therapies and vaccines that can mitigate or control this global health threat is among the most important challenges facing biomedical sciences. To construct a long-term strategy to fight both SARS-CoV-2 and other possible future threats from coronaviruses, it is critical to understand the molecular mechanisms underlying the virus action. The viral entry and associated infectivity stems from the formation of the SARS-CoV-2 spike protein complex with angiotensin-converting enzyme 2 (ACE2). The detection of putative allosteric sites on the viral spike protein molecule can be used to elucidate the molecular pathways that can be targeted with allosteric drugs to weaken the spike-ACE2 interaction and, thus, reduce viral infectivity. In this study, we present the results of the application of different computational methods aimed at detecting allosteric sites on the SARS-CoV-2 spike protein. The adopted tools consisted of the protein contact networks (PCNs), SEPAS (Affinity by Flexibility), and perturbation response scanning (PRS) based on elastic network modes. All of these methods were applied to the ACE2 complex with both the SARS-CoV2 and SARS-CoV spike proteins. All of the adopted analyses converged toward a specific region (allosteric modulation region [AMR]), present in both complexes and predicted to act as an allosteric site modulating the binding of the spike protein with ACE2. Preliminary results on hepcidin (a molecule with strong structural and sequence with AMR) indicated an inhibitory effect on the binding affinity of the spike protein toward the ACE2 protein.

Complex Stability is Encoded in Binding Patch Softness: a Monomer-Based Approach to Predict Inter-Subunit Affinity of Protein Dimers.

Knowledge about the structure and stability of protein-protein interactions is vital to decipher the behavior of protein systems. Prediction of protein complexes' stability is an interesting topic in the field of structural biology. There are some promising published computational approaches that predict the affinity between subunits of protein dimers using 3D structures of both subunits. In the current study, we classify protein complexes with experimentally measured affinities into distinct classes with different mean affinities. By predicting the mechanical stiffness of the protein binding patch (PBP) region on a single subunit, we successfully predict the assigned affinity class of the PBP in the classification step. Now to predict the experimentally measured affinity between protein monomers in solution, we just need the 3D structure of the suggested PBP on one subunit of the proposed dimer. We designed the SEPAS software and have made the software freely available for academic non-commercial research purposes at " http://biophysics.ir/affinity ". SEPAS predicts the stability of the intended dimer in a classwise manner by utilizing the computed mechanical stiffness of the introduced binding site on one subunit with the minimum accuracy of 0.72.

Soft regions of protein surface are potent for stable dimer formation.

By having knowledge about the characteristics of protein interaction interfaces, we will be able to manipulate protein complexes for therapies. Dimer state is considered as the primary alphabet of the most proteins' quaternary structure. The properties of binding interface between subunits and of noninterface region define the specificity and stability of the intended protein complex. Considering some topological properties and amino acids' affinity for binding in interfaces of protein dimers, we construct the interface-specific recurrence plots. The data obtained from recurrence quantitative analysis, and accessibility-related metrics help us to classify the protein dimers into four distinct classes. Some mechanical properties of binding interfaces are computed for each predefined class of the dimers. The computed mechanical characteristics of binding patch region are compared with those of nonbinding region of proteins. Our observations indicate that the mechanical properties of protein binding sites have a decisive impact on determining the dimer stability. We introduce a new concept in analyzing protein structure by considering mechanical properties of protein structure. We conclude that the interface region between subunits of stable dimers is usually mechanically softer than the interface of unstable protein dimers. AbbreviationsAABaverage affinity for bindingANManisotropic network modelAPCaffinity propagation clusteringASAaccessible surface areaCCDinter residues distanceCSCcomplex stability codeDMdistance matrixΔGdissPISA-computed dissociation free energyGNMGaussian normal mode analysisNMAnormal mode analysisPBPprotein binding patchPISAproteins, interfaces, structures and assembliesrASArelative accessible area in respect to unfolded state of residuesRMrecurrence matrixrPrelative protrusionRPrecurrence plotRQArecurrence quantitative analysisSEMstandard error of meanCommunicated by Ramaswamy H. Sarma.

Partner-Specific Prediction of Protein-Dimer Stability from Unbound Structure of Monomer.

Protein complexes play deterministic roles in live entities in sensing, compiling, controlling, and responding to external and internal stimuli. Thermodynamic stability is an important property of protein complexes; having knowledge about complex stability helps us to understand the basics of protein assembly-related diseases and the mechanism of protein assembly clearly. Enormous protein-protein interactions, detected by high-throughput methods, necessitate finding fast methods for predicting the stability of protein assemblies in a quantitative and qualitative manner. The existing methods of predicting complex stability need knowledge about the three-dimensional (3D) structure of the intended protein complex. Here, we introduce a new method for predicting dissociation free energy of subunits by analyzing the structural and topological properties of a protein binding patch on a single subunit of the desired protein complex. The method needs the 3D structure of just one subunit and the information about the position of the intended binding site on the surface of that subunit to predict dimer stability in a classwise manner. The patterns of structural and topological properties of a protein binding patch are decoded by recurrence quantification analysis. Nonparametric discrimination is then utilized to predict the stability class of the intended dimer with accuracy greater than 85%.

A hidden aggregation-prone structure in the heart of hypoxia inducible factor prolyl hydroxylase.

Prolyl hydroxylase domain-containing protein 2 (PHD2), as one of the most important regulators of angiogenesis and metastasis of cancer cells, is a promising target for cancer therapy drug design. Progressive studies imply that abnormality in PHD2 function may be due to misfolding. Therefore, study of the PHD2 unfolding pathway paves the way for a better understanding of the influence of PHD2 mutations and cancer cell metabolites on the protein folding pathway. We study the unfolding of the PHD2 catalytic domain using differential scanning calorimetry (DSC), fluorescence spectroscopy, and discrete molecular dynamics simulations (DMD). Using computational and experimental techniques, we find that PHD2 undergoes four transitions along the thermal unfolding pathway. To illustrate PHD2 unfolding events in atomic detail, we utilize DMD simulations. Analysis of computational results indicates an intermediate species in the PHD2 unfolding pathway that may enhance aggregation propensity, explaining mutation-independent PHD2 malfunction.

Journey of poly-nucleotides through OmpF porin.

OmpF is an abundant porin in many bacteria which attracts attention as a promising biological nanopore for DNA sequencing. We study the interactions of OmpF with pentameric poly-nucleotides (poly-Ns) in silico. The poly-N molecule is forced to translocate through the lumen of OmpF. Subsequently, the structural and dynamical effects of translocation steps on protein and poly-N molecules are explored in detail. The external loops of OmpF are introduced as the main region for discrimination of poly-Ns based on their organic bases. Structural network analyses of OmpF in the presence or absence of poly-Ns characterize special residues in the structural network of porin. These residues pave the way for engineering OmpF protein. The poly-N-specific pattern of OmpF's local conductance is detected in the current study. Computing the potential of mean force for translocation steps, we define the energetic barrier ahead of poly-N to move through OmpF's lumen. We suggest that fast translocation of the examined poly-N molecules through OmpF seems unattainable by small external driving forces. Our computational results suggest some abilities for OmpF porin like OmpF's potential for being used in poly-N sequencing.

Caseoperoxidase, mixed ß-casein-SDS-hemin-imidazole complex: a nano artificial enzyme.

A novel peroxidase-like artificial enzyme, named "caseoperoxidase", was biomimetically designed using a nano artificial amino acid apo-protein hydrophobic pocket. This four-component nano artificial enzyme containing heme-imidazole-ß-casein-SDS exhibited high activity growth and k(cat) performance toward the native horseradish peroxidase demonstrated by the steady state kinetics using UV-vis spectrophotometry. The hydrophobicity and secondary structure of the caseoperoxidase were studied by ANS fluorescence and circular dichroism spectroscopy. Camel ß-casein (Cß-casein) was selected as an appropriate apo-protein for the heme active site because of its innate flexibility and exalted hydrophobicity. This selection was confirmed by homology modeling method. Heme docking into the newly obtained Cß-casein structure indicated one heme was mainly incorporated with Cß-casein. The presence of a main electrostatic site for the active site in the Cß-casein was also confirmed by experimental methods through Wyman binding potential and isothermal titration calorimetry. The existence of Cß-casein protein in this biocatalyst lowered the suicide inactivation and provided a suitable protective role for the heme active-site. Additional experiments confirmed the retention of caseoperoxidase structure and function as an artificial enzyme.

Implication of disulfide bridge induced thermal reversibility, structural and functional stability for luciferase.

Firefly luciferase is a relatively unstable protein and commonly loses its activity at room temperature because of structural changes. The structural and functional stability of this protein is critical for its enzymatic applications. Different approaches are applied to increase the stability of this enzyme such as designing of covalent cross-links (disulfide bonds). In this study, luciferase mutants containing one or two disulfide bonds were compared to the native protein for their for their structural, thermodynamic, and functional properties. Mutant forms of P. Pyralis luciferase A²96C-A³²6C and A²96C-A³²6C/P45¹C-V469C were used. Thermodynamic and biophysical studies were carried out using UV-Vis, fluorescence, circular dichroism, luminescence spectroscopy and differential scanning calorimetry (DSC). We observed that both mutant forms of the protein were more stable than the wild-type enzyme. However, the single disulfide bond containing mutant was structurally and functionally more stable than the mutant protein containing two disulfide bonds. Furthermore, the enzymatic activity of the single disulfide bond containing mutant protein was 7-folds greater than the wild type and the double disulfide bond proteins. The A²96C-A³²6C mutation also increased the reversibility and disaggregation of the protein. The enhanced activity of the single disulfide bond mutant protein was contributed to the expansion of its active site cleft, which was confirmed by bioinformatics tools.

Application of OmpF nanochannel forming protein in polynucleotide sequence recognition.

Recognition of the sequence of human genome sequence is vital to address malfunctions occurring at molecular, cellular and tissue levels and requires a great deal of time, cost and efforts. Thus, various synthetic and natural pores were considered to fabricate high-throughput systems capable to fulfill the task in an efficient manner. Here, voltage gating OmpF nanochannel, whose structure is known at an atomic level, was used to recognize and differentiate between polynucleotide primers through voltage clamp technique. Our results showed that poly(T) occasionally blocked the channel at both polarities, while poly(C) and poly(G) obstructed it only at positive polarity. The channel was blocked at potential differences of as low as 80 mV in the presence of poly(T). The conductance of channel decreased in the presence of poly(C) and poly(G) by 61 and 5% respectively. Analysis of the number of events showed that poly(T) caused more closing events at higher voltages, while poly(G) and poly(C) induced it at lower voltages. Application of the hazard function as a statistical parameter and analysis of event closing times in various voltages demonstrated the most efficient differentiation at 60 mV. The results of practical and theoretical approaches presented here show that OmpF porin channel possesses the structural and dynamic characteristics required to be considered as a biosensor capable for continuous polynucleotide sequencing.

Inhibition study on insulin fibrillation and cytotoxicity by paclitaxel.

Alzheimer, a neurodegenerative disease, and a large variety of pathologic conditions are associated with a form of protein aggregation known as amyloid fibrils. Since fibrils and prefibrillar intermediates are cytotoxic, numerous attempts have been made to inhibit fibrillation process as a therapeutic strategy. Peptides, surfactants and aromatic small molecules have been used as fibrillation inhibitors. Here we studied the effects of paclitaxel, a polyphenol with a high tendency for interaction with proteins, on fibrillation of insulin as a model protein. The effects of paclitaxel on insulin fibrillation were determined by Thioflavin T fluorescence, Congo red absorbance, circular dichroism and atomic force microscopy. These studies indicated that paclitaxel considerably hindered nucleation, and therefore, fibrillation of insulin in a dose-dependant manner. The isothermal titration calorimetry studies showed that the interaction between paclitaxel and insulin was spontaneous. In addition, the van der Waal's interactions and hydrogen bonds were prominent forces contributing to this interaction. Computational results using molecular dynamic simulations and docking studies revealed that paclitaxel diminished the polarity of insulin dimer and electrostatic interactions by increasing the hydrophobicity of its dimer state. Furthermore, paclitaxel reduced disrupting effects of insulin fibrils on PC12 cell's neurite outgrowth and complexity, and enhanced their survival.

Interaction of 2-APB, dantrolene, and TDMT with IP3R and RyR modulates ER stress-induced programmed cell death I and II in neuron-like PC12 cells: an experimental and computational investigation.

Ca(2+) is an essential second messenger, playing a fundamental role in maintaining cell viability and neuronal activity. Two specific endoplasmic reticulum calcium channels, ryanodine receptors (RyRs) and inositol 1,4,5-trisphosphate receptors (IP3Rs) play an important role in Ca(2+) regulation. In the present study, we provided a 3D structure of RyR and IP3R by homology modeling, and we predicted their interactions with a known neuroprotective compound, 3-thiomethyl-5,6-(dimethoxyphenyl)-1,2,4-triazine (TDMT), as well as two inhibitors, dantrolene and 2-aminoethoxydiphenyl borate (2-APB). Interestingly, we found that dantrolene and 2-APB can bind to the IP3-binding domain of IP3R and RyR, while TDMT may directly block both channels by interacting with the putative resident domains in the pore. Cell culture experiments showed that these compounds could protect PC12 cells against H2O2-induced apoptosis and activate autophagic pathways. Collectively, our computational (in silico) and cell culture studies suggest that RyR and IP3R are novel and promising targets to be used against neurodegenerative diseases.

Thermal unfolding pathway of PHD2 catalytic domain in three different PHD2 species: computational approaches.

Prolyl hydroxylase domain 2 containing protein (PHD2) is a key protein in regulation of angiogenesis and metastasis. In normoxic condition, PHD2 triggers the degradation of hypoxia-inducible factor 1 (HIF-1α) that induces the expression of hypoxia response genes. Therefore the correct function of PHD2 would inhibit angiogenesis and consequent metastasis of tumor cells in normoxic condition. PHD2 mutations were reported in some common cancers. However, high levels of HIF-1α protein were observed even in normoxic metastatic tumors with normal expression of wild type PHD2. PHD2 malfunctions due to protein misfolding may be the underlying reason of metastasis and invasion in such cases. In this study, we scrutinize the unfolding pathways of the PHD2 catalytic domain's possible species and demonstrate the properties of their unfolding states by computational approaches. Our study introduces the possibility of aggregation disaster for the prominent species of PHD2 during its partial unfolding. This may justify PHD2 inability to regulate HIF-1α level in some normoxic tumor types.

Inhibition of mushroom tyrosinase by a newly synthesized ligand: inhibition kinetics and computational simulations.

Alterations in the synthesis of melanin contribute to a number of diseases; therefore, the design of new tyrosinase inhibitors is very important. Mushroom tyrosinase (MT) is a metalloenzyme, which plays an important role in melanin biosynthesis. In this study, the inhibitory effect of a novel designed compound, i.e. 2-((1Z)-(2-(2,4-dinitrophenyl)hydrazin-1-ylidene)methyl)phenol, as a specific ligand which can bind to the copper ion of MT, has been assessed. The ligand was found to competitively inhibit both the cresolase and catecholase activities of MT, with small inhibition constants of 2.8 and 2.6 µM, respectively. Intrinsic fluorescence studies were performed to gain more information on the binding constants. Docking results indicated that the ligand binds to copper ions in the active site of MT via the OH group of the ligand. The ligand makes four hydrogen bonds with aspartic acid and one hydrogen bond with the histidine residue in the active site. Molecular dynamics results show that ligand binds to the MT via both electrostatic and hydrophobic interactions with its different parts.

The inhibitory effect of ethylenediamine on mushroom tyrosinase.

The inhibitory effect of ethylenediamine on both activities of mushroom tyrosinase (MT) at 20 °C in a 10 mM phosphate buffer solution (pH 6.8), was studied. L-DOPA and L-tyrosine were used as substrates of catecholase and cresolase activities, respectively. The results showed that ethylenediamine competitively inhibits both activities of the enzyme with inhibition constants (K(i)) of 0.18±0.05 and 0.14±0.01 µM for catecholase and cresolase respectively, which are lower than the reported values for other MT inhibitors. For further insight a docking study between tyrosinase and ethylenediamine was performed. The docking simulation showed that ethylenediamine binds in the active site of the enzyme near the Cu atoms and makes 3 hydrogen bonds with two histidine residues of active site.

Disruption of tubulin polymerization and cell proliferation by 1-naphthylarsonic acid.

Arsenical compounds exhibit a differential toxicity to cancer cells. Microtubules are a primary target of a number of anticancer drugs, such as arsenical compounds. The interaction of 1-NAA (1-naphthylarsonic acid) has been investigated on microtubule polymerization under in vitro and cellular conditions. Microtubules were extracted from sheep brain. Transmission electron microscopy was used to show microtubule structure in the presence of 1-NAA. Computational docking method was applied for the discovery of ligand-binding sites on the microtubular proteins. Proliferation of HeLa cells and HF2 (human foreskin fibroblasts) was measured by the MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2H-tetrazolium bromide] assay method following their incubation with 1-NAA. Fluorescence microscopic labelling was done with the help of α-tubulin monoclonal antibody and Tunel kit was used to investigate the apoptotic effects of 1-NAA on the HeLa cells. 1-NAA inhibits the tubulin polymerization by the formation of abnormal polymers having high affinity to the inner cell wall.

Homology modeling, docking, molecular dynamics simulation, and structural analyses of coxsakievirus B3 2A protease: an enzyme involved in the pathogenesis of inflammatory myocarditis.

2A protease of the pathogenic coxsackievirus B3 is key to the pathogenesis of inflammatory myocarditis and, therefore, an attractive drug target. However lack of a crystal structure impedes design of inhibitors. Here we predict 3D structure of CVB3 2A(pro) based on sequence comparison and homology modeling with human rhinovirus 2A(pro). The two enzymes are remarkably similar in their core regions. However they have different conformations at the N-terminal. A large number of N-terminal hydrophobic residues reduce the thermal stability of CVB3 2A(pro), as we confirmed by fluorescence, western blot and turbidity measurement. Molecular dynamic simulation revealed that elevated temperature induces protein motion that results in frequent movement of the N-terminal coil. This may therefore induce successive active site changes and thus play an important role in destabilization of CVB3 2A(pro) structure.

A folding pathway-dependent score to recognize membrane proteins.

While various approaches exist to study protein localization, it is still a challenge to predict where proteins localize. Here, we consider a mechanistic viewpoint for membrane localization. Taking into account the steps for the folding pathway of α-helical membrane proteins and relating biophysical parameters to each of these steps, we create a score capable of predicting the propensity for membrane localization and call it FP(3)mem. This score is driven from the principal component analysis (PCA) of the biophysical parameters related to membrane localization. FP(3)mem allows us to rationalize the colocalization of a number of channel proteins with the Cav1.2 channel by their fewer propensities for membrane localization.