Design and analysis of polypeptide nanofiber using full atomistic Molecular Dynamic

Exploring How System Dimensions and Periodic Boundary Conditions Influence the Molecular Dynamics Simulation of A6H Peptide Self-Assembly Nanostructures

This work presents a study on the effects of periodic boundary conditions (PBC) on the energetic/structural properties and hydrogen bond dynamics (HB) using molecular dynamics (MD) simulations of peptide membranes composed of alanine and histidine. Our results highlight that simulations using small surface areas for the peptide membrane may result in nonconvergent values for membrane properties, which are only observed in regions simulated at a certain distance from the PBCs. Specifically, regarding hydrogen bonds, a property pervasive in peptide membranes, our findings indicate a significant increase in the lifetime of these interactions, reaching values ∼19% higher when observed in structures free from PBCs. For peptide mobility in these nanomembranes, our results compare regions simulated directly under the influence of PBCs with regions free from these conditions, emphasizing greater mobility of amino acid psi/phi angles in the latter model.

Polar Zipper on a Peptide Nanomembrane: A Characterization by Potential of Mean Force

In this work, nanomembranes formed by the I₃XGK (X = Q, S, or N) polar peptides are studied to characterize the average force and energy required to separate two neighboring β-sheets laterally joined by polar zippers. The results presented are obtained using a methodology (state of the art) involving the pulling umbrella method to generate the samples (umbrella sampling) and the potential of mean force (PMF) to evaluate the energetic variation evolved in the process of separating the polar zipper. It was observed that the maximum force required to separate the regions linked by polar zippers is 1.48 kJ/mol nm for the I₃NGK nanomembrane and 1.22 kJ/mol nm [1.30 kJ/mol nm] for the I₃QGK [I₃SGK] nanomembranes, emphasizing that polar zippers play an important role in the interaction that interconnects β-sheets in broad and robust two-dimensional structures (tapes and membranes), offering an agile route to the construction of distinct nanomaterials from β-sheets. Also, negative values were obtained for energy as a function of the reaction coordinate for the regions where the formation of polar zippers occurs, showing that the lateral union of neighboring β-sheets is energetically favorable, with a value up to -3.0 kJ/mol, in the case of I₃NGK nanomembranes.

Updating atomic charge parameters of aliphatic amino acids: a quest to improve the performance of molecular modeling via sequential molecular dynamics and DFT-GIAO-NMR calculations

In this work, we observe the behavior of the dipole moment, atomic charges, solute-solvent interactions and NMR spectroscopy of aliphatic amino acids in a water solution via the computational simulations of classical molecular dynamics and DFT quantum calculations. Our results indicate that the convergence of the atomic charge of the solute, from an iterative process, together with the dipole moment of the amino acid, alters the lifetime of hydrogen bonds present in the first solvation shell, resulting in the modification of its structure and dynamics. Using GIAO-DFT-NMR calculations, we assessed the impact of these structural solute-solvent modifications on the magnetic shielding constants of the solute carbon atoms. In this sense, we evaluate the importance of an update in parameters that describe atomic charges present in the CHARMM36 force field.