Computational Analysis of Interactions Between Drugs and Human Serum Albumin.

Drug molecules exist as complexed with serum proteins such as human serum albumin (HSA) and/or unbound free form in the blood circulation. Drugs can be effective only when they are free. Thus, it is important to understand aspects that are important for interaction between drugs and interacting proteins. In this study, interactions among 2990 FDA approved drugs and HSA were computational analyzed to unravel principles that are critical for drug-HSA interactions. Docking results showed that drugs have higher affinity toward cavity-1 (C1) than cavity-2 (C2). A total of 1131 drug molecules have docking score greater than 60 while 768 molecules have docking score greater than 60 when they are docked in C2. In addition, three solvent channels have potential to direct solvent to C1 cavity while C2 does not have any effective channel. The post MD analyses demonstrated that drugs are making polar interactions with basic amino acids in the binding cavities. Verbscoside and ceftazidime both have stable low RMSD values throughout MD simulation with 2 Å on average in C1 cavity. The ligand RMSD shows less stability for verbscoside, which is around 4 Å when it is in complex with HSA in C1. The individual contribution of the residues K192, K196, R215, and R254 to ceftazidime are -1.92 ± 0.18, -3.09 ± 0.09, -2.17 ± 0.17, and - 2.32 ± 0.098, respectively. These residues contribute the binding energy of the verbscoside by -6.06 ± 0.08, -2.10 ± 0.06, and - 1.57 ± 0.03 kcal/mol individually in C1 cavity. C2 is making polar interactions with drug via R469, K472, and K488 residues and their contribution to the two drugs are -3.13 ± 0.21 kcal/mol for R469, -1.94 ± 0.18 kcal/mol for K472, and -1.96 ± 0.11 kcal/mol for K488 to total binding energy of ceftazidime. The binding energy of verbscoside is 57.17 ± 7.00 kcal/mol and Arg-407 has the highest contribution this bind energy individually with -4.29 ± 0.12 kcal/mol. Drugs with hydrogen bond donor/acceptor chemical adducts such as verbscoside involve higher hydrogen bond formation in C1 pocket. Ceftazidime makes interaction with HSA toward hydrophobic residues, L384, L404, L487, and L488 in the C2 cavity.

Benchmarking ANI potentials as a rescoring function and screening FDA drugs for SARS-CoV-2 Mpro.

Here, we introduce the use of ANI-ML potentials as a rescoring function in the host-guest interaction in molecular docking. Our results show that the "docking power" of ANI potentials can compete with the current scoring functions at the same level of computational cost. Benchmarking studies on CASF-2016 dataset showed that ANI is ranked in the top 5 scoring functions among the other 34 tested. In particular, the ANI predicted interaction energies when used in conjunction with GOLD-PLP scoring function can boost the top ranked solution to be the closest to the x-ray structure. Rapid and accurate calculation of interaction energies between ligand and protein also enables screening of millions of drug candidates/docking poses. Using a unique protocol in which docking by GOLD-PLP, rescoring by ANI-ML potentials and extensive MD simulations along with end state free energy methods are combined, we have screened FDA approved drugs against the SARS-CoV-2 main protease (Mpro). The top six drug molecules suggested by the consensus of these free energy methods have already been in clinical trials or proposed as potential drug molecules in previous theoretical and experimental studies, approving the validity and the power of accuracy in our screening method.

Computational analysis of substrate recognition of Sars-Cov-2 Mpro main protease.

Mpro main protease takes an essential role in the Sars-Cov-2 viral life cycle by releasing the individual protein from the single poly-peptide chain via proteolytic cleavage in the beginning of the viral infection. Interfering with this step by inhibiting the protease with small compound-based inhibitors has been proven to be an effective strategy to treat the infection. Thus, understanding the substrate recognition mechanism of the Mpro main protease has gained great interest from the beginning of the pandemic. Here, we have studied the substrate recognition mechanism of the protease by means of the molecular dynamic methods. We have found that the glutamine residue at P1 has paramount effect in the interaction with the substrates as expected. In addition, we also have shown that for the first time, the arginine amino acid at the P3-P5 along with P4' can strengthen the interaction.

