Right- and left-handed three-helix proteins. II. Similarity and differences in mechanical unfolding of proteins

Molecular Dynamics Simulation Study of the Self-Assembly of Phenylalanine Peptide Nanotubes

In this paper, we propose and use a new approach for a relatively simple technique for conducting MD simulation (MDS) of various molecular nanostructures, determining the trajectory of the MD run and forming the final structure using external force actions. A molecular dynamics manipulator (MD manipulator) is a controlled MDS type. As an example, the applicability of the developed algorithm for assembling peptide nanotubes (PNT) from linear phenylalanine (F or Phe) chains of different chirality is presented. The most adequate regimes for the formation of nanotubes of right chirality D from the initial L-F and nanotubes of left chirality L of their initial dipeptides D-F modes were determined. We use the method of a mixed (vector–scalar) product of the vectors of the sequence of dipole moments of phenylalanine molecules located along the nanotube helix to calculate the magnitude and sign of chirality of self-assembled helical phenylalanine nanotubes, which shows the validity of the proposed approach. As result, all data obtained correspond to the regularity of the chirality sign change of the molecular structures with a hierarchical complication of their organization.

Role of Resultant Dipole Moment in Mechanical Dissociation of Biological Complexes

Protein-peptide interactions play essential roles in many cellular processes and their structural characterization is the major focus of current experimental and theoretical research. Two decades ago, it was proposed to employ the steered molecular dynamics (SMD) to assess the strength of protein-peptide interactions. The idea behind using SMD simulations is that the mechanical stability can be used as a promising and an efficient alternative to computationally highly demanding estimation of binding affinity. However, mechanical stability defined as a peak in force-extension profile depends on the choice of the pulling direction. Here we propose an uncommon choice of the pulling direction along resultant dipole moment (RDM) vector, which has not been explored in SMD simulations so far. Using explicit solvent all-atom MD simulations, we apply SMD technique to probe mechanical resistance of ligand-receptor system pulled along two different vectors. A novel pulling direction—when ligand unbinds along the RDM vector—results in stronger forces compared to commonly used ligand unbinding along center of masses vector. Our observation that RDM is one of the factors influencing the mechanical stability of protein-peptide complex can be used to improve the ranking of binding affinities by using mechanical stability as an effective scoring function.

Kinetics and mechanical stability of the fibril state control fibril formation time of polypeptide chains: A computational study

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aβ₄₀ and Aβ₄₂ and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics.

Comparative mechanical unfolding studies of spectrin domains R15, R16 and R17

Spectrins belong to repetitive three-helix bundle proteins that have vital functions in multicellular organisms and are of potential value in nanotechnology. To reveal the unique physical features of repeat proteins we have studied the structural and mechanical properties of three repeats of chicken brain α-spectrin (R15, R16 and R17) at the atomic level under stretching at constant velocities (0.01, 0.05 and 0.1 Å·ps^-1) and constant forces (700 and 900 pN) using molecular dynamics (MD) simulations at T = 300 K. 114 independent MD simulations were performed and their analysis has been done. Despite structural similarity of these domains we have found that R15 is less mechanically stable than R16, which is less stable than R17. This result is in agreement with the thermal unfolding rates. Moreover, we have observed the relationship between mechanical stability, flexibility of the domains and the number of aromatic residues involved in aromatic clusters.

What handedness and angles between helices has the studied three-helical protein domain?

We have created a new server FoldHandedness. Using this server it is possible: (i) to define the regions of helices from two issues (from the PDB file and using the last version of the DSSP program), (ii) to determine the handedness for any chosen three helices and (iii) to calculate the angle and sign between the chosen pairs of the helices for large proteins and complexes of proteins with DNA or RNA.