Conformational flexibility and receptor interaction

Dihedral Angle-Based Sampling of Natural Product Polyketide Conformations: Application to Permeability Prediction

Macrocycles pose challenges for computer-aided drug design due to their conformational complexity. One fundamental challenge is identifying all low-energy conformations of the macrocyclic ring, which is important for modeling target binding, passive membrane permeation, and other conformation-dependent properties. Macrocyclic polyketides are medically and biologically important natural products characterized by structural and functional diversity. Advances in synthetic biology and semisynthetic methods may enable creation of an even more diverse set of non-natural product polyketides for drug discovery and other applications. However, the conformational sampling of these flexible compounds remains demanding. We developed and optimized a dihedral angle-based macrocycle conformational sampling method for macrocycles of arbitrary structure, and here we apply it to diverse polyketide natural products. First, we evaluated its performance using a data set of 37 polyketides with available crystal structures, with 9-22 rotatable bonds in the macrocyclic ring. Our optimized protocol was able to reproduce the crystal structure of polyketides' aglycone backbone within 0.50 Å RMSD for 31 out of 37 polyketides. Consistent with prior structural studies, our analysis suggests that polyketides tend to have multiple distinct low-energy structures, including the bioactive (target-bound) conformation as well as others of unknown significance. For this reason, we also introduce a strategy to improve both efficiency and accuracy of the conformational search by utilizing torsional restraints derived from NMR vicinal proton couplings to restrict the conformational search. Finally, as a first application of the method, we made blinded predictions of the passive membrane permeability of a diverse set of polyketides, based on their predicted structures in low- and high-dielectric media.

Antibody recognition and conformational flexibility of a plaque-specific beta-amyloid epitope modulated by non-native peptide flanking regions

Here we report on the synthesis, antibody binding, and QSAR studies of a series of linear and cyclic peptides containing a beta-amyloid plaque-specific epitope (Abeta(4-10); FRHDSGY). In these constructs, two or three alpha- l-Ala, alpha- d-Ala, or beta-Ala residues were introduced at both N- and C-termini of the epitope as non-native flanking sequences. Cyclization of the linear Abeta(4-10) epitope peptide resulted in reduced antibody binding. However, the antibody binding could be fully compensated by insertion of alanine flanks into the corresponding cyclic peptides. These results indicate that the modification of a beta-amyloid plaque-specific epitope by combination of cyclization and flanking sequences could generate highly antigenic peptides compared to the native sequence. A novel 3D QSAR method, which explicitly handles conformational flexibility, was developed for the case of such molecular libraries. This method led to the prediction of the binding conformation for the common FRHDSGY sequence.

Comparative Analysis of Protein-Bound Ligand Conformations with Respect to Catalyst's Conformational Space Subsampling Algorithms

We examined the quality of Catalyst's conformational model generation algorithm via a large scale study based on the crystal structures of a sample of 510 pharmaceutically relevant protein-ligand complexes extracted from the Protein Data Bank (PDB). Our results show that the tested algorithms implemented within Catalyst are able to produce high quality conformers, which in most of the cases are well suited for in silico drug research. Catalyst-specific settings were analyzed, such as the method used for the conformational model generation (FAST vs BEST) and the maximum number of generated conformers. By setting these options for higher fitting quality, the average RMS values describing the similarity of experimental and simulated conformers were improved from an RMS of 1.06 with max. 50 FAST generated conformers to an RMS of 0.93 with max. 255 BEST generated conformers, which represents an improvement by 12%. Each method provides best fitting conformers with an RMS value<1.50 in more than 80% of all cases. We analyzed the computing time/quality ratio of various conformational model generation settings and examined ligands in high energy conformations. Furthermore, properties of the same ligands in various proteins were investigated, and the fitting qualities of experimental conformations from the PDB and the Cambridge Structural Database (CSD) were compared. One of the most important conclusions of former studies, the fact that bioactive conformers often have energy high above that of global minima, was confirmed.

Successful bridging from a peptide to a non peptide antagonist at the human tachykinin NK-2 receptor

Non peptide products have been found to show nanomolar binding and functional affinities at the human tachykinin NK-2 receptor. The new antagonists do not possess stereogenic centers and their thermal behaviour in solution is featured by a peculiar set of conformational stereoisomers. A macroscopic viewpoint is preferentially adopted to rationalize the obtained results.

Effect of protonation of the N-acetyl neuraminic acid residue of sialyl Lewis(X): a molecular orbital study with insights into its binding properties toward the carbohydrate recognition domain of E-selectin

Semiempirical molecular orbital (MO) calculations with an implicit treatment of the water environment were employed in order to assess whether the sialyl Lewis(X) (sLe(X)) tetrasaccharide binds to E-selectin in the anionic or neutral (i.e., protonated) state. The analysis of the frontier molecular orbitals, electrostatic potential surfaces, and conformational behavior of the sugar indicates that its neutral form possesses the necessary characteristics for binding. In particular, the LUMO level of the neutral sLe(X) molecule is localized both on the carboxylic group of the N-acetyl neuraminic acid (NeuNAc) residue and on the nearby glycosidic linkage. These two moieties interact with the Arg97 residue of E-selectin, as revealed by a recent crystal structure analysis of the E-selectin/sLe(X) complex. The energetics of this specific interaction was investigated with the aid of ab initio Hartree-Fock MO calculations, which resulted in a BSSE-corrected binding energy of 16.63 kcal/mol. Our observations could open up new perspectives in the design of sLe(X) mimics.

Classical QSAR and comparative molecular field analyses of the host-guest interaction of organic molecules with cyclodextrins

The application of classical QSAR and molecular modeling analysis using Comparative Molecular Field Analysis (CoMFA) to the complexation of some natural and modified cyclodextrins (CDs) with guest molecules was examined. For 1:1 complexation systems between natural beta-CD, modified alpha-, beta-, and gamma-CD that bear one p-(dimethylamino)benzoyl (DMAB) moiety (DMAB-alpha-, beta-, and gamma-CDs) and guest molecules of widely varying chemical structures and properties, the binding constants of the complexes were successfully fitted using multiple linear regression (MLR) with hydrophobic descriptor log P (the partition coefficient between 1-octanol and water phases) and molecular connectivity indices. A non-linear dependency of binding constants on the zero-th and/or first order molecular connectivity index as a measure of size becomes apparent. The modeling performance of the CoMFA models with steric/electrostatic fields to DMAB-alpha- and beta-CD systems was comparable to those of MLR models. However, statistically significant CoMFA models for gamma-CD systems which have higher conformational flexibility of the ring could not be obtained. The CoMFA models obtained for DMAB-alpha- and beta-CD systems showed that the predominant effects were steric for the DMAB-alpha-CD system and electrostatic for the DMAB-beta-CD system, respectively.