1
|
Behera S, Balasubramanian S. Lipase A from Bacillus subtilis: Substrate Binding, Conformational Dynamics, and Signatures of a Lid. J Chem Inf Model 2023; 63:7545-7556. [PMID: 37989487 DOI: 10.1021/acs.jcim.3c01681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
Protein-ligand binding studies are crucial for understanding the molecular basis of biological processes and for further advancing industrial biocatalysis and drug discovery. Using computational modeling and molecular dynamics simulations, we investigated the binding of a butyrate ester substrate to the lipase A (LipA) enzyme of Bacillus subtilis. Besides obtaining a close agreement of the binding free energy with the experimental value, the study reveals a remarkable reorganization of the catalytic triad upon substrate binding, leading to increased essential hydrogen bond populations. The investigation shows the distortion of the oxyanion hole in both the substrate-bound and unbound states of LipA and highlights the strengthening of the same in the tetrahedral intermediate complex. Principal component analysis of the unbound ensemble reveals the dominant motion in LipA to be the movement of Loop-1 (Tyr129-Arg142) between two states that cover and uncover the active site, mirroring that of a lid prevalent in several lipases. This lid-like motion of Loop-1 is also supported by its tendency to spontaneously open up at an oil-water interface. Overall, this study provides valuable insights into the impact of substrate binding on the structure, flexibility, and conformational dynamics of the LipA enzyme.
Collapse
Affiliation(s)
- Sudarshan Behera
- Chemistry and Physics of Materials Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore 560 064, India
| | - Sundaram Balasubramanian
- Chemistry and Physics of Materials Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore 560 064, India
| |
Collapse
|
2
|
Papadourakis M, Sinenka H, Matricon P, Hénin J, Brannigan G, Pérez-Benito L, Pande V, van Vlijmen H, de Graaf C, Deflorian F, Tresadern G, Cecchini M, Cournia Z. Alchemical Free Energy Calculations on Membrane-Associated Proteins. J Chem Theory Comput 2023; 19:7437-7458. [PMID: 37902715 PMCID: PMC11017255 DOI: 10.1021/acs.jctc.3c00365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Indexed: 10/31/2023]
Abstract
Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.
Collapse
Affiliation(s)
- Michail Papadourakis
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Hryhory Sinenka
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Pierre Matricon
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Jérôme Hénin
- Laboratoire
de Biochimie Théorique UPR 9080, CNRS and Université Paris Cité, 75005 Paris, France
| | - Grace Brannigan
- Center
for Computational and Integrative Biology, Rutgers University−Camden, Camden, New Jersey 08103, United States of America
- Department
of Physics, Rutgers University−Camden, Camden, New Jersey 08102, United States
of America
| | - Laura Pérez-Benito
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Vineet Pande
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Herman van Vlijmen
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Chris de Graaf
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Francesca Deflorian
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Gary Tresadern
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Marco Cecchini
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Zoe Cournia
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| |
Collapse
|
3
|
Gating Mechanism of the Voltage-Gated Proton Channel Studied by Molecular Dynamics Simulations. Molecules 2022; 27:molecules27072277. [PMID: 35408673 PMCID: PMC9000549 DOI: 10.3390/molecules27072277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Revised: 03/20/2022] [Accepted: 03/28/2022] [Indexed: 11/17/2022] Open
Abstract
The voltage-gated proton channel Hv1 has important roles in proton extrusion, pH homeostasis, sperm motility, and cancer progression. The Hv1 channel has also been found to be highly expressed in cell lines and tissue samples from patients with breast cancer. A high-resolution closed-state structure has been reported for the mouse Hv1 chimera channel (mHv1cc), solved by X-ray crystallography, but the open-state structure of Hv1 has not been solved. Since Hv1 is a promising drug target, various groups have proposed open conformations by molecular modeling and simulation studies. However, the gating mechanism and the open-state conformation under the membrane potential are still debate. Here, we present a molecular dynamics study considering membrane potential and pH conditions. The closed-state structure of mHv1cc was used to run molecular dynamics (MD) simulations with respect to electric field and pH conditions in order to investigate the mechanism of proton transfer. We observed a continuous hydrogen bond chain of water molecules called a water-wire to be formed through the channel pore in the channel opening, triggered by downward displacement of the S2 helix and upward movement of the S4 helix relative to other helices. Due to the movement of the S2 and S4 helices, the internal salt bridge network was rearranged, and the hydrophobic gating layers were destroyed. In line with previous experimental and simulation observations, our simulation results led us to propose a new gating mechanism for the Hv1 proton channel, and may provide valuable information for novel drug discovery.
Collapse
|
4
|
Wu Z, Biggin PC. Correction Schemes for Absolute Binding Free Energies Involving Lipid Bilayers. J Chem Theory Comput 2022; 18:2657-2672. [PMID: 35315270 PMCID: PMC9082507 DOI: 10.1021/acs.jctc.1c01251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Absolute
binding free-energy (ABFE) calculations are playing an
increasing role in drug design, especially as they can be performed
on a range of disparate compounds and direct comparisons between them
can be made. It is, however, especially important to ensure that they
are as accurate as possible, as unlike relative binding free-energy
(RBFE) calculations, one does not benefit as much from a cancellation
of errors during the calculations. In most modern implementations
of ABFE calculations, a particle mesh Ewald scheme is typically used
to treat the electrostatic contribution to the free energy. A central
requirement of such schemes is that the box preserves neutrality throughout
the calculation. There are many ways to deal with this problem that
have been discussed over the years ranging from a neutralizing plasma
with a post hoc correction term through to a simple co-alchemical
ion within the same box. The post hoc correction approach is the most
widespread. However, the vast majority of these studies have been
applied to a soluble protein in a homogeneous solvent (water or salt
solution). In this work, we explore which of the more common approaches
would be the most suitable for a simulation box with a lipid bilayer
within it. We further develop the idea of the so-called Rocklin correction
for lipid-bilayer systems and show how such a correction could work.
However, we also show that it will be difficult to make this generalizable
in a practical way and thus we conclude that the use of a “co-alchemical
ion” is the most useful approach for simulations involving
lipid membrane systems.
Collapse
Affiliation(s)
- Zhiyi Wu
- Department of Biochemistry, South Parks Road, Oxford OX1 3QU, U.K
| | - Philip C Biggin
- Department of Biochemistry, South Parks Road, Oxford OX1 3QU, U.K
| |
Collapse
|
5
|
Zhao C, Hong L, Galpin JD, Riahi S, Lim VT, Webster PD, Tobias DJ, Ahern CA, Tombola F. HIFs: New arginine mimic inhibitors of the Hv1 channel with improved VSD-ligand interactions. J Gen Physiol 2021; 153:212451. [PMID: 34228044 PMCID: PMC8263924 DOI: 10.1085/jgp.202012832] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Accepted: 06/14/2021] [Indexed: 12/14/2022] Open
Abstract
The human voltage-gated proton channel Hv1 is a drug target for cancer, ischemic stroke, and neuroinflammation. It resides on the plasma membrane and endocytic compartments of a variety of cell types, where it mediates outward proton movement and regulates the activity of NOX enzymes. Its voltage-sensing domain (VSD) contains a gated and proton-selective conduction pathway, which can be blocked by aromatic guanidine derivatives such as 2-guanidinobenzimidazole (2GBI). Mutation of Hv1 residue F150 to alanine (F150A) was previously found to increase 2GBI apparent binding affinity more than two orders of magnitude. Here, we explore the contribution of aromatic interactions between the inhibitor and the channel in the presence and absence of the F150A mutation, using a combination of electrophysiological recordings, classic mutagenesis, and site-specific incorporation of fluorinated phenylalanines via nonsense suppression methodology. Our data suggest that the increase in apparent binding affinity is due to a rearrangement of the binding site allowed by the smaller residue at position 150. We used this information to design new arginine mimics with improved affinity for the nonrearranged binding site of the wild-type channel. The new compounds, named “Hv1 Inhibitor Flexibles” (HIFs), consist of two “prongs,” an aminoimidazole ring, and an aromatic group connected by extended flexible linkers. Some HIF compounds display inhibitory properties that are superior to those of 2GBI, thus providing a promising scaffold for further development of high-affinity Hv1 inhibitors.
Collapse
Affiliation(s)
- Chang Zhao
- Department of Physiology and Biophysics, University of California, Irvine, Irvine, CA.,Chao Family Comprehensive Cancer Center, University of California, Irvine, Irvine, CA
| | - Liang Hong
- Department of Physiology and Biophysics, University of California, Irvine, Irvine, CA
| | - Jason D Galpin
- Department of Molecular Physiology and Biophysics, University of Iowa, Iowa City, IA
| | - Saleh Riahi
- Department of Chemistry, University of California, Irvine, Irvine, CA
| | - Victoria T Lim
- Department of Chemistry, University of California, Irvine, Irvine, CA
| | - Parker D Webster
- Department of Physiology and Biophysics, University of California, Irvine, Irvine, CA
| | - Douglas J Tobias
- Chao Family Comprehensive Cancer Center, University of California, Irvine, Irvine, CA.,Department of Chemistry, University of California, Irvine, Irvine, CA
| | - Christopher A Ahern
- Department of Molecular Physiology and Biophysics, University of Iowa, Iowa City, IA
| | - Francesco Tombola
- Department of Physiology and Biophysics, University of California, Irvine, Irvine, CA.,Chao Family Comprehensive Cancer Center, University of California, Irvine, Irvine, CA
| |
Collapse
|
6
|
Thermodynamics and Mechanism of the Membrane Permeation of Hv1 Channel Blockers. J Membr Biol 2020; 254:5-16. [PMID: 33196887 DOI: 10.1007/s00232-020-00149-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Accepted: 10/14/2020] [Indexed: 02/08/2023]
Abstract
The voltage-gated proton channel Hv1 mediates efflux of protons from the cell. Hv1 integrally contributes to various physiological processes including pH homeostasis and the respiratory burst of phagocytes. Inhibition of Hv1 may provide therapeutic avenues for the treatment of inflammatory diseases, breast cancer, and ischemic brain damage. In this work, we investigate two prototypical Hv1 inhibitors, 2-guanidinobenzimidazole (2GBI), and 5-chloro-2-guanidinobenzimidazole (GBIC), from an experimentally screened class of guanidine derivatives. Both compounds block proton conduction by binding the same site located on the intracellular side of the channel. However, when added to the extracellular medium, the compounds strongly differ in their ability to inhibit proton conduction, suggesting substantial differences in membrane permeability. Here, we compute the potential of mean force for each compound to permeate through the membrane using atomistic molecular dynamics simulations with the adaptive biasing force method. Our results rationalize the putative distinction between these two blockers with respect to their abilities to permeate the cellular membrane.
Collapse
|