Integrative modeling of guanylate binding protein dimers.

Guanylate-binding proteins (GBPs) are essential interferon-γ-activated large GTPases that play a crucial role in host defense against intracellular bacteria and parasites. While their protective functions rely on protein polymerization, our understanding of the structural intricacies of these multimerized states remains limited. To bridge this knowledge gap, we present dimer models for human GBP1 (hGBP1) and murine GBP2 and 7 (mGBP2 and mGBP7) using an integrative approach, incorporating the crystal structure of hGBP1's GTPase domain dimer, crosslinking mass spectrometry, small-angle X-ray scattering, protein-protein docking, and molecular dynamics simulations. Our investigation begins by comparing the protein dynamics of hGBP1, mGBP2, and mGBP7. We observe that the M/E domain in all three proteins exhibits significant mobility and hinge motion, with mGBP7 displaying a slightly less pronounced motion but greater flexibility in its GTPase domain. These dynamic distinctions can be attributed to variations in the sequences of mGBP7 and hGBP1/mGBP2, resulting in different dimerization modes. Unlike hGBP1 and its close ortholog mGBP2, which exclusively dimerize through their GTPase domains, we find that mGBP7 exhibits three equally probable alternative dimer structures. The GTPase domain of mGBP7 is only partially involved in its dimerization, primarily due to an accumulation of negative charge, allowing mGBP7 to dimerize independently of GTP. Instead, mGBP7 exhibits a strong tendency to dimerize in an antiparallel arrangement across its stalks. The results of this work go beyond the sequence-structure-function relationship, as the sequence differences in mGBP7 and mGBP2/hGBP1 do not lead to different structures, but to different protein dynamics and dimerization. The distinct GBP dimer structures are expected to encode specific functions crucial for disrupting pathogen membranes.

Domain motions, dimerization, and membrane interactions of the murine guanylate binding protein 2.

Guanylate-binding proteins (GBPs) are a group of GTPases that are induced by interferon-[Formula: see text] and are crucial components of cell-autonomous immunity against intracellular pathogens. Here, we examine murine GBP2 (mGBP2), which we have previously shown to be an essential effector protein for the control of Toxoplasma gondii replication, with its recruitment through the membrane of the parasitophorous vacuole and its involvement in the destruction of this membrane likely playing a role. The overall aim of our work is to provide a molecular-level understanding of the mutual influences of mGBP2 and the parasitophorous vacuole membrane. To this end, we performed lipid-binding assays which revealed that mGBP2 has a particular affinity for cardiolipin. This observation was confirmed by fluorescence microscopy using giant unilamellar vesicles of different lipid compositions. To obtain an understanding of the protein dynamics and how this is affected by GTP binding, mGBP2 dimerization, and membrane binding, assuming that each of these steps are relevant for the function of the protein, we carried out standard as well as replica exchange molecular dynamics simulations with an accumulated simulation time of more than 30 µs. The main findings from these simulations are that mGBP2 features a large-scale hinge motion in its M/E domain, which is present in each of the studied protein states. When bound to a cardiolipin-containing membrane, this hinge motion is particularly pronounced, leading to an up and down motion of the M/E domain on the membrane, which did not occur on a membrane without cardiolipin. Our prognosis is that this up and down motion has the potential to destroy the membrane following the formation of supramolecular mGBP2 complexes on the membrane surface.

SEC14-GOLD protein PATELLIN2 binds IRON-REGULATED TRANSPORTER1 linking root iron uptake to vitamin E.

