Energetics and dynamics of a light-driven sodium-pumping rhodopsin

The conversion of light energy into ion gradients across biological membranes is one of the most fundamental reactions in primary biological energy transduction. Recently, the structure of the first light-activated Na+ pump, Krokinobacter eikastus rhodopsin 2 (KR2), was resolved at atomic resolution [Kato HE, et al. (2015) Nature 521:48–53]. To elucidate its molecular mechanism for Na+ pumping, we perform here extensive classical and quantum molecular dynamics (MD) simulations of transient photocycle states. Our simulations show how the dynamics of key residues regulate water and ion access between the bulk and the buried light-triggered retinal site. We identify putative Na+ binding sites and show how protonation and conformational changes gate the ion through these sites toward the extracellular side. We further show by correlated ab initio quantum chemical calculations that the obtained putative photocycle intermediates are in close agreement with experimental transient optical spectroscopic data. The combined results of the ion translocation and gating mechanisms in KR2 may provide a basis for the rational design of novel light-driven ion pumps with optogenetic applications.

Download Full-text

Distinct effects of Q925 mutation on intracellular and extracellular Na+ and K+ binding to the Na+, K+-ATPase

Scientific Reports ◽

10.1038/s41598-019-50009-2 ◽

2019 ◽

Vol 9 (1) ◽

Cited By ~ 1

Author(s):

Hang N. Nielsen ◽

Kerri Spontarelli ◽

Rikke Holm ◽

Jens Peter Andersen ◽

Anja P. Einholm ◽

...

Keyword(s):

Crystal Structure ◽

Crystal Structures ◽

Binding Sites ◽

Functional Role ◽

Transmembrane Helix ◽

Binding Properties ◽

Ion Translocation ◽

Extracellular Side ◽

Insight Into

Abstract Three Na+ sites are defined in the Na+-bound crystal structure of Na+, K+-ATPase. Sites I and II overlap with two K+ sites in the K+-bound structure, whereas site III is unique and Na+ specific. A glutamine in transmembrane helix M8 (Q925) appears from the crystal structures to coordinate Na+ at site III, but does not contribute to K+ coordination at sites I and II. Here we address the functional role of Q925 in the various conformational states of Na+, K+-ATPase by examining the mutants Q925A/G/E/N/L/I/Y. We characterized these mutants both enzymatically and electrophysiologically, thereby revealing their Na+ and K+ binding properties. Remarkably, Q925 substitutions had minor effects on Na+ binding from the intracellular side of the membrane – in fact, mutations Q925A and Q925G increased the apparent Na+ affinity – but caused dramatic reductions of the binding of K+ as well as Na+ from the extracellular side of the membrane. These results provide insight into the changes taking place in the Na+-binding sites, when they are transformed from intracellular- to extracellular-facing orientation in relation to the ion translocation process, and demonstrate the interaction between sites III and I and a possible gating function of Q925 in the release of Na+ at the extracellular side.

Download Full-text

How does a small molecule bind at a cryptic binding site?

10.1101/2021.03.31.437917 ◽

2021 ◽

Author(s):

Yibing Shan ◽

Venkatesh P. Mysore ◽

Abba E. Leffler ◽

Eric T. Kim ◽

Shiori Sagawa ◽

...

Keyword(s):

Small Molecules ◽

Binding Site ◽

Small Molecule ◽

Protein Interactions ◽

Binding Sites ◽

Rational Design ◽

Structural Information ◽

Interleukin 2 ◽

Md Simulations ◽

Drug Binding

