Binding Free Energy Calculations Research Articles

Calculating Binding Free Energies in Model Host-Guest Systems with Unrestrained Advanced Sampling.

Host-guest interactions are important to the design of pharmaceuticals and, more broadly, to soft materials as they can enable targeted, strong, and specific interactions between molecules. The binding process between the host and guest may be classified as a "rare event" when viewing the system at atomic scales, such as those explored in molecular dynamics simulations. To obtain equilibrium binding conformations and dissociation constants from these simulations, it is essential to resolve these rare events. Advanced sampling methods such as the adaptive biasing force (ABF) promote the occurrence of less probable configurations in a system, therefore facilitating the sampling of essential collective variables that characterize the host-guest interactions. Here, we present the application of ABF to a rod-cavitand coarse-grained model of host-guest systems to acquire the potential of mean force. We show that the employment of ABF enables the computation of the configurational and thermodynamic properties of bound and unbound states, including the free energy landscape. Moreover, we identify important dynamic bottlenecks that limit sampling and discuss how these may be addressed in more general systems.

In-Silico study of Destabilizing Alzheimer’s Aβ42 Protofibrils with Curcumin.

Aggregation of amyloid beta (Aβ) peptides leads to formation of fibrilar, soluble oligomers, and their deposition is a key event in progression of Alzheimer’s disease (AD). Recent experimental studies of (Curcumin) showed\ significant Aβ aggregation inhibition, but its molecular mechanism is not yet clear. Hence, the present study aims at exploring the underlying mechanism of destabilization and inhibition of aggregation of the Aβ protofibril by Curcumin at the molecular level. Molecular docking analysis shows that Curcumin binds to chain E of the Aβ protofibril through hydrogen bonding interactions. Comparative molecular dynamics simulations depict the binding of Curcumin at the edge of chain E, and its partially inserted conformation at the hydrophobic core destabilizes the Aβ protofibril. Its binding causes loss of hydrophobic contacts, leading to a partial opening of tightly packed β-sheet protofibrils. The hydration effect of salt bridge between the amino group of Lys28 and the oxygen atom of Curcumin contributes in destabilization of Aβ protofibrils. Binding free energy calculations of Curcumin and the Aβ protofibril showed that van der Waals interactions are dominant over the others. Thus, our results revealed that Curcumin interacts mainly with the hydrophobic core along with positively charged residues of the Aβ protofibril for effective destabilization. Thus, this structural information could be useful to design new inhibitors to control the aggregation of Aβ protofibrils in AD patients.

Isolation, Virtual Screening, and Evaluation of Hazelnut-Derived Immunoactive Peptides for the Inhibition of SARS-CoV-2 Main Protease.

The aim of this study is to validate the activity of hazelnut (Corylus avellana L.)-derived immunoactive peptides inhibiting the main protease (Mpro) of SARS-CoV-2 and further unveil their interaction mechanism using in vitro assays, molecular dynamics (MD) simulations, and binding free energy calculations. In general, the enzymatic hydrolysis components, especially molecular weight < 3 kDa, possess good immune activity as measured by the proliferation ability of mouse splenic lymphocytes and phagocytic activity of mouse peritoneal macrophages. Over 866 unique peptide sequences were isolated, purified, and then identified by nanohigh-performance liquid chromatography/tandem mass spectrometry (NANO-HPLC-MS/MS) from hazelnut protein hydrolysates, but Trp-Trp-Asn-Leu-Asn (WWNLN) and Trp-Ala-Val-Leu-Lys (WAVLK) in particular are found to increase the cell viability and phagocytic capacity of RAW264.7 macrophages as well as promote the secretion of the cytokines nitric oxide (NO), tumor necrosis factor-α (TNF-α), and interleukin-1β (IL-1β). Fluorescence resonance energy transfer assay elucidated that WWNLN and WAVLK exhibit excellent inhibitory potency against Mpro, with IC50 values of 6.695 and 16.750 μM, respectively. Classical all-atom MD simulations show that hydrogen bonds play a pivotal role in stabilizing the complex conformation and protein-peptide interaction. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculation indicates that WWNLN has a lower binding free energy with Mpro than WAVLK. Furthermore, adsorption, distribution, metabolism, excretion, and toxicity (ADMET) predictions illustrate favorable drug-likeness and pharmacokinetic properties of WWNLN compared to WAVLK. This study provides a new understanding of the immunomodulatory activity of hazelnut hydrolysates and sheds light on peptide inhibitors targeting Mpro.