Exploring potential Plasmodium kinase inhibitors: a combined docking, MD and QSAR studies.

Malaria is a disease caused mostly by Plasmodium falciparum, affects millions of people each year. The kinases are validated targets for malaria infection. In this study, we investigate for real and hypothetical compounds that can inhibit cyclic guanosine monophosphate (CGMP)-dependent protein kinase using molecular docking via combined similarity analysis, molecular dynamics simulations, quantitative structure activity relationship (QSAR). Using Tanimoto similarity scores, ∼8.4 million compounds were screened. Compounds that have at least 70% similarity are used in further analysis. These compounds are assessed by means of docking, MMBPSA, MMGBSA and ANI_LIE. Based on consensus of different free energy methods and docking we revealed two potential inhibitors that can be useful for treatment of malaria. Apart from screening of real compounds, we have also selected the 10 most plausible hypothetical compounds by performing QSAR. By QSAR proposed pharmacophores, we generated over 247 hypothetical compounds and among them 19 molecules with lower QSAR predicted IC50 values and high docking scores were selected for further analysis. We selected the top 10 inhibitor candidates and performed MD simulations for free energy calculations like the protocol applied for real compounds. According to the free energy calculations, we suggest 2 real (C34H29F5N8O4S and C30H27F2N7O2S2, PubChem IDs: 140564801 and 89035196, respectively) and 2 hypothetical (C23H27FN6O2S, MOL3 and C23H25FN6O2S, MOL4) compounds that can be effective inhibitors against the protein kinase of Plasmodium falciparum.Communicated by Ramaswamy H. Sarma.

Revisiting MMPBSA by Adoption of MC-Based Surface Area/Volume, ANI-ML Potentials, and Two-Valued Interior Dielectric Constant.

Here, we report the accuracy improvements of molecular mechanics Poisson-Boltzmann surface area (MMPBSA) calculations by adoption of ANI-ML potentials in replacement of MM terms, the use of solvent-accessible surface area (SASA) and volume (SAV) values from the Monte Carlo sampling of the probe, and introducing two different interior dielectric constants for electrostatic interactions of protein-ligand (P-L) and polar solvation term in the MMPBSA calculations. Our results show that the Pearson correlation coefficients of MMPBSA-calculated values with respect to experimental binding free energies can be drastically improved from 0.48 to 0.90 by adoption of ANI-ML potentials in replacement of MM energy terms in the equation, referred to as ANI-PBSA. Moreover, we show that the SASA/SAV-combined equation in the scaled particle theory (SPT) can be a better choice to model nonpolar solvation term, reaching nearly the same accuracy by ANI-PBSA calculations. Finally, we introduce two different values of interior dielectric constants, which could be an alternative strategy between the single and variable constant definitions.

Accurate Binding Free Energy Method from End-State MD Simulations.

Herein, we introduce a new strategy to estimate binding free energies using end-state molecular dynamics simulation trajectories. The method is adopted from linear interaction energy (LIE) and ANI-2x neural network potentials (machine learning) for the atomic simulation environment (ASE). It predicts the single-point interaction energies between ligand-protein and ligand-solvent pairs at the accuracy of the wb97x/6-31G* level for the conformational space that is sampled by molecular dynamics (MD) simulations. Our results on 54 protein-ligand complexes show that the method can be accurate and have a correlation of R = 0.87-0.88 to the experimental binding free energies, outperforming current end-state methods with reduced computational cost. The method also allows us to compare BFEs of ligands with different scaffolds. The code is available free of charge (documentation and test files) at https://github.com/otayfuroglu/deepQM.

Molecular dynamics simulations reveal the plausible agonism/antagonism mechanism by steroids on androgen receptor mutations.

