Mutational Analysis of Neuropeptide Y Reveals Unusual Thermal Stability Linked to Higher-Order Self-Association

Tertiary Plasticity Drives the Efficiency of Enterocin 7B Interactions with Lipid Membranes

The ability of antimicrobial peptides to efficiently kill their bacterial targets depends on the efficiency of their binding to the microbial membrane. In the case of enterocins, there is a three-part interaction: initial binding, unpacking of helices on the membrane surface, and permeation of the lipid bilayer. Helical unpacking is driven by disruption of the peptide hydrophobic core when in contact with membranes. Enterocin 7B is a leaderless enterocin antimicrobial peptide produced from Enterococcus faecalis that functions alone, or with its cognate partner enterocin 7A, to efficiently kill a wide variety of Gram-stain positive bacteria. To better characterize the role that tertiary structural plasticity plays in the ability of enterocin 7B to interact with the membranes, a series of arginine single-site mutants were constructed that destabilize the hydrophobic core to varying degrees. A series of experimental measures of structure, stability, and function, including CD spectra, far UV CD melting profiles, minimal inhibitory concentrations analysis, and release kinetics of calcein, show that decreased stabilization of the hydrophobic core is correlated with increased efficiency of a peptide to permeate membranes and in killing bacteria. Finally, using the computational technique of adaptive steered molecular dynamics, we found that the atomistic/energetic landscape of peptide mechanical unfolding leads to free energy differences between the wild type and its mutants, whose trends correlate well with our experiment.

Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules

Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.

Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics

Long-time dynamical processes, such as those involving protein unfolding and ligand interactions, can be accelerated and realized through steered molecular dynamics (SMD). The challenge has been the extraction of information from such simulations that generalize for complex nonequilibrium processes. The use of Jarzynski's equality opened the possibility of determining the free energy along the steered coordinate, but sampling over the nonequilibrium trajectories is slow to converge. Adaptive steered molecular dynamics (ASMD) and other related techniques have been introduced to overcome this challenge through the use of stages. Here, we take advantage of these stages to address the numerical cost that arises from the required use of very large solvent boxes. We introduce telescoping box schemes within adaptive steered molecular dynamics (ASMD) in which we adjust the solvent box between stages and thereby vary (and optimize) the required number of solvent molecules. We have benchmarked the method on a relatively long α-helical peptide, Ala₃₀, with respect to the potential of mean force and hydrogen bonds. We show that the use of telescoping boxes introduces little numerical error while significantly reducing the computational cost.

Benchmarking Adaptive Steered Molecular Dynamics (ASMD) on CHARMM Force Fields

The potentials of mean force (PMFs) along the end‐to‐end distance of two different helical peptides have been obtained and benchmarked using the adaptive steered molecular dynamics (ASMD) method. The results depend strongly on the choice of force field driving the underlying all‐atom molecular dynamics, and are reported with respect to the three most popular CHARMM force field versions: c22, c27 and c36. Two small peptides, ALA10 and 1PEF, serve as the particular case studies. The comparisons between the versions of the CHARMM force fields provides both a qualitative and quantitative look at their performance in forced unfolding simulations in which peptides undergo large changes in structural conformations. We find that ASMD with the underlying c36 force field provides the most robust results for the selected benchmark peptides.

Development of a Therapeutic Peptide for Cachexia Suggests a Platform Approach for Drug-like Peptides

During the development of a melanocortin (MC) peptide drug to treat the condition of cachexia (a hypermetabolic state producing lean body mass wasting), we were confronted with the need for peptide transport across the blood-brain barrier (BBB): the MC-4 receptors (MC4Rs) for metabolic rate control are located in the hypothalamus, i.e., behind the BBB. Using the term "peptides with BBB transport", we screened the medical literature like a peptide library. This revealed numerous "hits"-peptides with BBB transport and/or oral activity. We noted several features common to most peptides in this class, including a dipeptide sequence of nonpolar residues, primary structure cyclization (whole or partial), and a Pro-aromatic motif usually within the cyclized region. Based on this, we designed an MC4R antagonist peptide, TCMCB07, that successfully treated many forms of cachexia. As part of our pharmacokinetic characterization of TCMCB07, we discovered that hepatobiliary extraction from blood accounted for a majority of the circulating peptide's excretion. Further screening of the literature revealed that TCMCB07 is a member of a long-forgotten peptide class, showing active transport by a multi-specific bile salt carrier. Bile salt transport peptides have predictable pharmacokinetics, including BBB transport, but rapid hepatic clearance inhibited their development as drugs. TCMCB07 shares the general characteristics of the bile salt peptide class but with a much longer half-life of hours, not minutes. A change in its C-terminal amino acid sequence slows hepatic clearance. This modification is transferable to other peptides in this class, suggesting a platform approach for producing drug-like peptides.

