pKa determination of histidine residues in α-conotoxin MII peptides by 1H NMR and constant pH molecular dynamics simulation

On the Use of the Discrete Constant pH Molecular Dynamics to Describe the Conformational Space of Peptides

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.

α-Conotoxin Peptidomimetics: Probing the Minimal Binding Motif for Effective Analgesia

Several analgesic α-conotoxins have been isolated from marine cone snails. Structural modification of native peptides has provided potent and selective analogues for two of its known biological targets-nicotinic acetylcholine and γ-aminobutyric acid (GABA) G protein-coupled (GABAB) receptors. Both of these molecular targets are implicated in pain pathways. Despite their small size, an incomplete understanding of the structure-activity relationship of α-conotoxins at each of these targets has hampered the development of therapeutic leads. This review scrutinises the N-terminal domain of the α-conotoxin family of peptides, a region defined by an invariant disulfide bridge, a turn-inducing proline residue and multiple polar sidechain residues, and focusses on structural features that provide analgesia through inhibition of high-voltage-activated Ca2+ channels. Elucidating the bioactive conformation of this region of these peptides may hold the key to discovering potent drugs for the unmet management of debilitating chronic pain associated with a wide range of medical conditions.

Snails In Silico: A Review of Computational Studies on the Conopeptides

Marine cone snails are carnivorous gastropods that use peptide toxins called conopeptides both as a defense mechanism and as a means to immobilize and kill their prey. These peptide toxins exhibit a large chemical diversity that enables exquisite specificity and potency for target receptor proteins. This diversity arises in terms of variations both in amino acid sequence and length, and in posttranslational modifications, particularly the formation of multiple disulfide linkages. Most of the functionally characterized conopeptides target ion channels of animal nervous systems, which has led to research on their therapeutic applications. Many facets of the underlying molecular mechanisms responsible for the specificity and virulence of conopeptides, however, remain poorly understood. In this review, we will explore the chemical diversity of conopeptides from a computational perspective. First, we discuss current approaches used for classifying conopeptides. Next, we review different computational strategies that have been applied to understanding and predicting their structure and function, from machine learning techniques for predictive classification to docking studies and molecular dynamics simulations for molecular-level understanding. We then review recent novel computational approaches for rapid high-throughput screening and chemical design of conopeptides for particular applications. We close with an assessment of the state of the field, emphasizing important questions for future lines of inquiry.

Mutagenesis of α-Conotoxins for Enhancing Activity and Selectivity for Nicotinic Acetylcholine Receptors

Nicotinic acetylcholine receptors (nAChRs) are found throughout the mammalian body and have been studied extensively because of their implication in a myriad of diseases. α-Conotoxins (α-CTxs) are peptide neurotoxins found in the venom of marine snails of genus Conus. α-CTxs are potent and selective antagonists for a variety of nAChR isoforms. Over the past 40 years, α-CTxs have proven to be valuable molecular probes capable of differentiating between closely related nAChR subtypes and have contributed greatly to understanding the physiological role of nAChRs in the mammalian nervous system. Here, we review the amino acid composition and structure of several α-CTxs that selectively target nAChR isoforms and explore strategies and outcomes for introducing mutations in native α-CTxs to direct selectivity and enhance binding affinity for specific nAChRs. This review will focus on structure-activity relationship studies involving native α-CTxs that have been rationally mutated and molecular interactions that underlie binding between ligand and nAChR isoform.

A Simple Add-and-Display Method for Immobilisation of Cancer Drug on His-tagged Virus-like Nanoparticles for Controlled Drug Delivery

pH-responsive virus-like nanoparticles (VLNPs) hold promising potential as drug delivery systems for cancer therapy. In the present study, hepatitis B virus (HBV) VLNPs harbouring His-tags were used to display doxorubicin (DOX) via nitrilotriacetic acid (NTA) conjugation. The His-tags served as pH-responsive nanojoints which released DOX from VLNPs in a controlled manner. The His-tagged VLNPs conjugated non-covalently with NTA-DOX, and cross-linked with folic acid (FA) were able to specifically target and deliver the DOX into ovarian cancer cells via folate receptor (FR)-mediated endocytosis. The cytotoxicity and cellular uptake results revealed that the His-tagged VLNPs significantly increased the accumulation of DOX in the ovarian cancer cells and enhanced the uptake of DOX, which improved anti-tumour effects. This study demonstrated that NTA-DOX can be easily displayed on His-tagged VLNPs by a simple Add-and-Display step with high coupling efficiency and the drug was only released at low pH in a controlled manner. This approach facilitates specific attachment of any drug molecule on His-tagged VLNPs at the very mild conditions without changing the biological structure and native conformation of the VLNPs.

The crystal structure of Ac-AChBP in complex with α-conotoxin LvIA reveals the mechanism of its selectivity towards different nAChR subtypes