Protein-protein interactions (PPIs) are ubiquitous biomolecular processes that are central to virtually all aspects of cellular function. Identifying small molecules that modulate specific disease-related PPIs is a strategy with enormous promise for drug discovery. The design of drugs to disrupt PPIs is challenging, however, because many potential drug-binding sites at PPI interfaces are "cryptic": When unoccupied by a ligand, cryptic sites are often flat and featureless, and thus not readily recognizable in crystal structures, with the geometric and chemical characteristics of typical small-molecule binding sites only emerging upon ligand binding. The rational design of small molecules to inhibit specific PPIs would benefit from a better understanding of how such molecules bind at PPI interfaces. To this end, we have conducted unbiased, all-atom MD simulations of the binding of four small-molecule inhibitors (SP4206 and three SP4206 analogs) to interleukin 2 (IL2)—which performs its function by forming a PPI with its receptor—without incorporating any prior structural information about the ligands' binding. In multiple binding events, a small molecule settled into a stable binding pose at the PPI interface of IL2, resulting in a protein–small-molecule binding site and pose virtually identical to that observed in an existing crystal structure of the IL2-SP4206 complex. Binding of the small molecule stabilized the IL2 binding groove, which when the small molecule was not bound emerged only transiently and incompletely. Moreover, free energy perturbation (FEP) calculations successfully distinguished between the native and non-native IL2–small-molecule binding poses found in the simulations, suggesting that binding simulations in combination with FEP may provide an effective tool for identifying cryptic binding sites and determining the binding poses of small molecules designed to disrupt PPI interfaces by binding to such sites.

Download Full-text

Mechanism and potential sites of potassium interaction with glutamate transporters

The Journal of General Physiology ◽

10.1085/jgp.202012577 ◽

2020 ◽

Vol 152 (10) ◽

Author(s):

Jiali Wang ◽

Kaiqi Zhang ◽

Puja Goyal ◽

Christof Grewer

Keyword(s):

Binding Site ◽

Binding Sites ◽

Glutamate Release ◽

Indirect Evidence ◽

Md Simulations ◽

Glutamate Transporters ◽

Site Directed Mutagenesis ◽

Regulatory Function ◽

Amino Acid Residues ◽

Extracellular Side

In the mammalian glutamate transporters, countertransported intracellular K+ is essential for relocating the glutamate binding site to the extracellular side of the membrane. This K+-dependent process is believed to be rate limiting for the transport cycle. In contrast, extracellular K+ induces glutamate release upon transporter reversal. Here, we analyzed potential K+ binding sites using molecular dynamics (MD) simulations and site-directed mutagenesis. Two candidate sites were identified by spontaneous K+ binding in MD simulations, one site (K1 site) overlapping with the Na1 Na+ binding site and the K2 site being localized under hairpin loop 2 (HP2). Mutations to conserved amino acid residues in these sites resulted in several transporters that were defective in K+-induced reverse transport and which bound K+ with reduced apparent affinity compared with the wild-type transporter. However, external K+ interaction was abolished in only one mutant transporter EAAC1D454A in the K1 site. Our results, for the first time, directly demonstrate effects of K1-site mutations on K+ binding, in contrast to previous reports on K+ binding sites based on indirect evidence. We propose that K+ binding to the K1 site is responsible for catalyzing the relocation step, whereas binding to the K2 site may have an as-of-yet unidentified regulatory function.

Download Full-text

Characterisation of plasmodial transketolases and identification of potential inhibitors: an in silico study

Malaria Journal ◽

10.1186/s12936-020-03512-1 ◽

2020 ◽

Vol 19 (1) ◽

Author(s):

Rita Afriyie Boateng ◽

Özlem Tastan Bishop ◽

Thommas Mutemi Musyoka

Keyword(s):

Molecular Docking ◽

Sequence Alignment ◽

Conformational Changes ◽

Active Site ◽

Motif Discovery ◽

Rational Design ◽

Binding Free Energy ◽

Md Simulations ◽

Starting Point ◽

Potential Inhibitors

