scholarly journals Protein pKa prediction with machine learning

Author(s):  
Zhitao Cai ◽  
Fangfang Luo ◽  
Yongxian Wang ◽  
Enling Li ◽  
Yandong Huang

Protein pKa prediction is essential for the investigation of pH-associated relationship between protein structure and function. In this work, we introduce a deep learning based protein pKa predictor DeepKa, which is trained and validated with the pKa values derived from continuous constant pH molecular dynamics (CpHMD) simulations of 279 soluble proteins. Here the CpHMD implemented in the Amber molecular dynamics package has been employed (Huang, Harris, and Shen J. Chem. Inf. Model. 2018, 58, 1372-1383). Notably, to avoid discontinuities at the boundary, grid charges are proposed to represent protein electrostatics. We show that the prediction accuracy by DeepKa is close to that by CpHMD benchmarking simulations, validating DeepKa as an efficient protein pKa predictor. In addition, the training and validation sets created in this study can be applied to the development of machine learning based protein pKa predictors in future. Finally, the grid charge representation is general and applicable to other topics, such as the protein-ligand binding affinity prediction.

2020 ◽  
Author(s):  
Stefan Hervø-Hansen ◽  
Casper Højgaard ◽  
Kristoffer Enøe Johansson ◽  
Yong Wang ◽  
Khadija Wahni ◽  
...  

ABSTRACTInteractions between charged residues are difficult to study because of the complex network of interactions found in most proteins. We have designed a purposely simple system to investigate this problem by systematically introducing individual and pairs of charged and titratable residues in a protein otherwise free of such residues. We used constant pH molecular dynamics simulations, NMR spectroscopy, and thermodynamic double mutant cycles to probe the structure and energetics of the interaction between the charged residues. We found that the partial burial of surface charges contributes to a shift in pKa value, causing an aspartate to titrate in the neutral pH range. Additionally, the interaction between pairs of residues was found to be highly context dependent, with some pairs having no apparent preferential interaction, while other pairs would engage in coupled titration forming a highly stabilized salt bridge. We find good agreement between experiments and simulations, and use the simulations to rationalize our observations and to provide a detailed mechanistic understanding of the electrostatic interactions.SignificanceElectrostatic forces are important for protein folding and are favored targets of protein engineering. However, despite the many advances in the field of protein electrostatics, the prediction of changes in protein structure and function upon introduction or removal of titratable residues is still complicated. In order to provide a basic understanding of protein electrostatics we here characterize a highly charge-depleted protein and its titratable variants by a combination of NMR spectroscopy and constant pH molecular dynamics simulations. Our investigations reveal how strongly interacting residues engaged in salt bridging, can be characterized. Furthermore, our study may also enrich and facilitate the understanding of dehydration of salt-bridges and its potential effect on protein stability.


2017 ◽  
Author(s):  
Jana Shen ◽  
Zhi Yue ◽  
Helen Zgurskaya ◽  
Wei Chen

AcrB is the inner-membrane transporter of E. coli AcrAB-TolC tripartite efflux complex, which plays a major role in the intrinsic resistance to clinically important antibiotics. AcrB pumps a wide range of toxic substrates by utilizing the proton gradient between periplasm and cytoplasm. Crystal structures of AcrB revealed three distinct conformational states of the transport cycle, substrate access, binding and extrusion, or loose (L), tight (T) and open (O) states. However, the specific residue(s) responsible for proton binding/release and the mechanism of proton-coupled conformational cycling remain controversial. Here we use the newly developed membrane hybrid-solvent continuous constant pH molecular dynamics technique to explore the protonation states and conformational dynamics of the transmembrane domain of AcrB. Simulations show that both Asp407 and Asp408 are deprotonated in the L/T states, while only Asp408 is protonated in the O state. Remarkably, release of a proton from Asp408 in the O state results in large conformational changes, such as the lateral and vertical movement of transmembrane helices as well as the salt-bridge formation between Asp408 and Lys940 and other sidechain rearrangements among essential residues.Consistent with the crystallographic differences between the O and L protomers, simulations offer dynamic details of how proton release drives the O-to-L transition in AcrB and address the controversy regarding the proton/drug stoichiometry. This work offers a significant step towards characterizing the complete cycle of proton-coupled drug transport in AcrB and further validates the membrane hybrid-solvent CpHMD technique for studies of proton-coupled transmembrane proteins which are currently poorly understood. <p><br></p>


