Search

Protein function is closely related to its structure and dynamics. Due to its large number of degrees of freedom, proteins adopt a large number of conformations, which describe a highly complex potential energy landscape. Considering the huge ensemble of conformations in dynamic equilibrium in solution, detailed investigation of proteins dynamics is extremely costly. Therefore, a significant number of different methods have emerged in order to improve the conformational sampling of biomolecules. One of these methods is Molecular Dynamics with excited Normal Modes (MDeNM) in which normal modes are used as collective variables in molecular dynamics. Here, we present a new implementation of the MDeNM method that allows a continuously controlled kinetic excitation energy in the normal mode space, while taking into account the natural constraints imposed either by the structure or the environment. These implementations prevent unphysical structural distortions. We tested the new approach on bacteriophage’s T4 lysozyme, Gallus gallus hen egg-white lysozyme and Staphylococcus aureus membrane-bound transglycosylase. Our results showed that the new approach outperformed free MD sampling and preserved the structural features comparatively to the original MDeNM approach. We also observed that by adaptively changing the excitation direction during calculations, proteins follow new transition paths preventing structural distortions.

We present the first structure of the human Kir2.1 channel containing both transmembrane domain (TMD) and cytoplasmic domain (CTD). Kir2.1 channels are strongly inward-rectifying potassium channels that play a key role in maintaining resting membrane potential. Their gating is modulated by phosphatidylinositol 4,5-bisphosphate (PIP 2 ). Genetically inherited defects in Kir2.1 channels are responsible for several rare human diseases, including Andersen’s syndrome. The structural analysis (cryo–electron microscopy), surface plasmon resonance, and electrophysiological experiments revealed a well-connected network of interactions between the PIP 2 -binding site and the G-loop through residues R312 and H221. In addition, molecular dynamics simulations and normal mode analysis showed the intrinsic tendency of the CTD to tether to the TMD and a movement of the secondary anionic binding site to the membrane even without PIP 2 . Our results revealed structural features unique to human Kir2.1 and provided insights into the connection between G-loop and gating and the pathological mechanisms associated with this channel.

Dimethyl phosphate (DMP- ) is a model for the phosphodiester backbone of DNA, RNA, and phospholipids. It is central for the binding of divalent cations and water along the backbone of nucleic acids. Significant polarization and charge-transfer contributions and nonadditivity come into play in the multimolecular complexes organized around phosphate. Prior to large-scale molecular dynamics (MD) with advanced polarizable potentials, it is essential to evaluate how well the values and trends of intermolecular interaction energies (ΔE) from ab initio quantum chemistry (QC) and their individual contributions are reproduced in a diversity of such complexes. These differ by the starting binding modes of a divalent cation, Zn(II), namely direct, bi- or mono-dentate to anionic and/or ester oxygens, versus through-water binding. We present first the results from automated refinements of the individual contributions of the SIBFA potential with respect to their QC counterparts using a Zn(II) or a water probe. This is followed by validations on eight relaxed multimolecular complexes of DMP- with Zn(II) or Mg(II) and seven waters, then on sixteen complexes of DMP- with Zn(II) and eight waters in arrangements extracted from MD or energy-minimization on a droplet of sixty-four waters. This monitors the compared evolutions of SIBFA and QC ΔE and their individual contributions in the competing arrangements. Some waters, bridging Zn(II) and DMP- , were found to have exceptionally large dipole moments, of up to 3.8 Debye. The perspectives of extension to a flexible phosphodiester backbone are discussed in the context of the SIBFA potential for DNA and RNA.

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.