Abstract Background Plasmodial transketolase (PTKT) enzyme is one of the novel pharmacological targets being explored as potential anti-malarial drug target due to its functional role and low sequence identity to the human enzyme. Despite this, features contributing to such have not been exploited for anti-malarial drug design. Additionally, there are no anti-malarial drugs targeting PTKTs whereas the broad activity of these inhibitors against PTKTs from other Plasmodium spp. is yet to be reported. This study characterises different PTKTs [Plasmodium falciparum (PfTKT), Plasmodium vivax (PvTKT), Plasmodium ovale (PoTKT), Plasmodium malariae (PmTKT) and Plasmodium knowlesi (PkTKT) and the human homolog (HsTKT)] to identify key sequence and structural based differences as well as the identification of selective potential inhibitors against PTKTs. Methods A sequence-based study was carried out using multiple sequence alignment, phylogenetic tree calculations and motif discovery analysis. Additionally, TKT models of PfTKT, PmTKT, PoTKT, PmTKT and PkTKT were modelled using the Saccharomyces cerevisiae TKT structure as template. Based on the modelled structures, molecular docking using 623 South African natural compounds was done. The stability, conformational changes and detailed interactions of selected compounds were accessed viz all-atom molecular dynamics (MD) simulations and binding free energy (BFE) calculations. Results Sequence alignment, evolutionary and motif analyses revealed key differences between plasmodial and the human TKTs. High quality homodimeric three-dimensional PTKTs structures were constructed. Molecular docking results identified three compounds (SANC00107, SANC00411 and SANC00620) which selectively bind in the active site of all PTKTs with the lowest (better) binding affinity ≤ − 8.5 kcal/mol. MD simulations of ligand-bound systems showed stable fluctuations upon ligand binding. In all systems, ligands bind stably throughout the simulation and form crucial interactions with key active site residues. Simulations of selected compounds in complex with human TKT showed that ligands exited their binding sites at different time steps. BFE of protein–ligand complexes showed key residues involved in binding. Conclusions This study highlights significant differences between plasmodial and human TKTs and may provide valuable information for the development of novel anti-malarial inhibitors. Identified compounds may provide a starting point in the rational design of PTKT inhibitors and analogues based on these scaffolds.

Download Full-text

Analysis of mutations leading to para-aminosalicylic acid resistance in Mycobacterium tuberculosis

Scientific Reports ◽

10.1038/s41598-019-48940-5 ◽

2019 ◽

Vol 9 (1) ◽

Cited By ~ 3

Author(s):

Bharati Pandey ◽

Sonam Grover ◽

Jagdeep Kaur ◽

Abhinav Grover

Keyword(s):

Mycobacterium Tuberculosis ◽

Conformational Changes ◽

Acid Resistance ◽

Principal Component ◽

Md Simulations ◽

Missense Mutations ◽

Aminosalicylic Acid ◽

Folate Pathway ◽

Key Enzyme ◽

Key Residues

Abstract Thymidylate synthase A (ThyA) is the key enzyme involved in the folate pathway in Mycobacterium tuberculosis. Mutation of key residues of ThyA enzyme which are involved in interaction with substrate 2′-deoxyuridine-5′-monophosphate (dUMP), cofactor 5,10-methylenetetrahydrofolate (MTHF), and catalytic site have caused para-aminosalicylic acid (PAS) resistance in TB patients. Focusing on R127L, L143P, C146R, L172P, A182P, and V261G mutations, including wild-type, we performed long molecular dynamics (MD) simulations in explicit solvent to investigate the molecular principles underlying PAS resistance due to missense mutations. We found that these mutations lead to (i) extensive changes in the dUMP and MTHF binding sites, (ii) weak interaction of ThyA enzyme with dUMP and MTHF by inducing conformational changes in the structure, (iii) loss of the hydrogen bond and other atomic interactions and (iv) enhanced movement of protein atoms indicated by principal component analysis (PCA). In this study, MD simulations framework has provided considerable insight into mutation induced conformational changes in the ThyA enzyme of Mycobacterium.

Download Full-text

Selectivity filter ion binding affinity determines inactivation in a potassium channel

Proceedings of the National Academy of Sciences ◽

10.1073/pnas.2009624117 ◽

2020 ◽

Vol 117 (47) ◽

