Thermodynamic resolution: how do errors in modeled protein structures affect binding affinity predictions?

A Study on the Binding Mechanism and the Impact of Key Residue Mutations between SND1 and MTDH Peptide through Molecular Dynamics Simulations

Metastasis of breast cancer is the main cause of death for patients with breast cancer. The interaction between metadherin (MTDH) and staphylococcal nuclease domain 1 (SND1) plays a pivotal role in promoting breast cancer development. However, the binding details between MTDH and SND1 remain unclear. In this study, we employed all-atom molecular dynamics simulations (MDs) and conducted binding energy calculations to investigate the binding details and the impact of key residue mutations on binding. The mutations in key residues have not significantly affected the overall stability of the structure and the fluctuation of residues near the binding site; they have exerted a substantial impact on the binding of SND1 and MTDH peptide. The electrostatic interactions and van der Waals interactions play an important role in the binding of SND1 and the MTDH peptide. The mutations in the key residues have a significant impact on electrostatic and van der Waals interactions, resulting in weakened binding. The energy contributions of key residues mainly come from the electrostatic energy and van der Waals interactions of the side chain. In addition, the key residues form an intricate and stable network of hydrogen bonds and salt-bridge interactions with the MTDH peptide. The mutations in key residues have directly disrupt the interactions formed between SND1 and MTDH peptide, consequently leading to changes in the binding mode of the MTDH peptide. These analyses unveil the detailed atomic-level interaction mechanism between SND1 and the MTDH peptide, providing a molecular foundation for the development of antibreast cancer drugs.

Exploring binding modes of the selected inhibitors to SND1 by all-atom molecular dynamics simulations

Breast cancer is the leading cause of cancer-related deaths in women. Previous studies have indicated that disrupting the interaction between Metadherin (MTDH) and Staphylococcal nuclease domain containing 1 (SND1) can inhibit breast cancer development. Understanding the binding mode of small molecule inhibitors with SND1 is of great significance for designing drugs targeting the MTDH-SND1 complex. In this study, we conducted all-atom molecular dynamics (MD) simulations in solution and performed binding energy calculations to gain insights into the binding mechanism of small molecules to SND1. The binding site of SND1 for small molecules is relatively rigid, and the binding of the small molecule and the mutation of key residues have little effect on the conformation of the binding site. SND1 binds more tightly to C26-A6 than to C26-A2, as C26-A2 undergoes a 180° directional change during the simulation process. The key residue mutations have a direct effect on the position and orientation of small molecule in the binding site. The key residues make primary contributions to the binding energy through van der Waals interaction and nonpolar solvation energy, although the contribution from nonpolar solvation is relatively minor. The key residue mutations also affect the formation of hydrogen bonds and ultimately the stability of the small molecule-SND1 complex.Communicated by Ramaswamy H. Sarma.

Revealing the binding mechanism of BACE1 inhibitors through molecular dynamics simulations

Alzheimer's disease is a debilitating neurodegenerative disorder, and the Beta-Site Amyloid Precursor Protein Cleaving Enzyme 1 (BACE1) is a key therapeutic target in its treatment. This study employs molecular dynamics simulations and binding energy analysis to investigate the binding interactions between BACE1 and four selected small molecules: CNP520, D9W, NB641, and NB360. The binding model analysis indicates that the binding of BACE1 with four molecules are stable, except the loop regions show significant fluctuation. The binding free energy analyses reveal that NB360 exhibits the highest binding affinity with BACE1, surpassing other molecules (CNP520, D9W, and NB641). Detailed energy component assessments highlight the critical roles of electrostatic interactions and van der Waals forces in the binding process. Furthermore, residue contribution analysis identifies key amino acids influencing the binding process across all systems. Hydrogen bond analysis reveals a limited number of bonds between BACE1 and each small molecule, highlighting the importance of structural modifications to enable more stable hydrogen bonds. This research provides valuable insights into the molecular mechanisms of potential Alzheimer's disease therapeutics, guiding the way for improved drug design and the development of effective treatments targeting BACE1.Communicated by Ramaswamy H. Sarma.

Structural characterization of bacterioferritin from Blastochloris viridis

Iron storage and elimination of toxic ferrous iron are the responsibility of bacterioferritins in bacterial species. Bacterioferritins are capable of oxidizing iron using molecular oxygen and import iron ions into the large central cavity of the protein, where they are stored in a mineralized form. We isolated, crystallized bacterioferritin from the microaerophilic/anaerobic, purple non-sulfur bacterium Blastochloris viridis and determined its amino acid sequence and X-ray structure. The structure and sequence revealed similarity to other purple bacterial species with substantial differences in the pore regions. Static 3- and 4-fold pores do not allow the passage of iron ions even though structural dynamics may assist the iron gating. On the other hand the B-pore is open to water and larger ions in its native state. In order to study the mechanism of iron import, multiple soaking experiments were performed. Upon Fe(II) and urea treatment the ferroxidase site undergoes reorganization as seen in bacterioferritin from Escherichia coli and Pseudomonas aeruginosa. When soaking with Fe(II) only, a closely bound small molecular ligand is observed close to Fe₁ and the coordination of Glu94 to Fe₂ changes from bidentate to monodentate. DFT calculations indicate that the bound ligand is most likely a water or a hydroxide molecule representing a product complex. On the other hand the different soaking treatments did not modify the conformation of other pore regions.