The α3* nAChRs, which are considered to be promising drug targets for problems such as pain, addiction, cardiovascular function, cognitive disorders etc., are found throughout the central and peripheral nervous system. The α-conotoxin (α-CTx) LvIA has been identified as the most selective inhibitor of α3β2 nAChRs known to date, and it can distinguish the α3β2 nAChR subtype from the α6/α3β2β3 and α3β4 nAChR subtypes. However, the mechanism of its selectivity towards α3β2, α6/α3β2β3, and α3β4 nAChRs remains elusive. Here we report the co-crystal structure of LvIA in complex with Aplysia californica acetylcholine binding protein (Ac-AChBP) at a resolution of 3.4 Å. Based on the structure of this complex, together with homology modeling based on other nAChR subtypes and binding affinity assays, we conclude that Asp-11 of LvIA plays an important role in the selectivity of LvIA towards α3β2 and α3/α6β2β3 nAChRs by making a salt bridge with Lys-155 of the rat α3 subunit. Asn-9 lies within a hydrophobic pocket that is formed by Met-36, Thr-59, and Phe-119 of the rat β2 subunit in the α3β2 nAChR model, revealing the reason for its more potent selectivity towards the α3β2 nAChR subtype. These results provide molecular insights that can be used to design ligands that selectively target α3β2 nAChRs, with significant implications for the design of new therapeutic α-CTxs.

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

Protein internal ionizable groups can exhibit large shifts in pK_a values. Although the environment and interaction changes have been extensively studied both experimentally and computationally, direct calculation of pK_a values of these internal ionizable groups in explicit water is challenging due to energy barriers in solvent interaction and in conformational transition. The virtual mixture of multiple states (VMMS) method is a new approach designed to study chemical state equilibrium. This method constructs a virtual mixture of multiple chemical states in order to sample the conformational space of all states simultaneously and to avoid crossing energy barriers related to state transition. By applying VMMS to 25 variants of staphylococcal nuclease with lysine residues at internal positions, we obtained the pK_a values of these lysine residues and investigated the physics underlining the pK_a shifts. Our calculation results agree reasonably well with experimental measurements, validating the VMMS method for pK_a calculation and providing molecular details of the protonation equilibrium for protein internal ionizable groups. Based on our analyses of protein conformation relaxation, lysine side chain flexibility, water penetration, and the microenvironment, we conclude that the hydrophobicity of the microenvironment around the lysine side chain (which affects water penetration differently for different protonation states) plays an important role in the pK_a shifts.

Assessing the Indirect Photochemical Transformation of Dissolved Combined Amino Acids through the Use of Systematically Designed Histidine-Containing Oligopeptides

Photooxidation is an important abiotic transformation pathway for amino acids (AAs) in sunlit waters. Although dissolved free AAs are well studied, the photooxidation of dissolved combined AAs (DCAAs) remains poorly investigated. This study is a systematic investigation of the effect of neighboring photostable AA residues (i.e., aliphatic, cationic, anionic, or aromatic residues) on the environmental indirect photochemical transformation of histidine (His) in His-containing oligopeptides. The pKa values of His residues in the studied oligopeptides were found to be between 4.3 and 8.1. Accordingly, the phototransformation rate constants of the His-containing oligopeptides were highly pH-dependent in an environmentally relevant pH range with higher reactivity for neutral His than for the protonated species. The photostable AA residues significantly modulated the photoreactivity of oligopeptides either through altering the accessibility of His to photochemically produced oxidants or through shifting the pKa values of His residues. In addition, the influence of neighboring photostable AA residues on the sorption-enhanced phototransformation of oligopeptides in solutions containing chromophoric dissolved organic matter (CDOM) was assessed. The constituent photostable AA residues promoted sorption of His-containing oligopeptides to CDOM macromolecules through electrostatic attraction, hydrophobic effects, and/or low-barrier hydrogen bonds, and subsequently increased the apparent phototransformation rate constants by up to 2 orders of magnitude.

The impact of active site protonation on substrate ring conformation in Melanocarpus albomyces cellobiohydrolase Cel7B

Understanding how the protonation state of active site residues impacts the enzyme's structure and substrate conformation is important for improving the efficiency and economic viability of the degradation of cellulosic materials as feedstock for liquid fuel and value-added chemicals.

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 pK_a 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.

Acetylcholine promotes binding of α-conotoxin MII at α3 β2 nicotinic acetylcholine receptors

α-Conotoxin MII (α-CTxMII) is a 16-residue peptide with the sequence GCCSNPVCHLEHSNLC, containing Cys2-Cys8 and Cys3-Cys16 disulfide bonds. This peptide, isolated from the venom of the marine cone snail Conus magus, is a potent and selective antagonist of neuronal nicotinic acetylcholine receptors (nAChRs). To evaluate the impact of channel-ligand interactions on ligand-binding affinity, homology models of the heteropentameric α3β2-nAChR were constructed. The models were created in MODELLER with the aid of experimentally characterized structures of the Torpedo marmorata-nAChR (Tm-nAChR, PDB ID: 2BG9) and the Aplysia californica-acetylcholine binding protein (Ac-AChBP, PDB ID: 2BR8) as templates for the α3- and β2-subunit isoforms derived from rat neuronal nAChR primary amino acid sequences. Molecular docking calculations were performed with AutoDock to evaluate interactions of the heteropentameric nAChR homology models with the ligands acetylcholine (ACh) and α-CTxMII. The nAChR homology models described here bind ACh with binding energies commensurate with those of previously reported systems, and identify critical interactions that facilitate both ACh and α-CTxMII ligand binding. The docking calculations revealed an increased binding affinity of the α3β2-nAChR for α-CTxMII with ACh bound to the receptor, and this was confirmed through two-electrode voltage clamp experiments on oocytes from Xenopus laevis. These findings provide insights into the inhibition and mechanism of electrostatically driven antagonist properties of the α-CTxMIIs on nAChRs.