Browse
Submit Data
Databases
API
Help

Dataset Information

17 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Changes in Protonation States of In-Pathway Residues can Alter Ligand Binding Pathways Obtained From Spontaneous Binding Molecular Dynamics Simulations.

ABSTRACT: Protein-ligand binding processes often involve changes in protonation states that can be key to recognize and orient the ligand in the binding site. The pathways through which (bio)molecules interplay to attain productively bound complexes are intricate and involve a series of interconnected intermediate and transition states. Molecular dynamics (MD) simulations and enhanced sampling techniques are commonly used to characterize the spontaneous binding of a ligand to its receptor. However, the effect of protonation state changes of in-pathway residues in spontaneous binding MD simulations remained mostly unexplored. Here, we used molecular dynamics simulations to reconstruct the trypsin-benzamidine binding pathway considering different protonation states of His57. This residue is part of the trypsin catalytic triad and is located more than 10 Å away from Asp189, which is responsible for benzamidine binding in the trypsin S1 pocket. Our MD simulations showed that the binding pathways that benzamidine follow to target the S1 binding site are critically dependent on the His57 protonation state. Binding of benzamidine frequently occurs when His57 is protonated in the delta nitrogen while the binding process is significantly less frequent when His57 is positively charged. Constant-pH MD simulations retrieved the equilibrium populations of His57 protonation states at trypsin active pH offering a clearer picture of benzamidine recognition and binding. These results indicate that properly accounting for protonation states of distal residues can be important in spontaneous binding MD simulations.

SUBMITTER: Girame H

PROVIDER: S-EPMC9289141 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

How Ligand Protonation State Controls Water in Protein-Ligand Binding.

Project description:The role of water in protein-ligand binding has been an intensely studied topic in recent years; however, how ligand protonation state change perturbs water has not been considered. Here we show that water dynamics and interactions can be controlled by the protonation state of ligand using continuous constant pH molecular dynamics simulations of two closely related model systems, β-secretase 1 and 2 (BACE1 and BACE2), in complex with a small-molecule inhibitor. Simulations revealed that, upon binding, the inhibitor pyrimidine ring remains deprotonated in BACE1 but becomes protonated in BACE2. Pyrimidine protonation results in water displacement, rigidification of the binding pocket, and shift in the ligand binding mode from water-mediated to direct hydrogen bonding. These findings not only support but also rationalize the most recent structure-selectivity data in BACE1 drug design. Binding-induced protonation state changes are likely common; our work offers a glimpse at how modeling protein-ligand binding while allowing ligand titration can further advance the understanding of water and structure-based drug design.

| S-EPMC6467052 | biostudies-literature

Protonation and pK changes in protein-ligand binding.

Project description:Formation of protein-ligand complexes causes various changes in both the receptor and the ligand. This review focuses on changes in pK and protonation states of ionizable groups that accompany protein-ligand binding. Physical origins of these effects are outlined, followed by a brief overview of the computational methods to predict them and the associated corrections to receptor-ligand binding affinities. Statistical prevalence, magnitude and spatial distribution of the pK and protonation state changes in protein-ligand binding are discussed in detail, based on both experimental and theoretical studies. While there is no doubt that these changes occur, they do not occur all the time; the estimated prevalence varies, both between individual complexes and by method. The changes occur not only in the immediate vicinity of the interface but also sometimes far away. When receptor-ligand binding is associated with protonation state change at particular pH, the binding becomes pH dependent: we review the interplay between sub-cellular characteristic pH and optimum pH of receptor-ligand binding. It is pointed out that there is a tendency for protonation state changes upon binding to be minimal at physiologically relevant pH for each complex (no net proton uptake/release), suggesting that native receptor-ligand interactions have evolved to reduce the energy cost associated with ionization changes. As a result, previously reported statistical prevalence of these changes - typically computed at the same pH for all complexes - may be higher than what may be expected at optimum pH specific to each complex. We also discuss whether proper account of protonation state changes appears to improve practical docking and scoring outcomes relevant to structure-based drug design. An overview of some of the existing challenges in the field is provided in conclusion.

| S-EPMC4437766 | biostudies-literature

Web application for studying the free energy of binding and protonation states of protein-ligand complexes based on HINT.