Theoretical study of the interactions between peptide tyrosine tyrosine [PYY (1-36)], a newly identified modulator in type 2 diabetes pathophysiology, with receptors NPY1R and NPY4R

PURPOSE

Diabetes mellitus is a common condition in the clinically obese. Bariatric surgery is one of the ways to put type 2 diabetes in remission. Recent findings propose the appetite-regulator peptide tyrosine tyrosine (PYY) as a therapeutic option for patients with type 2 diabetes. This novel gut hormone restores impaired insulin and glucagon secretion in pancreatic islets and is implicated in type 2 diabetes reversal after bariatric surgery. The current study elucidates the interactions between PYY and the NPY1R and NPY4R receptors using computational methods.

METHODS

Protein structure prediction, molecular docking simulation, and molecular dynamics (MD) simulation were performed to elucidate the interactions of PYY with NPY1R and NPY4R.

RESULTS

The predicted binding models of PYY-NPY receptors are in agreement with those described in the literature, although different interaction partners are presented for the C-terminal tail of PYY. Non-polar interactions are predicted to drive the formation of the protein complex. The calculated binding energies show that PYY has higher affinity for NPY4R (ΔG_GBSA = -65.08 and ΔG_PBSA = -87.62 kcal/mol) than for NPY1R (ΔG_GBSA = -23.11 and ΔG_PBSA = -50.56 kcal/mol).

CONCLUSIONS

Based on the constructed models, the binding conformations obtained from docking and MD simulation for both the PYY-NPY1R and PYY-NPY4R complexes provide a detailed map of possible interactions. The calculated binding energies show a higher affinity of PYY for NPY4R. These findings may help to understand the mechanisms behind the improvement of diabetes following bariatric surgery.

Energetics and structure of alanine-rich α-helices via adaptive steered molecular dynamics

The energetics and hydrogen bonding profiles of the helix-to-coil transition were found to be an additive property and to increase linearly with chain length, respectively, in alanine-rich α-helical peptides. A model system of polyalanine repeats was used to establish this hypothesis for the energetic trends and hydrogen bonding profiles. Numerical measurements of a synthesized polypeptide Ac-Y(AEAAKA)kF-NH2 and a natural α-helical peptide a2N (1-17) provide evidence of the hypothesis's generality. Adaptive steered molecular dynamics was employed to investigate the mechanical unfolding of all of these alanine-rich polypeptides. We found that the helix-to-coil transition is primarily dependent on the breaking of the intramolecular backbone hydrogen bonds and independent of specific side-chain interactions and chain length. The mechanical unfolding of the α-helical peptides results in a turnover mechanism in which a 310-helical structure forms during the unfolding, remaining at a near constant population and thereby maintaining additivity in the free energy. The intermediate partially unfolded structures exhibited polyproline II helical structure as previously seen by others. In summary, we found that the average force required to pull alanine-rich α-helical peptides in between the endpoints-namely the native structure and free coil-is nearly independent of the length or the specific primary structure.

Research on Correlation Between Psychological Factors, Mast Cells, and PAR-2 Signal Pathway in Irritable Bowel syndrome

PURPOSE

The study aimed to explore the level of psychological stress factors, mast cell (MC), and neuropeptide in the occurrence of irritable bowel syndrome (IBS) and the correlation among them, and to identify representative and effective indicators for the pathogenesis and clinical medication development of IBS.

SUBJECTS AND METHODS

Twenty-eight patients (20-64 years old) with diarrhea-predominant IBS (IBS-D) were included in the IBS-D group and 8 healthy subjects (35-63 years old) were enrolled in the control group. All subjects accepted the colonoscopic biopsies, self-rating depression scale (SDS) and self-rating anxiety scale (SAS) assessment. MC degranulation, neuropeptide S (NPS), neuropeptide Y (NPY), NPY receptor 2 (NPY2R) and Protease-activated receptor 2 (PAR-2) in colon tissues were performed by Strept Avidin-Biot complex (SABC) immunohistochemistry. Enzyme-linked immunosorbent assay (ELISA) detection was used to test the expression of NPS and NPY in peripheral blood plasma and colon tissues. Western blot was applied to examine the level of NPY2R and PAR-2.

RESULTS

The level of anxiety and depression of patients with IBS-D was more serious than that in the control. The expression of NPS, NPY and NPY2R was down-regulated in the IBS-D. The total MC and tryptase-positive MC increased significantly in the colon tissue of IBS-D and the expression level of PAR-2 was significantly up-regulated.

CONCLUSION

There has been a close connection among those indicators that the activated MC may up-regulate the function of PAR-2, resulting in the change of neuropeptide (NPS and NPY), successively leading to clinical symptoms and psychological negative changes in the IBS.

The relative stability of trpzip1 and its mutants determined by computation and experiment

The single-point mutations of tprzip1 are indicated at left, and their relative energetics are compared at right.