Selective Stabilization of Aspartic Acid Protonation State within a Given Protein Conformation Occurs via Specific "Molecular Association"

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Molecular dynamics simulations suggest Thiosemicarbazones can bind p53 cancer mutant R175H

Cancer pathologies are associated with the unfolding and aggregation of most recurring mutations in the DNA Binding Domain (DBD) of p53 that coordinate the destabilization of protein. Substitution at the 175th codon with arginine to histidine (R175H, a mutation of large to small side-chain amino acid) destabilizes the DBD by 3 kcal/mol and triggers breasts, lung cancer, etc. Stabilizing the p53 mutant by small molecules offers an attractive drug-targeted anti-cancer therapy. The thiosemicarbazone (TSC) molecules NPC and DPT are known to act as zinc-metallochaperones to reactivate p53R175H. Here, a combination of LESMD simulations for 10 TSC conformations with a p53R175H receptor, single ligand-protein conformation MD, and ensemble docking with multiple p53R175H conformations observed during simulations is suggested to identify the potential binding site of the target protein in light of their importance for the direct TSC - p53R175H binding. NPC binds mutant R175H in the loop region L2-L3, forming pivotal hydrogen bonds with HIS175, pi‑sulfur bonds with TYR163, and pi-alkyl linkages with ARG174 and PRO190, all of which are contiguous to the zinc-binding native site on p53DBD. DPT, on the other hand, was primarily targeting alternative binding sites such as the loop-helix L1/H2 region and the S8 strand. The similar structural characteristics of TSC-bound p53R175H complexes with wild-type p53DBD are thought to be attributable to involved interactions that favour binding free energy contributions of TSC ligands. Our findings may be useful in the identification of novel pockets with druggable properties.

CIPDB: A biological structure databank for studying cation and π interactions

