Inhibition of Mammalian 15-Lipoxygenase by Three Ebselen-like Drugs. A QM/MM and MM/PBSA Comparative Study

What Happens to a Pyrrole Hemithioindigo Photoswitch Trapped in a Fluorescent Protein?

Artificial molecular photoswitches that can be reversibly controlled into different configurations by external light stimulation have broad application prospects in various fields, such as materials and chemical biology. Among them, the pyrrole hemithioindigo (PHT) photoswitch has a geometric structure similar to that of the fluorescent protein chromophore. What happens when the chromophore is replaced by PHT, and does it achieve similar functions to the original one? To answer these questions, we carried out ONIOM(QM/MM) and classical molecular dynamics studies on the photoisomerization mechanism and spectroscopic properties of PHT in the fluorescent protein. The results showed that in the protein environment, the fate of excited PHT is governed by the competition between fluorescence emission and hula-twist isomerization. Due to the strong steric hindrance effects caused by the interlacing residues in the protein that restrict the PHT conformation transformation, the cis-to-trans isomerization process of PHT needs to overcome a barrier of at least 4.9 kcal/mol; thus, fluorescence emission is more dominant in competition. It is also found that the intermolecular interaction between LYS67 and the carbonyl oxygen of PHT has a significant effect on the fluorescence emission. These results revealed the photochemical reaction mechanism of a light-driven molecular switch in the fluorescent protein and provided further theoretical support for the field of chemical biology.

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.

An Updated Review on Developing Small Molecule Kinase Inhibitors Using Computer-Aided Drug Design Approaches

Small molecule kinase inhibitors (SMKIs) are of heightened interest in the field of drug research and development. There are 79 (as of July 2023) small molecule kinase inhibitors that have been approved by the FDA and hundreds of kinase inhibitor candidates in clinical trials that have shed light on the treatment of some major diseases. As an important strategy in drug design, computer-aided drug design (CADD) plays an indispensable role in the discovery of SMKIs. CADD methods such as docking, molecular dynamic, quantum mechanics/molecular mechanics, pharmacophore, virtual screening, and quantitative structure-activity relationship have been applied to the design and optimization of small molecule kinase inhibitors. In this review, we provide an overview of recent advances in CADD and SMKIs and the application of CADD in the discovery of SMKIs.

Inhibition mechanism and hot-spot prediction of nine potential drugs for SARS-CoV-2 Mpro by large-scale molecular dynamic simulations combined with accurate binding free energy calculations

Coronavirus disease 2019 (COVID-19), which is caused by a new coronavirus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is spreading around the world. However, a universally effective treatment regimen has not been developed to date. The main protease (M^pro), a key enzyme of SARS-CoV-2, plays a crucial role in the replication and transcription of this virus in cells and has become the ideal target for rational antiviral drug design. In this study, we performed molecular dynamics simulations three times for these complexes of M^pro (monomeric and dimeric) and nine potential drugs that have a certain effect on the treatment of COVID-19 to explore their binding mechanism. In addition, a total of 12 methods for calculating binding free energy were employed to determine the optimal drug. Ritonavir, Arbidol, and Chloroquine consistently showed an outstanding binding ability to monomeric M^pro under various methods. Ritonavir, Arbidol, and Saquinavir presented the best performance when binding to a dimer, which was independent of the protonated state of Hie41 (protonated at N_ε) and Hid41 (protonated at N_δ), and these findings suggest that Chloroquine may not effectively inhibit the activity of dimeric M^proin vivo. Furthermore, three common hot-spot residues of Met165, Hie41, and Gln189 of monomeric M^pro systems dominated the binding of Ritonavir, Arbidol, and Chloroquine. In dimeric M^pro, Gln189, Met165, and Met49 contributed significantly to binding with Ritonavir, Arbidol, and Saquinavir; therefore, Gln189 and Met165 might serve as the focus in the discovery and development of anti-COVID-19 drugs. In addition, the van der Waals interaction played a significant role in the binding process, and the benzene ring of the drugs showed an apparent inhibitory effect on the normal function of M^pro. The binding cavity had great flexibility to accommodate these different drugs. The results would be notably helpful for enabling a detailed understanding of the binding mechanisms for these important drug-M^pro interactions and provide valuable guidance for the design of potent inhibitors.

