Probing inhibition mechanisms of adenosine deaminase by using molecular dynamics simulations

In silico drug repurposing for coronavirus (COVID-19): screening known HCV drugs against the SARS-CoV-2 spike protein bound to angiotensin-converting enzyme 2 (ACE2) (6M0J)

In this study, FDA-approved HCV antiviral drugs and their structural analogues-several of them in clinical trials-were tested for their inhibitory properties toward the SARS-CoV-2 spike protein bound to angiotensin-converting enzyme 2 (6M0J) using a virtual screening approach and computational chemistry methods. The most stable structures and the corresponding binding affinities of thirteen such antiviral compounds were obtained. Frontier molecular orbital theory, global reactivity descriptors, molecular docking calculations and electrostatic potential analysis were used to hypothesize the bioactivity of these drugs against 6M0J. It is found that an increased affinity for the protein is shown by inhibitors with large compound volume, relatively higher electrophilicity index, aromatic rings and heteroatoms that participate in hydrogen bonding. Among the tested drugs, four compounds 10-13 showed excellent results-binding affinities - 11.2 to - 11.5 kcal mol-1. These four top scoring compounds may act as lead compounds for further experimental validation, clinical trials and even for the development of more potent antiviral agents against the SARS-CoV-2. Approved HCV drugs and analogues were tested for their bioactivity towards the SARS-CoV-2 (6M0J) using virtual screening, ESP and MD analysis.

What make malarial adenosine deaminase from PLASMODIUM VIVAX recognise adenosine and 5'-methylthioadenosine: simulation studies

Malaria is a life-threatening disease in humans caused by Plasmodium parasites. Plasmodium vivax (P. vivax) is one of the prevalent species found worldwide. An increase in an anti-malarial drug resistance suggests the urgent need for new drugs. Zn²⁺-containing adenosine deaminase (ADA) is a promising drug target because the ADA inhibition is fatal to the parasite. Malarial ADA accepts both adenosine (ADN) and 5'-methylthioadenosine (MTA) as substrates. The understanding of the substrate binding becomes crucial for an anti-malarial drug development. In this work, ADA from P. vivax (pvADA) is of interest due to its prevalence worldwide. The binding of ADN and MTA are studied here using Molecular Dynamics (MD) simulations. Upon binding, the open and closed states of pvADA are captured. The displacement of α7, linking loops of β3/α12, β4/α13, β5/α15, and α10/α11 is involved in the cavity closure and opening. Also, the inappropriate substrate orientation induces a failure in a complete cavity closure. Interactions with D46, D172, S280, D310, and D311 are important for ADN binding, whereas only hydrogen bonds with D172 and D311 are sufficient to anchor MTA inside the pocket. No Zn²⁺-coordinated histidine residues is acquired for substrate binding. D172 is found to play a role in ribose moiety recognition, while D311 is crucial for trapping the amine group of an adenine ring towards the Zn²⁺ site. Comparing between ADN and MTA, the additional interaction between D310 and an amine nitrogen on ADN supports a tighter fit that may facilitate the deamination.Communicated by Ramaswamy H. Sarma.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

Random acceleration and steered molecular dynamics simulations reveal the (un)binding tunnels in adenosine deaminase and critical residues in tunnels

Adenosine deaminase (ADA) is an important enzyme related to purine nucleoside metabolism in human serum and various tissues. Abnormal ADA levels are related to a wide variety of diseases such as rheumatoid arthritis, AIDS, anemia, lymphoma, and leukemia and ADA is considered as a useful target for various diseases. Currently, ADA can be divided into open conformation and closed conformation according to the inhibitors of binding. As a consequence, we chose two inhibitors, namely, 6-hydroxy-1,6-dihydro purine nucleoside (PRH) and N-[4,5-bis(4-hydroxyphenyl)-1,3-thiazol-2-yl]hexanamide (FRK) to bind to ADA in the closed conformation or open conformation respectively. In this study, we performed the random acceleration molecular dynamics (RAMD) method, steered molecular dynamics (SMD) simulations and adaptive basing force (ABF) simulation to explore the unbinding tunnels and tunnel characteristics of the two inhibitors in ADA. Our results showed that PRH and FRK escaped from ADA using three main tunnels (namely, T1, T2, and T3). Inhibitors (PRH and FRK) escape through T3 more frequently and more easily. The results from ABF simulations confirm that the free energy barrier in T1 or T2 is larger than that in T3 when inhibitors dissociate from the ADA and have potential mean of force (PMF) depth. Moreover, in the complexes (ADA-PRH, ADA-FRK), we also found that the most active residue that remarkably contributed to the binding affinity is W117 in T3, and the residue played an important role in the unbinding tunnel for inhibitor leaving. Our theoretical study provided insight into the ADA inhibitor passway mechanism and may be a clue for potent ADA inhibitor design.

Molecular mechanisms of metal ions in regulating the catalytic efficiency of D-psicose 3-epimerase revealed by multiple short molecular dynamic simulations and free energy predictions

Communicated by Ramaswamy H. Sarma.

Exploring the binding modes of cordycepin to human adenosine deaminase 1 (ADA1) compared to adenosine and 2'-deoxyadenosine

Cordycepin (3'-deoxyadenosine, abbreviated as COR) from Cordyceps shows a wide range of pharmacological activities, including antioxidant and anticancer effects, therefore representing a potential alternative medicine. However, COR has a short half-life in the human body, where it is metabolized by adenosine deaminase 1 (ADA1). ADA1 helps regulate adenosine levels by deaminating excess adenosine (ADE) and its derivatives, such as 2'-deoxyadenosine (DEO). Understanding binding mechanisms of ADA1 with COR in comparison with its other substrates will play a vital role in improving the bioactivity and lifetime of COR for commercial medicinal use. Recently, the first structure of human ADA1 in complex with DEO was solved. We therefore employed molecular dynamics (MD) simulations to predict structures and dynamics of ADA1 complexing with ADE, DEO, and COR in comparison to a ligand-free (LF) structure. Our data reveal that a large and highly water-exposed binding pocket of ADA1 is responsible for ligand translocation and reorientation. Two possible binding locations (site1 and site2) are identified. The binding affinities of the ligands are ADE > COR > DEO. Furthermore, the movements of two loop regions at the binding pocket entrance, residues 183-193 and 215-230, contribute to gating activity.

Enzymatic characterization and molecular mechanism of a novel aspartokinase mutant M372I/T379W from Corynebacterium pekinense

A novel aspartokinase mutant M372I/T379W from Corynebacterium pekinense was constructed by using site-directed mutagenesis. The enzyme was then purified, characterized, and its molecular mechanism was comprehensively analyzed. Compared with wild-type AK, the catalytic activity of M372I/T379W AK was 16.51 fold higher and the optimum temperature increased from 28 to 35 °C. The thermostability of M372I/T379W AK was significantly improved. Microscale thermophoresis analysis indicated that M372I/T379W AK not only weakened the inhibitory effect of Lys, but also had stronger binding force with Asp. Molecular dynamics simulation showed that mutations M372I and T379W could regulate the activity of CpAK through affecting the flexibility of Asp and ATP binding pocket residues and the hydrogen bond between CpAK and Asp. In addition, mutations could affect the relative position of protein domains. The width of the Asp binding pocket entrance gate Arg169-Ala60 of M372I/T379W AK was greater than that in wild-type AK and the CpAK switched from T-state to R-state, which promoted the binding of the enzyme to Asp and improving the catalytic efficiency of this enzyme. These results explain the molecular mechanism of M372I/T379W AK, which will greatly facilitate the rational design of more aspartokinase mutants, with have potential applications in aspartic acid metabolism.