A new intermolecular force field for the liquid-crystal-forming molecule 5CB (4-cyano, 4‘n-pentyl biphenyl) has been derived from two-body interaction energies, obtained by the fragmentation reconstruction method (FRM). The accuracy of this purely quantum mechanical approach has been verified by comparing FRM and directly ab initio computed interaction energies, obtaining a very satisfactory agreement with a maximum absolute error lower than 0.4 kcal/mol. The comparison was performed for a large number of internal geometries and dimer arrangements of two fragments of 5CB, namely, 4-cyano biphenyl (0CB) and n-pentyl benzene (5B). The potential energy surface of the 5CB dimer has then been used to parametrize a many-site model interaction potential suitable for computer simulations.

Due to its central role in immunosuppression and cell proliferation and due to its specific peptidyl-prolyl-isomerase (PPI) function, the FKBP protein family is at the crossroad of several important metabolic pathways. Members of this family, and notably FK506 binding protein (FKBP12), are thought to be involved in neurodegenerative diseases such as Alzheimer disease, Parkinson disease, multiple sclerosis, amyotrophic lateral sclerosis, as well as in proliferation disorders and cancer. Using an interdisciplinary approach based on computational, synthetic, and experimental techniques, we show that the best potential binders for FKBP proteins optimally expose the two contiguous carbonyl oxygen in the proline-mimetic chain for FKBP docking and are characterized by the abundance of rigid quasi-cyclic structures stabilized in aqueous solution by intraligand hydrophobic interactions mimicking the macrolide structure of the natural FKBP binders FK506 and Rapamycin. These peculiar structural and chemical-physical features define at the same time an ElteX compound and the minimal pharmacore in the FKBP family, shedding new light on the isomerization mechanism of the PPI domain. On the basis of the above hypothesis, we have successfully designed and synthesized several nanomolar ElteX FKBP12 ligands. Among these, ElteN378 is a new low atomic weight ligand with affinity comparable to that of the macrolide Rapamycin.

Understanding binding mechanisms between enzymes and potential inhibitors and quantifying protein-ligand affinities in terms of binding free energy is of primary importance in drug design studies. In this respect, several approaches based on molecular dynamics simulations, often combined with docking techniques, have been exploited to investigate the physicochemical properties of complexes of pharmaceutical interest. Even if the geometric properties of a modeled protein-ligand complex can be well predicted by computational methods, it is still challenging to rank with chemical accuracy a series of ligand analogues in a consistent way. In this article, we face this issue calculating relative binding free energies of a focal adhesion kinase, an important target for the development of anticancer drugs, with pyrrolopyrimidine-based ligands having different inhibitory power. To this aim, we employ steered molecular dynamics simulations combined with nonequilibrium work theorems for free energy calculations. This technique proves very powerful when a series of ligand analogues is considered, allowing one to tackle estimation of protein-ligand relative binding free energies in a reasonable time. In our cases, the calculated binding affinities are comparable with those recovered from experiments by exploiting the Michaelis-Menten mechanism with a competitive inhibitor.

Supporting Information. Compounds, synthesis of Elte421, computational details of REM simulations, chemical equilibria, fluorescence data and experimental details, NMR spectra, and references for Supporting Information. This material is available free of charge via the Internet at http://pubs.acs.org.

The conformational landscape of three FK506-related drugs with disparate inhibition constants is determined in bulk solution using a replica exchange simulation method with solute torsional tempering. Energetic fitness of important drug conformations with respect to the FKBP12 protein is evaluated by molecular mechanics. Results show that the experimental affinity toward peptidyl-prolyl cis-trans isomerase protein (FKBP12) of the analyzed ligands appears to be positively correlated to the observed population of specific chair structures of the drug piperidinic ring in bulk solution. This observation is rationalized on the basis that such structures, stabilized by stereospecific intramolecular hydrophobic interactions, allows the formation of a pair of protein-ligand hydrogen bonds upon binding.

We present a new computational strategy for calculating the absolute binding free energy for protein ligand association in the context of atomistic simulation in explicit solvent. The method is based on an appropriate definition of a solute tempering scheme enforced via Hamilton replica exchange method (HREM). The definition of "solute" includes both the ligand and the active site, with the remainder of the systems defined as "solvent". The hydrophilicity of the solute and the solute torsional plus nonbonded intrasolute interactions are increased and decreased, respectively, along the replica progression, thus favoring the extrusion of the drug form the active site in the scaled states of the generalized ensemble. The proposed technique, named "Energy Driven Undocking" (EDU-HREM), completely bypasses the need for defining and/or identifying the relevant reaction coordinates in a ligand receptor interactions and allows the calculation of the absolute binding free energy in one single generalized simulation of the drug-receptor system. The methodology is applied, with encouraging results, to the calculation of the absolute binding free energy of some FK506-related ligands of the peptidyl prolyl cis-trans isomerase protein (FKBP12) with known dissociation constants. Aspects of the binding/inhibition mechanism in FKBP12 are also analyzed and discussed.