We use state-of-the-art computer-aided drug design (CADD) techniques to identify prospective inhibitors of the main protease enzyme, 3CL
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 3CL
inhibitors. Importantly, some of the identified compounds had previously been reported to exhibit inhibitory activities against the 3CL
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 3CL
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 3CL
. 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 3CL
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.
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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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.
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.
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.
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.
Display omitted
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.
Computer simulations of protein-membrane systems Loschwitz, Jennifer; Olubiyi, Olujide O; Hub, Jochen S ...
Progress in molecular biology and translational science,
2020, Letnik:
170
Journal Article
Recenzirano
Odprti dostop
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.
Abstract
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.
The interaction between a plasma membrane root iron transporter and a SEC14-GOLD peripheral membrane protein is critical for root iron mobilization and antioxidant activities in response to iron.