Origin of pKa Shifts of Internal Lysine Residues in SNase Studied Via Equal-Molar VMMS Simulations in Explicit Water

Development and Recent Advances in Lysine and N-Terminal Bioconjugation for Peptides and Proteins

The demand for creation of protein diversity and regulation of protein function through native protein modification and post-translational modification has ignited the development of selective chemical modification methods for peptides and proteins. Chemical bioconjugation offers selective functionalization providing bioconjugates with desired properties and functions for diverse applications in chemical biology, medicine, and biomaterials. The amino group existing at the lysine residue and N-terminus of peptides and proteins has been extensively studied in bioconjugation because of its good nucleophilicity and high surface exposure. Herein, we review the development of chemical methods for modification of the amino groups on lysine residue and N-terminus featuring excellent selectivity, mild reaction conditions, short reaction time, high conversion, biocompatibility, and preservation of protein integrity. This review is organized based on the chemoselectivity and site-selectivity of the chemical bioconjugation reagents to the amino acid residues aiming to provide guidance for the selection of appropriate bioconjugation methods.

Protein pKa Prediction by Tree-Based Machine Learning

Protonation states of ionizable protein residues modulate many essential biological processes. For correct modeling and understanding of these processes, it is crucial to accurately determine their pKa values. Here, we present four tree-based machine learning models for protein pKa prediction. The four models, Random Forest, Extra Trees, eXtreme Gradient Boosting (XGBoost), and Light Gradient Boosting Machine (LightGBM), were trained on three experimental PDB and pKa datasets, two of which included a notable portion of internal residues. We observed similar performance among the four machine learning algorithms. The best model trained on the largest dataset performs 37% better than the widely used empirical pKa prediction tool PROPKA and 15% better than the published result from the pKa prediction method DelPhiPKa. The overall root-mean-square error (RMSE) for this model is 0.69, with surface and buried RMSE values being 0.56 and 0.78, respectively, considering six residue types (Asp, Glu, His, Lys, Cys, and Tyr), and 0.63 when considering Asp, Glu, His, and Lys only. We provide pKa predictions for proteins in human proteome from the AlphaFold Protein Structure Database and observed that 1% of Asp/Glu/Lys residues have highly shifted pKa values close to the physiological pH.

Towards solvent regulated self-activation of N-terminal disulfide bond oxidoreductase-D

N-terminal disulfide bond oxidoreductase-D (nDsbD), an essential redox enzyme in Gram-negative bacteria, consists of a single disulfide bond (Cys₁₀₃-Cys₁₀₉) in its active site. The enzymatic functions are believed to be regulated by an electron transfer mediated redox switching of the disulfide bond, which is vital in controlling bacterial virulence factors. In light of the disulfide bond's inclination towards nucleophilic cleavage, it is also plausible that an internal nucleophile could second the existing electron transfer mechanism in nDsbD. Using QM/MM MD metadynamics simulations, we explore different possibilities of generating an internal nucleophile near the nDsbD active site, which could serve as a fail-over mechanism in cleaving the disulfide bond. The simulations show the formation of the internal nucleophile Tyr₄₂O^- (F ≈ 9 kcal mol^-1) and its stabilization through the solvent medium. The static gas-phase calculations show that Tyr₄₂O^- could be a potential nucleophile for cleaving the S-S bond. Most strikingly, it is also seen that Tyr₄₂O^- and Asp₆₈OH communicate with each other through a proton-hole like water wire (F ≈ 12 kcal mol^-1), thus modulating the nucleophile formation. Accordingly, we propose the role of a solvent in regulating the internal nucleophilic reactions and the subsequent self-activation of nDsbD. We believe that this could be deterministic while designing enzyme-targeted inhibitor compounds.

pH-Dependent Conformational Changes Lead to a Highly Shifted pKa for a Buried Glutamic Acid Mutant of SNase

Ionizable residues are rarely present in the hydrophobic interior of proteins, but when they are, they play important roles in biological processes such as energy transduction and enzyme catalysis. Internal ionizable residues have anomalous experimental pK_a values with respect to their pK_a in bulk water. This work investigates the atomistic cause of the highly shifted pK_a of the internal Glu23 in the artificially mutated variant V23E of Staphylococcal Nuclease (SNase) using pH replica exchange molecular dynamics (pH-REMD) simulations. The pK_a of Glu23 obtained from our calculations is 6.55, which is elevated with respect to the glutamate pK_a of 4.40 in bulk water. The calculated value is close to the experimental pK_a of 7.10. Our simulations show that the highly shifted pK_a of Glu23 is the product of a pH-dependent conformational change, which has been observed experimentally and also seen in our simulations. We carry out an analysis of this pH-dependent conformational change in response to the protonation state change of Glu23. Using a four-state thermodynamic model, we estimate the two conformation-specific pK_a values of Glu23 and describe the coupling between the conformational and ionization equilibria.

Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor

The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.

Structure of the Cardiac Sodium Channel

