AMBER Force Field Parameters for Cobalt-Containing Biological Systems: A Systematic Derivation Study

Molecular dynamics simulations to explore the binding mode between the amyloid-β protein precursor (APP) and adaptor protein Mint2

Alzheimer's disease (AD) presents a significant challenge in neurodegenerative disease management, with limited therapeutic options available for its prevention and treatment. At the heart of AD pathogenesis is the amyloid-β (Aβ) protein precursor (APP), with the interaction between APP and the adaptor protein Mint2 being crucial. Despite previous explorations into the APP-Mint2 interaction, the dynamic regulatory mechanisms by which Mint2 modulates APP binding remain poorly understood. This study undertakes molecular dynamics simulations across four distinct systems-free Mint2, Mint2 bound to APP, a mutant form of Mint2, and the mutant form bound to APP-over an extensive 400 ns timeframe. Our findings reveal that the mutant Mint2 experiences significant secondary structural transformations, notably the formation of an α-helix in residues S55-K65 upon APP binding, within the 400 ns simulation period. Additionally, we observed a reduction in the active pocket size of the mutant Mint2 compared to its wild-type counterpart, enhancing its APP binding affinity. These insights hold promise for guiding the development of novel inhibitors targeting the Mints family, potentially paving the way for new therapeutic strategies in AD prevention and treatment.

Assign_v2: a novel bonded-force field parameterization software for square planar palladium molecular dynamics simulations

In the present study, a requisite bonded force field for conducting molecular dynamics simulations of palladium compounds was derived utilizing a novel technique termed the "numbering system," implemented through a program named assign.py. A significant challenge in deriving a bonded force field for square-planar palladium compounds emanates from their cis-trans characteristics, whereby the bond angles of identical atoms diverge across different compounds due to the presence of cis and trans isomers. To navigate this challenge, atoms surrounding the palladium were assigned numerical identifiers; consecutively numbered neighboring atoms were designated as "cis," while atoms alternately numbered were regarded as "trans" relative to each other. Through the employment of the newly developed assign.py program, atoms common to GAFF were automatically assigned, facilitating the seamless integration of GAFF with the newly derived force field for cis and trans configurations. The tested innovative force field demonstrated results in close proximity to experimental findings. This research introduced two novel elements: firstly, a comparatively extensive force field was derived for palladium compounds, enabling the simulation of the previously non-simulable square planar geometry of palladium. Secondly, new software, assign.py, was developed, capable of amalgamating "numbering system atom types" with GAFF atom types into a unified force field. Conclusive tests within the study affirmed that the innovative bonded force field and technique utilizing the assign.py program was successful, simulating the square-planar geometry of palladium in a manner highly congruent with the experimental geometry. The program, licensed under Apache2, can be accessed via the following link: https://github.com/bkurt00/palladiumFF2/blob/main/assign_v2/assign_v2.pyCommunicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulations Combined with Markov Model to Explore the Effect of Allosteric Inhibitor Binding on Bromodomain-Containing Protein 4

Bromodomain-Containing Protein 4 (BRD4) can play an important role in gene transcriptional regulation of tumor development and survival by participating in histone modification epigenetic mechanism. Although it has been reported that novel allosteric inhibitors such as ZL0590 have a high affinity with target protein BRD4 and good efficacy, their inhibitory mechanism has not been studied further. The aim of this study was to reveal the inhibition mechanism of allosteric inhibitor ZL0590 on Free-BRD4 and BRD4 binding MS436 (orthosteric inhibitor) by molecular dynamics simulation combined with a Markov model. Our results showed that BRD4-ZL0590 led to α-helices formation of 100-105 compared with Free-BRD4; the combination of MS436 caused residues 30-40 and 95-105 to form α-helices, while the combination of allosteric inhibitors untangled the α-helices formed by the MS436. The results of Markov flux analysis showed that the binding process of inhibitors mainly involved changes in the degree of α-helices at ZA loop. The binding of ZL0590 reduced the distance between ZA loop and BC loop, blocked the conformation at the active site, and inhibited the binding of MS436. After the allosteric inhibitor binding, the MS436 that could normally penetrate into the interior of the pocket was floating on the edge of the active pocket and did not continue to penetrate into the active pocket as expected. In summary, we provide a theoretical basis for the inhibition mechanism of ZL0590 against BRD4, which can be used as a reference for improving the development of drug targets for cancer therapy.