Molecular Mechanism of the Non-Covalent Orally Targeted SARS-CoV-2 Mpro Inhibitor S-217622 and Computational Assessment of Its Effectiveness against Mainstream Variants

4D-QSAR, ADMET properties, and molecular dynamics simulations for designing N-substituted urea/thioureas as human glutaminyl cyclase inhibitors

Human glutaminyl cyclase (hQC) inhibitors have great potential to be used as anti- Alzheimer's disease (AD) agents by reducing the toxic pyroform of β-amyloid in the brains of AD patients. The four-dimensional quantitative structure activity relationship (4D-QSAR) model of N-substituted urea/thioureas was established with satisfying predictive ability and statistical reliability (Q2 = 0.521, R2 = 0.933, R2prep = 0.619). By utilizing the developed 4D-QSAR model, a set of new N-substituted urea/thioureas was designed and evaluated for their Absorption Distribution Metabolism Excretion and Toxicity (ADMET) properties. The results of molecular dynamics (MD) simulations, Principal component analysis (PCA), free energy landscape (FEL), dynamic cross-correlation matrix (DCCM) and molecular mechanics generalized Born Poisson-Boltzmann surface area (MM-PBSA) free energy calculations, revealed that the designed compounds were remained stable in protein binding pocket and compounds b ∼ f (-35.1 to -44.55 kcal/mol) showed higher binding free energy than that of compound 14 (-33.51 kcal/mol). The findings of this work will be a theoretical foundation for further research and experimental validation of urea/thiourea derivatives as hQC inhibitors.

The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes

Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.

Advanced molecular mechanisms of modified DRV compounds in targeting HIV-1 protease mutations and interrupting monomer dimerization

HIV-1 protease (PR) plays a crucial role in the treatment of HIV as a key target. The global issue of emerging drug resistance is escalating, and PR mutations pose a substantial challenge to the effectiveness of inhibitors. HIV-1 PR is an ideal model for studying drug resistance to inhibitors. The inhibitor, darunavir (DRV), exhibits a high genetic barrier to viral resistance, but with mutations of residues in the PR, there is also some resistance to DRV. Inhibitors can impede PR in two ways: one involves binding to the active site of the dimerization protease, and the other involves binding to the PR monomer, thereby preventing dimerization. In this study, we aimed to investigate the inhibitory effect of DRV with a modified inhibitor on PR, comparing the differences between wild-type and mutated PR, using molecular dynamics simulations. The inhibitory effect of the inhibitors on PR monomers was subsequently investigated. And molecular mechanics Poisson-Boltzmann surface area evaluated the binding free energy. The energy contribution of individual residues in the complex was accurately calculated by the alanine scanning binding interaction entropy method. The results showed that these inhibitors had strong inhibitory effects against PR mutations, with GRL-142 exhibiting potent inhibition of both the PR monomer and dimer. Improved inhibitors could strengthen hydrogen bonds and interactions with PR, thereby boosting inhibition efficacy. The binding of the inhibitor and mutation of the PR affected the distance between D25 and I50, preventing their dimerization and the development of drug resistance. This study could accelerate research targeting HIV-1 PR inhibitors and help to further facilitate drug design targeting both mechanisms.

Binding Thermodynamics of Fourth-Generation EGFR Inhibitors Revealed by Absolute Binding Free Energy Calculations

The overexpression or mutation of the kinase domain of the epidermal growth factor receptor (EGFR) is strongly associated with non-small-cell lung cancer (NSCLC). EGFR tyrosine kinase inhibitors (TKIs) have proven to be effective in treating NSCLC patients. However, EGFR mutations can result in drug resistance. To elucidate the mechanisms underlying this resistance and inform future drug development, we examined the binding affinities of BLU-945, a recently reported fourth-generation TKI, to wild-type EGFR (EGFRWT) and its double-mutant (L858R/T790M; EGFRDM) and triple-mutant (L858R/T790M/C797S; EGFRTM) forms. We compared the binding affinities of BLU-945, BLU-945 analogues, CH7233163 (another fourth-generation TKI), and erlotinib (a first-generation TKI) using absolute binding free energy calculations. Our findings reveal that BLU-945 and CH7233163 exhibit binding affinities to both EGFRDM and EGFRTM stronger than those of erlotinib, corroborating experimental data. We identified K745 and T854 as the key residues in the binding of fourth-generation EGFR TKIs. Electrostatic forces were the predominant driving force for the binding of fourth-generation TKIs to EGFR mutants. Furthermore, we discovered that the incorporation of piperidinol and sulfone groups in BLU-945 substantially enhanced its binding capacity to EGFR mutants. Our study offers valuable theoretical insights for optimizing fourth-generation EGFR TKIs.

Entropy driven cooperativity effect in multi-site drug optimization targeting SARS-CoV-2 papain-like protease

Papain-like protease (PLpro), a non-structural protein encoded by SARS-CoV-2, is an important therapeutic target. Regions 1 and 5 of an existing drug, GRL0617, can be optimized to produce cooperativity with PLpro binding, resulting in stronger binding affinity. This work investigated the origin of the cooperativity using molecular dynamics simulations combined with the interaction entropy (IE) method. The regions' improvement exhibits cooperativity by calculating the binding free energies between the complex of PLpro-inhibitor. The thermodynamic integration method further verified the cooperativity generated in the drug improvement. To further determine the specific source of cooperativity, enthalpy and entropy in the complexes were calculated using molecular mechanics/generalized Born surface area and IE. The results show that the entropic change is an important contributor to the cooperativity. Our study also identified residues P248, Q269, and T301 that play a significant role in cooperativity. The optimization of the inhibitor stabilizes these residues and minimizes the entropic loss, and the cooperativity observed in the binding free energy can be attributed to the change in the entropic contribution of these residues. Based on our research, the application of cooperativity can facilitate drug optimization, and provide theoretical ideas for drug development that leverage cooperativity by reducing the contribution of entropy through multi-locus binding.

Antitarget, Anti-SARS-CoV-2 Leads, Drugs, and the Drug Discovery-Genetics Alliance Perspective

The most advanced antiviral molecules addressing major SARS-CoV-2 targets (Main protease, Spike protein, and RNA polymerase), compared with proteins of other human pathogenic coronaviruses, may have a short-lasting clinical efficacy. Accumulating knowledge on the mechanisms underlying the target structural basis, its mutational progression, and the related biological significance to virus replication allows envisaging the development of better-targeted therapies in the context of COVID-19 epidemic and future coronavirus outbreaks. The identification of evolutionary patterns based solely on sequence information analysis for those targets can provide meaningful insights into the molecular basis of host-pathogen interactions and adaptation, leading to drug resistance phenomena. Herein, we will explore how the study of observed and predicted mutations may offer valuable suggestions for the application of the so-called "synthetic lethal" strategy to SARS-CoV-2 Main protease and Spike protein. The synergy between genetics evidence and drug discovery may prioritize the development of novel long-lasting antiviral agents.