Organisms require micronutrients, and Arabidopsis (Arabidopsis thaliana) IRON-REGULATED TRANSPORTER1 (IRT1) is essential for iron (Fe2+) acquisition into root cells. Uptake of reactive Fe2+ exposes cells to the risk of membrane lipid peroxidation. Surprisingly little is known about how this is avoided. IRT1 activity is controlled by an intracellular variable region (IRT1vr) that acts as a regulatory protein interaction platform. Here, we describe that IRT1vr interacted with peripheral plasma membrane SEC14-Golgi dynamics (SEC14-GOLD) protein PATELLIN2 (PATL2). SEC14 proteins bind lipophilic substrates and transport or present them at the membrane. To date, no direct roles have been attributed to SEC14 proteins in Fe import. PATL2 affected root Fe acquisition responses, interacted with ROS response proteins in roots, and alleviated root lipid peroxidation. PATL2 had high affinity in vitro for the major lipophilic antioxidant vitamin E compound α-tocopherol. Molecular dynamics simulations provided insight into energetic constraints and the orientation and stability of the PATL2-ligand interaction in atomic detail. Hence, this work highlights a compelling mechanism connecting vitamin E with root metal ion transport at the plasma membrane with the participation of an IRT1-interacting and α-tocopherol-binding SEC14 protein.

Development of a First-in-Class Small-Molecule Inhibitor of the C-Terminal Hsp90 Dimerization.

Heat shock proteins 90 (Hsp90) are promising therapeutic targets due to their involvement in stabilizing several aberrantly expressed oncoproteins. In cancerous cells, Hsp90 expression is elevated, thereby exerting antiapoptotic effects, which is essential for the malignant transformation and tumor progression. Most of the Hsp90 inhibitors (Hsp90i) under investigation target the ATP binding site in the N-terminal domain of Hsp90. However, adverse effects, including induction of the prosurvival resistance mechanism (heat shock response or HSR) and associated dose-limiting toxicity, have so far precluded their clinical approval. In contrast, modulators that interfere with the C-terminal domain (CTD) of Hsp90 do not inflict HSR. Since the CTD dimerization of Hsp90 is essential for its chaperone activity, interfering with the dimerization process by small-molecule protein-protein interaction inhibitors is a promising strategy for anticancer drug research. We have developed a first-in-class small-molecule inhibitor (5b) targeting the Hsp90 CTD dimerization interface, based on a tripyrimidonamide scaffold through structure-based molecular design, chemical synthesis, binding mode model prediction, assessment of the biochemical affinity, and efficacy against therapy-resistant leukemia cells. 5b reduces xenotransplantation of leukemia cells in zebrafish models and induces apoptosis in BCR-ABL1+ (T315I) tyrosine kinase inhibitor-resistant leukemia cells, without inducing HSR.

Multiscale MD simulations of wild-type and sickle hemoglobin aggregation.

Sickle cell disease is a hemoglobinopathy resulting from a point mutation from glutamate to valine at position six of the ß-globin chains of hemoglobin. This mutation gives rise to pathological aggregation of the sickle hemoglobin and, as a result, impaired oxygen binding, misshapen and short-lived erythrocytes, and anemia. We aim to understand the structural effects caused by the single Glu6Val mutation leading to protein aggregation. To this end, we perform multiscale molecular dynamics simulations employing atomistic and coarse-grained models of both wild-type and sickle hemoglobin. We analyze the dynamics of hemoglobin monomers and dimers, study the aggregation of wild-type and sickle hemoglobin into decamers, and analyze the protein-protein interactions in the resulting aggregates. We find that the aggregation of sickle hemoglobin is driven by both hydrophobic and electrostatic protein-protein interactions involving the mutation site and surrounding residues, leading to an extended interaction area and thus stable aggregates. The wild-type protein can also self-assemble, which, however, results from isolated interprotein salt bridges that do not yield stable aggregates. This knowledge can be exploited for the development of sickle hemoglobin-aggregation inhibitors.

Disorder-to-order transition of the amyloid-ß peptide upon lipid binding.