Made available in DSpace on 2021-06-25T10:12:08Z (GMT). No. of bitstreams: 0 Previous issue date: 2020-12-01 Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) The activation process of phospholipase A2-like (PLA2-like) toxins is a key step in their molecular mechanism, which involves oligomeric changes leading to the exposure of specific sites. Few studies have focused on the characterization of allosteric activators and the features that distinguish them from inhibitors. Herein, a comprehensive study with the BthTX-I toxin from Bothrops jararacussu venom bound or unbound to α-tocopherol (αT) was carried out. The oligomerization state of BthTX-I bound or unbound to αT in solution was studied and indicated that the toxin is predominantly monomeric but tends to oligomerize when complexed with αT. In silico molecular simulations showed the toxin presents higher conformational changes in the absence of αT,which suggests that it is important to stabilize the structure of the toxin. The transition between the two states (active/inactive) was also studied, showing that only the unbound BthTX-I system could migrate to the inactive state. In contrast, the presence of αT induces the toxin to leave the inactive state, guiding it towards the active state, with more regions exposed to the solvent, particularly its active site. Finally, the structural determinants necessary for a molecule to be an inhibitor or activator were analyzed in light of the obtained results. Departamento de Biofísica e Farmacologia Instituto de Biociências Universidade Estadual Paulista (UNESP) Laboratoire de Biologie et de Pharmacologie Appliquée École Normale Supérieure Paris Saclay UMR 8113 CNRS, 4 Avenue des Sciences Departamento de Física Experimental Instituto de Física Universidade de São Paulo (USP) Departamento de Biotecnologia e Bioprocessos Faculdade de Ciências Agronômicas Universidade Estadual Paulista (UNESP) Instituto de Biotecnologia Universidade Estadual Paulista (UNESP) Departamento de Biofísica e Farmacologia Instituto de Biociências Universidade Estadual Paulista (UNESP) Departamento de Biotecnologia e Bioprocessos Faculdade de Ciências Agronômicas Universidade Estadual Paulista (UNESP) Instituto de Biotecnologia Universidade Estadual Paulista (UNESP) CNPq: 302883/2017-7 CNPq: 401190/2017-0 CAPES: 88881.134154/2016-01 CAPES: 88882.183567/2007-01

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. Macromolecules adopt many conformations in solution depending on their structure and shape, which determine their dynamics and function. In this context, several studies have reported important meta-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 can affect the rate of dopamine absorption. This work presents a protocol using hybrid simulation methods to study the conformational dynamics and stability of meta-states of hDAT under different phosphorylation states. With this protocol, we explored the conformational space of hDAT, identified the meta-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. The results suggest that the phosphorylation corresponding to NP-333, where (all sites [Ser/Thr](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4380854/) from residue 2 to 62 and residue 254 to 613 are phosphorylated, except residue 333) is the one that would affect dopamine transport the most, 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 meta-states necessary for dopamine transport. In this sense, the modification of phosphorylation influences protein movements and conformational preferences, affecting the stability of meta-states and the transition between them and, therefore, the transport.

International audience; ATP-sensitive potassium (K-ATP) channels are ubiquitously expressed on the plasma membrane of cells in several organs, including the heart, pancreas, and brain, and they govern a wide range of physiological processes. In pancreatic β-cells, K-ATP channels composed of Kir6.2 and SUR1 play a key role in coupling blood glucose and insulin secretion. A tryptophan residue located at the cytosolic end of the transmembrane helix is highly conserved in eukaryote and prokaryote Kir channels. Any mutation on this amino acid causes a gain of function and neonatal diabetes mellitus. In this study, we have investigated the effect of mutation on this highly conserved residue on a KirBac channel (prokaryotic homolog of mammalian Kir6.2). We provide the crystal structure of the mutant KirBac3.1 W46R (equivalent to W68R in Kir6.2) and its conformational flexibility properties using HDX-MS. In addition, the detailed dynamical view of the mutant during the gating was investigated using the in silico method. Finally, functional assays have been performed. A comparison of important structural determinants for the gating mechanism between the wild type KirBac and the mutant W46R suggests interesting structural and dynamical clues and a mechanism of action of the mutation that leads to the gain of function.

International audience; There is increasing support for water molecules playing a role in signal propagation through G protein-coupled receptors (GPCRs). However, exploration of the hydration features of GPCRs is still in its infancy. Here, we combined site-specific labeling with unnatural amino acids to molecular dynamics to delineate how local hydration of the ghrelin receptor growth hormone secretagogue receptor (GHSR) is rearranged upon activation. We found that GHSR is characterized by a specific hydration pattern that is selectively remodeled by pharmacologically distinct ligands and by the lipid environment. This process is directly related to the concerted movements of the transmembrane domains of the receptor. These results demonstrate that the conformational dynamics of GHSR are tightly coupled to the movements of internal water molecules, further enhancing our understanding of the molecular bases of GPCR-mediated signaling.