Impact of protein conformations on binding free energy calculations in the beta-secretase 1 system.

In binding free energy calculations, simulations must sample all relevant conformations of the system in order to obtain unbiased results. For instance, different ligands can bind to different metastable states of a protein, and if these protein conformational changes are not sampled in relative binding free energy calculations, the contribution of these states to binding is not accounted for and thus calculated binding free energies are inaccurate. In this work, we investigate the impact of different beta-sectretase 1 (BACE1) protein conformations obtained from x-ray crystallography on the binding of BACE1 inhibitors. We highlight how these conformational changes are not adequately sampled in typical molecular dynamics simulations. Furthermore, we show that insufficient sampling of relevant conformations induces substantial error in relative binding free energy calculations, as judged by a variation in calculated relative binding free energies up to 2kcal/mol depending on the starting protein conformation. These results emphasize the importance of protein conformational sampling and pose this BACE1 system as a challenge case for further method development in the area of enhanced protein conformational sampling, either in combination with binding calculations or as an endpoint correction.

MOLECULAR DOCKING DYNAMICS OF SELECTED BENZYLIDENE AMINO PHENYL ACETAMIDES AS TMK INHIBITORS USING HIGH THROUGHPUT VIRTUAL SCREENING (HTVS)

Objective: Thymidylate kinase (TMK) plays a crucial role in bacterial DNA synthesis by catalyzing the phosphorylation of deoxythymidine monophosphate (dTMP) to form deoxythymidine diphosphate (dTDP). Consequently, this enzyme emerges as a promising target for developing novel antibacterial drugs. However, no antibiotics were reported for this target, especially active against Staphylococcus aureus thymidylate kinase. Methods: Benzylidene acetamide-based ligands were examined for their potency using the in silico method. These novel ligand structures were built using ChemDraw software. The protein was retrieved from the Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB) website. The molecular docking and binding free energy calculation by prime Molecular Mechanics in Generalized Bond Surface Area (MM-GBSA) was performed for selected ligands. A 100 ns molecular dynamic simulation was also performed to assess the stability of the potential ligand as TMK inhibitors. Results: All ten molecules have shown good glide scores and hydrophobic and hydrogen hydrophobic hydrogen bonding interactions with Arg48, Arg36, and π-π stacking Phe66 in the TMK enzyme (PDB: 4HLC). Among them, N-(2-ethylphenyl)-2-(4-((4-nitrobenzylidene) amino) phenoxy) acetamide molecule had high XP-docking scores of-3.27 kcal/mol based on extra-precision data. Prime Molecular Mechanics in Generalized Bond Surface Area study (MM-GBSA) studies also showed promising binding affinities that are ΔBind (-65.80), ΔLipo (-28.55), and ΔVdW (-55.10). Phe66 amino acid residue maintained continuous connections with the ligand during MD simulation. This ligand showed promising binding affinity with the SaTMK target. Conclusion: The N-(2-ethylphenyl)-2-(4-((4-nitrobenzylidene) amino) phenoxy) acetamide ligand at the position of the benzene ring displayed nitrogen and oxygen group, thus indicating good potential activity as the inhibitor of TMK to treat antibacterial agents.

Evaluation of surfactant-aided polycyclic aromatic hydrocarbon biodegradation by molecular docking and molecular dynamic simulation in the marine environment