There is mounting evidence that Alzheimer's disease progression and severity are linked to neuronal membrane damage caused by aggregates of the amyloid-ß (Aß) peptide. However, the detailed mechanism behind the membrane damage is not well understood yet. Recently, the lipid-chaperone hypothesis has been put forward, based on which the formation of complexes between Aß and free lipids enables an easy insertion of Aß into membranes. In order to test this hypothesis, we performed numerous all-atom molecular dynamics simulations. We studied the complex formation between individual lipids, considering both POPC and DPPC, and Aß and examined whether the resulting complexes would be able to insert into lipid membranes. Complex formation at a one-to-one ratio was readily observed, yet with minimal effects on Aß's characteristics. Most importantly, the peptide remains largely disordered in 1:1 complexes, and the complex does not insert into the membrane; instead, it is adsorbed to the membrane surface. The results change considerably once Aß forms a complex with a POPC cluster composed of three lipid molecules. The hydrophobic interactions between Aß and the lipid tails cause the peptide to fold into either a helical or a ß-sheet structure. These observations provide atomic insight into the disorder-to-order transition that is needed for membrane insertion or amyloid aggregation to proceed.

Dataset of AMBER force field parameters of drugs, natural products and steroids for simulations using GROMACS.

We provide general AMBER force field (GAFF) parameters for 160 organic molecules including drugs, natural products, and steroids, which can be employed without further processing in molecular dynamics (MD) simulations using GROMACS. We determined these parameters based on quantum mechanical (QM) calculations involving geometry optimization at the HF6-31G* level of theory. For each molecule we provide a coordinate file of the three-dimensional molecular structure, the topology and the parameter file. The applicability of these parameters was demonstrated by MD simulations of these molecules bound to the active site of the main protease of the coronavirus SARS-CoV-2, 3CLpro, which is a main player during viral replication causing COVID-19.

Novel inhibitors of the main protease enzyme of SARS-CoV-2 identified via molecular dynamics simulation-guided in vitro assay.

For the COVID-19 pandemic caused by SARS-CoV-2, there are currently no effective drugs or vaccines to treat this coronavirus infection. In this study, we focus on the main protease enzyme of SARS-CoV-2, 3CLpro, which is critical for viral replication. We employ explicit solvent molecular dynamics simulations of about 150 compounds docked into 3CLpro's binding site and that had emerged as good main protease ligands from our previous in silico screening of over 1.2 million compounds. By incoporating protein dynamics and applying a range of structural descriptors, such as the ability to form specific contacts with the catalytic dyad residues of 3CLpro and the structural fluctuations of the ligands in the binding site, we are able to further refine our compound selection. Fourteen compounds including estradiol shown to be the most promising based on our calculations were procured and screened against recombinant 3CLpro in a fluorescence assay. Eight of these compounds have significant activity in inhibiting the SARS-CoV-2 main protease. Among these are corilagin, a gallotannin, and lurasidone, an antipsychotic drug, which emerged as the most promising natural product and drug, respectively, and might thus be candidates for drug repurposing for the treatment of COVID-19. In addition, we also tested the inhibitory activity of testosterone, and our results reveal testosterone as possessing moderate inhibitory potency against the 3CLpro enzyme, which may thus provide an explanation why older men are more severely affected by COVID-19.

High Throughput Virtual Screening to Discover Inhibitors of the Main Protease of the Coronavirus SARS-CoV-2.

We use state-of-the-art computer-aided drug design (CADD) techniques to identify prospective inhibitors of the main protease enzyme, 3CLpro of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing COVID-19. From our screening of over one million compounds including approved drugs, investigational drugs, natural products, and organic compounds, and a rescreening protocol incorporating enzyme dynamics via ensemble docking, we have been able to identify a range of prospective 3CLpro inhibitors. Importantly, some of the identified compounds had previously been reported to exhibit inhibitory activities against the 3CLpro enzyme of the closely related SARS-CoV virus. The top-ranking compounds are characterized by the presence of multiple bi- and monocyclic rings, many of them being heterocycles and aromatic, which are flexibly linked allowing the ligands to adapt to the geometry of the 3CLpro substrate site and involve a high amount of functional groups enabling hydrogen bond formation with surrounding amino acid residues, including the catalytic dyad residues H41 and C145. Among the top binding compounds we identified several tyrosine kinase inhibitors, which include a bioflavonoid, the group of natural products that binds best to 3CLpro. Another class of compounds that decently binds to the SARS-CoV-2 main protease are steroid hormones, which thus may be endogenous inhibitors and might provide an explanation for the age-dependent severity of COVID-19. Many of the compounds identified by our work show a considerably stronger binding than found for reference compounds with in vitro demonstrated 3CLpro inhibition and anticoronavirus activity. The compounds determined in this work thus represent a good starting point for the design of inhibitors of SARS-CoV-2 replication.

