Computational Investigation of BMAA and Its Carbamate Adducts as Potential GluR2 Modulators.

Beta-N-methylamino-l-alanine (BMAA) is a potential neurotoxic nonprotein amino acid, which can reach the human body through the food chain. When BMAA interacts with bicarbonate in the human body, carbamate adducts are produced, which share a high structural similarity with the neurotransmitter glutamate. It is believed that BMAA and its l-carbamate adducts bind in the glutamate binding site of ionotropic glutamate receptor 2 (GluR2). Chronic exposure to BMAA and its adducts could cause neurological illness such as neurodegenerative diseases. However, the mechanism of BMAA action and its carbamate adducts bound to GluR2 has not yet been elucidated. Here, we investigate the binding modes and the affinity of BMAA and its carbamate adducts to GluR2 in comparison to the natural agonist, glutamate, to understand whether these can act as GluR2 modulators. Initially, we perform molecular dynamics simulations of BMAA and its carbamate adducts bound to GluR2 to examine the stability of the ligands in the S1/S2 ligand-binding core of the receptor. In addition, we utilize alchemical free energy calculations to compute the difference in the free energy of binding of the beta-carbamate adduct of BMAA to GluR2 compared to that of glutamate. Our findings indicate that carbamate adducts of BMAA and glutamate remain stable in the binding site of the GluR2 compared to BMAA. Additionally, alchemical free energy results reveal that glutamate and the beta-carbamate adduct of BMAA have comparable binding affinity to the GluR2. These results provide a rationale that BMAA carbamate adducts may be, in fact, the modulators of GluR2 and not BMAA itself.

Natural Compounds for Bone Remodeling: A Computational and Experimental Approach Targeting Bone Metabolism-Related Proteins.

Osteoporosis, characterized by reduced bone density and increased fracture risk, affects over 200 million people worldwide, predominantly older adults and postmenopausal women. The disruption of the balance between bone-forming osteoblasts and bone-resorbing osteoclasts underlies osteoporosis pathophysiology. Standard treatment includes lifestyle modifications, calcium and vitamin D supplementation and specific drugs that either inhibit osteoclasts or stimulate osteoblasts. However, these treatments have limitations, including side effects and compliance issues. Natural products have emerged as potential osteoporosis therapeutics, but their mechanisms of action remain poorly understood. In this study, we investigate the efficacy of natural compounds in modulating molecular targets relevant to osteoporosis, focusing on the Mitogen-Activated Protein Kinase (MAPK) pathway and the gut microbiome's influence on bone homeostasis. Using an in silico and in vitro methodology, we have identified quercetin as a promising candidate in modulating MAPK activity, offering a potential therapeutic perspective for osteoporosis treatment.

Comparative Interaction Studies of Quercetin with 2-Hydroxyl-propyl-ß-cyclodextrin and 2,6-Methylated-ß-cyclodextrin.

Quercetin (QUE) is a well-known natural product that can exert beneficial properties on human health. However, due to its low solubility its bioavailability is limited. In the present study, we examine whether its formulation with two cyclodextrins (CDs) may enhance its pharmacological profile. Comparative interaction studies of quercetin with 2-hydroxyl-propyl-ß-cyclodextrin (2HP-ß-CD) and 2,6-methylated cyclodextrin (2,6Me-ß-CD) were performed using NMR spectroscopy, DFT calculations, and in silico molecular dynamics (MD) simulations. Using T1 relaxation experiments and 2D DOSY it was illustrated that both cyclodextrin vehicles can host quercetin. Quantum mechanical calculations showed the formation of hydrogen bonds between QUE with 2HP-ß-CD and 2,6Μe-ß-CD. Six hydrogen bonds are formed ranging between 2 to 2.8 Å with 2HP-ß-CD and four hydrogen bonds within 2.8 Å with 2,6Μe-ß-CD. Calculations of absolute binding free energies show that quercetin binds favorably to both 2,6Me-ß-CD and 2HP-ß-CD. MM/GBSA results show equally favorable binding of quercetin in the two CDs. Fluorescence spectroscopy shows moderate binding of quercetin in 2HP-ß-CD (520 M-1) and 2,6Me-ß-CD (770 M-1). Thus, we propose that both formulations (2HP-ß-CD:quercetin, 2,6Me-ß-CD:quercetin) could be further explored and exploited as small molecule carriers in biological studies.

FEPrepare: A Web-Based Tool for Automating the Setup of Relative Binding Free Energy Calculations.

Relative binding free energy calculations in drug design are becoming a useful tool in facilitating lead binding affinity optimization in a cost- and time-efficient manner. However, they have been limited by technical challenges such as the manual creation of large numbers of input files to set up, run, and analyze free energy simulations. In this Application Note, we describe FEPrepare, a novel web-based tool, which automates the setup procedure for relative binding FEP calculations for the dual-topology scheme of NAMD, one of the major MD engines, using OPLS-AA force field topology and parameter files. FEPrepare provides the user with all necessary files needed to run a FEP/MD simulation with NAMD. FEPrepare can be accessed and used at https://feprepare.vi-seem.eu/.

