Binding sites and hydrophobic pockets in Human Thioredoxin 1 determined by normal mode analysis

Effect of phosphorylation on the structural dynamics, thermal stability of human dopamine transporter: A simulation study using normal modes, molecular dynamics and Markov State Model

The Human Dopamine Transporter (hDAT) plays an essential role in modulating the Influx/Efflux of dopamine, and it is involved in the mechanism of certain neurodegenerative diseases such as Parkinson's disease. Several studies have reported important states for Dopamine transport: outward-facing open state (OFo), the outward-facing closed state (OFc), the holo-occluded state closed (holo), and the inward-facing open state (IFo). Furthermore, experimental assays have shown that different phosphorylation conditions in hDAT can affect the rate of dopamine absorption. We present a protocol using hybrid simulation methods to study the conformational dynamics and stability of states of hDAT under different phosphorylation sites. With this protocol, we explored the conformational space of hDAT, identified the states, and evaluated the free energy differences and the transition probabilities between them in each of the phosphorylation cases. We also presented the conformational changes and correlated them with those described in the literature. There is a thesis/hypothesis that the phosphorylation condition corresponding to NP-333 system (where all sites Ser/Thr from residue 2 to 62 and 254 to 613 are phosphorylated, except residue 333) would decrease the rate of dopamine transport from the extracellular medium to the intracellular medium by hDAT as previously described in the literature by Lin et al., 2003. Our results corroborated this thesis/hypothesis and the data reported. It is probably due to the affectation/changes/alteration of the conformational dynamics of this system that makes the intermediate states more likely and makes it difficult to initial states associated with the uptake of dopamine in the extracellular medium, corroborating the experimental results. Furthermore, our results showed that just single phosphorylation/dephosphorylation could alter intrinsic protein motions affecting the sampling of one or more states necessary for dopamine transport. In this sense, the modification of phosphorylation influences protein movements and conformational preferences, affecting the stability of states and the transition between them and, therefore, the transport.

HOPMA: Boosting Protein Functional Dynamics with Colored Contact Maps

In light of the recent very rapid progress in protein structure prediction, accessing the multitude of functional protein states is becoming more central than ever before. Indeed, proteins are flexible macromolecules, and they often perform their function by switching between different conformations. However, high-resolution experimental techniques such as X-ray crystallography and cryogenic electron microscopy can catch relatively few protein functional states. Many others are only accessible under physiological conditions in solution. Therefore, there is a pressing need to fill this gap with computational approaches. We present HOPMA, a novel method to predict protein functional states and transitions by using a modified elastic network model. The method exploits patterns in a protein contact map, taking its 3D structure as input, and excludes some disconnected patches from the elastic network. Combined with nonlinear normal mode analysis, this strategy boosts the protein conformational space exploration, especially when the input structure is highly constrained, as we demonstrate on a set of more than 400 transitions. Our results let us envision the discovery of new functional conformations, which were unreachable previously, starting from the experimentally known protein structures. The method is computationally efficient and available at https://github.com/elolaine/HOPMA and https://team.inria.fr/nano-d/software/nolb-normal-modes.

Towards gaining sight of multiscale events: utilizing network models and normal modes in hybrid methods

With the explosion of normal mode analyses (NMAs) based on elastic network models (ENMs) in the last decade, and the proven precision of MD simulations for visualizing interactions at atomic scale, many hybrid methods have been proposed in recent years. These aim at exploiting the best of both worlds: the atomic precision of MD that often fall short of exploring time and length scales of biological interest, and the capability of ENM-NMA to predict the cooperative and often functional rearrangements of large structures and assemblies, albeit at low resolution. We present an overview of recent progress in the field with examples of successful applications highlighting the utility of such hybrid methods and pointing to emerging future directions guided by advances in experimental characterization of biomolecular systems structure and dynamics.

Molecular, Biological and Structural Features of VL CDR-1 Rb44 Peptide, Which Targets the Microtubule Network in Melanoma Cells