The main protease of SARS-CoV-2 (called Mpro or 3CLpro) is essential for processing polyproteins encoded by viral RNA. Macromolecules adopt several favored conformations in solution depending on their structure and shape, determining their dynamics and function. Integrated methods combining the lowest-frequency movements obtained by Normal Mode Analysis (NMA), and the faster movements from Molecular Dynamics (MD), and data from biophysical techniques, are necessary to establish the correlation between complex structural dynamics of macromolecules and their function. In this article, we used a hybrid simulation method to sample the conformational space to characterize the structural dynamics and global motions of WT SARS-CoV-2 Mpro and 48 mutants, including several mutations that appear in P.1, B.1.1.7, B.1.351, B.1.525 and B.1.429+B.1.427 variants. Integrated Hybrid methods combining NMA and MD have been useful to study the correlation between the complex structural dynamics of macromolecules and their functioning mechanisms. Here, we applied this hybrid approach to elucidate the effects of mutation in the structural dynamics of SARS-CoV-2 Mpro, considering their flexibility, solvent accessible surface area analyses, global movements, and catalytic dyad distance. Furthermore, some mutants showed significant changes in their structural dynamics and conformation, which could lead to distinct functional properties.HighlightsSingle surface mutations lead to changes in Mpro structural dynamics.Mutants can be more stable than WT according to the structural dynamics properties.Mpromutants can present a distinct functionality in relation to the wild-type.Potential viral markers for more pathogenic or transmissible SARS-CoV-2 variants.

RalF is an Arf GEF from Legionella pneumophilia, the bacterium that causes severe pneumonia. In its crystal structure, RalF is in the autoinhibited form. A large-scale domain motion is expected to lift the autoinhibition, the mechanism of which is still unknown. Since RalF is activated in the presence of the membrane, its active structure and the structure of the RalF-Arf1 complex could not have been determined experimentally. On the simulation side, it has been proven that classical Molecular Dynamics (MD) alone is not efficient enough to map motions of such amplitude and determine the active conformation of RalF. In this article, using Molecular Dynamics with excited Normal Modes (MDeNM) combined with previous experimental findings we were able to determine the active RalF structure and the structure of the RalF-Arf1 complex in the presence of the membrane, bridging the gap between experiments and simulation.

Sulfotransferases (SULTs) are phase II drug-metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3′-phosphoadenosine 5′-phosphosulfate (PAPS) to a substrate. It has been previously suggested that a considerable shift of SULT structure caused by PAPS binding could control the capability of SULT to bind large substrates. We employed molecular dynamics (MD) simulations and the recently developed approach of MD with excited normal modes (MDeNM) to elucidate molecular mechanisms guiding the recognition of diverse substrates and inhibitors by SULT1A1. MDeNM allowed exploring an extended conformational space of PAPS-bound SULT1A1, which has not been achieved up to now by using classical MD. The generated ensembles combined with docking of 132 SULT1A1 ligands shed new light on substrate and inhibitor binding mechanisms. Unexpectedly, our simulations and analyses on binding of the substrates estradiol and fulvestrant demonstrated that large conformational changes of the PAPS-bound SULT1A1 could occur independently of the co-factor movements that could be sufficient to accommodate large substrates as fulvestrant. Such structural displacements detected by the MDeNM simulations in the presence of the co-factor suggest that a wider range of drugs could be recognized by PAPS-bound SULT1A1 and highlight the utility of including MDeNM in protein–ligand interactions studies where major rearrangements are expected.

In this study, we investigated the dynamics and functional characteristics of the KirBac3.1 S129R, a mutated bacterial potassium channel for which the inner pore-lining helix (TM2) was engineered so that the bundle crossing is trapped in an open conformation. The structure of this channel has been previously determined at high atomic resolution. We explored the dynamical characteristics of this open state channel using an in silico method MDeNM that combines molecular dynamics simulations and normal modes. We captured the global and local motions at the mutation level and compared these data with HDX-MS experiments. MDeNM provided also an estimation of the probability of the different opening states that are in agreement with our electrophysiological experiments. In the S129R mutant, the Arg129 mutation releases the two constriction points in the channel that existed in the wild type but interestingly creates another restriction point.