Project description:A public web server performing computational titration at the active site in a protein-ligand complex has been implemented. This calculation is based on the Hydropathic interaction noncovalent force field. From 3D coordinate data for the protein, ligand and bridging waters (if available), the server predicts the best combination of protonation states for each ionizable residue and/or ligand functional group as well as the Gibbs free energy of binding for the ionization-optimized protein-ligand complex. The 3D structure for the modified molecules is available as output. In addition, a graph depicting how this energy changes with acidity, i.e., as a function of added protons, can be obtained. This data may prove to be of use in preparing models for virtual screening and molecular docking. A few illustrative examples are presented. In beta secretase (2va7) computational titration flipped the amide groups of Gln12 and Asn37 and protonated a ligand amine yielding an improvement of 6.37 kcal mol(-1) in the protein-ligand binding score. Protonation of Glu139 in mutant HIV-1 reverse transcriptase (2opq) allows a water bridge between the protein and inhibitor that increases the protein-ligand interaction score by 0.16 kcal mol(-1). In human sialidase NEU2 complexed with an isobutyl ether mimetic inhibitor (2f11) computational titration suggested that protonating Glu218, deprotonating Arg237, flipping the amide bond on Tyr334, and optimizing the positions of several other polar protons would increase the protein-ligand interaction score by 0.71 kcal mol(-1).

| S-EPMC2730983 | biostudies-literature

Mannobiose Binding Induces Changes in Hydrogen Bonding and Protonation States of Acidic Residues in Concanavalin A As Revealed by Neutron Crystallography.

Project description:Plant lectins are carbohydrate-binding proteins with various biomedical applications. Concanavalin A (Con A) holds promise in treating cancerous tumors. To better understand the Con A carbohydrate binding specificity, we obtained a room-temperature neutron structure of this legume lectin in complex with a disaccharide Man?1-2Man, mannobiose. The neutron structure afforded direct visualization of the hydrogen bonding between the protein and ligand, showing that the ligand is able to alter both protonation states and interactions for residues located close to and distant from the binding site. An unprecedented low-barrier hydrogen bond was observed forming between the carboxylic side chains of Asp28 and Glu8, with the D atom positioned equidistant from the oxygen atoms having an O···D···O angle of 101.5°.

| S-EPMC5754002 | biostudies-literature

Absolute FKBP binding affinities obtained via nonequilibrium unbinding simulations.

Project description:We compute the absolute binding affinities for two ligands bound to the FKBP protein using nonequilibrium unbinding simulations. The methodology is straightforward requiring little or no modification to many modern molecular simulation packages. The approach makes use of a physical pathway, eliminating the need for complicated alchemical decoupling schemes. We compare our nonequilibrium results to those obtained via a fully equilibrium approach and to experiment. The results of this study suggest that to obtain accurate results using nonequilibrium approaches one should use the stiff-spring approximation with the second cumulant expansion. From this study we conclude that nonequilibrium simulation could provide a simple means to estimate protein-ligand binding affinities.

| S-EPMC2905451 | biostudies-literature

Elucidation of the protonation states of the catalytic residues in mtKasA: implications for inhibitor design.

Project description:KasA (?-ketoacyl ACP synthase I) is involved in the biosynthetic pathway of mycolic acids, an essential component of the cell wall in Mycobacterium tuberculosis. It was shown that KasA is essential for the survival of the pathogen and thus could serve as a new drug target for the treatment of tuberculosis. The active site of KasA was previously characterized by X-ray crystallography. However, questions regarding the protonation state of specific amino acids, the orientation of the histidine groups within the active site, and additional conformers being accessible at ambient temperatures remain open and have to be addressed prior to the design of new inhibitors. We investigate the active site of KasA in this work by means of structural motifs and relative energies. Molecular dynamics (MD) simulations, free energy perturbation computations, and calculations employing the hybrid quantum mechanics/molecular mechanics (QM/MM) method made it possible to determine the protonation status and reveal important details about the catalytic mechanism of KasA. Additionally, we can rationalize the molecular basis for the acyl-transfer activity in the H311A mutant. Our data strongly suggest that inhibitors should be able to inhibit different protonation states because the enzyme can switch easily between a zwitterionic and neutral state.

| S-EPMC4397578 | biostudies-literature

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations.