Voltage-gated sodium channel Na_v1.5 generates cardiac action potentials and initiates the heartbeat. Here, we report structures of Na_V1.5 at 3.2-3.5 Å resolution. Na_V1.5 is distinguished from other sodium channels by a unique glycosyl moiety and loss of disulfide-bonding capability at the Na_Vβ subunit-interaction sites. The antiarrhythmic drug flecainide specifically targets the central cavity of the pore. The voltage sensors are partially activated, and the fast-inactivation gate is partially closed. Activation of the voltage sensor of Domain III allows binding of the isoleucine-phenylalanine-methionine (IFM) motif to the inactivation-gate receptor. Asp and Ala, in the selectivity motif DEKA, line the walls of the ion-selectivity filter, whereas Glu and Lys are in positions to accept and release Na⁺ ions via a charge-delocalization network. Arrhythmia mutation sites undergo large translocations during gating, providing a potential mechanism for pathogenic effects. Our results provide detailed insights into Na_v1.5 structure, pharmacology, activation, inactivation, ion selectivity, and arrhythmias.

Threats to adhesive/dentin interfacial integrity and next generation bio-enabled multifunctional adhesives

Nearly 100 million of the 170 million composite and amalgam restorations placed annually in the United States are replacements for failed restorations. The primary reason both composite and amalgam restorations fail is recurrent decay, for which composite restorations experience a 2.0-3.5-fold increase compared to amalgam. Recurrent decay is a pernicious problem-the standard treatment is replacement of defective composites with larger restorations that will also fail, initiating a cycle of ever-larger restorations that can lead to root canals, and eventually, to tooth loss. Unlike amalgam, composite lacks the inherent capability to seal discrepancies at the restorative material/tooth interface. The low-viscosity adhesive that bonds the composite to the tooth is intended to seal the interface, but the adhesive degrades, which can breach the composite/tooth margin. Bacteria and bacterial by-products such as acids and enzymes infiltrate the marginal gaps and the composite's inability to increase the interfacial pH facilitates cariogenic and aciduric bacterial outgrowth. Together, these characteristics encourage recurrent decay, pulpal damage, and composite failure. This review article examines key biological and physicochemical interactions involved in the failure of composite restorations and discusses innovative strategies to mitigate the negative effects of pathogens at the adhesive/dentin interface. © 2019 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater 107B:2466-2475, 2019.

PKA17-A Coarse-Grain Grid-Based Methodology and Web-Based Software for Predicting Protein pK a Shifts

We have developed and tested PKA17, a coarse-grain grid-based model for predicting protein pK _a shifts. Our pK _a predictor is currently deployed via a website interface. We have carried out parameter fitting using 442 Asp, Glu, His, Lys, and Arg residues for which experimental results are available in the literature. PROPKA software has been used for benchmarking. The average unsigned error and root-mean-square deviation (RMSD) have been found to be 0.628 and 0.831 pH units, respectively, for PKA17. The corresponding results with PROPKA are 0.761 and 1.063 units. We have assessed the robustness of the developed PKA17 methodology with a number of tests and have also explored the possibility of using a combination of PROPKA and PKA17 calculations in order to improve the accuracy of predicted pK _a values for protein residues. We have also once again confirmed that protein acidity constants are influenced almost entirely by residues in the immediate spatial proximity of the ionizable amino acids. The resulting PKA17 software has been deployed online with a web-based interface at http://users.wpi.edu/~jpcvitkovic/pka_calc.html. © 2019 Wiley Periodicals, Inc.

Modulating pH through lysine integrated dental adhesives

OBJECTIVES

The objective of this study was to explore the effect of lysine integration to dental adhesives with respect to the polymerization kinetics, neutralization capacities in the acidic microenvironment, dynamic mechanical properties, and thermal properties.

MATERIALS AND METHOD

Lysine was incorporated into liquid resin formulations at 2.5 and 5.0wt % with additional water/ethanol co-solvents. The co-monomer system contained 2-hydroxyethyl-methacrylate (HEMA) and Bisphenol A glycerolate dimethacrylate (BisGMA) with a mass ratio of 45/55. The kinetics of photopolymerization, neutralization capacities, lysine-leaching, dynamic mechanical properties and thermal properties of the control and experimental adhesives were analyzed.

RESULTS

The degree of conversion of the experimental adhesive was increased substantially at 2.5wt% lysine as compared to the control. The experimental polymers provided acute neutralization of the acidic microenvironment. Approximately half of the lysine was released from the polymer network within one month. Under dry conditions and physiologic temperatures, the incorporation of lysine did not compromise the storage modulus. Comparison of the thermal properties suggests that the more compact structure of the control adhesive inhibits movement of the polymer chains resulting in increased T_g.

SIGNIFICANCE

Incorporating lysine in the adhesive formulations led to promising results regarding modulating pH, which may serve as one aspect of a multi-spectrum approach for enhancing the durability of composite restorations. The results provide insight and lay a foundation for incorporating amino acids or peptides into adhesive formulations for pH modulation or desired bioactivity at the interfacial margin between the composite and tooth.

Reservoir pH replica exchange