Inward rectifier potassium (Kir) channels play diverse and important roles in shaping action potentials in biological membranes. An increasing number of diseases are now known to be directly associated with abnormal Kir function. However, the gating of Kir still remains unknown. To increase our understanding of its gating mechanism, a dynamical view of the entire channel is essential. Here the gating activation was studied using a recent developped in silico method, MDeNM, which combines normal mode analysis and molecular dynamics simulations that showed for the very first time the importance of interrelated collective and localized conformational movements. In particular, we highlighted the role played by concerted movements of the different regions throughout the entire protein, such as the cytoplasmic and transmembrane domains and the slide helices. In addition, the HDX-MS analysis achieved in these studies provided a comprehensive and detailed view of the dynamics associated with open/closed transition of the Kir channel in coherence with the theoretical results. MDeNM gives access to the probability of the different opening states that are in agreement with our electrophysiological experiments. The investigations presented in this article are important to remedy dysfunctional channels and are of interest for designing new pharmacological compounds.

The activation process of phospholipase A

Zn-metalloproteins are a major class of targets for drug design. They constitute a demanding testing ground for polarizable molecular mechanics/dynamics aimed at extending the realm of quantum chemistry (QC) to very long-duration molecular dynamics (MD). The reliability of such procedures needs to be demonstrated upon comparing the relative stabilities of competing candidate complexes of inhibitors with the recognition site stabilized in the course of MD. This could be necessary when no information is available regarding the experimental structure of the inhibitor-protein complex. Thus, this study bears on the phosphomannose isomerase (PMI) enzyme, considered as a potential therapeutic target for the treatment of several bacterial and parasitic diseases. We consider its complexes with 5-phospho-d-arabinonohydroxamate and three analog ligands differing by the number and location of their hydroxyl groups. We evaluate the energy accuracy expectable from a polarizable molecular mechanics procedure, SIBFA. This is done by comparisons with ab initio quantum-chemistry (QC) calculations in the following cases: (a) the complexes of the four ligands in three distinct structures extracted from the entire PMI-ligand energy-minimized structures, and totaling up to 264 atoms; (b) the solvation energies of several energy-minimized complexes of each ligand with a shell of 64 water molecules; (c) the conformational energy differences of each ligand in different conformations characterized in the course of energy-minimizations; and (d) the continuum solvation energies of the ligands in different conformations. The agreements with the QC results appear convincing. On these bases, we discuss the prospects of applying the procedure to ligand-macromolecule recognition problems. © 2016 Wiley Periodicals, Inc.

Made available in DSpace on 2020-12-12T00:56:54Z (GMT). No. of bitstreams: 0 Previous issue date: 2019-12-01 NS3 is an important therapeutic target for direct-acting antiviral (DAA) drugs. However, many patients treated with DAAs have unsustained virologic response (UVR) due to the high mutation rate of HCV. The aim of this work was to shed some light on the puzzling molecular mechanisms of the virus's of patients who showed high viral loads even under treatment with DAA. Bioinformatics tools, molecular modelling analyses were employed to identify mutations associated with HCV resistance to boceprevir and possible structural features related to this phenomenon. We identified two mutations of NS3 that may be associated with HCV resistance: D168N and L153I. The substitution D168N was previously reported in the literature as related with drug failure. Additionally, we identified that its molecular resistance mechanism can be explained by the destabilization of receptor-ligand hydrogen bonds. For the L153I mutation, the resistance mechanism is different from previous models reported in the literature. The L153I substitution decreases the S139 deprotonation susceptibility, and consequently, this mutation impairs the covalent binding between the residue S139 from NS3 and the electrophilic trap on boceprevir, which can induce drug failure. These results were supported by the time course analysis of the mutations of the NS3 protease, which showed that boceprevir was designed for enzymes with an L residue at position 153; however, the sequences with I153 are predominant nowadays. The results presented here could be used to infer about resistance in others DAA, mainly protease inhibitors. Sao Paulo State University (UNESP) School of Agriculture Department of Bioprocess and Biotechnology, Avenue Universitária Sao Paulo State University (UNESP) Medical School Blood Center, Avenue Prof. Mário Rubens Guimarães Montenegro, s/n Max Planck Institut for Heart and Lung Research, Ludwigstraße 43 Sao Paulo State University (UNESP) Institute of Biosciences, Street Prof. Dr. Antônio Celso Wagner Zanin, 250 École Normale Supérieure Paris-Saclay Laboratory of Biology and Applied Pharmacology Sao Paulo State University (UNESP) School of Agriculture Department of Bioprocess and Biotechnology, Avenue Universitária Sao Paulo State University (UNESP) Medical School Blood Center, Avenue Prof. Mário Rubens Guimarães Montenegro, s/n Sao Paulo State University (UNESP) Institute of Biosciences, Street Prof. Dr. Antônio Celso Wagner Zanin, 250

