Anna Berteotti

Steered Molecular Dynamics Simulations for Studying Protein–Ligand Interaction in Cyclin-Dependent Kinase 5

In this study, we applied steered molecular dynamics (SMD) simulations to investigate the unbinding mechanism of nine inhibitors of the enzyme cyclin-dependent kinase 5 (CDK5). The study had two major objectives: (i) to create a correlation between the unbinding force profiles and the inhibition activities of these compounds expressed as IC50 values; (ii) to investigate the unbinding mechanism and to reveal atomistic insights, which could help identify accessory binding sites and transient interactions. Overall, we carried out 1.35 μs of cumulative SMD simulations. We showed that SMD could qualitatively discriminate binders from nonbinders, while it failed to properly rank series of inhibitors, particularly when IC50 values were too similar. From a mechanistic standpoint, SMD provided useful insights related to transient and dynamical interactions, which could complement static description obtained by X-ray crystallography experiments. In conclusion, the present study represents a further step toward a systematic exploitation of SMD and other dynamical approaches in structure-based drug design and computational medicinal chemistry.

Effect of urea on the β-hairpin conformational ensemble and protein denaturation mechanism.

Despite the daily use of urea to influence protein folding and stability, the molecular mechanism with which urea acts is still not well understood. Here the use of combined parallel tempering and metadynamics simulation allows us to study the free-energy landscape associated with the folding/unfolding of β-hairpin GB1 equilibrium in 8 M urea and pure water. The nature of the unfolded state in both solutions has been analyzed: in urea solution the addition of denaturants acts to expand the denatured state, while in pure water solution the unfolded state is noticeably more compact. For what concerns the mechanism by which urea acts as a denaturant, a preferential direct interaction between urea molecules and protein backbone has been found. However, the bias toward urea solvation is largest at intermediate values of the gyration radius.

The ligand binding mechanism to purine nucleoside phosphorylase elucidated via molecular dynamics and machine learning

The study of biomolecular interactions between a drug and its biological target is of paramount importance for the design of novel bioactive compounds. In this paper, we report on the use of molecular dynamics (MD) simulations and machine learning to study the binding mechanism of a transition state analogue (DADMe–immucillin-H) to the purine nucleoside phosphorylase (PNP) enzyme. Microsecond-long MD simulations allow us to observe several binding events, following different dynamical routes and reaching diverse binding configurations. These simulations are used to estimate kinetic and thermodynamic quantities, such as kon and binding free energy, obtaining a good agreement with available experimental data. In addition, we advance a hypothesis for the slow-onset inhibition mechanism of DADMe–immucillin-H against PNP. Combining extensive MD simulations with machine learning algorithms could therefore be a fruitful approach for capturing key aspects of drug–target recognition and binding.

Synthesis and Structure–Activity Relationship (SAR) of 2-Methyl-4-oxo-3-oxetanylcarbamic Acid Esters, a Class of Potent N-Acylethanolamine Acid Amidase (NAAA) Inhibitors

N-Acylethanolamine acid amidase (NAAA) is a lysosomal cysteine hydrolase involved in the degradation of saturated and monounsaturated fatty acid ethanolamides (FAEs), a family of endogenous lipid agonists of peroxisome proliferator-activated receptor-α, which include oleoylethanolamide (OEA) and palmitoylethanolamide (PEA). The β-lactone derivatives (S)-N-(2-oxo-3-oxetanyl)-3-phenylpropionamide (2) and (S)-N-(2-oxo-3-oxetanyl)-biphenyl-4-carboxamide (3) inhibit NAAA, prevent FAE hydrolysis in activated inflammatory cells, and reduce tissue reactions to pro-inflammatory stimuli. Recently, our group disclosed ARN077 (4), a potent NAAA inhibitor that is active in vivo by topical administration in rodent models of hyperalgesia and allodynia. In the present study, we investigated the structure-activity relationship (SAR) of threonine-derived β-lactone analogues of compound 4. The main results of this work were an enhancement of the inhibitory potency of β-lactone carbamate derivatives for NAAA and the identification of (4-phenylphenyl)-methyl-N-[(2S,3R)-2-methyl-4-oxo-oxetan-3-yl]carbamate (14q) as the first single-digit nanomolar inhibitor of intracellular NAAA activity (IC50 = 7 nM on both rat NAAA and human NAAA).