Marine oil spills directly cause polycyclic aromatic hydrocarbons (PAHs) pollution and affect marine organisms due to their toxic property. Chemical and bio-based dispersants composed of surfactants and solvents are considered effective oil spill-treating agents. Dispersants enhance oil biodegradation in the marine environment by rapidly increasing their solubility in the water column. However, the effect of dispersants, especially surfactants, on PAHs degradation by enzymes produced by microorganisms has not been studied at the molecular level. The role of the cytochrome P450 (CYP) enzyme in converting contaminants into reactive metabolites during the biodegradation process has been evidenced, but the activity in the presence of surfactants is still ambiguous. Thus, this study focused on the evaluation of the impact of chemical and bio-surfactants (i.e., Tween 80 (TWE) and Surfactin (SUC)) on the biodegradation of naphthalene (NAP), chrysene (CHR), and pyrene (PYR), the representative components of PAHs, with CYP enzyme from microalgae Parachlorella kessleri using molecular docking and molecular dynamics (MD) simulation. The molecular docking analysis revealed that PAHs bound to residues at the CYP active site through hydrophobic interactions for biodegradation. The MD simulation showed that the surfactant addition changed the enzyme conformation in the CYP-PAH complexes to provide more interactions between the enzyme and PAHs. This led to an increase in the enzyme's capability to degrade PAHs. Binding free energy (ΔG‌‌Bind) calculations confirmed that surfactant treatment could enhance PAHs degradation by the enzyme. The SUC gave a better result on NAP and PYR biodegradation based on ΔG‌‌Bind, while TWE facilitated the biodegradation of CHR. The research outputs could greatly facilitate evaluating the behaviors of oil spill-treating agents and oil spill response operations in the marine environment.

Phloretin Prolongs Lifespan of Caenorhabditis elegans via Inhibition of NDUFS1 and NDUFS6 at Mitochondrial Complex Ⅰ

Phloretin has been widely perceived as an antioxidant. However, the bioavailability of phloretin in vivo is generally far too low to elicit a direct antioxidant effect by scavenging reactive oxygen species (ROS). Here we showed that administration of phloretin of apple polyphenols extended lifespan of Caenorhabditis elegans and promoted fitness. Specially phloretin enhanced the survival rates of nematodes under oxidants in an inverted U-shaped dose-response manner. The lifespan-extending effects of phloretin were mediated by ROS via complex I inhibition. The increase of ROS stimulated p38 MAPK/PMK-1 as well as transcription factors of NRF2/SKN-1 and FOXO/DAF-16. Consistent with the involvement of NRF2/SKN-1 and FOXO/DAF-16 in lifespan-extending effects, activities of SOD and CAT were enhanced by phloretin. The exogenous application of antioxidants BHA and NAC abolished the increase of ROS, the enhancement of SOD and CAT activities, and the lifespan extending effects of phloretin. Meanwhile, with the inhibition of mitochondrial complex I, ATP was instantly decreased. Both energy sensors of AMPK/AAK-2 and SIRT1/SIR-2.1 were involved in the lifespan extension by phloretin. Transcriptomic, real-time qPCR and molecular docking analyses demonstrated that the binding of phloretin at complex I located at NDUFS1/NUO-5, NDUFS2/GAS-1, and NDUFS6/NDUF-6. The molecular dynamic simulation and binding free energy calculations showed that phloretin had high binding affinities towards NDUFS1 (-7.21 kcal/mol) and NDUFS6 (-7.02 kcal/mol). Collectively, our findings suggested phloretin had effects of life expectancy enhancement and fitness promotion via redox regulations in vivo. NDUFS1/NUO-5 and NDUFS6/NDUF-6 might be new targets in the lifespan and wellness regulations.

In silico analysis and characterization of potential inhibitors of MmaA3, a methoxy mycolic acid synthase from Mycobacterium tuberculosis