Many receptors elicit signal transduction by activating multiple intracellular pathways. This transduction may be triggered by a nonspecific ligand, which simultaneously activates all receptor’s signaling paths. In addition, the binding of a biased ligand preferentially elicits one path over another, in a process called biased signaling. Identifying the functional motions related to each one of these distinct pathways have direct impact in the development of new efficient and specific drugs. Here we show how to detect functional motions considering the case of the NGF/TrkA-Ig2 complex. TrkA receptors activation mediated by NGF depends on specific structural motions that trigger neuronal growth, development and survival in nervous system. R221W mutation in the ngf gene impairs the nociceptive signaling. Here we show the wide spread structural effects of this mutation in the NGF/TrkA-Ig2 complex, leading to the deletion of collective motions important to TrkA activation of nociceptive signaling. Our results suggest that subtle changes in the interaction network within neurotrophic factors due to the point mutation are sufficient to inhibit necessary conformational changes for TrkA receptors activation. The methodological approach presented in this article based conjunctively on normal mode analysis and the experimentally observed functional alterations due to point mutations provides an essential tool for unveiling the structural changes and motions related to the disease, that in turn could be important for a drug design study.

BackgroundPlacebo-controlled clinical trials have led to concern over possible increased risk of suicide-related events in some populations exposed to antidepressants.AimsTo evaluate the risk of suicide attempts by antidepressant drug class and the presence or absence of depression.MethodA retrospective propensity-matched new-user cohort study was used to compare participants with incident depression classified by antidepressant treatment with each other and with the general population.ResultsAmong the treated group, the suicide attempt rate peaked in the month prior to diagnosis then decreased steadily over the next 6 months. Among the pharmacologically untreated group, the highest rate was seen in the second month after diagnosis. Cohorts with depression had significantly higher suicide attempt risk than the general population, but the treated group did not differ significantly from the untreated group.ConclusionsPatients on antidepressants did not have significantly higher risk compared with untreated patients. No significant differences were observed for patients treated with individual serotonin–noradrenaline reuptake inhibitors (SNRIs) or selective serotonin reuptake inhibitors (SSRIs) or by class (SSRI v. SNRI cohorts).