The SAMPL6 SAMPLing challenge: assessing the reliability and efficiency of binding free energy calculations.

Approaches for computing small molecule binding free energies based on molecular simulations are now regularly being employed by academic and industry practitioners to study receptor-ligand systems and prioritize the synthesis of small molecules for ligand design. Given the variety of methods and implementations available, it is natural to ask how the convergence rates and final predictions of these methods compare. In this study, we describe the concept and results for the SAMPL6 SAMPLing challenge, the first challenge from the SAMPL series focusing on the assessment of convergence properties and reproducibility of binding free energy methodologies. We provided parameter files, partial charges, and multiple initial geometries for two octa-acid (OA) and one cucurbit[8]uril (CB8) host-guest systems. Participants submitted binding free energy predictions as a function of the number of force and energy evaluations for seven different alchemical and physical-pathway (i.e., potential of mean force and weighted ensemble of trajectories) methodologies implemented with the GROMACS, AMBER, NAMD, or OpenMM simulation engines. To rank the methods, we developed an efficiency statistic based on bias and variance of the free energy estimates. For the two small OA binders, the free energy estimates computed with alchemical and potential of mean force approaches show relatively similar variance and bias as a function of the number of energy/force evaluations, with the attach-pull-release (APR), GROMACS expanded ensemble, and NAMD double decoupling submissions obtaining the greatest efficiency. The differences between the methods increase when analyzing the CB8-quinine system, where both the guest size and correlation times for system dynamics are greater. For this system, nonequilibrium switching (GROMACS/NS-DS/SB) obtained the overall highest efficiency. Surprisingly, the results suggest that specifying force field parameters and partial charges is insufficient to generally ensure reproducibility, and we observe differences between seemingly converged predictions ranging approximately from 0.3 to 1.0 kcal/mol, even with almost identical simulations parameters and system setup (e.g., Lennard-Jones cutoff, ionic composition). Further work will be required to completely identify the exact source of these discrepancies. Among the conclusions emerging from the data, we found that Hamiltonian replica exchange-while displaying very small variance-can be affected by a slowly-decaying bias that depends on the initial population of the replicas, that bidirectional estimators are significantly more efficient than unidirectional estimators for nonequilibrium free energy calculations for systems considered, and that the Berendsen barostat introduces non-negligible artifacts in expanded ensemble simulations.

Blinded predictions of standard binding free energies: lessons learned from the SAMPL6 challenge.

In the context of the SAMPL6 challenges, series of blinded predictions of standard binding free energies were made with the SOMD software for a dataset of 27 host-guest systems featuring two octa-acids hosts (OA and TEMOA) and a cucurbituril ring (CB8) host. Three different models were used, ModelA computes the free energy of binding based on a double annihilation technique; ModelB additionally takes into account long-range dispersion and standard state corrections; ModelC additionally introduces an empirical correction term derived from a regression analysis of SAMPL5 predictions previously made with SOMD. The performance of each model was evaluated with two different setups; buffer explicitly matches the ionic strength from the binding assays, whereas no-buffer merely neutralizes the host-guest net charge with counter-ions. ModelC/no-buffer shows the lowest mean-unsigned error for the overall dataset (MUE 1.29 < 1.39 < 1.50 kcal mol-1, 95% CI), while explicit modelling of the buffer improves significantly results for the CB8 host only. Correlation with experimental data ranges from excellent for the host TEMOA (R2 0.91 < 0.94 < 0.96), to poor for CB8 (R2 0.04 < 0.12 < 0.23). Further investigations indicate a pronounced dependence of the binding free energies on the modelled ionic strength, and variable reproducibility of the binding free energies between different simulation packages.

Alchemical Free Energy Calculations on Membrane-Associated Proteins.

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Energetics of a protein disorder-order transition in small molecule recognition.

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered 'lid' region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the apparent dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

The effect of S427F mutation on RXRα activity depends on its dimeric partner.

RXRs are nuclear receptors acting as transcription regulators that control key cellular processes in all tissues. All type II nuclear receptors require RXRs for transcriptional activity by forming heterodimeric complexes. Recent whole-exome sequencing studies have identified the RXRα S427F hotspot mutation in 5% of the bladder cancer patients, which is always located at the interface of RXRα with its obligatory dimerization partners. Here, we show that mutation of S427 deregulates transcriptional activity of RXRα dimers, albeit with diverse allosteric mechanisms of action depending on its dimeric partner. S427F acts by allosteric mechanisms, which range from inducing the collapse of the binding pocket to allosteric stabilization of active co-activator competent RXRα states. Unexpectedly, RXR S427F heterodimerization leads to either loss- or gain-of-function complexes, in both cases likely compromising its tumor suppressor activity. This is the first report of a cancer-associated single amino acid substitution that affects the function of the mutant protein variably depending on its dimerization partner.