The emergence of the multi-and extensively drug-resistant (MDR and XDR) strains of Mycobacterium tuberculosis (M.tb), necessitates paradigm-shifting therapeutic approaches. The impermeable waxy lipid layer, primarily composed of mycolic acids, is a key factor in conferring resistance to conventional drugs. This study introduces a novel strategy to combat drug resistance by targeting Methoxy mycolic acid synthase 3 (MmaA3), a critical enzyme in the mycolic acid biosynthesis pathway. MmaA3 is responsible for the O-methylation of hydroxymycolate precursors and emerges as a promising therapeutic target. Through homology-based modeling, we generated a three-dimensional structure of MmaA3, providing crucial insights into its structural characteristics. High throughput virtual screening was performed against the MmaA3 model, using diverse sources: knowledge-based, FDA-approved Drugbank, and Asinex-Elite libraries. Through rigorous computational analyses, including binding affinity assessments, molecular interactions analysis, and binding free energy calculations, potential inhibitors of MmaA3 have been identified. Subsequent validation studies evaluated the stability of top protein-ligand complexes, and free energy calculations using molecular dynamics simulations. The stability of complexes within the catalytic site was confirmed through RMSD and RMSF profile analyses. Furthermore, binding free energy calculations using the MM-GBSA approach revealed significant binding affinity of identified ligands for MmaA3 target protein, comparable to its substrate/cofactors. These findings underscore the potential of the proposed molecules as candidates for further experimental exploration, offering promising avenues for the development of effective inhibitors against M.tb. Overall, our research contributes to significantly advancing the formulation of progressive therapeutic strategies in combating drug-resistant tuberculosis. Communicated by Ramaswamy H. Sarma

In silico evaluation of anti-colorectal cancer inhibitors by Resveratrol derivatives targeting Armadillo repeats domain of APC: molecular docking and molecular dynamics simulation.

Colorectal cancer is the second leading cause of cancer-related deaths. In 2018, there were an estimated 1.8 million cases, and this number is expected to increase to 2.2 million by 2030. Despite its prevalence, the current therapeutic option has a lot of side effects and limitations. Therefore, this study was designed to employ a computational approach for the identification of anti-cancer inhibitors against colorectal cancer using Resveratrol derivatives. Initially, the pass prediction spectrum of 50 derivatives was conducted and selected top seven compounds based on the maximum pass prediction score. After that, a comprehensive analysis, including Lipinski Rule, pharmacokinetics, ADMET profile study, molecular orbitals analysis, molecular docking, molecular dynamic simulations, and MM-PBSA binding free energy calculations. The reported binding affinity ranges of Resveratrol derivatives from molecular docking were -6.1 kcal/mol to -7.9 kcal/mol against the targeted receptor of human armadillo repeats domain of adenomatous polyposis coli (APC) (PDB ID: 3NMW). Specifically, our findings reported that two compounds [(03) Resveratrol 3-beta-mono-D-glucoside, and (29) Resveratrol 3-Glucoside] displayed the highest level of effectiveness compared to all other derivatives (-7.7 kcal/mol and -7.9 kcal/mol), and favorable drug-likeness, and exceptional safety profiles. Importantly, almost all the molecules were reported as free from toxic effects. Subsequently, molecular dynamic simulations conducted over 100ns confirmed the stability of the top two ligand-protein complexes. These findings suggest that Resveratrol derivatives may be effective drug candidate to manage the colorectal cancer. However, further experimental research, such as in vitro/in vivo studies, is essential to validate these computational findings and confirm their practical value.

A structure-based pharmacophore modelling approach to identify and design new neprilysin (NEP) inhibitors: An in silico-based investigation

Neutral endopeptidase or neprilysin (NEP) cleaves the natriuretic peptides, bradykinin, endothelin, angiotensin II, amyloid β protein, substance P, etc., thus modulating their effects on heart, kidney, and other organs. NEP has a proven role in hypertension, heart disease, renal disease, Alzheimer's, diabetes, and some cancers. NEP inhibitor development has been in focus since the US FDA approved a combination therapy of angiotensin II type 1 receptor inhibitor (valsartan) and NEP inhibitor (sacubitril) for use in heart failure. Considering the importance of NEP inhibitors the present work focuses on the designing of a potential lead for NEP inhibition. A structure-based pharmacophore modelling approach was employed to identify NEP inhibitors from the pool of 1140 chemical entities obtained from the ZINC database. Based on the docking score and pivotal interactions, ten molecules were selected and subjected to binding free energy calculations and ADMET predictions. The top two compounds were studied further by molecular dynamics simulations to determine the stability of the ligand-receptor complex. ZINC0000004684268, a phenylalanine derivative, showed affinity and complex stability comparable to sacubitril. However, in silico studies indicated that it may have poor pharmacokinetic parameters. Therefore, the molecule was optimized using bioisosteric replacements, keeping the phenylalanine moiety intact, to obtain five potential lead molecules with an acceptable pharmacokinetic profile. The works thus open up the scope to further corroborate the present in silico findings with the biological analysis.