Androgen receptors (AR) are the primary drug target in prostate cancer (PCa). There are several drugs developed against its activity for prostate cancer treatment, but cancer cells revive AR signaling against those drugs by using alternative steroids such as glucocorticoids. In addition, antagonists become agonists due to emergence of mutations in AR gene. The mechanism by which antagonists are converted into agonists and how AR signaling is recovered by other steroids has yet to be fully elucidated. In this study, we interrogated the role of bicalutamide conformation in its antagonist function and how glucocorticoids such as prednisolone and dexamethasone revive AR signaling at the molecular level by means of molecular dynamics. We found that the ''closed'' conformation of bicalutamide is essential for its antagonist function and W741 residue is forcing it into this conformation. Moreover, we show that prednisolone and dexamethasone behave like natural agonist DHT which confirm the experimental results that show their role in the reviving AR signaling in the case of ARL701H mutation.

Cysteine Disulfide Traps Reveal Distinct Conformational Ensembles in Dengue Virus NS2B-NS3 Protease.

The dengue virus protease (NS2B-NS3pro) plays a critical role in the dengue viral life cycle, making it an attractive drug target for dengue-related pathologies, including dengue hemorrhagic fever. A number of studies indicate that NS2B-NS3pro undergoes a transition between two widely different conformational states: an "open" (inactive) conformation and a "closed" (active) conformation. For the past several years, the equilibrium between these states and the resting conformation of NS2B-NS3pro have been debated, although a strong consensus is emerging. To investigate the importance of such conformational states, we developed versions of NS2B-NS3pro that allow us to trap the enzyme in various distinct conformations. Our data from these variants suggest that the enzymatic activity appears to be dependent on the movement of NS2B and may rely on the flexibility of the protease core. Locking the enzyme into the "closed" conformation dramatically increased activity, strongly suggesting that the "closed" conformation is the active conformation. The observed resting state of the enzyme depends largely on the construct used to express the NS2B-NS3pro complex. In an "unlinked" construct, in which the NS2B and NS3 regions exist as independent, co-expressed polypeptides, the enzyme rests predominantly in a "closed", active conformation. In contrast, in a "linked" construct, in which NS2B and NS3 are attached by a nine-amino acid linker, NS2B-NS3pro adopts a more relaxed, alternative conformation. Nevertheless, even the unlinked construct samples both the "closed" and other alternative conformations. Given our findings, and the more realistic resemblance of NS2B-NS3pro to the native enzyme, these data strongly suggest that studies should focus on the "unlinked" constructs moving forward. Additionally, the results from these studies provide a more detailed understanding of the various poses of the dengue virus NS2B-NS3 protease and should help guide future drug discovery efforts aimed at this enzyme.

Molecular Dynamics Studies of Histo-Blood Group Antigen Blocking Human Immunoglobulin A Antibody and Escape Mechanism in Noroviruses Upon Mutation.

Norovirus is the causing agent of acute gastroenteritis disease globally. Efforts in developing therapeutics against virus infection mostly fail due to emergence of drug resistance that is a consequence of presence of high mutation rates in virus genome during virus' life cycle. In this study, we computationally analyzed the affinity of a drug target, wild type VP1 envelope protein and its three variants to a therapeutic antibody FAB5I2. We have found that mutations break important hydrogen bonds and cause high fluctuations in residues that form VP1-FAB5I2 complex interface. In addition to changes in dynamics, we also revealed that the affinity of FAB5I2 to VP1 protein drops significantly upon mutations in terms of relative binding free energy.

Docking, molecular dynamics and free energy studies on aspartoacylase mutations involved in Canavan disease.

The disruption of aspartoacylase enzyme's catalytic activity causes fatal neurodegenerative Canavan disease. By molecular dynamics and docking methods, here we studied two deleterious mutations that have been identified in the Canavan patients' genotype E285A, F295S, and revealed the possible cause for the enzyme inhibition due to the drastic changes in active site dynamics, loss of interactions among Arg 71, Arg 168 and the substrate and pKa value of critical Glu178 residue. In addition to changes in the enzyme dynamics, free energy calculations show that the binding energy of substrate decreases dramatically up on mutations.