Computational analysis of binding free energies, hotspots and the binding mechanism of Bcl-xL/Bcl-2 binding to Bad/Bax

Hierarchical clustering tree of residues providing contributions to system binding based on the binding free energy of specific residues for (A) Bcl-xL systems (B) Bcl-2 systems.

Entropic effect and residue specific entropic contribution to the cooperativity in streptavidin-biotin binding

Molecular dynamics (MD) simulations were performed employing the polarized protein-specific charge (PPC) to explore the origin of the cooperativity in streptavidin-biotin systems (wild type, two single mutations and one double-mutation). The results of the experiment found that the existence of cooperativity is mainly the result of the entropic effect. In this study, the entropic contribution to the binding free energy was calculated using the recently developed interaction entropy (IE) method, and computational results are in excellent agreement with the experimental observations and are further verified by the calculation of the thermodynamic integration. Comparison of different force fields in terms of predicted binding strength ordering, cooperativity of energy and the stability of hydrogen bonding suggests that the PPC force field combined IE method is a suitable choice. In addition, the IE method enables us to obtain the residue-specific entropic contributions to the streptavidin-biotin binding affinity and identify ten hot-spot residues providing the dominant contribution to the cooperative binding. Importantly, the overall cooperativity obtained from the ten residues also comes mainly from the entropic effect in our study. The calculation of the potential of mean force shows that the unbinding of streptavidin-biotin is a multi-step process, and each step corresponds to the formation and rupture of the hydrogen bond network. And S45A mutation may increase the rigidity of the linker region, making the flap region relatively difficult to open. The present study provides significant molecular insight into the binding cooperativity of the streptavidin-biotin complex.

Dengzhan Shengmai capsules and their active component scutellarin prevent cognitive decline in APP/PS1 mice by accelerating Aβ aggregation and reducing oligomers formation

There is currently no effective treatment to prevent the progress of Alzheimer's disease (AD). The traditional Chinese herbs Dengzhan Shengmai (DZSM) capsules and their active component scutellarin possess multiple effects and are clinically used for the treatment of cerebrovascular diseases. Scutellarin has been reported to affect Aβ aggregation. However, the effects of DZSM capsules on AD remain unknown. Through in vivo experiments, our study proved that the alleviating effects of DZSM capsules on cognitive deficits of AD mice were due to the role of scutellarin, which up-regulated low toxic amyloid plaques and down-regulated highly toxic soluble Aβ42 and Aβ40 levels in cortex. In vitro, we confirmed scutellarin's role in accelerating transforming Aβ42 monomers into high-molecular-mass aggregates by biochemical assays, which supported the results observed in drug-treated APP/PS1 mice. In detail, the 1:10 ratio of scutellarin/Aβ42 mixtures promoted production of large β-sheet-rich fibrils whereas the 1:1 ratio promoted production of protofibrils. In addition, the binding between scutellarin and Aβ monomers was quantified by microscale thermophoresis test and the apparent dissociation constant (Kd) was 1284.4 ± 238.8 μM. What's more, binding regions between scutellarin and Aβ fibrils were predicted by computational docking models and scutellarin might bind parallel to the long axis of Aβ42 fibrils targeting hydrophobic grooves at residues 35-36 or 39. In conclusion, DZSM capsules protected against cognitive defects of AD through scutellarin-mediated acceleration of Aβ aggregation into fibrils or protofibrils and reduction of soluble Aβ oligomers, thus suggesting potential clinical applications of DZSM capsules and scutellarin in the treatment of AD.