QM, molecular docking and molecular dynamics investigation on acidic phospholipase A2 2 protein and acidic phospholipase A2 3 protein with silane dimethyl

Phytochemicals and non-nutritive plant substances can protect humans from various diseases by possessing anti-inflammatory or anti-disease properties. The Andrographis Paniculata species has been the primary focus of extensive pharmacological and phytochemical research. QM concepts are being applied scientifically in chemical and biological problems, with molecular docking and molecular dynamics techniques being crucial in-silico drug development. Absorptions in Andrographis Paniculata come from charge transfer due to visible light's energy emissions, with GCMS– results indicating its chemical components possess antibacterial properties. The PES concept aids in the intuitive understanding of mathematical concepts, while RDG recreates the complete distribution of electronic pairings, including interatomic bonding regions, lone pairs, and cores. Molecular dynamics is crucial for protein-ligand docking due to its accurate binding free energy calculations and profound understanding of the dynamic nature of protein-ligand interactions.

Structural insights into the potential binding sites of Cathepsin D using molecular modelling techniques

Alzheimer’s disease (AD) is the most prevalent type of dementia caused by the accumulation of amyloid beta (Aβ) peptides. The extracellular deposition of Aβ peptides in human AD brain causes neuronal death. Therefore, it has been found that Aβ peptide degradation is a possible therapeutic target for AD. CathD has been known to breakdown amyloid beta peptides. However, the structural role of CathD is not yet clear. Hence, for the purpose of gaining a deeper comprehension of the structure of CathD, the present computational investigation was performed using virtual screening technique to predict CathD's active site residues and substrate binding mode. Ligand-based virtual screening was implemented on small molecules from ZINC database against crystal structure of CathD. Further, molecular docking was utilised to investigate the binding mechanism of CathD with substrates and virtually screened inhibitors. Localised compounds obtained through screening performed by PyRx and AutoDock 4.2 with CathD receptor and the compounds having highest binding affinities were picked as; ZINC00601317, ZINC04214975 and ZINCC12500925 as our top choices. The hydrophobic residues Viz. Gly35, Val31, Thr34, Gly128, Ile124 and Ala13 help stabilising the CathD-ligand complexes, which in turn emphasises substrate and inhibitor selectivity. Further, MM-GBSA approach has been used to calculate binding free energy between CathD and selected compounds. Therefore, it would be beneficial to understand the active site pocket of CathD with the assistance of these discoveries. Thus, the present study would be helpful to identify active site pocket of CathD, which could be beneficial to develop novel therapeutic strategies for the AD.

Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning.

Bromodomain 4 and 9 (BRD4 and BRD9) have been regarded as important targets of drug designs in regard to the treatment of multiple diseases. In our current study, molecular dynamics (MD) simulations, deep learning (DL) and binding free energy calculations are integrated to probe the binding modes of three inhibitors (H1B, JQ1 and TVU) to BRD4 and BRD9. The MD trajectory-based DL successfully identify significant functional function domains, such as BC-loop and ZA-loop. The information from the post-processing analysis of MD simulations indicates that inhibitor binding highly influences the structural flexibility and dynamic behavior of BRD4 and BRD9. The results of the MM-GBSA calculations not only suggest that the binding ability of H1B, JQ1 and TVU to BRD9 are stronger than to BRD4, but they also verify that van der Walls interactions are the primary forces responsible for inhibitor binding. The hot spots of BRD4 and BRD9 revealed by residue-based free energy estimation provide target sites of drug design in regard to BRD4 and BRD9. This work is anticipated to provide useful theoretical aids for the development of selective inhibitors over BRD family members.

On Analytical Corrections for Restraints in Absolute Binding Free Energy Calculations.