Computational insights into the protonation states of catalytic dyad in BACE1-acyl guanidine based inhibitor complex.

Developing small compound based drugs targeting the ß-secretase (BACE) enzyme is one of the most promising strategies in treatment of the Alzheimer's disease. As the enzyme shows the activity based on the acid-base reaction at a very narrow pH range, the protonation state of aspartic acids with the residue number 32 and 228 (Asp32 and Asp228), which forms the active site dyad, along with the protonation state of the ligand (substrate or inhibitor) play very critical role in interactions between the ligand and enzyme. Thus, understanding the nature of the protonation state of both enzyme's active site dyad and ligand is crucial for drug design in Alzheimer's disease field. Here we have investigated the protonation state of the Asp32 and Asp228 residues in the presence of a highly potent beta secretase inhibitor, containing acyl guanidine warhead that have recently been devised but not extensively studied. Our Quantum Mechanical, Molecular Dynamics and Docking studies on all the possible protonation states have suggested that the dyad residues are in di-deprotonated states in the presence of protonated inhibitor.

The impact of metabolic syndrome on serum total testosterone level in patients with erectile dysfunction.

OBJECTIVES: To determine the association between metabolic syndrome (MetS) and serum testosterone levels (TT) in patients with erectile dysfunction (ED). METHODS: This study included 280 ED patients above 40-years-of-age. Participants were divided into two groups according to 2005 criteria of International Diabetes Federation. The severity of ED was determined according to the International Index of Erectile Function-EF (IIEF-EF score; 0-10 severe ED, 11-25 mild to moderate ED). The severity of ED, serum TT levels and other MetS components were compared between the groups. RESULTS: The mean age of the patients was 55.7 ± 8.2 years. One hundred eighteen patients (%42.1) had MetS. Sixty-eight patients with MetS (57.6%) and 71 patients without MetS (43.8%) had severe ED (p = 0.031). A total of 46 (16.4%) patients had hypogonadism. Hypogonadism was seen more prevalent in patients with MetS (22.9% vs. 11.7%, p = 0.013). Logistic regression analyses for ED risk factors demonstrated that abnormal FBG increased the relative risk of severe ED up to 10.7-fold (p < 0.001) but not presence of hypogonadism (p = 0.706). CONCLUSION: Metabolic syndrome was seen in almost half of the patients with ED. ED was more severe among MetS patients. Hypogonadism alone is a not risk factor for severe ED.

Allosteric inhibition of the NS2B-NS3 protease from dengue virus.

Dengue virus is the flavivirus that causes dengue fever, dengue hemorrhagic disease, and dengue shock syndrome, which are currently increasing in incidence worldwide. Dengue virus protease (NS2B-NS3pro) is essential for dengue virus infection and is thus a target of therapeutic interest. To date, attention has focused on developing active-site inhibitors of NS2B-NS3pro. The flat and charged nature of the NS2B-NS3pro active site may contribute to difficulties in developing inhibitors and suggests that a strategy of identifying allosteric sites may be useful. We report an approach that allowed us to scan the NS2B-NS3pro surface by cysteine mutagenesis and use cysteine reactive probes to identify regions of the protein that are susceptible to allosteric inhibition. This method identified a new allosteric site utilizing a circumscribed panel of just eight cysteine variants and only five cysteine reactive probes. The allosterically sensitive site is centered at Ala125, between the 120s loop and the 150s loop. The crystal structures of WT and modified NS2B-NS3pro demonstrate that the 120s loop is flexible. Our work suggests that binding at this site prevents a conformational rearrangement of the NS2B region of the protein, which is required for activation. Preventing this movement locks the protein into the open, inactive conformation, suggesting that this site may be useful in the future development of therapeutic allosteric inhibitors.