Molecular Insights into the Impact of Mutations on the Binding Affinity of Targeted Covalent Inhibitors of BTK.

Targeted covalent inhibitors (TCIs) have witnessed a significant resurgence in recent years, particularly in the kinase drug discovery field for treating diverse clinical indications. The inhibition of Bruton's tyrosine kinase (BTK) for treating B-cell cancers is a classic example where TCIs such as ibrutinib have had breakthroughs in targeted therapy. However, selectivity remains challenging, and the emergence of resistance mutations is a critical concern for clinical efficacy. Computational methods that can accurately predict the impact of mutations on inhibitor binding affinity could prove helpful in informing targeted approaches─providing insights into drug resistance mechanisms. In addition, such systems could help guide the systematic evaluation and impact of mutations in disease models for optimal experimental design. Here, we have employed in silico physics-based methods to understand the effects of mutations on the binding affinity and conformational dynamics of select TCIs of BTK. The TCIs studied include ibrutinib, acalabrutinib, and zanubrutinib─all of which are FDA-approved drugs for treating multiple forms of leukemia and lymphoma. Our results offer useful molecular insights into the structural determinants, thermodynamics, and conformational energies that impact ligand binding for this biological target of clinical relevance.

Synthesis, Characterization, and DFT Calculations of a New Sulfamethoxazole Schiff Base and Its Metal Complexes.

