Understanding the Conformational Properties of Fluorinated Polypeptides: Molecular Modelling of Unguisin A

Rational Discovery of Microtubule-Stabilizing Peptides

Microtubule (MT) stabilization is an attractive pharmacological strategy to hamper the progress of neurodegenerative diseases. In this regard, seeking peptides with MT-stabilizing properties has awoken great interest. This work reports the rational discovery of two structurally related MT-stabilizing octapeptides using a combination of protein-peptide docking, conventional molecular dynamics, Gaussian accelerated molecular dynamics (GaMD), and tubulin polymerization assays. FASTA sequences for ∼1000 peptides were crafted from single and double mutants of davunetide (NAP) and docked against the Taxol (TX) site on an octameric MT model representing a portion of the MT wall. Docked peptides were rescored after MM minimization and binding free energy refinement through single-point MM/GBSA calculations. The 60 best-ranked peptides were subjected to 50 ns MD simulations on peptide-MT complexes at the terminal TX site in the octameric Tau-MT model resulting in 11 complexes with occupancies greater than 99% and peptide-protein binding free energies less than -40 kcal/mol. Selected peptides were then examined through 300 ns GaMD simulations in complexes containing two identical ligands at the terminal and intermediate TX sites in the Tau-MT model to account for the differential association of MT-binding peptides to different regions of the MT structure. Six candidates showed a favorable MT-binding potential based on the analysis of interaction frequencies and relative mobilities of the complex components, suggesting a pivotal role of Arg278, Gln281, and Arg369 residues for peptides recognition. Four candidates were predicted to preserve an adequate balance of longitudinal and lateral interactions between tubulin dimers in peptide-MT complexes such that MT-stabilizing effects could be expected. MT polymerization experiments confirmed that four peptides (HAPVSIHQ, NYPVSIHQ, NWPVSIWQ, HAPVSIIQ) exhibit MT-stabilizing activity in vitro with NWPVSIWQ (P43) and HAPVSIIQ (P52) being the most active. Tryptophan quenching assays verified that P43 and P52 bind to nonpolymeric tubulin, whereas viability experiments on HEK cells confirmed their safety to pursue future pharmacological studies. The results herein presented are valuable to making progress in the rational design of MT-stabilizing peptides.

Toward Reliable and Insightful Entropy Calculations on Flexible Molecules

The absolute entropy of a flexible molecule can be approximated by the sum of a rigid-rotor-harmonic-oscillator (RRHO) entropy and a Gibbs-Shannon entropy associated to the Boltzmann distribution for the occupation of the conformational energy levels. Herein, we show that such partitioning, which has received renewed interest, leads to accurate entropies of single molecules of increasing size provided that the conformational part is estimated by means of a set of discretization and expansion techniques that are able to capture the significant correlation effects among the torsional motions. To ensure a reliable entropy estimation, we rely on extensive sampling as that produced by classical molecular dynamics simulations on the microsecond time scale, which is currently affordable for small- and medium-sized molecules. According to test calculations, the gas-phase entropy of simple organic molecules is predicted with a mean unsigned error of 0.9 cal/(mol K) when the RRHO entropies are computed at the B3LYP-D3/cc-pVTZ level. Remarkably, the same protocol gives small errors [<1 cal/(mol K)] for the extremely flexible linear alkane molecules (C_nH_2n+2, n = 14, 16, and 18). Similarly, we obtain well-converged entropies for a more challenging test of drug molecules, which exhibit more pronounced correlation effects. We also perform equivalent entropy calculations on a 76 amino acid protein, ubiquitin, by taking advantage of the cutoff-dependent formulation of an expansion technique (correlation-consistent multibody local approximation, CC-MLA), which incorporates genuine correlation effects among the neighboring dihedral angles. Moreover, we show that insightful descriptors of the coupled torsional motions can be obtained with the CC-MLA approach.