As major forces for modulating protein folding and molecular recognition, cation and π interactions are extensively identified in protein structures. They are even more competitive than hydrogen bonds in molecular recognition, thus, are vital in numerous biological processes. In this review, we introduce the methods for the identification and quantification of cation and π interactions, provide insights into the characteristics of cation and π interactions in the natural state, and reveal their biological function together with our developed database (Cation and π Interaction in Protein Data Bank; CIPDB; http://chemyang.ccnu.edu.cn/ccb/database/CIPDB). This review lays the foundation for the in-depth study of cation and π interactions and will guide the use of molecular design for drug discovery.

Water Network in the Binding Pocket of Fluorinated BPTI-Trypsin Complexes─Insights from Simulation and Experiment

Structural waters in the S1 binding pocket of β-trypsin are critical for the stabilization of the complex of β-trypsin with its inhibitor bovine pancreatic trypsin inhibitor (BPTI). The inhibitor strength of BPTI can be modulated by replacing the critical lysine residue at the P1 position by non-natural amino acids. We study BPTI variants in which the critical Lys15 in BPTI has been replaced by α-aminobutyric acid (Abu) and its fluorinated derivatives monofluoroethylglycine (MfeGly), difluoroethylglycine (DfeGly), and trifluoroethylglycine (TfeGly). We investigate the hypothesis that additional water molecules in the binding pocket can form specific noncovalent interactions with the fluorinated side chains and thereby act as an extension of the inhibitors. We report potentials of mean force (PMF) of the unbinding process for all four complexes and enzyme activity inhibition assays. Additionally, we report the protein crystal structure of the Lys15MfeGly-BPTI-β-trypsin complex (pdb: 7PH1). Both experimental and computational data show a stepwise increase in inhibitor strength with increasing fluorination of the Abu side chain. The PMF additionally shows a minimum for the encounter complex and an intermediate state just before the bound state. In the bound state, the computational analysis of the structure and dynamics of the water molecules in the S1 pocket shows a highly dynamic network of water molecules that does not indicate a rigidification or stabilizing trend in regard to energetic properties that could explain the increase in inhibitor strength. The analysis of the energy and the entropy of the water molecules in the S1 binding pocket using grid inhomogeneous solvation theory confirms this result. Overall, fluorination systematically changes the binding affinity, but the effect cannot be explained by a persistent water network in the binding pocket. Other effects, such as the hydrophobicity of fluorinated amino acids and the stability of the encounter complex as well as the additional minimum in the potential of mean force in the bound state, likely influence the affinity more directly.

Conformational Stability and Denaturation Processes of Proteins Investigated by Electrophoresis under Extreme Conditions

The functional structure of proteins results from marginally stable folded conformations. Reversible unfolding, irreversible denaturation, and deterioration can be caused by chemical and physical agents due to changes in the physicochemical conditions of pH, ionic strength, temperature, pressure, and electric field or due to the presence of a cosolvent that perturbs the delicate balance between stabilizing and destabilizing interactions and eventually induces chemical modifications. For most proteins, denaturation is a complex process involving transient intermediates in several reversible and eventually irreversible steps. Knowledge of protein stability and denaturation processes is mandatory for the development of enzymes as industrial catalysts, biopharmaceuticals, analytical and medical bioreagents, and safe industrial food. Electrophoresis techniques operating under extreme conditions are convenient tools for analyzing unfolding transitions, trapping transient intermediates, and gaining insight into the mechanisms of denaturation processes. Moreover, quantitative analysis of electrophoretic mobility transition curves allows the estimation of the conformational stability of proteins. These approaches include polyacrylamide gel electrophoresis and capillary zone electrophoresis under cold, heat, and hydrostatic pressure and in the presence of non-ionic denaturing agents or stabilizers such as polyols and heavy water. Lastly, after exposure to extremes of physical conditions, electrophoresis under standard conditions provides information on irreversible processes, slow conformational drifts, and slow renaturation processes. The impressive developments of enzyme technology with multiple applications in fine chemistry, biopharmaceutics, and nanomedicine prompted us to revisit the potentialities of these electrophoretic approaches. This feature review is illustrated with published and unpublished results obtained by the authors on cholinesterases and paraoxonase, two physiologically and toxicologically important enzymes.

3D Interaction Homology: Computational Titration of Aspartic Acid, Glutamic Acid and Histidine Can Create pH-Tunable Hydropathic Environment Maps

Aspartic acid, glutamic acid and histidine are ionizable residues occupying various protein environments and perform many different functions in structures. Their roles are tied to their acid/base equilibria, solvent exposure, and backbone conformations. We propose that the number of unique environments for ASP, GLU and HIS is quite limited. We generated maps of these residue's environments using a hydropathic scoring function to record the type and magnitude of interactions for each residue in a 2703-protein structural dataset. These maps are backbone-dependent and suggest the existence of new structural motifs for each residue type. Additionally, we developed an algorithm for tuning these maps to any pH, a potentially useful element for protein design and structure building. Here, we elucidate the complex interplay between secondary structure, relative solvent accessibility, and residue ionization states: the degree of protonation for ionizable residues increases with solvent accessibility, which in turn is notably dependent on backbone structure.

A practical guide to large-scale docking

Structure-based docking screens of large compound libraries have become common in early drug and probe discovery. As computer efficiency has improved and compound libraries have grown, the ability to screen hundreds of millions, and even billions, of compounds has become feasible for modest-sized computer clusters. This allows the rapid and cost-effective exploration and categorization of vast chemical space into a subset enriched with potential hits for a given target. To accomplish this goal at speed, approximations are used that result in undersampling of possible configurations and inaccurate predictions of absolute binding energies. Accordingly, it is important to establish controls, as are common in other fields, to enhance the likelihood of success in spite of these challenges. Here we outline best practices and control docking calculations that help evaluate docking parameters for a given target prior to undertaking a large-scale prospective screen, with exemplification in one particular target, the melatonin receptor, where following this procedure led to direct docking hits with activities in the subnanomolar range. Additional controls are suggested to ensure specific activity for experimentally validated hit compounds. These guidelines should be useful regardless of the docking software used. Docking software described in the outlined protocol (DOCK3.7) is made freely available for academic research to explore new hits for a range of targets.

DeepCys: Structure-based multiple cysteine function prediction method trained on deep neural network: Case study on domains of unknown functions belonging to COX2 domains

Cysteine (Cys) is the most reactive amino acid participating in a wide range of biological functions. In-silico predictions complement the experiments to meet the need of functional characterization. Multiple Cys function prediction algorithm is scarce, in contrast to specific function prediction algorithms. Here we present a deep neural network-based multiple Cys function prediction, available on web-server (DeepCys) (https://deepcys.herokuapp.com/). DeepCys model was trained and tested on two independent datasets curated from protein crystal structures. This prediction method requires three inputs, namely, PDB identifier (ID), chain ID and residue ID for a given Cys and outputs the probabilities of four cysteine functions, namely, disulphide, metal-binding, thioether and sulphenylation and predicts the most probable Cys function. The algorithm exploits the local and global protein properties, like, sequence and secondary structure motifs, buried fractions, microenvironments and protein/enzyme class. DeepCys outperformed most of the multiple and specific Cys function algorithms. This method can predict maximum number of cysteine functions. Moreover, for the first time, explicitly predicts thioether function. This tool was used to elucidate the cysteine functions on domains of unknown functions belonging to cytochrome C oxidase subunit-II like transmembrane domains. Apart from the web-server, a standalone program is also available on GitHub (https://github.com/vam-sin/deepcys).