Comparative QM/MM Molecular Dynamics and Umbrella Sampling Simulations: Interaction of the Zinc-Bound Intermediate Gem-Diolate Trapoxin A Inhibitor and Acetyl-l-lysine Substrate with Histone Deacetylase 8

Molecular Dynamics Simulation of an Iron(III) Binding Site on the Fc Domain of IgG1 Relevant for Visible Light-Induced Protein Fragmentation

Molecular dynamics simulations were employed to investigate the interaction between Fe(III) and an iron-binding site composed of THR259, ASP252, and GLU261 on the Fc domain of an IgG1. The goal was to provide microscopic mechanistic information for the photochemical, iron-dependent site-specific oxidative fragmentation of IgG1 at THR259 reported in our previous paper. The distance between Fe(III) and residues of interest as well as the binding pocket size was examined for both protonated and deprotonated THR259. The Fe(III) binding free energy (ΔG) was estimated by using an umbrella sampling approach. The pKa shift of the THR259 hydroxyl group caused by the presence of nearby Fe(III) was estimated based on a thermodynamic cycle. The simulation results show that Fe(III) resides inside the proposed binding pocket and profoundly changes the pocket configuration. The ΔG values indicate that the pocket possesses a strong binding affinity for Fe(III). Furthermore, Fe(III) profoundly lowers the pKa value of the THR259 hydroxyl group by 5.4 pKa units.

Lysine Deacetylase Substrate Selectivity: Distinct Interaction Surfaces Drive Positive and Negative Selection for Residues Following Acetyllysine

Lysine acetylation is a post-translational modification that is reversed by lysine deacetylases (KDACs). The goal of this work was to identify determinants of substrate specificity for KDACs, focusing on short-range interactions occurring with residues immediately following the acetyllysine. Using a fluorescence-based in vitro assay, we determined the activity for each enzyme with a limited panel of derivative substrate peptides, revealing a distinct reactivity profile for each enzyme. We mapped the interaction surface for KDAC6, KDAC8, and KDAC1 with the +1 and +2 substrate residues (with respect to acetyllysine) based on enzyme-substrate interaction pairs observed in molecular dynamics simulations. Characteristic residues in each KDAC interact preferentially with particular substrate residues and correlate with either enhanced or inhibited activity. Although nonpolar aromatic residues generally enhanced activity with all KDACs, the manner in which each enzyme interacted with these residues is distinct. Furthermore, each KDAC has distinctive interactions that correlate with lower activity, primarily ionic in nature. KDAC8 exhibited the most diverse and widest range of effects, while KDAC6 was sensitive only to the +1 position and KDAC1 selectivity was primarily driven by negative selection. The substrate preferences were validated for KDAC6 and KDAC8 using a set of peptides derived from known acetylated proteins. Overall, we determined how KDAC6, KDAC8, and KDAC1 achieve substrate specificity with residues following the acetyllysine. These new insights into KDAC specificity will be critical for identifying novel substrates of particular KDACs, designing KDAC-specific inhibitors, and demonstrate a general framework for understanding substrate specificity for other enzyme classes.

Mechanism of Chiral-Selective Aminoacylation of an RNA Minihelix Explored by QM/MM Free-Energy Simulations

Aminoacylation of a primordial RNA minihelix composed of D-ribose shows L-amino acid preference over D-amino acid without any ribozymes or enzymes. This preference in the amino acylation reaction likely plays an important role in the establishment of homochirality in L-amino acid in modern proteins. However, molecular mechanisms of the chiral selective reaction remain unsolved mainly because of difficulty in direct observation of the reaction at the molecular scale by experiments. For seeking a possible mechanism of the chiral selectivity, quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations of the aminoacylation reactions in a modeled RNA were performed to investigate differences in their free-energy profiles along the reactions for L- and D-alanine and its physicochemical origin. The reaction is initiated by approaching a 3′-oxygen of the RNA minihelix to the carbonyl carbon of an aminoacyl phosphate oligonucleotide. The QM/MM umbrella sampling MD calculations showed that the height of the free-energy barrier for L-alanine aminoacylation reaction was 17 kcal/mol, which was 9 kcal/mol lower than that for the D-alanine system. At the transition state, the distance between the negatively charged 3′-oxygen and the positively charged amino group of L-alanine was shorter than that of D-alanine, which was caused by the chirality difference of the amino acid. These results indicate that the transition state for L-alanine is more electrostatically stabilized than that for D-alanine, which would be a plausible mechanism previously unexplained for chiral selectivity in the RNA minihelix aminoacylation.

Mechanism of the Early Catalytic Events in the Collagenolysis by Matrix Metalloproteinase-1

Metalloproteinase-1 (MMP-1) catalyzed collagen degradation is essential for a wide variety of normal physiological processes, while at the same time contributing to several diseases in humans. Therefore, a comprehensive understanding of this process is of great importance. Although crystallographic and spectroscopic studies provided fundamental information about the structure and function of MMP-1, the precise mechanism of collagen degradation especially considering the complex and flexible structure of the substrate, remains poorly understood. In addition, how the protein environment dynamically reorganizes at the atomic scale into a catalytically active state capable of collagen hydrolysis remains unknown. In this study, we applied experimentally-guided multiscale molecular modeling methods including classical molecular dynamics (MD), well-tempered (WT) classical metadynamics (MetD), combined quantum mechanics/molecular mechanics (QM/MM) MD and QM/MM MetD simulations to explore and characterize the early catalytic events of MMP-1 collagenolysis. Importantly the study provided a complete atomic and dynamic description of the transition from the open to the closed form of the MMP-1•THP complex. Notably, the formation of catalytically active Michaelis complex competent for collagen cleavage was characterized. The study identified the changes in the coordination state of the catalytic zinc(II) associated with the conformational transformation and the formation of catalytically productive ES complex. Our results confirm the essential role of the MMP-1 catalytic domain's α-helices (hA, hB and hC) and the linker region in the transition to the catalytically competent ES complex. Overall, the results provide unique mechanistic insight into the conformational transformations and associated changes in the coordination state of the catalytic zinc(II) that would be important for the design of effective MMP-1 inhibitors.

Exploring the Role of Chemical Reactions in the Selectivity of Tyrosine Kinase Inhibitors

A variety of diseases are associated with tyrosine kinase enzymes that activate many proteins via signal transduction cascades. The similar ATP-binding pockets of these tyrosine kinases make it extremely difficult to design selective covalent inhibitors. The present study explores the contribution of the chemical reaction steps to the selectivity of the commercialized inhibitor acalabrutinib over the Bruton's tyrosine kinase (BTK) and the interleukin-2-inducible T-cell kinase (ITK). Ab initio and empirical valence bond (EVB) simulations of the two kinases indicate that the most favorable reaction path involves a water-assisted mechanism of the 2-butynamide reactive group of acalabrutinib. BTK reacts with acalabrutinib with a substantially lower barrier than ITK, according to our calculated free-energy profile and kinetic simulations. Such a difference is due to the microenvironment of the active site, as further supported by a sequence-based analysis of specificity determinants for several commercialized inhibitors. Our study involves a new approach of simulating directly the IC50 and inactivation efficiency k_eff, instead of using the standard formulas. This new strategy is particularly important in studies of covalent inhibitors with a very exothermic bonding step. Overall, our results demonstrate the importance of understanding the chemical reaction steps in designing selective covalent inhibitors for tyrosine kinases.