Computer simulations of protein-membrane systems.

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Biochemical and structural characterization of murine GBP7, a guanylate binding protein with an elongated C-terminal tail.

Guanylate-binding proteins (GBPs) constitute a family of interferon-inducible guanosine triphosphatases (GTPases) that are key players in host defense against intracellular pathogens ranging from protozoa to bacteria and viruses. So far, human GBP1 and GBP5 as well as murine GBP2 (mGBP2) have been biochemically characterized in detail. Here, with murine GBP7 (mGBP7), a GBP family member with an unconventional and elongated C-terminus is analyzed. The present study demonstrates that mGBP7 exhibits a concentration-dependent GTPase activity and an apparent GTP turnover number of 20 min-1. In addition, fluorescence spectroscopy analyses reveal that mGBP7 binds GTP with high affinity (KD = 0.22 µM) and GTPase activity assays indicate that mGBP7 hydrolyzes GTP to GDP and GMP. The mGBP7 GTPase activity is inhibited by incubation with γ-phosphate analogs and a K51A mutation interfering with GTP binding. SEC-MALS analyses give evidence that mGBP7 forms transient dimers and that this oligomerization pattern is not influenced by the presence of nucleotides. Moreover, a structural model for mGBP7 is provided by homology modeling, which shows that the GTPase possesses an elongated C-terminal (CT) tail compared with the CaaX motif-containing mGBP2 and human GBP1. Molecular dynamics simulations indicate that this tail has transmembrane characteristics and, interestingly, confocal microscopy analyses reveal that the CT tail is required for recruitment of mGBP7 to the parasitophorous vacuole of Toxoplasma gondii.

Large-scale, dynamin-like motions of the human guanylate binding protein 1 revealed by multi-resolution simulations.

Guanylate binding proteins (GBPs) belong to the dynamin-related superfamily and exhibit various functions in the fight against infections. The functions of the human guanylate binding protein 1 (hGBP1) are tightly coupled to GTP hydrolysis and dimerization. Despite known crystal structures of the hGBP1 monomer and GTPase domain dimer, little is known about the dynamics of hGBP1. To gain a mechanistic understanding of hGBP1, we performed sub-millisecond multi-resolution molecular dynamics simulations of both the hGBP1 monomer and dimer. We found that hGBP1 is a highly flexible protein that undergoes a hinge motion similar to the movements observed for other dynamin-like proteins. Another large-scale motion was observed for the C-terminal helix α13, providing a molecular view for the α13-α13 distances previously reported for the hGBP1 dimer. Most of the loops of the GTPase domain were found to be flexible, revealing why GTP binding is needed for hGBP1 dimerization to occur.

Interference with Amyloid-ß Nucleation by Transient Ligand Interaction.

Amyloid-ß peptide (Aß) is an intrinsically disordered protein (IDP) associated with Alzheimer's disease. The structural flexibility and aggregation propensity of Aß pose major challenges for elucidating the interaction between Aß monomers and ligands. All-D-peptides consisting solely of D-enantiomeric amino acid residues are interesting drug candidates that combine high binding specificity with high metabolic stability. Here we characterized the interaction between the 12-residue all-D-peptide D3 and Aß42 monomers, and how the interaction influences Aß42 aggregation. We demonstrate for the first time that D3 binds to Aß42 monomers with submicromolar affinities. These two highly unstructured molecules are able to form complexes with 1:1 and other stoichiometries. Further, D3 at substoichiometric concentrations effectively slows down the ß-sheet formation and Aß42 fibrillation by modulating the nucleation process. The study provides new insights into the molecular mechanism of how D3 affects Aß assemblies and contributes to our knowledge on the interaction between two IDPs.