We present the reservoir pH replica exchange (R-pH-REM) method for constant pH simulations. The R-pH-REM method consists of a two-step procedure; the first step involves generation of one or more reservoirs of conformations. Each reservoir is obtained from a standard or enhanced molecular dynamics simulation with a constrained (fixed) protonation state. In the second step, fixed charge constraints are relaxed, as the structures from one or more reservoirs are periodically injected into a constant pH or a pH-replica exchange (pH-REM) simulation. The benefit of this two-step process is that the computationally intensive part of conformational search can be decoupled from constant pH simulations, and various techniques for enhanced conformational sampling can be applied without the need to integrate such techniques into the pH-REM framework. Simulations on blocked Lys, KK, and KAAE peptides were used to demonstrate an agreement between pH-REM and R-pH-REM simulations. While the reservoir simulations are not needed for these small test systems, the real need arises in cases when ionizable molecules can sample two or more conformations separated by a large energy barrier, such that adequate sampling is not achieved on a time scale of standard constant pH simulations. Such problems might be encountered in protein systems that exploit conformational transitions for function. A hypothetical case is studied, a small molecule with a large torsional barrier; while results of pH-REM simulations depend on the starting structure, R-pH-REM calculations on this model system are in excellent agreement with a theoretical model.

A fast and accurate computational method for the linear-combination-based isotropic periodic sum

An isotropic periodic sum (IPS) is a powerful technique to reasonably calculate intermolecular interactions for wide range of molecular systems under periodic boundary conditions. A linear-combination-based IPS (LIPS) has been developed to attain computational accuracy close to an exact lattice sum, such as the Ewald sum. The algorithm of the original LIPS method has a high computational cost because it needs long-range interaction calculations in real space. This becomes a performance bottleneck for long-time molecular simulations. In this work, the combination of an LIPS and fast Fourier transform (FFT) was developed, and evaluated on homogeneous and heterogeneous molecular systems. This combinational approach of LIPS/FFT attained computational efficiency close to that of a smooth particle mesh Ewald while maintaining the same high accuracy as the original LIPS. We concluded that LIPS/FFT has great potential to extend the capability of IPS techniques for the fast and accurate computation of many types of molecular systems.

Hydronium Ions Accompanying Buried Acidic Residues Lead to High Apparent Dielectric Constants in the Interior of Proteins

Internal ionizable groups are known to play important roles in protein functions. A mystery that has attracted decades of extensive experimental and theoretical studies is the apparent dielectric constants experienced by buried ionizable groups, which are much higher than values expected for protein interiors. Many interpretations have been proposed, such as water penetration, conformational relaxation, local unfolding, protein intrinsic backbone fluctuations, etc. However, these interpretations conflict with many experimental observations. The virtual mixture of multiple states (VMMS) simulation method developed in our lab provides a direct approach for studying the equilibrium of multiple chemical states and can monitor p K_a values along simulation trajectories. Through VMMS simulations of staphylococcal nuclease (SNase) variants with internal Asp or Glu residues, we discovered that cations were attracted to buried deprotonated acidic groups and the presence of the nearby cations were essential to reproduce experimentally measured p K_a values. This finding, combined with structural analysis and validation simulations, suggests that the proton released from a deprotonation process stays near the deprotonated group inside proteins, possibly in the form of a hydronium ion. The existence of a proton near a buried charge has many implications in our understanding of protein functions.

Critical test of isotropic periodic sum techniques with group-based cut-off schemes

Truncation is still chosen for many long-range intermolecular interaction calculations to efficiently compute free-boundary systems, macromolecular systems and net-charge molecular systems, for example. Advanced truncation methods have been developed for long-range intermolecular interactions. Every truncation method can be implemented as one of two basic cut-off schemes, namely either an atom-based or a group-based cut-off scheme. The former computes interactions of “atoms” inside the cut-off radius, whereas the latter computes interactions of “molecules” inside the cut-off radius. In this work, the effect of group-based cut-off is investigated for isotropic periodic sum (IPS) techniques, which are promising cut-off treatments to attain advanced accuracy for many types of molecular system. The effect of group-based cut-off is clearly different from that of atom-based cut-off, and severe artefacts are observed in some cases. However, no severe discrepancy from the Ewald sum is observed with the extended IPS techniques.

Anomalous Properties of Lys Residues Buried in the Hydrophobic Interior of a Protein Revealed with 15N-Detect NMR Spectroscopy

Ionizable residues buried in hydrophobic environments in proteins are essential for many fundamental biochemical processes. These residues titrate with anomalous pK_a values that are challenging to reproduce with structure-based calculations owing to the conformational reorganization coupled to their ionization. Detailed characterization of this conformational reorganization is of interest; unfortunately, the properties of buried Lys residues are difficult to study experimentally. Here we demonstrate the utility of ¹⁵N NMR spectroscopy to gain insight into the protonation state, state of hydration and conformational dynamics of the Nζ amino group of buried Lys residues. The experiments were applied to five variants of staphylococcal nuclease, with internal Lys residues that titrate with pK_a values ranging from 6.2 to 8.1. Direct detection of buried Lys residues with these NMR spectroscopy methods will enable correlation between thermodynamic and structural data as well as unprecedented examination of how conformational transitions coupled to their ionization affect their pK_a values.