Double decoupling absolute binding free energy simulations require an intermediate state at which the ligand is held solely by restraints in a position and orientation resembling the bound state. One possible choice consists of one distance, two angle, and three dihedral angle restraints. Here, I demonstrate that in practically all cases the analytical correction derived under the rigid rotator harmonic oscillator approximation is sufficient to account for the free energy of the restraints.

Effect of quercetin on the protein-substrate interactions in SIRT6: Insight from MD simulations

SIRT6 is of interest for its promising effect in the treatment of aging-related diseases. Studies have shown quercetin (QUE) and its derivatives have varying degrees of effect on the catalytic effect of SIRT6. In the research, the effect of QUE on the protein-substrate interaction in the SIRT6-mediated mono-ADP ribosylation system was investigated by conventional molecular dynamics (MD) simulations combined with MM/PBSA binding free energy calculations. The results show that QUE can bind stably to SIRT6 with the binding energy of −22.8 kcal/mol and further affect the atomic interaction between SIRT6 and NAD+ (or H3K9), resulting in an increased affinity between SIRT6-NAD+ and decreased SIRT6-H3K9 binding capacity. At the same time, the binding of QUE can also alter some structural characteristics of the protein, with large shifts occurring in the residue regions involving the N-terminal (residues 1–27), Rossmann fold regions (residues 55–92), and ZBD (residues 164–179). Thus, QUE shows great potential as a scaffold for the design of novel potent SIRT6 modulators.

Investigating the role of functional mutations in leucine binding to Sestrin2 in aging and age-associated degenerative pathologies using structural and molecular simulation approaches

Leucine is the native known ligand of Sestrin2 (Sesn2) and its interaction with Sesn2 is particularly important, as it influences the activity of mTOR in aging and its associated pathologies. It is important to find out how leucine interacts with Sesn2 and how mutations in the binding pocket of leucine affect the binding of leucine. Therefore, this study was committed to investigating the impact of non-synonymous mutations by incorporating a broad spectrum of simulation techniques, from molecular dynamics to free energy calculations. Our study was designed to model the atomic-scale interactions between leucine and mutant forms of Sesn2. Our results demonstrated that the interaction paradigm for the mutants has been altered thus showing a significant decline in the hydrogen bonding network. Moreover, these mutations compromised the dynamic stability by altering the conformational flexibility, sampling time, and leucine-induced structural constraints that consequently caused variation in the binding and structural stability. Molecular dynamics-based flexibility analysis revealed that the regions 217–339 and 371–380 demonstrated a higher fluctuation. Noteworthy, these regions correspond to a linker (217–339) and a loop (371–380) that cover the leucine binding cavity that is critical for the ‘latch’ mechanism in the N-terminal, which is essential for leucine binding. Further validation of reduced binding and modified internal motions caused by the mutants was obtained through binding free energy calculations, principal components analysis (PCA), and free energy landscape (FEL) analysis. By unraveling the molecular intricacies of Sesn2-leucine interactions and their mutations, we hope to pave the way for innovative strategies to combat the inevitable tide of aging and its associated diseases. Communicated by Ramaswamy H. Sarma

Long-timescale atomistic simulations uncover loss-of-function mechanisms of uncharacterized Angiogenin mutants associated with ALS

Amyotrophic Lateral Sclerosis (ALS) is a devastating neurodegenerative disease characterized by progressive degeneration of motor neurons, resulting in respiratory failure and mortality within 3–5 years. Mutations in the Angiogenin (ANG) cause loss of ribonucleolytic and nuclear translocation activities, contributing to ALS pathogenesis. This study focused on investigating two uncharacterized ANG mutations, T11S and R122H, newly identified in the Project Mine consortium. Using extensive computational analysis, including structural modeling and microsecond-timescale molecular dynamics (MD) simulations, we observed conformational changes in the catalytic residue His114 of ANG induced by T11S and R122H mutations. These alterations impaired ribonucleolytic activity, as inferred through molecular docking and binding free energy calculations. Gibbs free energy landscape and residue-residue interaction network analysis further supported our findings, revealing the energetic states and allosteric pathway from the mutated site to His114. Additionally, we assessed the binding of NCI-65828, an inhibitor of ribonucleolytic activity of ANG, and found reduced effectiveness in binding to T11S and R122H mutants when His114 assumed a non-native conformation. This highlights the crucial role of His114 and its association with ALS. Elucidating the relationship between physical structure and functional dynamics of frequently mutated ANG mutants is essential for understanding ALS pathogenesis and developing more effective therapeutic interventions.