Microtubules are important drug targets in tumor cells, owing to their role in supporting and determining the cell shape, organelle movement and cell division. The complementarity-determining regions (CDRs) of immunoglobulins have been reported to be a source of anti-tumor peptide sequences, independently of the original antibody specificity for a given antigen. We found that, the anti-Lewis B mAb light-chain CDR1 synthetic peptide Rb44, interacted with microtubules and induced depolymerization, with subsequent degradation of actin filaments, leading to depolarization of mitochondrial membrane-potential, increase of ROS, cell cycle arrest at G2/M, cleavage of caspase-9, caspase-3 and PARP, upregulation of Bax and downregulation of Bcl-2, altogether resulting in intrinsic apoptosis of melanoma cells. The in vitro inhibition of angiogenesis was also an Rb44 effect. Peritumoral injection of Rb44L1 delayed growth of subcutaneously grafted melanoma cells in a syngeneic mouse model. L1-CDRs from immunoglobulins and their interactions with tubulin-dimers were explored to interpret effects on microtubule stability. The opening motion of tubulin monomers allowed for efficient L1-CDR docking, impairment of dimer formation and microtubule dissociation. We conclude that Rb44 V_L-CDR1 is a novel peptide that acts on melanoma microtubule network causing cell apoptosis in vitro and melanoma growth inhibition in vivo.

Structural Dynamics of DPP-4 and Its Influence on the Projection of Bioactive Ligands

Dipeptidyl peptidase-4 (DPP-4) is a target to treat type II diabetes mellitus. Therefore, it is important to understand the structural aspects of this enzyme and its interaction with drug candidates. This study involved molecular dynamics simulations, normal mode analysis, binding site detection and analysis of molecular interactions to understand the protein dynamics. We identified some DPP-4 functional motions contributing to the exposure of the binding sites and twist movements revealing how the two enzyme chains are interconnected in their bioactive form, which are defined as chains A (residues 40-767) and B (residues 40-767). By understanding the enzyme structure, its motions and the regions of its binding sites, it will be possible to contribute to the design of new DPP-4 inhibitors as drug candidates to treat diabetes.

Molecular cloning and in silico characterization of knottin peptide, U2-SCRTX-Lit2, from brown spider (Loxosceles intermedia) venom glands

Inhibitor cystine knots (ICKs) are a family of structural peptides with a large number of cysteine residues that form intramolecular disulfide bonds, resulting in a knot. These peptides are involved in a variety of biological functions including predation and defense, and are found in various species, such as spiders, scorpions, sea anemones, and plants. The Loxosceles intermedia venom gland transcriptome identified five groups of ICK peptides that represent more than 50 % of toxin-coding transcripts. Here, we describe the molecular cloning of U2-Sicaritoxin-Lit2 (U2-SCRTX-Lit2), bioinformatic characterization, structure prediction, and molecular dynamic analysis. The sequence of U2-SCRTX-Lit2 obtained from the transcriptome is similar to that of μ-Hexatoxin-Mg2, a peptide that inhibits the insect Nav channel. Bioinformatic analysis of sequences classified as ICK family members also showed a conservation of cysteine residues among ICKs from different spiders, with the three dimensional molecular model of U2-SCRTX-Lit2 similar in structure to the hexatoxin from μ-hexatoxin-Mg2a. Molecular docking experiments showed the interaction of U2-SCRTX-Lit2 to its predictable target-the Spodoptera litura voltage-gated sodium channel (SlNaVSC). After 200 ns of molecular dynamic simulation, the final structure of the complex showed stability in agreement with the experimental data. The above analysis corroborates the existence of a peptide toxin with insecticidal activity from a novel ICK family in L. intermedia venom and demonstrates that this peptide targets Nav channels.

Binding of phenothiazines into allosteric hydrophobic pocket of human thioredoxin 1