Made available in DSpace on 2018-12-11T17:20:01Z (GMT). No. of bitstreams: 0 Previous issue date: 2018-06-02 Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Porcine circovirus 2 (PCV2) is an icosahedral, non-enveloped, and single-stranded circular DNA virus that belongs to the family Circoviridae, genus Circovirus, and is responsible for a complex of different diseases defined as porcine circovirus diseases (PCVDs). These diseases – including postweaning multisystemic wasting syndrome (PMWS), enteric disease, respiratory disease, porcine dermatitis and nephropathy syndrome (PDNS), and reproductive failure – are responsible for large economic losses in the pig industry. After serial passages in swine testicle (ST) cells of a wild-type virus isolated from an animal with PMWS, we identified three PCV2b viruses with capsid protein (known as Cap protein) cumulative mutations, including two novel mutants. The mutant viruses were introduced into new ST cell cultures for reisolation and showed, in comparison to the wild-type PCV2b, remarkable viral replication efficiency (> 1011 DNA copies/ml) and cell death via necrosis, which were clearly related to the accretion of capsid protein mutations. The analysis of a Cap protein/capsid model showed that the mutated residues were located in solvent-accessible positions on the external PCV2b surface. Additionally, the mutated residues were found in linear epitopes and participated in pockets on the capsid surface, indicating that these residues could also be involved in antibody recognition. Taking into account the likely natural emergence of PCV2b variants, it is possible to consider that the results of this work increase knowledge of Circovirus biology and could help to prevent future serious cases of vaccine failure that could lead to heavy losses to the swine industry. Departamento de Microbiologia e Imunologia Instituto de Biociências (IB) Universidade Estadual Paulista (Unesp) Instituto de Biotecnologia (IBTEC) Universidade Estadual Paulista (Unesp) Departamento de Bioprocessos e Biotecnologia Faculdade de Ciências Agronômicas (FCA) Universidade Estadual Paulista (Unesp) Departamento de Medicina Veterinária Preventiva e Saúde Animal Universidade de São Paulo (USP) Departamento de Salud Animal Universidad de Caldas Departamento de Bioestatística Instituto de Biociências (IB) Universidade Estadual Paulista (Unesp) Laboratoire de Biologie et de Pharmacologie Appliquée ENS Cachan/Université Paris-Saclay Complexo Educacional Faculdades Metropolitanas Unidas (FMU/HOVET) São Paulo São Paulo Departamento de Microbiologia e Imunologia Instituto de Biociências (IB) Universidade Estadual Paulista (Unesp) Instituto de Biotecnologia (IBTEC) Universidade Estadual Paulista (Unesp) Departamento de Bioprocessos e Biotecnologia Faculdade de Ciências Agronômicas (FCA) Universidade Estadual Paulista (Unesp) Departamento de Bioestatística Instituto de Biociências (IB) Universidade Estadual Paulista (Unesp) FAPESP: 2006/57976-6 FAPESP: 2006/59002-9 FAPESP: 2013/14530-1

Proteins are found in solution as ensembles of conformations in dynamic equilibrium. Exploration of functional motions occurring on micro- to millisecond time scales by molecular dynamics (MD) simulations still remains computationally challenging. Alternatively, normal mode (NM) analysis is a well-suited method to characterize intrinsic slow collective motions, often associated with protein function, but the absence of anharmonic effects preclude a proper characterization of conformational distributions in a multidimensional NM space. Using both methods jointly appears to be an attractive approach that allows an extended sampling of the conformational space. In line with this view, the MDeNM (molecular dynamics with excited normal modes) method presented here consists of multiple-replica short MD simulations in which motions described by a given subset of low-frequency NMs are kinetically excited. This is achieved by adding additional atomic velocities along several randomly determined linear combinations of NM vectors, thus allowing an efficient coupling between slow and fast motions. The relatively high-energy conformations generated with MDeNM are further relaxed with standard MD simulations, enabling free energy landscapes to be determined. Two widely studied proteins were selected as examples: hen egg lysozyme and HIV-1 protease. In both cases, MDeNM provides a larger extent of sampling in a few nanoseconds, outperforming long standard MD simulations. A high degree of correlation with motions inferred from experimental sources (X-ray, EPR, and NMR) and with free energy estimations obtained by metadynamics was observed. Finally, the large sets of conformations obtained with MDeNM can be used to better characterize relevant dynamical populations, allowing for a better interpretation of experimental data such as SAXS curves and NMR spectra.