pp. 29968-29978

Author(s):

Céline Boiteux ◽

David J. Posson ◽

Toby W. Allen ◽

Crina M. Nimigean

Keyword(s):

Conformational Changes ◽

Binding Sites ◽

Bk Channels ◽

Selectivity Filter ◽

Ion Binding ◽

X Ray Crystallography ◽

Extracellular Side ◽

Opening And Closing ◽

Low K ◽

Conductive State

Potassium channels can become nonconducting via inactivation at a gate inside the highly conserved selectivity filter (SF) region near the extracellular side of the membrane. In certain ligand-gated channels, such as BK channels and MthK, a Ca2+-activated K+channel fromMethanobacterium thermoautotrophicum, the SF has been proposed to play a role in opening and closing rather than inactivation, although the underlying conformational changes are unknown. Using X-ray crystallography, identical conductive MthK structures were obtained in wide-ranging K+concentrations (6 to 150 mM), unlike KcsA, whose SF collapses at low permeant ion concentrations. Surprisingly, three of the SF’s four binding sites remained almost fully occupied throughout this range, indicating high affinities (likely submillimolar), while only the central S2 site titrated, losing its ion at 6 mM, indicating low K+affinity (∼50 mM). Molecular simulations showed that the MthK SF can also collapse in the absence of K+, similar to KcsA, but that even a single K+binding at any of the SF sites, except S4, can rescue the conductive state. The uneven titration across binding sites differs from KcsA, where SF sites display a uniform decrease in occupancy with K+concentration, in the low millimolar range, leading to SF collapse. We found that ions were disfavored in MthK’s S2 site due to weaker coordination by carbonyl groups, arising from different interactions with the pore helix and water behind the SF. We conclude that these differences in interactions endow the seemingly identical SFs of KcsA and MthK with strikingly different inactivating phenotypes.

Download Full-text

Biochemical study of fibrinolytic protease from Euphausia superba possessing multifunctional serine protease activity

Protein and Peptide Letters ◽

10.2174/0929866527666201112123714 ◽

2020 ◽

Vol 27 ◽

Author(s):

Guo-Ying Qian ◽

Gyutae Lim ◽

Shang-Jun Yin ◽

Jun-Mo Yang ◽

Jinhyuk Lee ◽

...

Keyword(s):

Defense Mechanisms ◽

Fluorescence Spectra ◽

Biochemical Study ◽

Euphausia Superba ◽

Md Simulations ◽

Time Interval ◽

Serine Residue ◽

Catalytic Function ◽

Serine Protease Activity ◽

Key Residues

Background: Background: Fibrinolytic protease from Euphausia superba (EFP) was isolated. Objective: Biochemical distinctions, regulation of the catalytic function, and the key residues of EFP were investigated. Methods: The serial inhibition kinetic evaluations coupled with measurements of fluorescence spectra in the presence of 4- (2-aminoethyl) benzene sulfonyl fluoride hydrochloride (AEBSF) was conducted. The computational molecular dynamics (MD) simulations were also applied for a comparative study. Results: The enzyme behaved as a monomeric protein with a molecular mass of about 28.6 kD with Km BApNA = 0.629 ± 0.02 mM and kcat/Km BApNA = 7.08 s-1 /mM. The real-time interval measurements revealed that the inactivation was a first-order reaction, with the kinetic processes shifting from a monophase to a biphase. Measurements of fluorescence spectra showed that serine residue modification by AEBSF directly caused conspicuous changes of the tertiary structures and exposed hydrophobic surfaces. Some osmolytes were applied to find protective roles. These results confirmed that the active region of EFP is more flexible than the overall enzyme molecule and serine, as the key residue, is associated with the regional unfolding of EFP in addition to its catalytic role. The MD simulations were supportive to the kinetics data. Conclusion: Our study indicated that EFP has an essential serine residue for its catalyst function and associated folding behaviors. Also, the functional role of osmolytes such as proline and glycine that may play a role in defense mechanisms from environmental adaptation in a krill’s body was suggested.

