Assessing protein-ligand binding modes with computational tools: the case of PDE4B.

In a first step in the discovery of novel potent inhibitor structures for the PDE4B family with limited side effects, we present a protocol to rank newly designed molecules through the estimation of their IC[Formula: see text] values. Our protocol is based on reproducing the linear relationship between the logarithm of experimental IC[Formula: see text] values [[Formula: see text](IC[Formula: see text])] and their calculated binding free energies ([Formula: see text]). From 13 known PDE4B inhibitors, we show here that (1) binding free energies obtained after a docking process by AutoDock are not accurate enough to reproduce this linear relationship; (2) MM-GB/SA post-processing of molecular dynamics (MD) trajectories of the top ranked AutoDock pose improves the linear relationship; (3) by taking into account all representative structures obtained by AutoDock and by averaging MM-GB/SA computations on a series of 40 independent MD trajectories, a linear relationship between [Formula: see text](IC[Formula: see text]) and the lowest [Formula: see text] is achieved with [Formula: see text].

Investigation of allosteric coupling in human ß2-adrenergic receptor in the presence of intracellular loop 3.

BACKGROUND: This study investigates the allosteric coupling that exists between the intra- and extracellular parts of human ß2-adrenergic receptor (ß2-AR), in the presence of the intracellular loop 3 (ICL3), which is missing in all crystallographic experiments and most of the simulation studies reported so far. Our recent 1 µs long MD run has revealed a transition to the so-called very inactive state of the receptor, in which ICL3 packed under the G protein's binding cavity and completely blocked its accessibility to G protein. Simultaneously, an outward tilt of transmembrane helix 5 (TM5) caused an expansion of the extracellular ligand-binding site. In the current study, we performed independent runs with a total duration of 4 µs to further investigate the very inactive state with packed ICL3 and the allosteric coupling event (three unrestrained runs and five runs with bond restraints at the ligand-binding site). RESULTS: In all three independent unrestrained runs (each 500 ns long), ICL3 preserved its initially packed/closed conformation within the studied time frame, suggesting an inhibition of the receptor's activity. Specific bond restraints were later imposed between some key residues at the ligand-binding site, which have been experimentally determined to interact with the ligand. Restraining the binding site region to an open state facilitated ICL3 closure, whereas a relatively constrained/closed binding site hindered ICL3 packing. However, the reverse operation, i.e. opening of the packed ICL3, could not be realized by restraining the binding site region to a closed state. Thus, any attempt failed to free the ICL3 from its locked state due to the presence of persistent hydrogen bonds. CONCLUSIONS: Overall, our simulations indicated that starting with very inactive states, the receptor stayed almost irreversibly inhibited, which in turn decreased the overall mobility of the receptor. Bond restraints which represented the geometric restrictions caused by ligands of various sizes when bound at the ligand-binding site, induced the expected conformational changes in TM5, TM6 and consequently, ICL3. Still, once ICL3 was packed, the allosteric coupling became ineffective due to strong hydrogen bonds connecting ICL3 to the core of the receptor.

Identification of Alternative Allosteric Sites in Glycolytic Enzymes for Potential Use as Species-Specific Drug Targets.

Three allosteric glycolytic enzymes, phosphofructokinase, glyceraldehyde-3 phosphate dehydrogenase and pyruvate kinase, associated with bacterial, parasitic and human species, were explored to identify potential allosteric sites that would be used as prime targets for species-specific drug design purposes using a newly developed approach which incorporates solvent mapping, elastic network modeling, sequence and structural alignments. The majority of binding sites detected by solvent mapping overlapped with the interface regions connecting the subunits, thus appeared as promising target sites for allosteric regulation. Each binding site was then evaluated by its ability to alter the global dynamics of the receptor defined by the percentage change in the frequencies of the lowest-frequency modes most significantly and as anticipated, the most effective ones were detected in the vicinity of the well-reported catalytic and allosteric sites. Furthermore, some of our proposed regions intersected with experimentally resolved sites which are known to be critical for activity regulation, which further validated our approach. Despite the high degree of structural conservation encountered between bacterial/parasitic and human glycolytic enzymes, the majority of the newly presented allosteric sites exhibited a low degree of sequence conservation which further increased their likelihood to be used as species-specific target regions for drug design studies.

A docking study using atomistic conformers generated via elastic network model for cyclosporin A/cyclophilin A complex.

Anisotropic network model is used to generate a set of distinct conformations for cylophilin A (CypA). The native structure is deformed to different extents along each of the lowest-frequency modes (first 7 modes) both in negative and positive directions. Each node of the elastic network represents either a single atom in the high-resolution model or a single residue in the low-resolution model. Realistic conformations with energies close to or lower than the crystal structure and with satisfactory internal geometry are recovered by energy minimization using implicit solvation model. These conformations are then used for ensemble docking to the ligand cyclosporin A for both a further test of accuracy of generated conformers and exploration of different binding modes. Higher number of correctly docked ligands are obtained for conformations with low deformation factors as a result of lower root mean square distances with respect to crystal structure. Yet, surprisingly, the lowest binding energy is obtained for one of the highly deformed conformations as a result of its special contact with arginine side chain oriented towards binding site. Considering the fact that the cyclic ligand's backbone and protein's side chains are held rigid during docking, the conformers generated by high- and low-resolution elastic network models are almost equally successful in providing the correct binding mode. The shape of the binding pocket that incorporates crucial interaction sites for hydrogen bond formation is found to be another important determining factor for the success of the dock. Also, the small backbone variations of a few Angstroms in magnitude at the loop regions surrounding the binding pocket can cause amino acids' side chains to be displaced by magnitudes of up to 10 A and therefore have a strong influence on the efficiency of the conformational search during docking.

