1
|
Castelli M, Marchetti F, Osuna S, F. Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations. J Am Chem Soc 2024; 146:901-919. [PMID: 38116743 PMCID: PMC10785808 DOI: 10.1021/jacs.3c11396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 12/06/2023] [Accepted: 12/07/2023] [Indexed: 12/21/2023]
Abstract
Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.
Collapse
Affiliation(s)
- Matteo Castelli
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| | - Filippo Marchetti
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
- INSTM, via G. Giusti 9, 50121 Florence, Italy
- E4
Computer Engineering, via Martiri delle libertà 66, 42019 Scandiano (RE), Italy
| | - Sílvia Osuna
- Institut
de Química Computacional i Catàlisi (IQCC) and Departament
de Química, Universitat de Girona, Girona, Catalonia E-17071, Spain
- ICREA, Barcelona, Catalonia E-08010, Spain
| | - A. Sofia F. Oliveira
- Centre for
Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Adrian J. Mulholland
- Centre for
Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Stefano A. Serapian
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| | - Giorgio Colombo
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| |
Collapse
|
2
|
Castelli M, Magni A, Bonollo G, Pavoni S, Frigerio F, Oliveira ASF, Cinquini F, Serapian SA, Colombo G. Molecular mechanisms of chaperone-directed protein folding: Insights from atomistic simulations. Protein Sci 2023; 33:e4880. [PMID: 38145386 PMCID: PMC10895457 DOI: 10.1002/pro.4880] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 12/06/2023] [Accepted: 12/19/2023] [Indexed: 12/26/2023]
Abstract
Molecular chaperones, a family of proteins of which Hsp90 and Hsp70 are integral members, form an essential machinery to maintain healthy proteomes by controlling the folding and activation of a plethora of substrate client proteins. This is achieved through cycles in which Hsp90 and Hsp70, regulated by task-specific co-chaperones, process ATP and become part of a complex network that undergoes extensive compositional and conformational variations. Despite impressive advances in structural knowledge, the mechanisms that regulate the dynamics of functional assemblies, their response to nucleotides, and their relevance for client remodeling are still elusive. Here, we focus on the glucocorticoid receptor (GR):Hsp90:Hsp70:co-chaperone Hop client-loading and the GR:Hsp90:co-chaperone p23 client-maturation complexes, key assemblies in the folding cycle of glucocorticoid receptor (GR), a client strictly dependent upon Hsp90/Hsp70 for activity. Using a combination of molecular dynamics simulation approaches, we unveil with unprecedented detail the mechanisms that underpin function in these chaperone machineries. Specifically, we dissect the processes by which the nucleotide-encoded message is relayed to the client and how the distinct partners of the assemblies cooperate to (pre)organize partially folded GR during Loading and Maturation. We show how different ligand states determine distinct dynamic profiles for the functional interfaces defining the interactions in the complexes and modulate their overall flexibility to facilitate progress along the chaperone cycle. Finally, we also show that the GR regions engaged by the chaperone machinery display peculiar energetic signatures in the folded state, which enhance the probability of partial unfolding fluctuations. From these results, we propose a model where a dynamic cross-talk emerges between the chaperone dynamics states and remodeling of client-interacting regions. This factor, coupled to the highly dynamic nature of the assemblies and the conformational heterogeneity of their interactions, provides the basis for regulating the functions of distinct assemblies during the chaperoning cycle.
Collapse
Affiliation(s)
| | - Andrea Magni
- Dipartimento di Chimica, Università di Pavia, Pavia, Italy
| | | | - Silvia Pavoni
- Department of Physical Chemistry, R&D Eni SpA, San Donato Milanese, Italy
| | - Francesco Frigerio
- Department of Physical Chemistry, R&D Eni SpA, San Donato Milanese, Italy
| | - A Sofia F Oliveira
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| | - Fabrizio Cinquini
- Upstream & Technical Services - TECS/STES - Eni Spa, San Donato Milanese, Italy
| | | | | |
Collapse
|
3
|
Daniels DE, Ferrer-Vicens I, Hawksworth J, Andrienko TN, Finnie EM, Bretherton NS, Ferguson DCJ, Oliveira ASF, Szeto JYA, Wilson MC, Brewin JN, Frayne J. Human cellular model systems of β-thalassemia enable in-depth analysis of disease phenotype. Nat Commun 2023; 14:6260. [PMID: 37803026 PMCID: PMC10558456 DOI: 10.1038/s41467-023-41961-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2022] [Accepted: 09/26/2023] [Indexed: 10/08/2023] Open
Abstract
β-thalassemia is a prevalent genetic disorder causing severe anemia due to defective erythropoiesis, with few treatment options. Studying the underlying molecular defects is impeded by paucity of suitable patient material. In this study we create human disease cellular model systems for β-thalassemia by gene editing the erythroid line BEL-A, which accurately recapitulate the phenotype of patient erythroid cells. We also develop a high throughput compatible fluorometric-based assay for evaluating severity of disease phenotype and utilize the assay to demonstrate that the lines respond appropriately to verified reagents. We next use the lines to perform extensive analysis of the altered molecular mechanisms in β-thalassemia erythroid cells, revealing upregulation of a wide range of biological pathways and processes along with potential novel targets for therapeutic investigation. Overall, the lines provide a sustainable supply of disease cells as research tools for identifying therapeutic targets and as screening platforms for new drugs and reagents.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | - Jenn-Yeu A Szeto
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK
| | | | - John N Brewin
- Haematology Department, King's college Hospital NHS Foundation, London, SE5 9RS, UK
- Red Cell Biology Group, Kings College London, London, SE5 9NU, UK
| | - Jan Frayne
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK.
| |
Collapse
|
4
|
Oliveira ASF, Shoemark DK, Davidson AD, Berger I, Schaffitzel C, Mulholland AJ. SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid. J Mol Cell Biol 2023; 15:mjad021. [PMID: 36990513 PMCID: PMC10563148 DOI: 10.1093/jmcb/mjad021] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 11/07/2022] [Accepted: 03/28/2023] [Indexed: 03/31/2023] Open
Abstract
The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol BS8 1TS, UK
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | | | - Andrew D Davidson
- School of Cellular and Molecular Medicine, University of Bristol, University Walk, Bristol BS8 1TD, UK
| | - Imre Berger
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
- School of Chemistry, Max Planck Bristol Centre for Minimal Biology, Bristol BS8 1TS, UK
| | | | - Adrian J Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol BS8 1TS, UK
| |
Collapse
|
5
|
Hutchins GH, Noble CEM, Bunzel HA, Williams C, Dubiel P, Yadav SKN, Molinaro PM, Barringer R, Blackburn H, Hardy BJ, Parnell AE, Landau C, Race PR, Oliver TAA, Koder RL, Crump MP, Schaffitzel C, Oliveira ASF, Mulholland AJ, Anderson JLR. An expandable, modular de novo protein platform for precision redox engineering. Proc Natl Acad Sci U S A 2023; 120:e2306046120. [PMID: 37487099 PMCID: PMC10400981 DOI: 10.1073/pnas.2306046120] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Accepted: 06/20/2023] [Indexed: 07/26/2023] Open
Abstract
The electron-conducting circuitry of life represents an as-yet untapped resource of exquisite, nanoscale biomolecular engineering. Here, we report the characterization and structure of a de novo diheme "maquette" protein, 4D2, which we subsequently use to create an expanded, modular platform for heme protein design. A well-folded monoheme variant was created by computational redesign, which was then utilized for the experimental validation of continuum electrostatic redox potential calculations. This demonstrates how fundamental biophysical properties can be predicted and fine-tuned. 4D2 was then extended into a tetraheme helical bundle, representing a 7 nm molecular wire. Despite a molecular weight of only 24 kDa, electron cryomicroscopy illustrated a remarkable level of detail, indicating the positioning of the secondary structure and the heme cofactors. This robust, expressible, highly thermostable and readily designable modular platform presents a valuable resource for redox protein design and the future construction of artificial electron-conducting circuitry.
Collapse
Affiliation(s)
- George H. Hutchins
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - Claire E. M. Noble
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, BristolBS8 1TQ, United Kingdom
| | - H. Adrian Bunzel
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | | | - Paulina Dubiel
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - Sathish K. N. Yadav
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - Paul M. Molinaro
- Department of Physics, The City College of New York, New York, NY10031
- Graduate Programs of Physics, Biology, Chemistry and Biochemistry, The Graduate Center of The City University of New York, New York, NY10016
| | - Rob Barringer
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - Hector Blackburn
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - Benjamin J. Hardy
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - Alice E. Parnell
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, BristolBS8 1TQ, United Kingdom
| | - Charles Landau
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - Paul R. Race
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, BristolBS8 1TQ, United Kingdom
| | | | - Ronald L. Koder
- Department of Physics, The City College of New York, New York, NY10031
- Graduate Programs of Physics, Biology, Chemistry and Biochemistry, The Graduate Center of The City University of New York, New York, NY10016
| | - Matthew P. Crump
- School of Chemistry, University of Bristol, BristolBS8 1TS, United Kingdom
| | - Christiane Schaffitzel
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
| | - A. Sofia F. Oliveira
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
- School of Chemistry, University of Bristol, BristolBS8 1TS, United Kingdom
| | - Adrian J. Mulholland
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, BristolBS8 1TQ, United Kingdom
- School of Chemistry, University of Bristol, BristolBS8 1TS, United Kingdom
| | - J. L. Ross Anderson
- School of Biochemistry, University of Bristol, University Walk, BristolBS8 1TD, United Kingdom
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, BristolBS8 1TQ, United Kingdom
| |
Collapse
|
6
|
Freeman SL, Oliveira ASF, Gallio AE, Rosa A, Simitakou MK, Arthur CJ, Mulholland AJ, Cherepanov P, Raven EL. Heme binding to the SARS-CoV-2 spike glycoprotein. J Biol Chem 2023; 299:105014. [PMID: 37414149 PMCID: PMC10416065 DOI: 10.1016/j.jbc.2023.105014] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Revised: 06/12/2023] [Accepted: 07/01/2023] [Indexed: 07/08/2023] Open
Abstract
The target for humoral immunity, SARS-CoV-2 spike glycoprotein, has become the focus of vaccine research and development. Previous work demonstrated that the N-terminal domain (NTD) of SARS-CoV-2 spike binds biliverdin-a product of heme catabolism-causing a strong allosteric effect on the activity of a subset of neutralizing antibodies. Herein, we show that the spike glycoprotein is also able to bind heme (KD = 0.5 ± 0.2 μM). Molecular modeling indicated that the heme group fits well within the same pocket on the SARS-CoV-2 spike NTD. Lined by aromatic and hydrophobic residues (W104, V126, I129, F192, F194, I203, and L226), the pocket provides a suitable environment to stabilize the hydrophobic heme. Mutagenesis of N121 has a substantive effect on heme binding (KD = 3000 ± 220 μM), confirming the pocket as a major heme binding location of the viral glycoprotein. Coupled oxidation experiments in the presence of ascorbate indicated that the SARS-CoV-2 glycoprotein can catalyze the slow conversion of heme to biliverdin. The heme trapping and oxidation activities of the spike may allow the virus to reduce levels of free heme during infection to facilitate evasion of the adaptive and innate immunity.
Collapse
Affiliation(s)
- Samuel L Freeman
- School of Chemistry, Cantock's Close, University of Bristol, Bristol, United Kingdom
| | - A Sofia F Oliveira
- School of Chemistry, Cantock's Close, University of Bristol, Bristol, United Kingdom
| | - Andrea E Gallio
- School of Chemistry, Cantock's Close, University of Bristol, Bristol, United Kingdom
| | - Annachiara Rosa
- Chromatin Structure and Mobile DNA Laboratory, The Francis Crick Institute, London, United Kingdom
| | - Maria K Simitakou
- Chromatin Structure and Mobile DNA Laboratory, The Francis Crick Institute, London, United Kingdom
| | - Christopher J Arthur
- School of Chemistry, Cantock's Close, University of Bristol, Bristol, United Kingdom
| | - Adrian J Mulholland
- School of Chemistry, Cantock's Close, University of Bristol, Bristol, United Kingdom
| | - Peter Cherepanov
- Chromatin Structure and Mobile DNA Laboratory, The Francis Crick Institute, London, United Kingdom; Department of Infectious Disease, St-Mary's Campus, Imperial College London, United Kingdom.
| | - Emma L Raven
- School of Chemistry, Cantock's Close, University of Bristol, Bristol, United Kingdom.
| |
Collapse
|
7
|
Chan HT, Oliveira ASF, Schofield CJ, Mulholland AJ, Duarte F. Dynamical Nonequilibrium Molecular Dynamics Simulations Identify Allosteric Sites and Positions Associated with Drug Resistance in the SARS-CoV-2 Main Protease. JACS Au 2023; 3:1767-1774. [PMID: 37384148 PMCID: PMC10262681 DOI: 10.1021/jacsau.3c00185] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Revised: 05/19/2023] [Accepted: 05/19/2023] [Indexed: 06/30/2023]
Abstract
The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.
Collapse
Affiliation(s)
- H. T.
Henry Chan
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| | - A. Sofia F. Oliveira
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, UK
- School
of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | - Christopher J. Schofield
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| | - Adrian J. Mulholland
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, UK
| | - Fernanda Duarte
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| |
Collapse
|
8
|
Elseginy SA, Oliveira ASF, Shoemark DK, Sessions RB. Identification and validation of novel microtubule suppressors with an imidazopyridine scaffold through structure-based virtual screening and docking. RSC Med Chem 2022; 13:929-943. [PMID: 36092142 PMCID: PMC9384815 DOI: 10.1039/d1md00392e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Accepted: 05/13/2022] [Indexed: 11/29/2022] Open
Abstract
Targeting the colchicine binding site of α/β tubulin microtubules can lead to suppression of microtubule dynamics, cell cycle arrest and apoptosis. Therefore, the development of microtubule (MT) inhibitors is considered a promising route to anticancer agents. Our approach to identify novel scaffolds as MT inhibitors depends on a 3D-structure-based pharmacophore approach and docking using three programs MOE, Autodock and BUDE (Bristol University Docking Engine) to screen a library of virtual compounds. From this work we identified the compound 7-(3-hydroxy-4-methoxy-phenyl)-3-(3-trifluoromethyl-phenyl)-6,7-dihydro-3H-imidazo[4,5-b]pyridin-5-ol (6) as a novel inhibitor scaffold. This compound inhibited several types of cancer cell proliferation at low micromolar concentrations with low toxicity. Compound 6 caused cell cycle arrest in the G2/M phase and blocked tubulin polymerization at low micromolar concentration (IC50 = 6.1 ±0.1 μM), inducing apoptosis via activation of caspase 9, increasing the level of the pro-apoptotic protein Bax and decreasing the level of the anti-apoptotic protein Bcl2. In summary, our approach identified a lead compound with potential antimitotic and antiproliferative activity.
Collapse
Affiliation(s)
- Samia A. Elseginy
- School of Biochemistry, University of BristolBiomedical Sciences Building, University WalkBristolBS8 1TDUK+44 (0)1173312146,Green Chemistry Department, Chemical Industries Research Division, National Research Centre12622Egypt
| | - A. Sofia F. Oliveira
- School of Biochemistry, University of BristolBiomedical Sciences Building, University WalkBristolBS8 1TDUK+44 (0)1173312146,School of Chemistry, University of BristolBristolBS8 1TSUK
| | - Deborah K. Shoemark
- School of Biochemistry, University of BristolBiomedical Sciences Building, University WalkBristolBS8 1TDUK+44 (0)1173312146
| | - Richard B. Sessions
- School of Biochemistry, University of BristolBiomedical Sciences Building, University WalkBristolBS8 1TDUK+44 (0)1173312146
| |
Collapse
|
9
|
Allen WJ, Corey RA, Watkins DW, Oliveira ASF, Hards K, Cook GM, Collinson I. Rate-limiting transport of positively charged arginine residues through the Sec-machinery is integral to the mechanism of protein secretion. eLife 2022; 11:77586. [PMID: 35486093 PMCID: PMC9110029 DOI: 10.7554/elife.77586] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Accepted: 04/29/2022] [Indexed: 11/24/2022] Open
Abstract
Transport of proteins across and into membranes is a fundamental biological process with the vast majority being conducted by the ubiquitous Sec machinery. In bacteria, this is usually achieved when the SecY-complex engages the cytosolic ATPase SecA (secretion) or translating ribosomes (insertion). Great strides have been made towards understanding the mechanism of protein translocation. Yet, important questions remain – notably, the nature of the individual steps that constitute transport, and how the proton-motive force (PMF) across the plasma membrane contributes. Here, we apply a recently developed high-resolution protein transport assay to explore these questions. We find that pre-protein transport is limited primarily by the diffusion of arginine residues across the membrane, particularly in the context of bulky hydrophobic sequences. This specific effect of arginine, caused by its positive charge, is mitigated for lysine which can be deprotonated and transported across the membrane in its neutral form. These observations have interesting implications for the mechanism of protein secretion, suggesting a simple mechanism through which the PMF can aid transport by enabling a 'proton ratchet', wherein re-protonation of exiting lysine residues prevents channel re-entry, biasing transport in the outward direction.
Collapse
Affiliation(s)
- William J Allen
- School of Biochemistry, University of Bristol, Bristol, United Kingdom
| | - Robin A Corey
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - Daniel W Watkins
- School of Biochemistry, University of Bristol, Bristol, United Kingdom
| | | | - Kiel Hards
- Department of Microbiology and Immunology, University of Otago, Dunedin, New Zealand
| | - Gregory M Cook
- Department of Microbiology and Immunology, University of Otago, Duneding, New Zealand
| | - Ian Collinson
- School of Biochemistry, University of Bristol, Bristol, United Kingdom
| |
Collapse
|
10
|
Sofia F. Oliveira A, Shoemark DK, Avila Ibarra A, Davidson AD, Berger I, Schaffitzel C, Mulholland AJ. The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour. Comput Struct Biotechnol J 2021; 20:139-147. [PMID: 34934478 PMCID: PMC8670790 DOI: 10.1016/j.csbj.2021.12.011] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 12/08/2021] [Accepted: 12/08/2021] [Indexed: 12/12/2022] Open
Abstract
The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.
Collapse
Key Words
- ACE2, angiotensin-converting 2 enzyme
- CD, connector domain
- CH, central helix
- FA, fatty acid
- FP, fusion peptide
- FPPR, fusion-peptide proximal region
- HR1, heptad repeat 1
- LA, Linoleic acid
- MD, Molecular dynamics
- MERS, middle east respiratory syndrome
- NTD, N-terminal domain
- RBD, receptor binding domain
- RBM, receptor-binding motif
- RMB, receptor binding motif
- SARS, severe acute respiratory syndrome
- SARS-CoV-2, severe acute respiratory syndrome 2
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
| | - Deborah K. Shoemark
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | - Amaurys Avila Ibarra
- Research Software Engineering, Advanced Computing Research Centre, University of Bristol, Bristol BS1 5QD, UK
| | - Andrew D. Davidson
- School of Cellular and Molecular Medicine, University of Bristol, University Walk, Bristol BS8 1TD, UK
| | - Imre Berger
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
- Max Planck Bristol Centre for Minimal Biology, Cantock's Close, Bristol BS8 1TS, UK
| | - Christiane Schaffitzel
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| |
Collapse
|
11
|
Abreu B, Cruz C, Oliveira ASF, Soares CM. ATP hydrolysis and nucleotide exit enhance maltose translocation in the MalFGK 2E importer. Sci Rep 2021; 11:10591. [PMID: 34012037 PMCID: PMC8134467 DOI: 10.1038/s41598-021-89556-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2020] [Accepted: 04/27/2021] [Indexed: 02/03/2023] Open
Abstract
ATP binding cassette (ABC) transporters employ ATP hydrolysis to harness substrate translocation across membranes. The Escherichia coli MalFGK2E maltose importer is an example of a type I ABC importer and a model system for this class of ABC transporters. The MalFGK2E importer is responsible for the intake of malto-oligossacharides in E.coli. Despite being extensively studied, little is known about the effect of ATP hydrolysis and nucleotide exit on substrate transport. In this work, we studied this phenomenon using extensive molecular dynamics simulations (MD) along with potential of mean force calculations of maltose transport across the pore, in the pre-hydrolysis, post-hydrolysis and nucleotide-free states. We concluded that ATP hydrolysis and nucleotide exit trigger conformational changes that result in the decrease of energetic barriers to maltose translocation towards the cytoplasm, with a concomitant increase of the energy barrier in the periplasmic side of the pore, contributing for the irreversibility of the process. We also identified key residues that aid in positioning and orientation of maltose, as well as a novel binding pocket for maltose in MalG. Additionally, ATP hydrolysis leads to conformations similar to the nucleotide-free state. This study shows the contribution of ATP hydrolysis and nucleotide exit in the transport cycle, shedding light on ABC type I importer mechanisms.
Collapse
Affiliation(s)
- Bárbara Abreu
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - Carlos Cruz
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - A. Sofia F. Oliveira
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal ,grid.5337.20000 0004 1936 7603School of Biochemistry and Centre for Computational Chemistry, University of Bristol, Bristol, UK
| | - Cláudio M. Soares
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| |
Collapse
|
12
|
Oliveira ASF, Ibarra AA, Bermudez I, Casalino L, Gaieb Z, Shoemark DK, Gallagher T, Sessions RB, Amaro RE, Mulholland AJ. A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors. Biophys J 2021; 120:983-993. [PMID: 33609494 PMCID: PMC7889469 DOI: 10.1016/j.bpj.2021.01.037] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 01/11/2021] [Accepted: 01/13/2021] [Indexed: 01/08/2023] Open
Abstract
Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4β2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom; Bristol Synthetic Biology Centre, BrisSynBio, Bristol, United Kingdom
| | - Amaurys Avila Ibarra
- Research Software Engineering, Advanced Computing Research Centre, University of Bristol, Bristol, United Kingdom
| | - Isabel Bermudez
- Department of Biological and Medical Sciences, Oxford Brookes University, Oxford, United Kingdom
| | - Lorenzo Casalino
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Deborah K Shoemark
- School of Biochemistry, University of Bristol, Bristol, United Kingdom; Bristol Synthetic Biology Centre, BrisSynBio, Bristol, United Kingdom
| | - Timothy Gallagher
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| | | | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom.
| |
Collapse
|
13
|
Oliveira ASF, Ciccotti G, Haider S, Mulholland AJ. Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems. Eur Phys J B 2021; 94:144. [PMID: 34720710 PMCID: PMC8549953 DOI: 10.1140/epjb/s10051-021-00157-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 07/05/2021] [Indexed: 05/03/2023]
Abstract
ABSTRACT A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in β -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol, BS8 1TQ UK
| | - Giovanni Ciccotti
- Institute for Applied Computing “Mauro Picone” (IAC), CNR, Via dei Taurini 19, 00185 Rome, Italy
- School of Physics, University College of Dublin, UCD-Belfield, Dublin 4, Ireland
- Università di Roma La Sapienza, Ple. A. Moro 5, 00185 Rome, Italy
| | - Shozeb Haider
- School of Pharmacy, University College London, London, WC1N 1AX UK
| | - Adrian J. Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
| |
Collapse
|
14
|
Oliveira ASF, Ibarra AA, Bermudez I, Casalino L, Gaieb Z, Shoemark DK, Gallagher T, Sessions RB, Amaro RE, Mulholland AJ. Simulations support the interaction of the SARS-CoV-2 spike protein with nicotinic acetylcholine receptors. bioRxiv 2020:2020.07.16.206680. [PMID: 32743575 PMCID: PMC7386492 DOI: 10.1101/2020.07.16.206680] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Changeux et al. recently suggested that the SARS-CoV-2 spike (S) protein may interact with nicotinic acetylcholine receptors (nAChRs). Such interactions may be involved in pathology and infectivity. Here, we use molecular simulations of validated atomically detailed structures of nAChRs, and of the S protein, to investigate this 'nicotinic hypothesis'. We examine the binding of the Y674-R685 loop of the S protein to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results indicate that Y674-R685 has affinity for nAChRs and the region responsible for binding contains the PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. In particular, R682 has a key role in the stabilisation of the complexes as it forms interactions with loops A, B and C in the receptor's binding pocket. The conformational behaviour of the bound Y674-R685 region is highly dependent on the receptor subtype, adopting extended conformations in the α4β2 and α7 complexes and more compact ones when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket where it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1 and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of the simulations of the full-length S protein show that the Y674-R685 region is accessible for binding, and suggest a potential binding orientation of the S protein with nAChRs.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Amaurys Avila Ibarra
- Research Software Engineering, Advanced Computing Research Centre, University of Bristol, Bristol BS1 5QD, UK
| | - Isabel Bermudez
- Department of Biological and Medical Sciences, Oxford Brookes University, Oxford OX30BP, UK
| | - Lorenzo Casalino
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, 92093 USA
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, 92093 USA
| | | | - Timothy Gallagher
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | | | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, 92093 USA
| | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| |
Collapse
|
15
|
Lythell E, Suardíaz R, Hinchliffe P, Hanpaibool C, Visitsatthawong S, Oliveira ASF, Lang EJM, Surawatanawong P, Lee VS, Rungrotmongkol T, Fey N, Spencer J, Mulholland AJ. Resistance to the “last resort” antibiotic colistin: a single-zinc mechanism for phosphointermediate formation in MCR enzymes. Chem Commun (Camb) 2020; 56:6874-6877. [DOI: 10.1039/d0cc02520h] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Simulations show the mono-zinc form of MCR to be stable and competent for covalent phospho(ethanolamine) intermediate formation.
Collapse
|
16
|
Abreu B, Lopes EF, Oliveira ASF, Soares CM. F508del disturbs the dynamics of the nucleotide binding domains of CFTR before and after ATP hydrolysis. Proteins 2019; 88:113-126. [PMID: 31298435 DOI: 10.1002/prot.25776] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2019] [Revised: 06/17/2019] [Accepted: 07/06/2019] [Indexed: 12/20/2022]
Abstract
The cystic fibrosis transmembrane conductance regulator (CFTR) channel is an ion channel responsible for chloride transport in epithelia and it belongs to the class of ABC transporters. The deletion of phenylalanine 508 (F508del) in CFTR is the most common mutation responsible for cystic fibrosis. Little is known about the effect of the mutation in the isolated nucleotide binding domains (NBDs), on dimer dynamics, ATP hydrolysis and even on nucleotide binding. Using molecular dynamics simulations of the human CFTR NBD dimer, we showed that F508del increases, in the prehydrolysis state, the inter-motif distance in both ATP binding sites (ABP) when ATP is bound. Additionally, a decrease in the number of catalytically competent conformations was observed in the presence of F508del. We used the subtraction technique to study the first 300 ps after ATP hydrolysis in the catalytic competent site and found that the F508del dimer evidences lower conformational changes than the wild type. Using longer simulation times, the magnitude of the conformational changes in both forms increases. Nonetheless, the F508del dimer shows lower C-α RMS values in comparison to the wild-type, on the F508del loop, on the residues surrounding the catalytic site and the portion of NBD2 adjacent to ABP1. These results provide evidence that F508del interferes with the NBD dynamics before and after ATP hydrolysis. These findings shed a new light on the effect of F508del on NBD dynamics and reveal a novel mechanism for the influence of F508del on CFTR.
Collapse
Affiliation(s)
- Bárbara Abreu
- Protein Modelling Lab, ITQB-NOVA, Universidade Nova de Lisboa, Oeiras, Portugal
| | - Emanuel F Lopes
- Protein Modelling Lab, ITQB-NOVA, Universidade Nova de Lisboa, Oeiras, Portugal
| | - A S F Oliveira
- Protein Modelling Lab, ITQB-NOVA, Universidade Nova de Lisboa, Oeiras, Portugal.,School of Biochemistry & Center for Computational Chemistry, University of Bristol, Bristol, UK
| | - Cláudio M Soares
- Protein Modelling Lab, ITQB-NOVA, Universidade Nova de Lisboa, Oeiras, Portugal
| |
Collapse
|
17
|
Oliveira ASF, Shoemark DK, Campello HR, Wonnacott S, Gallagher T, Sessions RB, Mulholland AJ. Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations. Structure 2019; 27:1171-1183.e3. [PMID: 31130483 DOI: 10.1016/j.str.2019.04.008] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 01/28/2019] [Accepted: 04/10/2019] [Indexed: 02/02/2023]
Abstract
Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- School of Biochemistry, University of Bristol, Bristol BS8 1DT, UK; Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | | | - Hugo Rego Campello
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Susan Wonnacott
- Department of Biology and Biochemistry, University of Bath, Bath BA2 7AY, UK
| | - Timothy Gallagher
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | | | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK.
| |
Collapse
|
18
|
Abstract
Type-II NADH:quinone oxidoreductases (NDH-2s) are membrane proteins involved in respiratory chains and the only enzymes with NADH:quinone oxidoreductase activity expressed in Staphylococcus aureus (S. aureus), one of the most common causes of clinical infections. NDH-2s are members of the two-Dinucleotide Binding Domains Flavoprotein (tDBDF) superfamily, having a flavin adenine dinucleotide, FAD, as prosthetic group and NAD(P)H as substrate. The establishment of a Charge-Transfer Complex (CTC) between the isoalloxazine ring of the reduced flavin and the nicotinamide ring of NAD+ in NDH-2 was described, and in this work we explored its role in the kinetic mechanism using different electron donors and electron acceptors. We observed that CTC slows down the rate of the second half reaction (quinone reduction) and determines the effect of HQNO, an inhibitor. Also, protonation equilibrium simulations clearly indicate that the protonation probability of an important residue for proton transfer to the active site (D302) is influenced by the presence of the CTC. We propose that CTC is critical for the overall mechanism of NDH-2 and possibly relevant to keep a low quinol/quinone ratio and avoid excessive ROS production in vivo.
Collapse
Affiliation(s)
- Filipa V Sena
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - Filipe M Sousa
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - A Sofia F Oliveira
- School of Biochemistry, University of Bristol, University Walk, Bristol BS8 1TD, UK
| | - Cláudio M Soares
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - Teresa Catarino
- Departamento de Química, Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, 2829-516 Caparica, Portugal
| | - Manuela M Pereira
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal; University of Lisbon, Faculty of Sciences, BioISI - Biosystems & Integrative Sciences Institute, Campo Grande, C8, 1749-016 Lisboa, Portugal.
| |
Collapse
|
19
|
Frazao C, Pereira L, Saraiva IH, Oliveira ASF, Soares CM, Louro RO. Structure of FoxE, the Rhodobacter ferrooxidans SW2 putative iron oxidase. Acta Crystallogr A Found Adv 2017. [DOI: 10.1107/s2053273317083991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
|
20
|
Marreiros BC, Sena FV, Sousa FM, Oliveira ASF, Soares CM, Batista AP, Pereira MM. Structural and Functional insights into the catalytic mechanism of the Type II NADH:quinone oxidoreductase family. Sci Rep 2017; 7:42303. [PMID: 28181562 PMCID: PMC5299459 DOI: 10.1038/srep42303] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2016] [Accepted: 01/05/2017] [Indexed: 12/29/2022] Open
Abstract
Type II NADH:quinone oxidoreductases (NDH-2s) are membrane proteins involved in respiratory chains. These proteins contribute indirectly to the establishment of the transmembrane difference of electrochemical potential by catalyzing the reduction of quinone by oxidation of NAD(P)H. NDH-2s are widespread enzymes being present in the three domains of life. In this work, we explored the catalytic mechanism of NDH-2 by investigating the common elements of all NDH-2s, based on the rationale that conservation of such elements reflects their structural/functional importance. We observed conserved sequence motifs and structural elements among 1762 NDH-2s. We identified two proton pathways possibly involved in the protonation of the quinone. Our results led us to propose the first catalytic mechanism for NDH-2 family, in which a conserved glutamate residue, E172 (in NDH-2 from Staphylococcus aureus) plays a key role in proton transfer to the quinone pocket. This catalytic mechanism may also be extended to the other members of the two-Dinucleotide Binding Domains Flavoprotein (tDBDF) superfamily, such as sulfide:quinone oxidoreductases.
Collapse
Affiliation(s)
- Bruno C Marreiros
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - Filipa V Sena
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - Filipe M Sousa
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - A Sofia F Oliveira
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - Cláudio M Soares
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - Ana P Batista
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| | - Manuela M Pereira
- Instituto de Tecnologia Química e Biológica - António Xavier, Universidade Nova de Lisboa, Av. da Republica EAN, 2780-157 Oeiras, Portugal
| |
Collapse
|
21
|
Magalhães PR, Oliveira ASF, Campos SRR, Soares CM, Baptista AM. Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study. J Chem Inf Model 2017; 57:256-266. [PMID: 28095694 DOI: 10.1021/acs.jcim.6b00575] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Cytochrome c oxidase (CcO) couples the reduction of dioxygen to water with transmembrane proton pumping, which leads to the generation of an electrochemical gradient. In this study we analyze how one of the components of the electrochemical gradient, the difference in pH across the membrane, or ΔpH, influences the protonation states of residues in CcO. We modified our continuum electrostatics/Monte Carlo (CE/MC) method in order to include the ΔpH and applied it to the study of CcO, in what is, to our best knowledge, the first CE/MC study of CcO in the presence of a pH gradient. The inclusion of a transmembrane pH gradient allows for the identification of residues whose titration behavior depends on the pH on both sides of the membrane. Among the several residues with unusual titration profiles, three are well-known key residues in the proton transfer process of CcO: E286I, Y288I, and K362I. All three residues have been previously identified as being critical for the catalytic or proton pumping functions of CcO. Our results suggest that when the pH gradient increases, these residues may be part of a regulatory mechanism to stem the proton flow.
Collapse
Affiliation(s)
- Pedro R Magalhães
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa , 2781-901 Oeiras, Portugal
| | - A Sofia F Oliveira
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa , 2781-901 Oeiras, Portugal
| | - Sara R R Campos
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa , 2781-901 Oeiras, Portugal
| | - Cláudio M Soares
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa , 2781-901 Oeiras, Portugal
| | - António M Baptista
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa , 2781-901 Oeiras, Portugal
| |
Collapse
|
22
|
Oliveira ASF, Campos SRR, Baptista AM, Soares CM. Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations. Biochim Biophys Acta 2016; 1857:759-71. [PMID: 27033303 DOI: 10.1016/j.bbabio.2016.03.024] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/17/2015] [Revised: 03/18/2016] [Accepted: 03/23/2016] [Indexed: 12/11/2022]
Abstract
Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- ITQB, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - Sara R R Campos
- ITQB, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - António M Baptista
- ITQB, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal.
| | - Cláudio M Soares
- ITQB, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal.
| |
Collapse
|
23
|
Oliveira ASF, Damas JM, Baptista AM, Soares CM. Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site. PLoS Comput Biol 2014; 10:e1004010. [PMID: 25474152 PMCID: PMC4256069 DOI: 10.1371/journal.pcbi.1004010] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2014] [Accepted: 10/29/2014] [Indexed: 12/04/2022] Open
Abstract
Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O2) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O2 to travel from the solvent/membrane to the heme a3-CuB binuclear center (BNC). Moreover, the identification of all possible O2 channels as well as the atomic details of O2 diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O2 distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O2 channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O2 molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O2 inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O2 to reach the BNC. Both alternative O2 channels start in the membrane spanning region and terminate close to Y288I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O2 to reach the BNC are the alternative channels, rather than the X-ray inferred pathway. Cytochrome c oxidases (Ccoxs), the terminal enzymes of the respiratory electron transport chain in eukaryotes and many prokaryotes, are key enzymes in aerobic respiration. These proteins couple the reduction of molecular dioxygen to water with the creation of a transmembrane electrochemical proton gradient. Over the last decades, most of the Ccoxs research focused on the mechanisms and energetics of reduction and/or proton pumping, and little emphasis has been given to the pathways used by dioxygen to reach the binuclear center, where dioxygen reduction takes place. In particular, the existence and the characteristics of the channel(s) used by O2 to travel from the solvent/membrane to the binuclear site are still unclear. In this work, we combine all-atom molecular dynamics simulations and implicit ligand sampling calculations in order to identify and characterize the O2 delivery channels in the Ccox from Rhodobacter sphaeroides. Altogether, our results suggest that, in this Ccox, O2 can diffuse via three well-defined channels that start in membrane region (where O2 solubility is higher than in the water). One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O2 to reach the binuclear center.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- ITQB - Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - João M. Damas
- ITQB - Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - António M. Baptista
- ITQB - Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - Cláudio M. Soares
- ITQB - Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
- * E-mail:
| |
Collapse
|
24
|
Campos E, Asandei A, McVey CE, Dias JC, Oliveira ASF, Soares CM, Luchian T, Astier Y. The role of Lys147 in the interaction between MPSA-gold nanoparticles and the α-hemolysin nanopore. Langmuir 2012; 28:15643-15650. [PMID: 23046444 DOI: 10.1021/la302613g] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Single channel recordings were used to determine the effect of direct electrostatic interactions between sulfonate-coated gold nanoparticles and the constriction of the Staphylococcus aureus α-hemolysin protein channel on the ionic current amplitude. We provide evidence that Lys147 of α-hemolysin can interact with the sulfonate groups at the nanoparticle surface, and these interactions can reversibly block 100% of the residual ionic current. Lys147 is normally involved in a salt bridge with Glu111. The capture of a nanoparticle leads to a partial current block at neutral pH values, but protonation of Glu111 at pH 2.8 results in a full current block when the nanoparticle is captured. At pH 2.8, we suggest that Lys147 is free to engage in electrostatic interactions with sulfonates at the nanoparticle surface. To verify our results, we engineered a mutation in the α-hemolysin protein, where Glu111 is substituted by Ala (E111A), thus removing Glu111-Lys147 interactions and facilitating Lys147-sulfonate electrostatic interactions. This mutation leads to a 100% current block at pH 2.8 and a 92% block at pH 8.0, showing that electrostatic interactions are formed between the nanopore and the nanoparticle surface. Besides demonstrating the effect of electrostatic interactions on cross channel ionic current, this work offers a novel approach to controlling open and closed states of the α-hemolysin nanopore as a function of external gears.
Collapse
Affiliation(s)
- Elisa Campos
- Single Molecule Processes Laboratory, Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Oeiras, Portugal
| | | | | | | | | | | | | | | |
Collapse
|
25
|
Damas JM, Oliveira ASF, Baptista AM, Soares CM. Structural consequences of ATP hydrolysis on the ABC transporter NBD dimer: molecular dynamics studies of HlyB. Protein Sci 2011; 20:1220-30. [PMID: 21563222 DOI: 10.1002/pro.650] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
ABC transporters are a large and important family of membrane proteins involved in substrate transport across the membrane. The transported substrates are quite diverse, ranging from monatomic ions to large biomolecules. Consequently, some ABC transporters are involved in biomedically relevant situations, from genetic diseases to multidrug resistance. The most conserved domains in ABC transporters are the nucleotide binding domains (NBDs), which form a dimer responsible for the binding and hydrolysis of ATP, concomitantly with substrate translocation. To elucidate how ATP hydrolysis structurally affects the NBD dimer, and consequently the transporter, we performed a molecular dynamics study on the NBD dimer of the HlyB ABC exporter. We have observed a change in the contact surface between the monomers after hydrolysis, even though we have not seen dimer opening in any of the five 100 ns simulations. We have also identified specific regions that respond to ATP hydrolysis, in particular the X-loop motif of ABC exporters, which has been shown to be in contact with the coupling helices of the transmembrane domains (TMDs). We propose that this motif is an important part of the NBD-TMD communication in ABC exporters. Through nonequilibrium analysis, we have also identified gradual conformational changes within a short time scale after ATP hydrolysis.
Collapse
Affiliation(s)
- João M Damas
- Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Oeiras, Portugal
| | | | | | | |
Collapse
|
26
|
Oliveira ASF, Baptista AM, Soares CM. Insights into the molecular mechanism of an ABC transporter: conformational changes in the NBD dimer of MJ0796. J Phys Chem B 2010; 114:5486-96. [PMID: 20369870 DOI: 10.1021/jp905735y] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Despite the rapid advances in the study of ABC transporters, many fundamental questions linked to ATP binding/hydrolysis and its relation to the transport cycle remain unanswered. In particular, it is still neither clear nor consensual how the ATP energy is used by the nucleotide binding domains (NBDs) to produce mechanical work and drive the substrate translocation. The major conformational changes in the NBDs following ATP hydrolysis during the transport cycle and the role played by the conserved family motifs in harnessing the energy associated with nucleotide hydrolysis are yet unknown. Additionally, the way energy is transmitted from the catalytic to the membrane domains, in order to drive substrate translocation, is also a fundamental question that remains unanswered. Due to the high structure similarities of the NBD architecture throughout the whole ABC family, it is likely that the mechanism of ATP binding, hydrolysis, and communication with the transmembrane domains is similar in all family members, independently of the nature of the transported substrate. In this work, we focused our attention on the consequences of ATP hydrolysis in the NBDs, especially on the structural changes that occur during this process. For that, we use molecular dynamics simulation techniques taking as a starting point the X-ray structure of the MJ0796 dimer from Methanococcus jannaschii. Several potential intermediate states of the ATP hydrolytic cycle are investigated, each consisting of different combinations of nucleotide-bound forms. The results obtained allowed us to identify the conformational rearrangements induced by hydrolysis on the catalytic subunits, as well as the residues involved in this reorganization. The major changes are localized at specific regions of the protein, namely, involving segments 11-19 and 93-124. Additionally, our results together with the knowledge of complete ABC transporter X-ray structures suggest a possible NBD:TMD signal transmission interface.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, 2780-157 Oeiras, Portugal
| | | | | |
Collapse
|
27
|
Valente FMA, Oliveira ASF, Gnadt N, Pacheco I, Coelho AV, Xavier AV, Teixeira M, Soares CM, Pereira IAC. Hydrogenases in Desulfovibrio vulgaris Hildenborough: structural and physiologic characterisation of the membrane-bound [NiFeSe] hydrogenase. J Biol Inorg Chem 2005; 10:667-82. [PMID: 16187073 DOI: 10.1007/s00775-005-0022-4] [Citation(s) in RCA: 64] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2005] [Accepted: 08/11/2005] [Indexed: 10/25/2022]
Abstract
The genome of Desulfovibrio vulgaris Hildenborough (DvH) encodes for six hydrogenases (Hases), making it an interesting organism to study the role of these proteins in sulphate respiration. In this work we address the role of the [NiFeSe] Hase, found to be the major Hase associated with the cytoplasmic membrane. The purified enzyme displays interesting catalytic properties, such as a very high H(2) production activity, which is dependent on the presence of phospholipids or detergent, and resistance to oxygen inactivation since it is isolated aerobically in a Ni(II) oxidation state. Evidence was obtained that the [NiFeSe] Hase is post-translationally modified to include a hydrophobic group bound to the N-terminal, which is responsible for its membrane association. Cleavage of this group originates a soluble, less active form of the enzyme. Sequence analysis shows that [NiFeSe] Hases from Desulfovibrionacae form a separate family from the [NiFe] enzymes of these organisms, and are more closely related to [NiFe] Hases from more distant bacterial species that have a medial [4Fe4S](2+/1+) cluster, but not a selenocysteine. The interaction of the [NiFeSe] Hase with periplasmic cytochromes was investigated and is similar to the [NiFe](1) Hase, with the Type I cytochrome c (3) as the preferred electron acceptor. A model of the DvH [NiFeSe] Hase was generated based on the structure of the Desulfomicrobium baculatum enzyme. The structures of the two [NiFeSe] Hases are compared with the structures of [NiFe] Hases, to evaluate the consensual structural differences between the two families. Several conserved residues close to the redox centres were identified, which may be relevant to the higher activity displayed by [NiFeSe] Hases.
Collapse
Affiliation(s)
- Filipa M A Valente
- Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Oeiras, Portugal
| | | | | | | | | | | | | | | | | |
Collapse
|
28
|
Abstract
Molecular dynamics simulation (MD) constitutes an alternative to time-consuming experiments for studying conformational changes. We apply MD on a redox system where experimental information exists for the fully oxidized and fully reduced states: tetraheme cytochrome c3. Instead of doing one simulation for each state, we apply 10 4-ns replicas for both states, which provides robust statistics to characterize the redox changes. Besides these long simulations, we perform 120 short ones (50 ps), where an equilibrated oxidized state is perturbed to a reduced state. This allows the application of a nonequilibrium method, the subtraction technique, which makes it possible to characterize the different timescales of conformational changes. Reduction induces conformational changes in the N-terminus and on the loops spanning residues 36-42 and 88-93, which correlate very well with experiments, demonstrating the applicability of this methodology. We also analyze the effect of reduction on hydrogen bonds, solvent accessible surface and bound water, the changes being found to involve the hemes and propionate groups. Redox-induced protonation is also investigated, by protonating the propionates D from hemes I and IV. Although this change in the former does not have major conformational consequences, it induces in the latter conformational changes beyond the ones obtained with reduction.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, 2781-901 Oeiras, Portugal
| | | | | | | |
Collapse
|
29
|
Teixeira VH, Cunha CA, Machuqueiro M, Oliveira ASF, Victor BL, Soares CM, Baptista AM. On the Use of Different Dielectric Constants for Computing Individual and Pairwise Terms in Poisson−Boltzmann Studies of Protein Ionization Equilibrium. J Phys Chem B 2005; 109:14691-706. [PMID: 16852854 DOI: 10.1021/jp052259f] [Citation(s) in RCA: 76] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Poisson-Boltzmann (PB) models are a fast and common tool for studying electrostatic processes in proteins, particularly their ionization equilibrium (protonation and/or reduction), often yielding quite good results when compared with more detailed models. Yet, they are conceptually very simple and necessarily approximate, their empirical character being most evident when it comes to the choice of the dielectric constant assigned to the protein region. The present study analyzes several factors affecting the ability of PB-based methods to model protein ionization equilibrium. We give particular attention to a suggestion made by Warshel and co-workers (e.g., Sham et al. J. Phys. Chem. B 1997, 101, 4458) of using different protein dielectric constants for computing the individual (site) and the pairwise (site-site) terms of the ionization free energies. Our prediction of pK(a) values for several proteins indicates that no advantage is obtained by such a procedure, even for sites that are buried and/or display large pK(a) shifts relative to the solution values. In particular, the present methodology gives the best predictions using a dielectric constant around 20, for shifted/buried and nonshifted/exposed sites alike. The similarities and differences between the PB model and Warshel's PDLD/S model are discussed, as well as the reasons behind their apparently discrepant results. The present PB model is shown to predict also good reduction potentials in redox proteins.
Collapse
Affiliation(s)
- Vitor H Teixeira
- Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Av. da República, EAN, Apartado 127, 2781-901 Oeiras, Portugal
| | | | | | | | | | | | | |
Collapse
|