2014 ◽  
Vol 40 (10-11) ◽  
pp. 830-838 ◽  
Author(s):  
Wei Chen ◽  
Brian H. Morrow ◽  
Chuanyin Shi ◽  
Jana K. Shen

Polymers ◽  
2020 ◽  
Vol 13 (1) ◽  
pp. 99
Author(s):  
Cristian Privat ◽  
Sergio Madurga ◽  
Francesc Mas ◽  
Jaime Rubio-Martínez

Solvent pH is an important property that defines the protonation state of the amino acids and, therefore, modulates the interactions and the conformational space of the biochemical systems. Generally, this thermodynamic variable is poorly considered in Molecular Dynamics (MD) simulations. Fortunately, this lack has been overcome by means of the Constant pH Molecular Dynamics (CPHMD) methods in the recent decades. Several studies have reported promising results from these approaches that include pH in simulations but focus on the prediction of the effective pKa of the amino acids. In this work, we want to shed some light on the CPHMD method and its implementation in the AMBER suitcase from a conformational point of view. To achieve this goal, we performed CPHMD and conventional MD (CMD) simulations of six protonatable amino acids in a blocked tripeptide structure to compare the conformational sampling and energy distributions of both methods. The results reveal strengths and weaknesses of the CPHMD method in the implementation of AMBER18 version. The change of the protonation state according to the chemical environment is presumably an improvement in the accuracy of the simulations. However, the simulations of the deprotonated forms are not consistent, which is related to an inaccurate assignment of the partial charges of the backbone atoms in the CPHMD residues. Therefore, we recommend the CPHMD methods of AMBER program but pointing out the need to compare structural properties with experimental data to bring reliability to the conformational sampling of the simulations.


Biomolecules ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 334
Author(s):  
Shih-Ting Hong ◽  
Yu-Cheng Su ◽  
Yu-Jen Wang ◽  
Tian-Lu Cheng ◽  
Yeng-Tseng Wang

Humira is a monoclonal antibody that binds to TNF alpha, inactivates TNF alpha receptors, and inhibits inflammation. Neonatal Fc receptors can mediate the transcytosis of Humira–TNF alpha complex structures and process them toward degradation pathways, which reduces the therapeutic effect of Humira. Allowing the Humira–TNF alpha complex structures to dissociate to Humira and soluble TNF alpha in the early endosome to enable Humira recycling is crucial. We used the cytoplasmic pH (7.4), the early endosomal pH (6.0), and pKa of histidine side chains (6.0–6.4) to mutate the residues of complementarity-determining regions with histidine. Our engineered Humira (W1-Humira) can bind to TNF alpha in plasma at neutral pH and dissociate from the TNF alpha in the endosome at acidic pH. We used the constant-pH molecular dynamics, Gaussian accelerated molecular dynamics, two-dimensional potential mean force profiles, and in vitro methods to investigate the characteristics of W1-Humira. Our results revealed that the proposed Humira can bind TNF alpha with pH-dependent affinity in vitro. The W1-Humira was weaker than wild-type Humira at neutral pH in vitro, and our prediction results were close to the in vitro results. Furthermore, our approach displayed a high accuracy in antibody pH-dependent binding characteristics prediction, which may facilitate antibody drug design. Advancements in computational methods and computing power may further aid in addressing the challenges in antibody drug design.


2014 ◽  
Vol 82 (7) ◽  
pp. 1319-1331 ◽  
Author(s):  
Garrett B. Goh ◽  
Benjamin S. Hulbert ◽  
Huiqing Zhou ◽  
Charles L. Brooks

2018 ◽  
Vol 20 (5) ◽  
pp. 3523-3530 ◽  
Author(s):  
Zhi Yue ◽  
Jana Shen

Constant pH molecular dynamics simulations of BBL reveals negligible folding free energy barrier that is pH dependent and a sparsely populated dry molten globule state.


Sign in / Sign up

Export Citation Format

Share Document