Project description:Accurate prediction of binding free energies is critical to streamlining the drug development and protein design process. With the advent of GPU acceleration, absolute alchemical methods, which simulate the removal of ligand electrostatics and van der Waals interactions with the protein, have become routinely accessible and provide a physically rigorous approach that enables full consideration of flexibility and solvent interaction. However, standard explicit solvent simulations are unable to model protonation or electronic polarization changes upon ligand transfer from water to the protein interior, leading to inaccurate prediction of binding affinities for charged molecules. Here, we perform extensive simulation totaling ∼540 μs to benchmark the impact of modeling conditions on predictive accuracy for absolute alchemical simulations. Binding to urokinase plasminogen activator (UPA), a protein frequently overexpressed in metastatic tumors, is evaluated for a set of 10 inhibitors with extended flexibility, highly charged character, and titratable properties. We demonstrate that the alchemical simulations can be adapted to utilize the MBAR/PBSA method to improve the accuracy upon incorporating electronic polarization, highlighting the importance of polarization in alchemical simulations of binding affinities. Comparison of binding energy prediction at various protonation states indicates that proper electrostatic setup is also crucial in binding affinity prediction of charged systems, prompting us to propose an alternative binding mode with protonated ligand phenol and Hid-46 at the binding site, a testable hypothesis for future experimental validation.

| S-EPMC8254375 | biostudies-literature

Protoss: a holistic approach to predict tautomers and protonation states in protein-ligand complexes.

Project description:The calculation of hydrogen positions is a common preprocessing step when working with crystal structures of protein-ligand complexes. An explicit description of hydrogen atoms is generally needed in order to analyze the binding mode of particular ligands or to calculate the associated binding energies. Due to the large number of degrees of freedom resulting from different chemical moieties and the high degree of mutual dependence this problem is anything but trivial. In addition to an efficient algorithm to take care of the complexity resulting from complicated hydrogen bonding networks, a robust chemical model is needed to describe effects such as tautomerism and ionization consistently. We present a novel method for the placement of hydrogen coordinates in protein-ligand complexes which takes tautomers and protonation states of both protein and ligand into account. Our method generates the most probable hydrogen positions on the basis of an optimal hydrogen bonding network using an empirical scoring function. The high quality of our results could be verified by comparison to the manually adjusted Astex diverse set and a remarkably low rate of undesirable hydrogen contacts compared to other tools.

| S-EPMC4019353 | biostudies-literature

The influence of amino acid protonation states on molecular dynamics simulations of the bacterial porin OmpF.

Project description:Several groups, including our own, have found molecular dynamics (MD) calculations to result in the size of the pore of an outer membrane bacterial porin, OmpF, to be reduced relative to its size in the x-ray crystal structure. At the narrowest portion of its pore, loop L3 was found to move toward the opposite face of the pore, resulting in decreasing the cross-section area by a factor of approximately 2. In an earlier work, we computed the protonation states of titratable residues for this system and obtained values different from those that had been used in previous MD simulations. Here, we show that MD simulations carried out with these recently computed protonation states accurately reproduce the cross-sectional area profile of the channel lumen in agreement with the x-ray structure. Our calculations include the investigation of the effect of assigning different protonation states to the one residue, D(127), whose protonation state could not be modeled in our earlier calculation. We found that both assumptions of charge states for D(127) reproduced the lumen size profile of the x-ray structure. We also found that the charged state of D(127) had a higher degree of hydration and it induced greater mobility of polar side chains in its vicinity, indicating that the apparent polarizability of the D(127) microenvironment is a function of the D(127) protonation state.

| S-EPMC1367011 | biostudies-literature

Phosphite binding by the HtxB periplasmic binding protein depends on the protonation state of the ligand.

Project description:Phosphorus acquisition is critical for life. In low phosphate conditions, some species of bacteria have evolved mechanisms to import reduced phosphorus compounds, such as phosphite and hypophosphite, as alternative phosphorus sources. Uptake is facilitated by high-affinity periplasmic binding proteins (PBPs) that bind cargo in the periplasm and shuttle it to an ATP-binding cassette (ABC)-transporter in the bacterial inner membrane. PtxB and HtxB are the PBPs responsible for binding phosphite and hypophosphite, respectively. They recognize the P-H bond of phosphite/hypophosphite via a conserved P-H...π interaction, which confers nanomolar dissociation constants for their respective ligands. PtxB also has a low-level binding affinity for phosphate and hypophosphite, whilst HtxB can facilitate phosphite uptake in vivo. However, HtxB does not bind phosphate, thus the HtxBCDE transporter has recently been successfully exploited for biocontainment of genetically modified organisms by phosphite-dependent growth. Here we use a combination of X-ray crystallography, NMR and Microscale Thermophoresis to show that phosphite binding to HtxB depends on the protonation state of the ligand, suggesting that pH may effect the efficiency of phosphite uptake by HtxB in biotechnology applications.

| S-EPMC6629693 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data