Understanding peptide competitive inhibition of botulinum neurotoxin a binding to SV2 protein via molecular dynamics simulations

Botulinum neurotoxins (BoNTs) are known as the most toxic natural substances. Synaptic vesicle protein 2 (SV2) has been proposed to be a protein receptor for BoNT/A. Recently, two short peptides (BoNT/A‐A2 and SV2C‐A3) were designed to inhibit complex formation between the BoNT/A receptor‐binding domain (BoNT/A‐RBD) and the synaptic vesicle protein 2C luminal domain (SV2C‐LD). In this article, the two peptide complex systems are studied by molecular dynamics (MD) simulations. The structural stability analysis indicates that BoNT/A‐A2 system is more stable than SV2C‐A3 system. The conformational analysis implies that the β‐sheet in BoNT/A‐A2 system maintains its secondary structure but the two β‐strands in SV2C‐A3 system have remarkable conformational changes. Based on the calculation of hydrogen bonds, hydrophobic interactions and cation‐π interactions, it is found that the internal hydrogen bonds play crucial roles in the structural stability of the peptides. Because of the stable secondary structure, the β‐sheet in BoNT/A‐A2 system establishes effective interactions at the interface and inhibits BoNT/A‐RBD binding to SV2C‐LD. In contrast, without other β‐strands forming internal hydrogen bonds, the two isolated β‐strands in SV2C‐A3 system become the random coil. This conformational change breaks important hydrogen bonds and weakens cation‐π interaction in the interface, so the complex formation is only partially inhibited by the two β‐strands. These results are consistent with experimental studies and may be helpful in understanding the inhibition mechanisms of peptide inhibitors. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 597–608, 2015.     

Molecular dynamics simulations of a double unit cell in a protein crystal: volume relaxation at constant pressure and correlation of motions between the two unit cells

Eight molecular dynamics simulations of a double crystal unit cell of ubiquitin were performed to investigate the effects of simulating at constant pressure and of simulating two unit cells compared to a single unit cell. To examine the influence of different simulation conditions, the constant-pressure and constant-volume simulations were each performed with and without counterions and using two different treatments of the long-range electrostatic interactions (lattice-sum and reaction-field methods). The constant-pressure simulations were analyzed in terms of unit cell deformation and accompanying protein deformations. Energetic and structural properties of the proteins in the simulations of the double unit cell were compared to the results of previously reported one-unit-cell simulations. Correlation between the two unit cells was also investigated based on relative translational and rotational movements of the proteins and on dipole fluctuations. The box in the constant-pressure simulations is found to deform slowly to reach convergence only after 5-10 ns. This deformation does not result from a distortion in the structure of the proteins but rather from changes in protein packing within the unit cell. The results of the double-unit-cell simulations are closely similar to the results of the single-unit-cell simulations, and little motional correlation is found between the two unit cells.     

Molecular dynamics simulation study of AG10 and tafamidis binding to the Val122Ile transthyretin variant

Molecular dynamics (MD) simulations were used to investigate the binding of four ligands to the Val122Ile mutant of the protein transthyretin. Dissociation, misfolding, and subsequent aggregation of mutated transthyretin proteins are associated with the disease Familial Amyloidal Cardiomyopathy. The ligands investigated were the drug candidate AG10 and its decarboxy and N-methyl derivatives along with the drug tafamidis. These ligands bound to the receptor in two halogen binding pockets (HBP) designated AB and A’B’. Inter-ligand distances, solvent accessible surface areas, root mean squared deviation measurements, and extracted structures showed very little change in the AG10 ligands' conformations or locations within the HBP during the MD simulation. In addition, the AG10 ligands experienced stable, two-point interactions with the protein by forming hydrogen bonds with Ser-117 residues in both the AB and A’B’ binding pockets and Lysine-15 residues found near the surface of the receptor. Distance measurements showed these H-bonds formed simultaneously during the MD simulation. Removal of the AG10 carboxylate functional group to form decarboxy-AG10 disrupted this two-point interaction causing the ligand in the AB pocket to undergo a conformational change during the MD simulation. Likewise, addition of a methyl group to the AG10 hydrazone functional group also disrupted the two-point interaction by decreasing hydrogen bonding interactions with the receptor. Finally, MD simulations showed that the tafamidis ligands experienced fewer hydrogen bonding interactions than AG10 with the protein receptor. The tafamidis ligand in pocket A’B’ was also found to move deeper into the HBP during the MD simulation.