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.

Dissecting the sequence determinants for dephosphorylation by the catalytic subunits of phosphatases PP1 and PP2A.

The phosphatases PP1 and PP2A are responsible for the majority of dephosphorylation reactions on phosphoserine (pSer) and phosphothreonine (pThr), and are involved in virtually all cellular processes and numerous diseases. The catalytic subunits exist in cells in form of holoenzymes, which impart substrate specificity. The contribution of the catalytic subunits to the recognition of substrates is unclear. By developing a phosphopeptide library approach and a phosphoproteomic assay, we demonstrate that the specificity of PP1 and PP2A holoenzymes towards pThr and of PP1 for basic motifs adjacent to the phosphorylation site are due to intrinsic properties of the catalytic subunits. Thus, we dissect this amino acid specificity of the catalytic subunits from the contribution of regulatory proteins. Furthermore, our approach enables discovering a role for PP1 as regulator of the GRB-associated-binding protein 2 (GAB2)/14-3-3 complex. Beyond this, we expect that this approach is broadly applicable to detect enzyme-substrate recognition preferences.

Protein conformational transitions: the closure mechanism of a kinase explored by atomistic simulations.

Kinase large-scale conformational rearrangement is an issue of enormous biological and pharmacological relevance. Atomistic simulations able to capture the dynamics and the energetics of kinase large-scale motions are still in their infancy. Here, we present a computational study in which the atomistic dynamics of the "open-to-closed" movement of the cyclin-dependent kinase 5 (CDK5) have been simulated. Simulations were carried out using a new sampling method that is able to find the lowest free-energy channel between an initial state and a final state. This large-scale movement has a two-step mechanism: first, the alphaC-helix rotates by approximately 45 degrees , allowing the interaction between Glu51 and Arg149; then the CDK5 activation loop refolds to assume the closed conformation. We have also estimated the free-energy profile associated with the global motion and identified a CDK5 intermediate, which could be exploited for drug-design purposes. Our new sampling method turned out to be well-suited for investigating at an atomistic level the energetics and dynamics of kinase large-scale conformational motions.

Mutational Analysis of a Conserved Glutamate Reveals Unique Mechanistic and Structural Features of the Phosphatase PRL-3.

Phosphatase of regenerating liver (PRL)-3 (PTP4A3) has gained much attention in cancer research due to its involvement in tumor promoting and metastatic processes. It belongs to the protein tyrosine phosphatase (PTP) superfamily and is thought to follow the catalytic mechanism shared by this family, which aside from the conserved active-site amino acids includes a conserved glutamic acid residue that is usually required for the integrity of the active site in PTPs. We noted that in structures of PRL-3, PRL-1, and PTEN these residues do not clearly align and therefore we sought to investigate if the glutamic acid residue fulfills its usual function in these proteins. Although this residue was essential for PTEN's catalytic activity, it was nonessential for PRL-1 and PRL-3. Surprisingly, the mutation E50R increased PRL-3 activity against all tested in vitro substrates and also enhanced PRL-3-promoted cell adhesion and migration. We show that the introduction of Arg50 leads to an enhancement of substrate turnover for both PRL-3 and, to a lesser extent, PRL-1, and that the stronger gain in activity correlates with a higher structural flexibility of PRL-3, likely allowing for conformational adaptation during catalysis. Thus, in contrast to its crucial functions in other PTPs, this conserved glutamic acid can be replaced in PRL-3 without impairing the structural integrity. The variant with enhanced activity might serve as a tool to study PRL-3 in the future.

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.

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.

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.

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.

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).

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.