Thioredoxins are multifunctional oxidoreductase proteins implicated in the antioxidant cellular apparatus and oxidative stress. They are involved in several pathologies and are promising anticancer targets. Identification of noncatalytic binding sites is of great interest for designing new allosteric inhibitors of thioredoxin. In a recent work, we predicted normal mode motions of human thioredoxin 1 and identified two major putative hydrophobic binding sites. In this work we investigated noncovalent interactions of human thioredoxin 1 with three phenotiazinic drugs acting as prooxidant compounds by using molecular docking and circular dichroism spectrometry to probe ligand binding into the previously predicted allosteric hydrophobic pockets. Our in silico and CD spectrometry experiments suggested one preferred allosteric binding site involving helix 3 and adopting the best druggable conformation identified by NMA. The CD spectra showed binding of thioridazine into thioredoxin 1 and suggested partial helix unfolding, which most probably concerns helix 3. Taken together, these data support the strategy to design thioredoxin inhibitors targeting a druggable allosteric binding site.

Oxidation of structural cysteine residues in thioredoxin 1 by aromatic arsenicals enhances cancer cell cytotoxicity caused by the inhibition of thioredoxin reductase 1

Thioredoxin systems, composed of thioredoxin reductase (TrxR), thioredoxin (Trx) and NADPH, play important roles in maintaining cellular redox homeostasis and redox signaling. Recently the cytosolic Trx1 system has been shown to be a cellular target of arsenic containing compounds. To elucidate the relationship of the structure of arsenic compounds with their ability of inhibiting TrxR1 and Trx1, and cytotoxicity, we have investigated the reaction of Trx1 system with seven arsenic trithiolates: As(Cys)3, As(GS)3, As(Penicillamine)3, As(Mercaptoethanesulfonate)3, As(Mercaptopurine)3, As(2-mercaptopyridine)3 and As(2-mercaptopyridine N-oxide)3. The cytotoxicity of these arsenicals was consistent with their ability to inhibit TrxR1 in vitro and in cells. Unlike other arsenicals, As(Mercaptopurine)3 which did not show inhibitory effects on TrxR1 had very weak cytotoxicity, indicating that TrxR1 is a reliable drug target for arsenicals. Moreover, the two aromatic compounds As(2-mercaptopyridine)3 and As(2-mercaptopyridine N-oxide)3 showed stronger cytotoxicity than the others. As(2-mercaptopyridine)3 which selectively oxidized two structural cysteines (Cys62 and Cys69) in Trx1 showed mild improvement in cytotoxicity. As(2-mercaptopyridine N-oxide)3 oxidized all the Cys residues in Trx1, exhibiting the strongest cytotoxicity. Oxidation of Trx1 by As(2-mercaptopyridine)3 and As(2-mercaptopyridine N-oxide)3 affected electron transfer from NADPH and TrxR1 to peroxiredoxin 1 (Prx1), which could result in the reactive oxygen species elevation and trigger cell death process. These results suggest that oxidation of structural cysteine residues in Trx1 by aromatic group in TrxR1-targeting drugs may sensitize tumor cells to cell death, providing a novel approach to regulate cellular redox signaling and also a basis for rational design of new anticancer agents.

New insights into the posttranslational regulation of human cytosolic thioredoxin by S-palmitoylation

High level of palmitate is associated with metabolic disorders. We recently showed that enhanced level of S-palmitoylated cytosolic thioredoxin (Trx1) in mouse liver was new characteristic feature of insulin resistance. However, our understanding of the effect of S-palmitoylation on Trx1 is limited, and the tissue specificity of Trx1 S-palmitoylation is unclear. Here we show that S-palmitoylation also occurs at Cys73 of Trx1 in living endothelial cells, and the level of S-palmitoylated Trx1 undergoes regulation by insulin signaling. Trx1 prefers thiol-thioester exchange with palmitoyl-CoA to acetyl-CoA. S-palmitoylation alters conformation or secondary structure of Trx1, as well as decreases the ability of Trx1 to transfer electrons from thioredoxin reductase to S-nitrosylated protein-tyrosine phosphatase 1B and S-nitroso-glutathione. Our results demonstrate that S-palmitoylation is an important post-translational modification of human Trx1.