Intrinsic Dynamics and Causality in Correlated Motions Unraveled in Two Distinct Inactive States of Human ß2-Adrenergic Receptor.

The alternative inactive state of the human ß2-adrenergic receptor originally exposed in molecular dynamics simulations was investigated using various analysis tools to evaluate causality between correlated residue-pair fluctuations and suggest allosteric communication pathways. A major conformational shift observed in the third intracellular loop (ICL3) displayed a novel inactive state, featuring an inaccessible G protein binding site blocked by ICL3 and an expanded orthosteric ligand binding site. Residue-based mean-square fluctuation and stiffness calculations revealed a significant mobility decrease in ICL3, which induced a mobility increase in the remaining loop regions. This indicates conformational entropy loss in one mobile region being compensated by residual intermolecular motions in other mobile regions. Moreover, the extent of significantly correlated motions decreased, and correlations that once existed between transmembrane helices shifted toward regions with increased mobility. Conditional time-delayed cross-correlation analysis identified distinct driver-follower relationship profiles. Prior to its packing, freely moving ICL3 was markedly driven by transmembrane helix-8 whereas once packed, ICL3 controlled future fluctuations of nearby helices. Moreover, two transmembrane helices, (H5 and H6), started to control future fluctuations of a remote site, the extracellular loop, ECL2. This clearly suggests that allosteric coupling between extra- and intracellular parts intensified, in agreement with the receptor's well recognized feature, which is the inverse proportionality between activity and the degree of coupling.

Structural analysis of peptide fragments following the hydrolysis of bovine serum albumin by trypsin and chymotrypsin.

Peptide bond hydrolysis of bovine serum albumin (BSA) by chymotrypsin and trypsin was investigated by employing time-resolved fluorescence spectroscopy. As a fluorescent cross-linking reagent, N-(1-pyrenyl) maleimide (PM) was attached to BSA, through all free amine groups of arginine, lysine, and/or single free thiol (Cys34). Time-resolved fluorescence spectroscopy was used to monitor fluorescence decays analyzed by exponential series method to obtain the changes in lifetime distributions. After the exposure of synthesized protein substrate PM-BSA to chymotrypsin and trypsin, it is observed that each protease produced a distinct change in the lifetime distribution profile, which was attributed to distinct chemical environments created by short peptide fragments in each hydrolysate. The persistence of excimer emission at longer lifetime regions for chymotrypsin, as opposed to trypsin, suggested the presence of small-scale hydrophobic clusters that might prevent some excimers from being completely quenched. It is most likely that the formation of these clusters is due to hydrophobic end groups of peptide fragments in chymotrypsin hydrolysate. A similar hydrophobic shield was not suggested for trypsin hydrolysis, as the end groups of peptide fragments would be either arginine or lysine. Overall, in case the target protein's 3D structure is known, the structural analysis of possible excimer formation presented here can be used as a tool to explain the differences in activity between two proteases, i.e. the peak's intensity and location in the profile. Furthermore, this structural evaluation might be helpful in obtaining the optimum experimental conditions in order to generate the highest amount of PM-BSA complexes.

Molecular Docking Study Based on Pharmacophore Modeling for Novel PhosphodiesteraseIV Inhibitors.

In this study, pharmacophore modelling was carried out for novel PhosphodiesteraseIV (PDEIV) inhibitors. A pharmacophore-based virtual screening, which resulted in 1959 hit compounds was performed with six chemical databases. The pharmacophore screening was proven to be successful in discriminating active and inactive inhibitors using a set of compounds with known activity obtained from ChEMBL database. Furthermore, the Lipinski's rule of five was applied for physicochemical filtering of the hit molecules and this yielded 1840 compounds. Three docking software tools, AutoDock 4.0, AutoDock Vina, and Gold v5.1 were used for the docking process. All 1840 compounds and the known selective inhibitor, rolipram, were docked into the active site of the target protein. A total of 234 compounds with all three scoring values higher than those of rolipram were determined with the three docking tools. The interaction maps of 14 potent inhibitors complexed with PDEIV B and D isoforms have been analyzed and seven key residues (Asn 395, Gln 443, Tyr 233, Ile 410, Phe 446, Asp 392, Thr 407) were found to interact with more than 80 % of the potent inhibitors. For each one of the 234 hit compounds, using the bound conformation with the highest AutoDock score, the interacting residues were determined. 117 out of 234 compounds are found to interact with at least five of the seven key residues and these were selected for further evaluation. The conformation with the highest AutoDock score for each 117 compounds were rescored using the DSX scoring function. This yielded a total of 101 compounds with better score values than the natural ligand rolipram. For ADME/TOX calculations, the FAF-Drugs2 server was used and 32 out of 101 compounds were found to be non-toxic.