In-silico screening and identification of glycomimetic as potential human sodium-glucose co-transporter 2 inhibitor

Sodium-glucose co-transporter 2 (SGLT2) is one of the important targets against type II diabetes mellitus. A typical SGLT2 inhibitor acts by inhibiting glucose reabsorption, thus lowering the blood glucose level. Unlike SGLT1, SGLT2 is responsible for almost 90% glucose reabsorption from glomerular filtrate. The current SGLT2 inhibitors include gliflozins, often prescribed as second or third-line agents in diabetes mellitus. The SGLT2 inhibitors also benefit patients with heart and kidney disease. Due to instability issues with the natural O-aryl glycoside analogues C-glycoside analogues were developed and showed improved stability. Despite enhanced bioavailability and selectivity of newer derivatives, some serious side effects are associated with gliflozin analogues. At the present study, we applied in-silico approaches to find new glycomimetic compounds as potent SGLT2 inhibitors that could show improvement in side effects associated with current analogues. This work applied both ligand-based and structure-based drug approaches to find potential compounds. We developed a 3D-QSAR method to screen potential inhibitors from a library of ten thousand compounds and performed docking studies. The compounds were ranked based on predicted pIC50 and docking score. An initial screening of five thousand compounds was conducted, and the subsequently selected top 12 compounds were based on binding free energy calculations. These selected compounds were subjected to molecular dynamics (MD) simulations. Remarkably, our simulations identified nine compounds that exhibited significant and sustained binding affinity compared to the co-crystallized Empagliflozin. Collectively, considering the anticipated pharmacokinetic profiles and toxicity assessments, several of these compounds emerged as promising candidates for further in-depth evaluation.

An Ensemble Docking Approach for Analyzing and Designing Aptamer Heterodimers Targeting VEGF165.

Vascular endothelial growth factor 165 (VEGF165) is a prominent isoform of the VEGF-A protein that plays a crucial role in various angiogenesis-related diseases. It is homodimeric, and each of its monomers is composed of two domains connected by a flexible linker. DNA aptamers, which have emerged as potent therapeutic molecules for many proteins with high specificity and affinity, can also work for VEGF165. A DNA aptamer heterodimer composed of monomers of V7t1 and del5-1 connected by a flexible linker (V7t1:del5-1) exhibits a greater binding affinity with VEGF165 compared to either of the two monomers alone. Although the structure of the complex formed between the aptamer heterodimer and VEGF165 is unknown due to the highly flexible linkers, gaining structural information will still be valuable for future developments. Toward this end of accessing structural information, we adopt an ensemble docking approach here. We first obtain an ensemble of structures for both VEGF165 and the aptamer heterodimer by considering both small- and large-scale motions. We then proceed through an extraction process based on ensemble docking, molecular dynamics simulations, and binding free energy calculations to predict the structures of the VEGF165/V7t1:del5-1 complex. Through the same procedures, we reach a new aptamer heterodimer that bears a locked nucleic acid-modified counterpart of V7t1, namely RNV66:del5-1, which also binds well with VEGF165. We apply the same protocol to the monomeric units V7t1, RNV66, and del5-1 to target VEGF165. We observe that V7t1:del5-1 and RNV66:del5-1 show higher binding affinities with VEGF165 than any of the monomers, consistent with experiments that support the notion that aptamer heterodimers are more effective anti-VEGF165 aptamers than monomeric aptamers. Among the five different aptamers studied here, the newly designed RNV66:del5-1 shows the highest binding affinity with VEGF165. We expect that our ensemble docking approach can help in de novo designs of homo/heterodimeric anti-angiogenic drugs to target the homodimeric VEGF165.

De novo design of drug-binding proteins with predictable binding energy and specificity.

The de novo design of small molecule-binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase-1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule-binding proteins with tuned interaction energies is feasible entirely from computation.