Insight Into the Binding Mechanism of p53/pDIQ-MDMX/MDM2 With the Interaction Entropy Method

The study of the p53-MDMX/MDM2 binding sites is a research hotspot for tumor drug design. The inhibition of p53-targeted MDMX/MDM2 has become an effective approach in anti-tumor drug development. In this paper, a theoretically rigorous and computationally accurate method, namely, the interaction entropy (IE) method, combined with the polarized protein-specific charge (PPC) force field, is used to explore the difference in the binding mechanism between p53-MDMX and p53-MDM2. The interaction of a 12mer peptide inhibitor (pDIQ), which is similar to p53 in structure, with MDMX/MDM2 is also studied. The results demonstrate that p53/pDIQ with MDM2 generates a stronger interaction than with MDMX. Compared to p53, pDIQ has larger binding free energies with MDMX and MDM2. According to the calculated binding free energies, the differences in the binding free energy among the four complexes that are obtained from the combination of PPC and IE are more consistent with the experimental values than with the results from the combination of the non-polarizable AMBER force field and IE. In addition, according to the decomposition of the binding free energy, the van der Waals (vdW) interactions are the main driving force for the binding of the four complexes. They are also the main source of the weaker binding affinity of p53/pDIQ-MDMX relative to p53/pDIQ-MDM2. Compared with p53-MDMX/MDM2, according to the analysis of the residue decomposition, the predicated total residue contributions are higher in pDIQ-MDMX/MDM2 than in p53-MDMX/MDM2, which explains why pDIQ has higher binding affinity than p53 with MDMX/MDM2. The current study provides theoretical guidance for understanding the binding mechanisms and designing a potent dual inhibitor that is targeted to MDMX/MDM2.

The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction

At present, the calculated binding free energy obtained using the molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MM/PB(GB)SA) method is overestimated due to the lack of knowledge of suitable interior dielectric constants in the simulation on the interaction of Human Immunodeficiency Virus (HIV-1) protease systems with inhibitors. Therefore, the impact of different values of the interior dielectric constant and the entropic contribution when using the MM/PB(GB)SA method to calculate the binding free energy was systemically evaluated. Our results show that the use of higher interior dielectric constants (1.4-2.0) can clearly improve the predictive accuracy of the MM/PBSA and MM/GBSA methods, and computational errors are significantly reduced by including the effects of electronic polarization and using a new highly efficient interaction entropy (IE) method to calculate the entropic contribution. The suitable range for the interior dielectric constant is 1.4-1.6 for the MM/PBSA method; within this range, the correlation coefficient fluctuates around 0.84, and the mean absolute error fluctuates around 2 kcal/mol. Similarly, an interior dielectric constant of 1.8-2.0 produces a correlation coefficient of approximately 0.76 when using the MM/GBSA method. In addition, the entropic contribution of each individual residue was further calculated using the IE method to predict hot-spot residues, and the detailed binding mechanisms underlying the interactions of the HIV-1 protease, its inhibitors, and bridging water molecules were investigated. In this study, the use of a higher interior dielectric constant and the IE method can improve the calculation accuracy of the HIV-1 system.

Protein-Ligand Interaction Energy-Based Entropy Calculations: Fundamental Challenges For Flexible Systems

Entropy calculations represent one of the most challenging steps in obtaining the binding free energy in biomolecular systems. A novel computationally effective approach (IE) was recently proposed to calculate the entropy based on the computation of protein-ligand interaction energy directly from molecular dynamics (MD) simulations. We present a study focused on the application of this method to flexible molecular systems and compare its performance with well-established normal mode (NM) and quasiharmonic (QH) entropy calculation approaches. Our results demonstrated that the IE method is intended for calculating entropy change for binding partners in fixed conformations, as by the original definition of IE, and is not applicable to the molecular complexes in which the interacting partners undergo significant conformational changes during the binding process.