Impact of L718Q mutation on the inhibitory activity of osimertinib on EGFR revealed by free-energy simulations.

The duration of drug efficacy in vivo is a key aspect primarily addressed during the lead optimization phase of drug discovery. Hence, the availability of robust computational approaches that can predict the residence time of a compound at its target would accelerate candidate selection. Nowadays the theoretical prediction of this parameter is still very challenging. Starting from methods reported in the literature, we set up and validated a new metadynamics (META-D)-based protocol that was used to rank the experimental residence times of 10 arylpyrazole cyclin-dependent kinase 8 (CDK8) inhibitors for which target-bound X-ray structures are available. The application of reported methods based on the detection of the escape from the first free energy well gave a poor correlation with the experimental values. Our protocol evaluates the energetics of the whole unbinding process, accounting for multiple intermediates and transition states. Using seven collective variables (CVs) encoding both roto-translational and conformational motions of the ligand, a history-dependent biasing potential is deposited as a sum of constant-height Gaussian functions until the ligand reaches an unbound state. The time required to achieve this state is proportional to the integral of the deposited potential over the CV hyperspace. Average values of this time, for replicated META-D simulations, provided an accurate classification of CDK8 inhibitors spanning short, medium, and long residence times.

The Eph receptor-ephrin system is an emerging target for the development of novel anti-angiogenic therapies. Research programs aimed at developing small-molecule antagonists of the Eph receptors are still in their initial stage as available compounds suffer from pharmacological drawbacks, limiting their application in vitro and in vivo. In the present work, we report the design, synthesis and evaluation of structure-activity relationships of a class of Δ(5)-cholenoyl-amino acid conjugates as Eph-ephrin antagonists. As a major achievement of our exploration, we identified N-(3β-hydroxy-Δ(5)-cholen-24-oyl)-L-tryptophan (UniPR1331) as the first small molecule antagonist of the Eph-ephrin system effective as an anti-angiogenic agent in endothelial cells, bioavailable in mice by the oral route and devoid of biological activity on G protein-coupled and nuclear receptors targeted by bile acid derivatives.

No abstract

Metadynamics (META-D) is emerging as a powerful method for the computation of the multidimensional free-energy surface (FES) describing the protein-ligand binding process. Herein, the FES of unbinding of the antagonist N-(3α-hydroxy-5β-cholan-24-oyl)-l-β-homotryptophan (UniPR129) from its EphA2 receptor was reconstructed by META-D simulations. The characterization of the free-energy minima identified on this FES proposes a binding mode fully consistent with previously reported and new structure-activity relationship data. To validate this binding mode, new N-(3α-hydroxy-5β-cholan-24-oyl)-l-β-homotryptophan derivatives were designed, synthesized, and tested for their ability to displace ephrin-A1 from the EphA2 receptor. Among them, two antagonists, namely compounds 21 and 22, displayed high affinity versus the EphA2 receptor and resulted endowed with better physicochemical and pharmacokinetic properties than the parent compound. These findings highlight the importance of free-energy calculations in drug design, confirming that META-D simulations can be used to successfully design novel bioactive compounds.

The UPC ratio in dogs was minimally affected in urine samples collected by free catch, thus allowing correct grading of proteinuria with this method. The high reliability of the UPC ratio in free-catch urine samples coupled with the ease of collection should increase the use of this value for assessment of proteinuria.

The free-energy surface (FES) of protein-ligand binding contains information useful for drug design. Here we show how to exploit a free-energy minimum of a protein-ligand complex identified by metadynamics simulations to design a new EphA2 antagonist with improved inhibitory potency.

The N-terminal cysteine hydrolase N-acylethanolamine acid amidase (NAAA) catalyzes the hydrolytic deactivation of the lipid messenger palmitoylethanolamide (PEA), with optimal activity at acidic pH. Using the crystal structure of human NAAA as a starting point, we investigated the catalytic mechanism of PEA hydrolysis with a multiscale approach based on classic molecular dynamics (MD) and quantum mechanical/molecular mechanics (QM/MM) simulations coupled with enhanced sampling and path-collective variables (PCVs). The proton configuration of the catalytic nucleophile, Cys126, and of the surrounding carboxylates was critical to preserve the active site architecture. A stable Michaelis complex was then used to reconstruct the free-energy surfaces of NAAA acylation and deacylation during PEA hydrolysis. Acylation emerged as the critical step, with Cys126 acting both as an acid, to protonate the ethanolamine leaving group, and as a nucleophile, to attack the PEA carbonyl carbon. The ethanol fragment of PEA did not appear to play an indispensable role in acylation, a result further supported by kinetic experiments showing that NAAA hydrolyzes palmitoyl methyl amide (PMA) with high catalytic efficiency. Our multiscale approach identified a distinctive protonation state and catalytic mechanism for NAAA which accounts for pH-dependent activity, mutagenesis data, and mechanism of covalent inhibitors.