Trypsin inhibition by a peptide hormone: crystal structure of trypsin-vasopressin complex

Unsaponifiable Compounds and Phenols Content, Antioxidant and Antitrypsin Activities of Prunus persica Kernel Oil

Although peach kernels are rich in oil, there is a lack of information about its chemical and biological properties. Therefore, the purpose of this study was to determine the lipid profile, antioxidant capacity, and trypsin inhibitory propriety of peach oil extracted from two varieties (sweet cap and O'Henry) cultivated in Tunisia. The investigated peach kernel oil contains significant amount of unsaponifiable (2.1±0.5-2.8±0.2% of oil) and phenolic compounds (45.8±0.92-74.6±1.3 mg GAE/g of oil). Its n-alkane profile was characterized by the predominance of tetracosane n-C24 (47.24%) followed by tricosane n-C23 (34.43%). An important total tocopherol content (1192.83±3.1 mg/kg oil) has been found in sweet cap cultivar. Although rich in polyphenols and tocopherols, the tested oil did not display an inhibitory effect on trypsin. However, all peach oil samples showed effective antioxidant capacity and the highest values (86.34±1.3% and 603.50±2.6 μmol TE/g oil for DPPH test and ORAC assay, respectively) were observed for sweet cap oil. Peach oil has an excellent potential for application in the food and pharmaceutical industries as source of naturally-occurring bioactive substances.

Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT

The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.

Significance of D- tryptophan in Contryphan-Ar1131 Conus peptide: Oxidative folding, trypsin binding, and photostabilization activity

Distinct differences have been observed between L-tryptophan and D-tryptophan containing contryphan-Ar1131 in oxidative folding, trypsin binding, and photostabilization activity on avobenzone. [W⁵] contryphan-Ar1131 and [w⁵] contryphan-Ar1131 were chemically synthesized and characterized using RP-HPLC and mass spectrometry. Structural differences due to the change of configuration of tryptophan were evident from the optimized structures of contryphan-Ar1131 using density functional theory (DFT). The comparison of early events of oxidative folding has revealed the role of D-tryptophan in accelerating the formation of a disulfide bond. The optimized structures of the reduced form of peptides revealed the occurrence of aromatic-aromatic and aromatic-proline interactions in [w⁵] contryphan-Ar1131 which may be critical in aiding the oxidative folding reaction. The presence of the Lys6-Pro7 peptide bond indicates that contryphan-Ar1131 is resistant but may bind to trypsin allowing to assign the binding affinity of peptides to the protein surface. Competitive binding studies and molecular docking along with molecular dynamic (MD) simulations have revealed that [w⁵] contryphan-Ar1131 has more affinity for the active site of trypsin. Given tryptophan is a photostabilizer of FDA-approved chemical UV-A filter avobenzone, the report has compared the photostabilization activity of [W⁵]/ [w⁵] contryphan-Ar1131 on avobenzone under natural sunlight. [w⁵] contryphan-Ar1131 has better photostabilization activity than that of [W⁵] contryphan-Ar1131 and also individual D-tryptophan and L-tryptophan amino acids. These biochemical studies have highlighted the significance of D-tryptophan in contryphan-Ar1131 and its photostabilization activity on avobenzone may find applications in cosmetics.

Conformational selection of vasopressin upon V1a receptor binding

The neuropeptide vasopressin (VP) and its three G protein-coupled receptors (V_1aR, V_1bR and V₂R) are of high interest in a wide array of drug discovery programs. V_1aR is of particular importance due to its cardiovascular functions and diverse roles in the central nervous system. The structure–activity relationships underpinning ligand-receptor interactions remain however largely unclear, hindering rational drug design. This is not least due to the high structural flexibility of VP in its free as well as receptor-bound states. In this work, we developed a novel approach to reveal features of conformational selectivity upon VP-V_1aR complex formation. We employed virtual screening strategies to probe VP’s conformational space for transiently adopted structures that favor binding to V_1aR. To this end, we dissected the VP conformational space into three sub-ensembles, each containing distinct structural sets for VP’s three-residue C-terminal tail. We validated the computational results with experimental nuclear magnetic resonance (NMR) data and docked each sub-ensemble to V_1aR. We observed that the conformation of VP’s three-residue tail significantly modulated the complex dissociation constants. Solvent-exposed and proline trans-configured VP tail conformations bound to the receptor with three-fold enhanced affinities compared to compacted or cis-configured conformations. The solvent-exposed and more flexible structures facilitated unique interaction patterns between VP and V_1aR transmembrane helices 3, 4, and 6 which led to high binding energies. The presented “virtual conformational space screening” approach, integrated with NMR spectroscopy, thus enabled identification and characterization of a conformational selection-type complex formation mechanism that confers novel perspectives on targeting the VP-V_1aR interactions at the level of the encounter complex – an aspect that opens novel research avenues for understanding the functionality of the evolutionary selected conformational properties of VP, as well as guidance for ligand design strategies to provide more potent and selective VP analogues.

Cryo-electron microscopy structure of the antidiuretic hormone arginine-vasopressin V2 receptor signaling complex

The antidiuretic hormone arginine-vasopressin (AVP) forms a signaling complex with the V2 receptor (V2R) and the G_s protein, promoting kidney water reabsorption. Molecular mechanisms underlying activation of this critical G protein-coupled receptor (GPCR) signaling system are still unknown. To fill this gap of knowledge, we report here the cryo-electron microscopy structure of the AVP-V2R-G_s complex. Single-particle analysis revealed the presence of three different states. The two best maps were combined with computational and nuclear magnetic resonance spectroscopy constraints to reconstruct two structures of the ternary complex. These structures differ in AVP and G_s binding modes. They reveal an original receptor-G_s interface in which the Gα_s subunit penetrates deep into the active V2R. The structures help to explain how V2R R137H or R137L/C variants can lead to two severe genetic diseases. Our study provides important structural insights into the function of this clinically relevant GPCR signaling complex.

Aggregation of Aβ40/42 chains in the presence of cyclic neuropeptides investigated by molecular dynamics simulations

Alzheimer’s disease is associated with the formation of toxic aggregates of amyloid beta (Aβ) peptides. Despite tremendous efforts, our understanding of the molecular mechanisms of aggregation, as well as cofactors that might influence it, remains incomplete. The small cyclic neuropeptide somatostatin-14 (SST₁₄) was recently found to be the most selectively enriched protein in human frontal lobe extracts that binds Aβ₄₂ aggregates. Furthermore, SST₁₄’s presence was also found to promote the formation of toxic Aβ₄₂ oligomers in vitro. In order to elucidate how SST₁₄ influences the onset of Aβ oligomerization, we performed all-atom molecular dynamics simulations of model mixtures of Aβ₄₂ or Aβ₄₀ peptides with SST₁₄ molecules and analyzed the structure and dynamics of early-stage aggregates. For comparison we also analyzed the aggregation of Aβ₄₂ in the presence of arginine vasopressin (AVP), a different cyclic neuropeptide. We observed the formation of self-assembled aggregates containing the Aβ chains and small cyclic peptides in all mixtures of Aβ₄₂–SST₁₄, Aβ₄₂–AVP, and Aβ₄₀–SST₁₄. The Aβ₄₂–SST₁₄ mixtures were found to develop compact, dynamically stable, but small aggregates with the highest exposure of hydrophobic residues to the solvent. Differences in the morphology and dynamics of aggregates that comprise SST₁₄ or AVP appear to reflect distinct (1) regions of the Aβ chains they interact with; (2) propensities to engage in hydrogen bonds with Aβ peptides; and (3) solvent exposures of hydrophilic and hydrophobic groups. The presence of SST₁₄ was found to impede aggregation in the Aβ₄₂–SST₁₄ system despite a high hydrophobicity, producing a stronger “sticky surface” effect in the aggregates at the onset of Aβ₄₂–SST₁₄ oligomerization.

Improper folding of proteins causes disorders known as protein misfolding diseases. Under normal conditions most proteins adopt particular folds, which allow them functioning properly. However, for reasons that are not yet fully understood, proteins may misfold and aggregate, forming deposits known as amyloid fibrils, which accumulate in the brain or other tissues. This process affects functioning of the nervous system, gradually causing loss of cognitive abilities. Alzheimer’s disease is one of the most common diseases from this group. A better understanding of the aggregation of peptides implicated in Alzheimer’s disease, known as amyloid beta (Aβ) peptides, may facilitate the development of treatments that ameliorate or prevent the disease. We use detailed molecular dynamics simulations to investigate the influence of somatostatin-14 (SST₁₄), a cyclic neuropeptide that might be involved in the amyloidogenic aggregation of Aβ, on molecular processes occurring at the onset of Aβ aggregation. Results of these simulations explain how the presence of SST₁₄ might alter pathways of aggregation of Aβ, shedding light upon the possible role of extrinsic factors in the aggregation at a molecular level.

Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations

Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.

Mechanistic Insights into an Unusual Side-Chain-Mediated N-Cα Bond Cleavage under Collision-Induced Dissociation Conditions in the Disulfide-Containing Peptide Conopressin

Conopressin, a nonapeptide disulfide CFIRNCPKG amide present in cone snail venom, undergoes a facile cleavage at the Cys6-Pro7 peptide bond to yield a disulfide bridged b₆ ion. Analysis of the mass spectral fragmentation pattern reveals the presence of a major fragment ion, which is unambiguously assigned as the tripeptide sequence IRN amide. The sequence dependence of this unusual fragmentation process has been investigated by comparing it with the fragmentation patterns of related peptides, oxytocin (CYIQNCPLG amide), Lys-vasopressin (CYFQNCPKG amide), and a series of synthetic analogues. The results establish the role of the Arg4 residue in facilitating the unusual N-C^α bond cleavage at Cys6. Structures are proposed for a modified disulfide bridged fragment containing the Cys1 and Cys6 residues. Gas-phase molecular dynamics simulations provide evidence for the occurrence of conformational states that permit close approach of the Arg4 side chain to the Cys6 C^β methylene protons.

Understanding Russell's viper venom factor V activator's substrate specificity by surface plasmon resonance and in-silico studies

Blood coagulation factor V (FV) is activated either by Factor X or thrombin, cleaving at three different sites viz., Site I (Arg709-Ser710), site II (Arg1018-Thr1019), and site III (Arg1545-Ser1546). Russell's viper venom factor V activator (RVV-V) is a thrombin-like serine proteinase that activates FV with selective, single cleavage at site III. A long lasting effort is being pending in understanding the 'selective' binding specificity of the RVV-V towards site III. Here, we present the binding kinetic study of RVV-V with two designed peptides corresponding to the regions from site I (Gln699-Asn713) and site II (1008Lys-Pro1022), respectively, that include 15 amino acids. Our investigation for justifying the binding efficacy and kinetics of peptides includes SPR method, protein-peptide docking, molecular dynamics simulation, and principal component analysis (PCA). Surprisingly, the SPR experiment disclosed that the Peptide II showed a lower binding affinity with KD of 2.775 mM while the Peptide I showed none. Docking and simulation of both the peptides with RVV-V engaged either rooted or shallow binding for Peptide II and Peptide I respectively. The peptide binding resulted in global conformational changes in the native fold of RVV-V, whereas the similar studies for thrombin failed to make major changes in the native fold. In support, the PCA analysis for RVV-V showed the dislocation of catalytic triad upon binding both the peptides. Hence, RVV-V, a serine protease, is incompetent in cleaving these two sites. This study suggests a transition in RVV-V from the native rigid to the distorted flexible structure and paves a way to design a new peptide substrate/inhibitor.

Structural Elements in the Gαs and Gαq C Termini That Mediate Selective G Protein-coupled Receptor (GPCR) Signaling

Although the importance of the C terminus of the α subunit of the heterotrimeric G protein in G protein-coupled receptor (GPCR)-G protein pairing is well established, the structural basis of selective interactions remains unknown. Here, we combine live cell FRET-based measurements and molecular dynamics simulations of the interaction between the GPCR and a peptide derived from the C terminus of the Gα subunit (Gα peptide) to dissect the molecular mechanisms of G protein selectivity. We observe a direct link between Gα peptide binding and stabilization of the GPCR conformational ensemble. We find that cognate and non-cognate Gα peptides show deep and shallow binding, respectively, and in distinct orientations within the GPCR. Binding of the cognate Gα peptide stabilizes the agonist-bound GPCR conformational ensemble resulting in favorable binding energy and lower flexibility of the agonist-GPCR pair. We identify three hot spot residues (Gαs/Gαq-Gln-384/Leu-349, Gln-390/Glu-355, and Glu-392/Asn-357) that contribute to selective interactions between the β2-adrenergic receptor (β2-AR)-Gαs and V1A receptor (V1AR)-Gαq The Gαs and Gαq peptides adopt different orientations in β2-AR and V1AR, respectively. The β2-AR/Gαs peptide interface is dominated by electrostatic interactions, whereas the V1AR/Gαq peptide interactions are predominantly hydrophobic. Interestingly, our study reveals a role for both favorable and unfavorable interactions in G protein selection. Residue Glu-355 in Gαq prevents this peptide from interacting strongly with β2-AR. Mutagenesis to the Gαs counterpart (E355Q) imparts a cognate-like interaction. Overall, our study highlights the synergy in molecular dynamics and FRET-based approaches to dissect the structural basis of selective G protein interactions.

Probing the Binding Interfaces of Protein Complexes Using Gas-Phase H/D Exchange Mass Spectrometry

Fast gas-phase hydrogen/deuterium exchange mediated by ND3 gas and measured by mass spectrometry (gas-phase HDX-MS) is a largely unharnessed, fast, and sensitive method for probing primary- and higher-order polypeptide structure. Labeling of heteroatom-bound non-amide hydrogens in a sub-millisecond time span after electrospray ionization by ND3 gas can provide structural insights into protein conformers present in solution. Here, we have explored the use of gas-phase HDX-MS for probing the higher-order structure and binding interfaces of protein complexes originating from native solution conditions. Lysozyme ions bound by an oligosaccharide incorporated less deuterium than the unbound ion. Similarly, trypsin ions showed reduced deuterium uptake when bound by the peptide ligand vasopressin. Our results are in good agreement with crystal structures of the native protein complexes, and illustrate that gas-phase HDX-MS can provide a sensitive and simple approach to measure the number of heteroatom-bound non-amide side-chain hydrogens involved in the binding interface of biologically relevant protein complexes.

Conservation weighting functions enable covariance analyses to detect functionally important amino acids

The explosive growth in the number of protein sequences gives rise to the possibility of using the natural variation in sequences of homologous proteins to find residues that control different protein phenotypes. Because in many cases different phenotypes are each controlled by a group of residues, the mutations that separate one version of a phenotype from another will be correlated. Here we incorporate biological knowledge about protein phenotypes and their variability in the sequence alignment of interest into algorithms that detect correlated mutations, improving their ability to detect the residues that control those phenotypes. We demonstrate the power of this approach using simulations and recent experimental data. Applying these principles to the protein families encoded by Dscam and Protocadherin allows us to make testable predictions about the residues that dictate the specificity of molecular interactions.

Disordered structural ensembles of vasopressin and oxytocin and their mutants

Vasopressin and oxytocin are intrinsically disordered cyclic nonapeptides belonging to a family of neurohypophysial hormones. Although unique in their functions, these peptides differ only by two residues and both feature a tocin ring formed by the disulfide bridge between first and sixth cysteine residues. This sequence and structural similarity are experimentally linked to oxytocin agonism at vasopressin receptors and vasopressin antagonism at oxytocin receptors. Yet single- or double-residue mutations in both peptides have been shown to have drastic impacts on their activities at either receptor, and possibly the ability to bind to their neurophysin carrier protein. In this study we perform molecular dynamics simulations of the unbound native and mutant sequences of the oxytocin and vasopressin hormones to characterize their structural ensembles. We classify the subpopulations of these structural ensembles on the basis of the distributions of radius of gyration and secondary structure and hydrogen-bonding features of the canonical tocin ring and disordered tail region. We then relate the structural changes observed in the unbound form of the different hormone sequences to experimental information about peptide receptor binding, and more indirectly, carrier protein binding affinity, receptor activity, and protease degradation. This study supports the hypothesis that the structural characteristics of the unbound form of an IDP can be used to predict structural or functional preferences of its functional bound form.

C4ORF48, a gene from the Wolf-Hirschhorn syndrome critical region, encodes a putative neuropeptide and is expressed during neocortex and cerebellar development

In order to identify novel genes involved in mental retardation/intellectual disability, we focused on a microdeletion reported in a patient with a mild form of Wolf-Hirschhorn syndrome. This patient presented with attention-deficit hyperactivity disorder, some learning and fine motor deficits as well as facial abnormalities. The deleted region included three genes. Here, we report the first characterization of one of these genes, C4ORF48. C4ORF48 encodes a short (139 aa) evolutionarily conserved protein with a predicted signal peptide and two potential dibasic convertase cleavage sites. In mice, we demonstrated expression of the corresponding protein exclusively in brain tissue using an anti-mouse C4Orf48 polyclonal antibody. Detailed RNA in situ hybridization experiments revealed expression of C4Orf48 in different zones during cortical and cerebellar development, as well as in almost all cortical and subcortical regions of the adult mouse brain. Based on the present data, we propose that C4Orf48 probably encodes a novel neuropeptide, which, if hemizygously deleted, may be involved in the observed intellectual and fine motor disabilities and thus in the overall neurological aspects of Wolf-Hirschhorn syndrome.

Structural basis of specificity of a peptidyl urokinase inhibitor, upain-1

Urokinase-type plasminogen activator (uPA) plays a crucial role in the regulation of plasminogen activation, tumor cell adhesion and migration. The inhibition of uPA activity is a promising mechanism for anti-cancer therapy. A cyclic peptidyl inhibitor, upain-1, CSWRGLENHRMC, was identified recently as a competitive and highly specific uPA inhibitor. We determined the crystal structure of uPA in complex with upain-1 at 2.15 A. The structure reveals that the cyclic peptide adopts a rigid conformation stabilized by a disulfide bond (residues 1-12) and three tight beta turns (residues 3-6, 6-9, 9-12). The Glu7 residue of upain-1 forms hydrogen bonds with the main chain nitrogen atoms of residues 4, 5, and 6 of upain-1, and is also critical for maintaining the active conformation of upain-1. The Arg4 of upain-1 is inserted into the uPA's specific S1 pocket. The Ser2 residue of upain-1 locates close to the S1beta pocket of uPA. The Gly5 and Glu7 residues of upain-1 occupy the S2 pocket and the oxyanion hole of uPA, respectively. Furthermore, the Asn8 residue of upain-1 binds to the 37- and 60-loops of uPA and renders the specificity of upain-1 for uPA. Based on this structure, a new pharmacophore for the design of highly specific uPA inhibitors was proposed.

A vasopressin/oxytocin-related conopeptide with gamma-carboxyglutamate at position 8

Vasopressins and oxytocins are homologous, ubiquitous and multifunctional peptides present in animals. Conopressins are vasopressin/oxytocin-related peptides that have been found in the venom of cone snails, a genus of marine predatory molluscs that envenom their prey with a complex mixture of neuroactive peptides. In the present paper, we report the purification and characterization of a unique conopressin isolated from the venom of Conus villepinii, a vermivorous cone snail species from the western Atlantic Ocean. This novel peptide, designated gamma-conopressin-vil, has the sequence CLIQDCPgammaG* (gamma is gamma-carboxyglutamate and * is C-terminal amidation). The unique feature of this vasopressin/oxytocin-like peptide is that the eighth residue is gamma-carboxyglutamate instead of a neutral or basic residue; therefore it could not be directly classified into either the vasopressin or the oxytocin peptide families. Nano-NMR spectroscopy of the peptide isolated directly from the cone snails revealed that the native gamma-conopressin-vil undergoes structural changes in the presence of calcium. This suggests that the peptide binds calcium, and the calcium-binding process is mediated by the gamma-carboxyglutamate residue. However, the negatively charged residues in the sequence of gamma-conopressin-vil may mediate calcium binding by a novel mechanism not observed in other peptides of this family.