Chemical system biology approach to identify multi-targeting FDA inhibitors for treating COVID-19 and associated health complications.

In view of many European countries and the USA leading to the second wave of COVID-19 pandemic, winter season, the evolution of new mutations in the spike protein, and no registered drugs and vaccines for COVID-19 treatment, the discovery of effective and novel therapeutic agents is urgently required. The degrees and frequencies of COVID-19 clinical complications are related to uncontrolled immune responses, secondary bacterial infections, diabetes, cardiovascular disease, hypertension, and chronic pulmonary diseases. It is essential to recognize that the drug repurposing strategy so far remains the only means to manage the disease burden of COVID-19. Despite some success of using single-target drugs in treating the disease, it is beyond suspicion that the virus will acquire drug resistance by acquiring mutations in the drug target. The possible synergistic inhibition of drug efficacy due to drug-drug interaction cannot be avoided while treating COVID-19 and allied clinical complications. Hence, to avoid the unintended development drug resistance and loss of efficacy due to drug-drug interaction, multi-target drugs can be promising tools for the most challenging disease. In the present work, we have carried out molecular docking studies of compounds from the FDA approved drug library, and the FDA approved and passed phase -1 drug libraries with ten therapeutic targets of COVID-19. Results showed that known drugs, including nine anti-inflammatory compounds, four antibiotics, six antidiabetic compounds, and one cardioprotective compound, could effectively inhibit multiple therapeutic targets of COVID-19. Further in-vitro, in vivo, and clinical studies will guide these drugs' proper allocation to treat COVID-19.Communicated by Ramaswamy H. Sarma.

Computational Investigation on the p53-MDM2 Interaction Using the Potential of Mean Force Study.

Murine double minute 2 (MDM2) proteins are found to be overproduced by many human tumors in order to inhibit the functioning of p53 molecules, a tumor suppressor protein. Thus, reactivating p53 functioning in cancer cells by disrupting p53-MDM2 interactions may offer a significant approach in cancer treatment. However, the structural characterization of the p53-MDM2 complex at the atomistic level and the mechanism of binding/unbinding of the p53-MDM2 complex still remain unclear. Therefore, we demonstrate here the probable binding (unbinding) pathway of transactivation domain 1 of p53 during the formation (dissociation) of the p53-MDM2 complex in terms of free energy as a function of reaction coordinate from the potential of mean force (PMF) study using two different force fields: ff99SB and ff99SB-ILDN. From the PMF plot, we noticed the PMF to have a minimum value at a p53-MDM2 separation of 12 Å, with a dissociation energy of 30 kcal mol-1. We also analyzed the conformational dynamics and stability of p53 as a function of its distance of separation from MDM2. The secondary structure content (helix and turns) in p53 was found to vary with its distance of separation from MDM2. The p53-MDM2 complex structure with lowest potential energy was isolated from the ensemble at the reaction coordinate corresponding to the minimum PMF value and subjected to molecular dynamics simulation to identify the interface surface area, interacting residues at the interface, and the stability of the complex. The simulation results highlight the importance of hydrogen bonds and the salt bridge between Lys94 of MDM2 and Glu17 of p53 in the stability of the p53-MDM2 complex. We also carried out the binding free energy calculations and the per residue energy decomposition analyses of the interface residues of the p53-MDM2 complex. We found that the binding affinity between MDM2 and p53 is indeed high [ΔG bind = -7.29 kcal mol-1 from molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) and ΔG bind = -53.29 kcal mol-1 from molecular mechanics/generalized borne surface area]. The total binding energy obtained using the MM/PBSA method was noticed to be closer to the experimental values (-6.4 to -9.0 kcal mol-1). The p53-MDM2 complex binding profile was observed to follow the same trend even in the duplicate simulation run and also in the simulation carried out with different force fields. We found that Lys51, Leu54, Tyr100, and Tyr104 from MDM2 and the residues Phe19, Trp23, and Leu26 from p53 provide the highest energy contributions for the p53-MDM2 interaction. Our findings highlight the prominent structural and binding characteristics of the p53-MDM2 complex that may be useful in designing potential inhibitors to disrupt the p53-MDM2 interactions.

Characterizing the Binding Interactions between DNA-Binding Proteins, XPA and XPE: A Molecular Dynamics Approach.

The scaffold nature of Xeroderma pigmentosum complementation group A (XPA) protein makes it an important member of nucleotide excision repair (NER) that removes bulky DNA lesions with the help of various protein-protein interactions (PPI) and DNA-protein interactions. However, many structural insights of XPA's interaction and the binding patterns with other NER proteins are yet to be understood. Here, we have studied one such crucial PPI of XPA with another NER protein, Xeroderma pigmentosum complementation group A (XPE), by using the previously identified binding site of XPA (residues 185-226) in the Assisted Model Building With Energy Refinement force-field-mediated dynamic system. We studied the relationship between XPA185-226-XPE complex using three different docked models. The major residues observed in all of the models that were responsible for the PPI of this complex were Arg20, Arg47, Asp51, and Leu57 from XPE and the residues Leu191, Gln192, Val193, Trp194, Glu198, Glu202, Glu205, Arg207, Glu209, Gln216, and Phe219 from XPE185-226. During the simulation study, the orientation of XPA was also noted to be changed by almost 180° in models 1 and 3, which remain unchanged in model 2, indicating that XPA interacts with XPE with its N-terminal end facing downward and C-terminal end facing upward. The same was concurrent with the binding of DNA-binding domain region of XPA (aa98-239) with XPE. The N-terminal of XPE was stretched for accommodating XPA. Using the per-residue energy decomposition analysis for the interface residues of all models, the binding affinity between these proteins were found to be dependent on R20, R47, and L57 of XPE and the residues L191, V193, W194, E198, E202, E205, R207, and F219 of XPA. The net binding free energy of the XPA185-226-XPE protein complex was found to be -48.3718 kcal mol-1 for model 1, -49.09 kcal mol-1 for model 2, and -56.51 kcal mol-1 for model 3.