A new Schiff base, 4-((1E,2E)-3-(furan-2-yl)allylidene)amino)-N-(5-methylisoxazol-3-yl) benzene-sulfonamide (L), was synthesized by thermal condensation of 3-(2-furyl)acrolein and sulfamethoxazole (SMX), and the furan Schiff base (L) was converted to a phenol Schiff base (L') according to the Diels-Alder [4 + 2] cycloaddition reaction and studied experimentally. The structural and spectroscopic properties of the Schiff base were also corroborated by utilizing density functional theory (DFT) calculations. Furthermore, a series of lanthanide and transition metal complexes of the Schiff base were synthesized from the nitrate salts of Gd, Sm, Nd, and Zn (L1, L2, L3, and L4), respectively. Various spectroscopic studies confirmed the chemical structures of the Schiff-base ligand and its complexes. Based on the spectral studies, a nine-coordinated geometry was assigned to the lanthanide complexes and a six-coordinated geometry to the zinc complex. The elemental analysis data confirmed the suggested structure of the metal complexes, and the TGA studies confirmed the presence of one coordinated water molecule in the lanthanide complexes and one crystalline water molecule in the zinc complex; in addition, the conductivity showed the neutral nature of the complexes. Therefore, it is suggested that the ligand acts as a bidentate through coordinates to each metal atom by the isoxazole nitrogen and oxygen atoms of the sulfur dioxide moiety of the SMX based on FTIR studies. The ligand and its complexes were tested for their anti-inflammatory, anti-hemolytic, and antioxidant activities by various colorimetric methods. These complexes were found to exhibit potential effects of the selected biological activities.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA.

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.

Accelerating the discovery of the beyond rule of five compounds that have high affinities toward SARS-CoV-2 spike RBD.

The battle against SARS-CoV-2 coronavirus is the focal point for the global pandemic that has affected millions of lives worldwide. The need for effective and selective therapeutics for the treatment of the disease caused by SARS-CoV-2 is critical. Herein, we performed a hierarchical computational approach incorporating molecular docking studies, molecular dynamics simulations, absolute binding energy calculations, and steered molecular dynamics simulations for the discovery of potential compounds with high affinity towards SARS-CoV-2 spike RBD. By leveraging ZINC15 database, a total of 1282 in-clinical and FDA approved drugs were filtered out from nearly 0.5 million protomers of relatively large compounds (MW > 500, and LogP ≤ 5). Our results depict plausible mechanistic aspects related to the blockage of SARS-CoV-2 spike RBD by the top hits discovered. We found that the most promising candidates, namely, ZINC95628821, ZINC95617623, ZINC3979524, and ZINC261494658, strongly bind to the spike RBD and interfere with the human ACE2 receptor. These findings accelerate the rational design of selective inhibitors targeting the spike RBD protein of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

High-throughput virtual screening of drug databanks for potential inhibitors of SARS-CoV-2 spike glycoprotein.

COVID-19, which is caused by a novel coronavirus known as SARS-CoV-2, has spread rapidly around the world, and it has infected more than 29 million individuals as recorded on 16 September 2020. Much effort has been made to stop the virus from spreading, and there are currently no approved pharmaceutical products to treat COVID-19. Here, we apply an in silico approach to investigate more than 3800 FDA approved drugs on the viral RBD S1-ACE2 interface as a target. The compounds were investigated through flexible ligand docking, ADME property calculations and protein-ligand interaction maps. Molecular dynamics (MD) simulations were also performed on eleven compounds to study the stability and the interactions of the protein-ligand complexes. The MD simulations show that bagrosin, chidamide, ebastine, indacaterol, regorafenib, salazosulfadimidine, silodosin and tasosartan are relatively stable near the C terminal domain (CTD1) of the S1 subunit of the viral S protein. The relative MMGBSA binding energies show that silodosin has the best binding to the target. The constant velocity steered molecular dynamics (SMD) simulations show that silodosin preferentially interacts with the RBD S1 and has potential to act as an interfering compound between viral spike-host ACE2 interactions. Communicated by Ramaswamy H. Sarma.

Covalent and non-covalent binding free energy calculations for peptidomimetic inhibitors of SARS-CoV-2 main protease.

COVID-19, the disease caused by the newly discovered coronavirus-SARS-CoV-2, has created a global health, social, and economic crisis. As of mid-January 2021, there are over 90 million confirmed cases and more than 2 million reported deaths due to COVID-19. Currently, there are very limited therapeutics for the treatment or prevention of COVID-19. For this reason, it is important to find drug targets that will lead to the development of safe and effective therapeutics against the disease. The main protease (Mpro) of the virus is an attractive target for the development of effective antiviral therapeutics because it is required for proteolytic cleavage of viral polyproteins. Furthermore, the Mpro has no human homologues, so drugs designed to bind to this target directly have less risk for off-target effects. Recently, several high-resolution crystallographic structures of the Mpro in complex with inhibitors have been determined-to guide drug development and to spur efforts in structure-based drug design. One of the primary objectives of modern structure-based drug design is the accurate prediction of receptor-ligand binding affinities for rational drug design and discovery. Here, we perform rigorous alchemical absolute binding free energy calculations and QM/MM calculations to give insight into the total binding energy of two recently crystallized inhibitors of SARS-CoV-2 Mpro, namely, N3 and α-ketoamide 13b. The total binding energy consists of both covalent and non-covalent binding components since both compounds are covalent inhibitors of the Mpro. Our results indicate that the covalent and non-covalent binding free energy contributions of both inhibitors to the Mpro target differ significantly. The N3 inhibitor has more favourable non-covalent interactions, particularly hydrogen bonding, in the binding site of the Mpro than the α-ketoamide inhibitor. Also, the Gibbs energy of reaction for the Mpro-N3 covalent adduct is greater than the Gibbs reaction energy for the Mpro-α-ketoamide covalent adduct. These differences in the covalent and non-covalent binding free energy contributions for both inhibitors could be a plausible explanation for their in vitro differences in antiviral activity. Our findings are consistent with the reversible and irreversible character of both inhibitors as reported by experiment and highlight the importance of both covalent and non-covalent binding free energy contributions to the absolute binding affinity of a covalent inhibitor towards its target. This information could prove useful in the rational design, discovery, and evaluation of potent SARS-CoV-2 Mpro inhibitors for targeted antiviral therapy.

Discovery of potent inhibitors for SARS-CoV-2's main protease by ligand-based/structure-based virtual screening, MD simulations, and binding energy calculations.

COVID-19 has caused lockdowns all over the world in early 2020, as a global pandemic. Both theoretical and experimental efforts are seeking to find an effective treatment to suppress the virus. In silico drug design can play a vital role in identifying promising drug candidates against COVID-19. Herein, we focused on the main protease of SARS-CoV-2 that has crucial biological functions in the virus. We performed a ligand-based virtual screening followed by a docking screening for testing approved drugs and bioactive compounds listed in the DrugBank and ChEMBL databases. The top 8 docking results were advanced to all-atom MD simulations to study the relative stability of the protein-ligand interactions. MD simulations support that the catalytic residue, His41, has a neutral side chain with a protonated delta position. An absolute binding energy (ΔG) of -42 kJ mol-1 for the protein-ligand (Mpro-N3) complex has been calculated using the potential-of-mean-force (geometrical) approach. Furthermore, the relative binding energies were computed for the top docking results. Our results suggest several promising approved and bioactive inhibitors of SARS-CoV-2 Mpro as follows: a bioactive compound, ChEMBL275592, which has the best MM/GBSA binding energy; the second-best compound, montelukast, is an approved drug used in the treatment of asthma and allergic rhinitis; the third-best compound, ChEMBL288347, is a bioactive compound. Bromocriptine and saquinavir are other approved drugs that also demonstrate stability in the active site of Mpro, albeit their relative binding energies are low compared to the N3 inhibitor. This study provides useful insights into de novo protein design and novel inhibitor development, which could reduce the cost and time required for the discovery of a potent drug to combat SARS-CoV-2.

Role of Asp190 in the Phosphorylation of the Antibiotic Kanamycin Catalyzed by the Aminoglycoside Phosphotransferase Enzyme: A Combined QM:QM and MD Study.

The aminoglycoside phosphotransferase (APH(3')-IIIa) kinases form a clinically central group of antibiotic-resistant enzymes. Computationally, we have studied the catalytic mechanism of the APH(3')-IIIa enzyme at the atomic-level. The proposed reaction mechanism involves protonation of Asp190 by the kanamycin 3'-hydroxyl group mediated through an explicit neighboring water molecule, which leads to a simultaneous nucleophilic attack on the γ-phosphate of the ATP by the deprotonated kanamycin 3'-hydroxyl group. The second step is a proton abstraction from the protonated Asp190 to the phosphate group of the phosphorylated kanamycin mediated by an explicit water molecule. The calculated Gibbs energy of activation (ΔG⧧) of the rate-determining step for the phosphorylation reaction is 77 kJ mol-1 at the M06-2X/6-311++G(2df,p)//ONIOM(M06-2X/6-31+G(d):HF/6-31G(d)) level of theory. This study has provided a new understanding of the APH(3')-IIIa catalytic mechanism that agrees with the available experimental data (ΔG⧧ = 75 ± 4 kJ mol-1) and could provide a starting point for the rational design of mechanism-based inhibitors of aminoglycoside modifying enzyme to circumvent antibiotic resistance.