Download Full-text

Closing the gap: an atomistic structural and functional perspective of S. mansoni universal stress G4LZI3 protein in complex with phenolic compounds

Current Drug Discovery Technologies ◽

10.2174/1570163817666201001122436 ◽

2020 ◽

Vol 17 ◽

Author(s):

Priscilla Masamba ◽

Geraldene Munsamy ◽

Abidemi Paul Kappo

Keyword(s):

Phenolic Compounds ◽

Conformational Changes ◽

Stress Protein ◽

Md Simulations ◽

Binding Energies ◽

Developmental Cycle ◽

Structure Based Drug Design ◽

Bioinformatic Tools ◽

Drug Of Choice ◽

Functional Perspective

Background: For decades, Praziquantel has been the undisputed drug of choice for all schistosome infections, but rising concerns due to the unelucidated mechanism of action of the drug and unavoidable reports of emerging drug resistant strains has necessitated the need for alternative treatment drug. Moreover, current apprehension has been reinforced by total dependence on the drug for treatment hence, the search for novel and effective anti-schistosomal drugs. Uses: This study made use of bioinformatic tools to determine the structural binding of the Universal G4LZI3 stress protein (USP) in complex with ten polyphenol compounds, thereby highlighting the effectiveness of these recently identified ‘lead’ molecules in the design of novel therapeutics targeted against schistosomiasis. Upregulation of the G4LZI3 USP throughout the schistosome multifaceted developmental cycle sparks interest in its potential role as a druggable target. The integration of in silico tools provides an atomistic perspective into the binding of potential inhibitors to target proteins. Conclusion: This study therefore, implemented the use of molecular dynamic (MD) simulations to provide functional and structural insight into key conformational changes upon binding of G4ZLI3 to these key phenolic compounds. Post-MD analyses revealed unique structural and conformational changes in the G4LZI3 protein in complex with curcumin and catechin respectively. These systems exhibited the highest binding energies, while the major interacting residues conserved in all the complexes provides a route map for structure-based drug design of novel compounds with enhanced inhibitory potency against the G4LZI3 protein. This study suggests an alternative approach for the development of anti-schistosomal drugs using natural compounds.

Download Full-text

Correlation of membrane protein conformational and functional dynamics

Nature Communications ◽

10.1038/s41467-021-24660-1 ◽

2021 ◽

Vol 12 (1) ◽

Author(s):

Raghavendar Reddy Sanganna Gari ◽

Joel José Montalvo‐Acosta ◽

George R. Heath ◽

Yining Jiang ◽

Xiaolong Gao ◽

...

Keyword(s):

Membrane Protein ◽

Conformational Changes ◽

Single Channel ◽

Conformational Dynamics ◽

Md Simulations ◽

High Temporal Resolution ◽

Dependent Manner ◽

Functional Dynamics ◽

Ph Dependent ◽

Open States

AbstractConformational changes in ion channels lead to gating of an ion-conductive pore. Ion flux has been measured with high temporal resolution by single-channel electrophysiology for decades. However, correlation between functional and conformational dynamics remained difficult, lacking experimental techniques to monitor sub-millisecond conformational changes. Here, we use the outer membrane protein G (OmpG) as a model system where loop-6 opens and closes the β-barrel pore like a lid in a pH-dependent manner. Functionally, single-channel electrophysiology shows that while closed states are favored at acidic pH and open states are favored at physiological pH, both states coexist and rapidly interchange in all conditions. Using HS-AFM height spectroscopy (HS-AFM-HS), we monitor sub-millisecond loop-6 conformational dynamics, and compare them to the functional dynamics from single-channel recordings, while MD simulations provide atomistic details and energy landscapes of the pH-dependent loop-6 fluctuations. HS-AFM-HS offers new opportunities to analyze conformational dynamics at timescales of domain and loop fluctuations.

Download Full-text