Evidence for a high pKa of an aspartic acid residue in the active site of CALB by a fully atomistic multiscale approach.

Candida antarctica Lipase B (CALB) is a paradigm for the family of lipases. At pH 7, the optimal pH for catalysis, the protonation state of an aspartic acid of the active site (Asp134) could not be conclusively assigned. In fact, the pKa estimate provided by a widely used computational tool, namely PropKa, that predicts pKa values of ionizable groups in proteins based on the crystallographic structure, is only slightly above 7 (pKa = 7.25). This, along with the lack of an experimental evaluation, makes the assignment of its protonation state at neutral pH challenging. Here, we calculate the pKa of Asp134 by means of a fully atomistic multiscale computational approach based on classical molecular dynamics (MD) simulation and the perturbed matrix method (PMM), namely the MD-PMM approach. MD-PMM is able to take into account the dynamics of the system and, at the same time, to treat the deprotonation step at the quantum level. The calculations provide a pKa value of 8.9 ± 1.1, hence suggesting that Asp134 in CALB should be protonated at neutral, and even at slightly basic, pH.Communicated by Ramaswamy H. Sarma.

Origins of the pH-Responsive Photoluminescence of Peptide-Functionalized Au Nanoclusters.

Ultrasmall peptide-protected gold nanoclusters are a promising class of bioresponsive material exhibiting pH-sensitive photoluminescence. We present a theoretical insight into the effect peptide-ligand environment has on pH-responsive fluorescence, with the aim of enhancing the rational design of gold nanoclusters for bioapplications. Employing a hybrid quantum/classical computational methodology, we systematically calculate deprotonation free energies of N-terminal cysteine amine groups in proximity to the inherently fluorescent core of Au25(Peptide)18 nanoclusters. We find that subtle changes in hexapeptide sequence alter the electrostatic environment and significantly shift the conventional N-terminal amine pKa expected for amino acids free-in-solution. Our findings provide an insight into how the deprotonation equilibrium of N-terminal amine and side chain carboxyl groups cooperatively respond to solution pH changes, explaining the experimentally observed, yet elusive, pH-responsive fluorescence of peptide-functionalized Au25 clusters.

In silico design of novel proton-pump inhibitors with reduced adverse effects.

The development of new proton-pump inhibitors (PPIs) with less adverse effects by lowering the pKa values of nitrogen atoms in pyrimidine rings has been previously suggested by our group. In this work, we proposed that new PPIs should have the following features: (1) number of ring II = number of ring I + 1; (2) preferably five, six, or seven-membered heteroatomic ring for stability; and (3) 1 < pKa1 < 4. Six molecular scaffolds based on the aforementioned criteria were constructed, and R groups were extracted from compounds in extensive data sources. A virtual molecule dataset was established, and the pKa values of specific atoms on the molecules in the dataset were calculated to select the molecules with required pKa values. Drug-likeness screening was further conducted to obtain the candidates that significantly reduced the adverse effects of long-term PPI use. This study provided insights and tools for designing targeted molecules in silico that are suitable for practical applications.

In silico design of novel proton-pump inhibitors with reduced adverse effects / 医学前沿

The development of new proton-pump inhibitors (PPIs) with less adverse effects by lowering the pKa values of nitrogen atoms in pyrimidine rings has been previously suggested by our group. In this work, we proposed that new PPIs should have the following features: (1) number of ring II = number of ring I + 1; (2) preferably five, six, or seven-membered heteroatomic ring for stability; and (3) 1 < pKa1 < 4. Six molecular scaffolds based on the aforementioned criteria were constructed, and R groups were extracted from compounds in extensive data sources. A virtual molecule dataset was established, and the pKa values of specific atoms on the molecules in the dataset were calculated to select the molecules with required pKa values. Drug-likeness screening was further conducted to obtain the candidates that significantly reduced the adverse effects of long-term PPI use. This study provided insights and tools for designing targeted molecules in silico that are suitable for practical applications.

Recent development and application of constant pH molecular dynamics.

Solution pH is a critical environmental factor for chemical and biological processes. Over the last decade, significant efforts have been made in the development of constant pH molecular dynamics (pHMD) techniques for gaining detailed insights into pH-coupled dynamical phenomena. In this article we review the advancement of this field in the past five years, placing a special emphasis on the development of the all-atom continuous pHMD technique. We discuss various applications, including the prediction of pKa shifts for proteins, nucleic acids and surfactant assemblies, elucidation of pH-dependent population shifts, protein-protein and protein-RNA binding, as well as the mechanisms of pH-dependent self-assembly and phase transitions of surfactants and peptides. We also discuss future directions for the further improvement of the pHMD techniques.

A detailed representation of electrostatic energy in prediction of sequence and pH dependence of protein stability.

A molecular mechanics model, previously validated in applications to structure prediction, is shown to reproduce experiment in predictions of protein ionization state, and in predictions of sequence and pH dependence of protein stability. Over a large dataset, 1876 values of ΔΔG of folding, the RMSD is 1.34 kcal/mol. Using an alternative measure of accuracy, either the sign of the calculated ΔΔG agrees with experiment or the absolute value of the deviation is less than 1.0 kcal/mol, 1660 of 1876 data points (88.5%) pass the condition. Relative to models used previously in computer-aided protein design, the concept, we propose, most responsible for the performance of our model, and for the extensibility to non-neutral values of pH, is the treatment of electrostatic energy. The electronic structure of the protein is modeled using distributed atomic multipoles. The structured liquid state of the solvent is modeled using a dielectric continuum. A modification to the energetics of the reaction field, induced by the protein in the dielectric continuum, attempts to account for preformed multipoles of solvent water molecules and ions. An adjustable weight (with optimal value.141) applied to the total vacuum energy accounts implicitly for electronic polarization. A threshold distance, beyond which pairwise atomic interactions are neglected, is not used. In searches through subspaces of sequences and conformations, efficiency remains acceptable for useful applications.