Our previous studies showed that (-)-epigallocatechin-3-gallate (EGCG) inhibits signal transducer activator of transcription 1 (STAT1) activation. Since EGCG may be a promising lead compound for new anti-STAT1 drug design, 15 synthetic catechins, characterized by the (-)-gallocatechin-3-gallate stereochemistry, were studied in the human mammary MDA-MB-231 cell line to identify the minimal structural features that preserve the anti-STAT1 activity. We demonstrate that the presence of three hydroxyl groups of B ring and one hydroxyl group in D ring is essential to preserve their inhibitory action. Moreover, a possible molecular target of these compounds in the STAT1 pathway was investigated. Our results demonstrate a direct interaction between STAT1 protein and catechins displaying anti-STAT1 activity. In particular, surface plasmon resonance (SPR) analysis and molecular modeling indicate the presence of two putative binding sites (a and b) with different affinity. Based on docking data, site-directed mutagenesis was performed, and interaction of the most active catechins with STAT1 was studied with SPR to test whether Gln518 on site a and His568 on site b could be important for the catechin-STAT1 interaction. Data indicate that site b has higher affinity for catechins than site a as the highest affinity constant disappears in the H568A-STAT1 mutant. Furthermore, Janus kinase 2 (JAK2) kinase assay data suggest that the contemporary presence in vitro of STAT1 and catechins inhibits JAK2-elicited STAT1 phosphorylation. The very tight catechin-STAT1 interaction prevents STAT1 phosphorylation and represents a novel, specific and efficient molecular mechanism for the inhibition of STAT1 activation.

We describe a new method and a software, called GANM (Genetic Algorithm with Normal Modes), that combines genetic algorithm and normal mode analysis in order to achieve protein–ligand docking simulations that take into account global and local conformational adjustments. The energetical and structural optimization during the ligand binding process is achieved with GA that uses a rotamer library for changing the side chain conformations within and around the binding site. The methodology uses normal modes to account for the global conformational changes of the protein. This software is implemented in distinct versions using different levels of approximations. In this article, we describe the methodology and discuss the results and the performance obtained with the first two versions of GANM. The protein–ligand complexes considered for the test calculations are: (1) Dihydrofolate Reductase (DHFR) and MTX as ligands, (2) HIV-1 Aspartyl Protease with ritonavir, nelfinavir and DMP323 as ligands, and (3) CDK-2 and DT-5 as ligands.

Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials & Methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results & discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility.[Formula: see text]

The antiretroviral chemotherapy helps to reduce the mortality of HIVs infected patients. However, RNA dependant virus replication has a high mutation rate. Human immunodeficiency virus Type 1 protease plays an essential role in viral replication cycle. This protein is an important target for therapy with viral protein inhibitors. There are few works using normal mode analysis to investigate this problem from the structural changes viewpoint. The investigation of protein flexibility may be important for the study of processes associated with conformational changes and state transitions. The normal mode analysis allowed us to investigate structural changes in the protease (such as flexibility) in a straightforward way and try to associate these changes with the increase of fitness for each positively selected HIV-1 mutant protease of patients treated with several protease inhibitors (saquinavir, indinavir, ritonavir, nelfinavir, lopinavir, fosamprenavir, atazanavir, darunavir, and tripanavir) in combination or separately. These positively selected mutations introduce significant flexibility in important regions such as the active site cavity and flaps. These mutations were also able to cause changes in accessible solvent area. This study showed that the majority of HIV-1 protease mutants can be grouped into two main classes of protein flexibility behavior. We presented a new approach to study structural changes caused by positively selected mutations in a pathogen protein, for instance the HIV-1 protease and their relationship with their resistance mechanism against known inhibitors. The method can be applied to any pharmaceutically relevant pathogen proteins and could be very useful to understand the effects of positively selected mutations in the context of structural changes.

Molecular dynamics (MD) simulation is an established method for studying the conformational changes that are important for protein function. Recent advances in hardware and software have allowed MD simulations over the same timescales as experiment, improving the agreement between theory and experiment to a large extent. However, running such simulations are costly, in terms of resources, storage, and trajectory analysis. There is still a place for techniques that involve short MD simulations. In order to overcome the sampling paucity of short time-scales, hybrid methods that include some form of MD simulation can exploit certain features of the system of interest, often combining experimental information in surprising ways. Here, we review some recent hybrid approaches to the simulation of proteins.

International audience; Describing biological macromolecular energetics from computer simulations can pose major challenges, and often necessitates enhanced conformational sampling. We describe the calculation of conformational free-energy profiles along carefully chosen collective coordinates: “consensus” normal modes, developed recently as robust alternatives to conventional normal modes. In an application to the HIV-1 protease, we obtain efficient sampling of significant flap opening movements governing inhibitor binding from relatively short simulations, in close correspondence with experimental results.

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.