Synthesis, structure-activity, and structure-stability relationships of 2-substituted-N-(4-oxo-3-oxetanyl) N-acylethanolamine acid amidase (NAAA) inhibitors.

N‐Acylethanolamine acid amidase (NAAA) is a cysteine amidase that preferentially hydrolyzes saturated or monounsaturated fatty acid ethanolamides (FAEs), such as palmitoylethanolamide (PEA) and oleoylethanolamide (OEA), which are endogenous agonists of nuclear peroxisome proliferator‐activated receptor‐α (PPAR‐α). Compounds that feature an α‐amino‐β‐lactone ring have been identified as potent and selective NAAA inhibitors and have been shown to exert marked anti‐inflammatory effects that are mediated through FAE‐dependent activation of PPAR‐α. We synthesized and tested a series of racemic, diastereomerically pure β‐substituted α‐amino‐β‐lactones, as either carbamate or amide derivatives, investigating the structure–activity and structure–stability relationships (SAR and SSR) following changes in β‐substituent size, relative stereochemistry at the α‐ and β‐positions, and α‐amino functionality. Substituted carbamate derivatives emerged as more active and stable than amide analogues, with the cis configuration being generally preferred for stability. Increased steric bulk at the β‐position negatively affected NAAA inhibitory potency, while improving both chemical and plasma stability.

Predicting the Reactivity of Nitrile-Carrying Compounds with Cysteine: A Combined Computational and Experimental Study

Here, we report on a mechanistic investigation based on DFT calculations and kinetic measures aimed at determining the energetics related to the cysteine nucleophilic attack on nitrile-carrying compounds. Activation energies were found to correlate well with experimental kinetic measures of reactivity with cysteine in phosphate buffer. The agreement between computations and experiments points to this DFT-based approach as a tool for predicting both nitrile reactivity toward cysteines and the toxicity of nitriles as electrophile agents.

Cyclin-dependent kinases: bridging their structure and function through computations

Cyclin-dependent kinases (CDKs) are one of the most promising target families for drug discovery for several diseases, such as cancer and neurodegenerative disorders. Over the years, structural insights on CDKs have demonstrated high protein plasticity, with several cases where two or more structures of the same protein adopt different conformations. This has generated a great deal of interest in understanding the relationship between CDK structure and function. Here, we highlight how computer simulations have recently contributed in characterizing some key rare and transient events in CDKs, such as the reaction transition state and activation loop movement. Although not yet fully defined, we can now portray the enzymatic mechanism and plasticity of CDKs at high spatial and temporal resolution. These theoretical studies bridge with experiments and highlight structural determinants that could help in designing specific CDK inhibitors.

Synthesis, biological evaluation, and 3D QSAR study of 2-methyl-4-oxo-3-oxetanylcarbamic acid esters as N-acylethanolamine acid amidase (NAAA) inhibitors.

N-(2-Oxo-3-oxetanyl)carbamic acid esters have recently been reported to be noncompetitive inhibitors of the N-acylethanolamine acid amidase (NAAA) potentially useful for the treatment of pain and inflammation. In the present study, we further explored the structure-activity relationships of the carbamic acid ester side chain of 2-methyl-4-oxo-3-oxetanylcarbamic acid ester derivatives. Additional favorable features in the design of potent NAAA inhibitors have been found together with the identification of a single digit nanomolar inhibitor. In addition, we devised a 3D QSAR using the atomic property field method. The model turned out to be able to account for the structural variability and was prospectively validated by designing, synthesizing, and testing novel inhibitors. The fairly good agreement between predictions and experimental potency values points to this 3D QSAR model as the first example of quantitative structure-activity relationships in the field of NAAA inhibitors.

Protein Conformational Plasticity: the “off‐on” Switching Movement in Cdk5

Cyclin‐dependent kinases (CDKs) are mostly known for their role in the cell cycle regulation. The activation mechanism of all CDKs involves the association with a regulatory protein, generally a cyclin, that binds to the kinase unit and stabilizes a catalytically active conformation. Active and inactive conformations of CDKs are characterized by the different spatial localization of two typical elements, namely the activation loop and an □‐helix, whose amino‐acid composition varies throughout the family.