Ribosome structure: binding site of macrolides studied by photoaffinity labeling

The macrolide antibiotics carbomycin A, niddamycin, and tylosin have been radioactively labeled by reducing their aldehyde group at the C-18 position. Dihydro derivatives with specific activities around 2.5 Ci/mmol can be obtained that, although partially affected in their activity, still bind to the ribosomes with high affinity. The presence in the chemical structure of these antibiotics of alpha-beta-unsaturated ketone groups makes them photochemically reactive, and by irradiation above 300 nm, covalent incorporation of the radioactive dihydro derivatives into ribosomes has been achieved. The covalent binding seems to take place at the specific binding sites for macrolides as deduced from binding saturation studies and competition experiments with unmodified drugs. Analysis of the ribosomal components labeled by the drugs indicated that most radioactivity is associated with the proteins L27, L2, and L28 when 50S subunits are labeled, and with L27, L2, L32/33, S9, and S12 in the case of 70S ribosomes. These results agree well with a model of macrolides' mode of action that assumes an interaction of the drug at the peptidyl transferase P site that would block the exit channel for the growing peptide chain.     

A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTA LIGAND , we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density‐95/Dlg/ZO‐1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.     

3D-QSAR and docking studies of aldehyde inhibitors of human cathepsin K

In order to better understand the structural and chemical features of human cathepsin K (CatK), which is an important cysteine protease in the pathogenesis of osteoporosis, the 3D-QSAR (CoMFA) studies were conducted on recently explored aldehyde compounds with known CatK inhibitory activities. The genetic algorithm of GOLD2.2 has been employed to position 59 aldehyde compounds into the active sites of CatK to determine the probable binding conformation. Good correlations between the predicted binding free energies and the experimental inhibitory activities suggested that the identified binding conformations of these potential inhibitors are reliable. The docking results also provided a reliable conformational alignment scheme for 3D-QSAR model. Based on the docking conformations, highly predictive comparative molecular field analysis (CoMFA) was performed with q2 value of 0.723. The predictive ability was validated by some compounds that were not included in the training set. Furthermore, the CoMFA model was mapped back to the binding sites of CatK, to get a better understanding of vital interactions between the aldehyde compounds and the protease. The CoMFA field distributions are in good agreement with the structural characteristics of the binding groove of the CatK, which suggested that the n-Bu in R4 position is the favor group substitute at P1 and moderate groups in R2 group are required on P2 substitute. In addition, 3D-QSAR results also demonstrated that aldehyde is an important pharmacophore because of electrostatic effect. These results, together with the good correlations between the inhibitory activities and the binding free energies predicted by GOLD2.2, demonstrated the power of combining docking/QSAR approach to explore the probable binding conformations of compounds at the active sites of the protein target, and further provided useful information in understanding the structural and chemical features of CatK in designing and finding new potential inhibitors.     

Recognition of the gluconeogenic enzyme,Pck1, via the Gid4 E3 ligase: An in silico perspective

Gluconeogenesis, the reverse process of glycolysis, is a favorable mechanism at conditions of glucose deprivation. Pck1 is a rate‐limiting gluconeogenic enzyme, where its deficiency or mutation contributes to serious clinical situations as neonatal hypoglycemia and liver failure. A recent report confirms that Pck1 is a target for proteasomal degradation through its proline residue at the penultimate position, recognized by Gid4 E3 ligase, but with a lack of informative structural details. In this study, we delineate the localized sequence motif, degron, that specifically interact with Gid4 ligase and unravel the binding mode of Pck1 to the Gid4 ligase by using molecular docking and molecular dynamics. The peptide/protein docking HPEPDOCK web server along with molecular dynamic simulations are applied to demonstrate the binding mode and interactions of a Pck1 wild type (SPSK) and mutant (K4V) with the recently solved structure of Gid4 ligase. Results unveil a distinct binding mode of the mutated peptide compared with the wild type despite having comparable binding affinities to Gid4. Moreover, the four‐residue peptide is found insufficient for Gid4 binding, while the seven‐residue peptide suffices for binding to Gid4. The amino acids S134, K135, and N137 in the loop L1 (between β1 and β2) of the Gid4 are essential for the stabilization of the seven‐residue peptide in the binding site of the ligase. The presence of Val⁴ instead of Lys⁴ smashes the H‐bonds that are formed between Lys⁴ and Gid4 in the wild type peptide, making the peptide prone to bind with the other side of the binding pocket (L4 loop of Gid4). The dynamics of Gid4 L3 loop is affected dramatically once K4V mutant Pck1 peptide is introduced. This opens the door to explore the mutation effects on the binding mode and smooth the path to target protein degradation by design competitive and non‐competitive inhibitors.     

Nanoenzymology of the 20S proteasome: proteasomal actions are controlled by the allosteric transition

The proteasome is a major cytosolic proteolytic assembly, essential for the physiology of eukaryotic cells. Both the architecture and enzymatic properties of the 20S proteasome are relatively well understood. However, despite longstanding interest, the integration of structural and functional properties of the proteasome into a coherent model explaining the mechanism of its enzymatic actions has been difficult. Recently, we used tapping mode atomic force microscopy (AFM) in liquid to demonstrate that the alpha-rings of the proteasome imaged in a top-view position repeatedly switched between their open and closed conformations, apparently to control access to the central channel. Here, we show with AFM that the molecules in a side-view position acquired two stable conformations. The overall shapes of the 20S particles were classified as either barrel-like or cylinder-like. The relative abundance of the two conformers depended on the nature of their interactions with ligands. Similarly to the closed molecules in top view, the barrels predominated in control or inhibited molecules. The cylinders and open molecules prevailed when the proteasome was observed in the presence of peptide substrates. Based on these data, we developed the two-state model of allosteric transitions to explain the dynamics of proteasomal structure. This model helps to better understand the observed properties of the 20S molecule, and sets foundations for further studies of the structural dynamics of the proteasome.     

Analysis of binding modes of ligands to multiple conformations of CYP3A4

Cytochromes P450 (CYPs) are extremely versatile enzymes capable of catalyzing a vast number of compounds, and CYP3A4 is no exception metabolizing approximately half of the currently marketed drugs, besides endogenous compounds. To metabolize such a variety of compounds, CYP3A4 has to be extremely flexible, which makes interaction studies difficult. We employ a multi-conformational docking setup where conformations are generated by several molecular dynamics simulations to analyze the binding modes of various ligands, and the docking is considered successful if the ligand site of catalysis (SOC) is within 6.0 Å of the haem Fe. While docking with the X-ray structure proved unsuccessful with all ligands, the multi-conformational docking achieved successful binding of each ligand to at least one protein conformation. Analysis of the docked solutions highlights residues in the active site cavity that may have an important role in access, binding and stabilization of the ligand.     

Binding of hydrophobic peptides to several non-catalytic sites promotes peptide hydrolysis by all active sites of 20 S proteasomes. Evidence for peptide-induced channel opening in the alpha-rings

The eukaryotic 20 S proteasome contains the following 6 active sites: 2 chymotrypsin-like, 2 trypsin-like, and 2 caspase-like. We previously showed that hydrophobic peptide substrates of the chymotrypsin-like sites allosterically stimulate peptide hydrolysis by the caspase-like sites and their own cleavage. More thorough analysis revealed that these peptides also stimulate peptide hydrolysis by the trypsin-like site. This general activation by hydrophobic peptides occurred even if the chymotrypsin-like sites were occupied by a covalent inhibitor and was highly cooperative, with an average Hill coefficient of 7. Therefore, this stimulation of peptide hydrolysis at all active sites occurs upon binding of hydrophobic peptides to several non-catalytic sites. The stimulation by hydrophobic peptides was not observed in the yeast Delta N alpha 3 mutant 20 S proteasomes, in 20 S-PA26 complexes, or SDS-activated proteasomes and was significantly lower in 26 S proteasomes, all of which appear to have the gated channel in the alpha-rings in an open conformation and hydrolyze peptides at much faster rates than 20 S proteasomes. Also the hydrophobic peptides altered K(m), V(max) of active sites in a similar fashion as PA26 and the Delta N alpha 3 mutation. The activation by hydrophobic peptides was decreased in K(+)-containing buffer, which favors the closed state of the channels. Therefore, hydrophobic peptides stimulate peptide hydrolysis most likely by promoting the opening of the channels in the alpha-rings. During protein breakdown, this peptide-induced channel opening may function to facilitate the release of products from the proteasome.     

Identification of a RON tyrosine kinase receptor binding peptide using phage display technique and computational modeling of its binding mode

RON (Recepteur d’Origine Nantais) tyrosine kinase receptor is a promising target for therapeutic intervention in cancer therapy. The aim of this work was identification of RON-binding peptides using phage display and computational modeling their mode of binding. A 12-mer peptide phage library was utilized to perform biopanning against RON. The RON-binding ability of the selected peptide-displaying phage and their possible binding sites were examined by ELISA. Binding modes and affinities were also predicted by docking and molecular dynamics (MD) simulation. The results of ELISA experiment showed that P6 peptide displaying phage has higher affinity for RON compared to others and its binding site is located out of ligand binding site. Docking and MD simulation results also indicated higher affinity of P6 to RON as well as its exosite-binding feature. Taken together, our data suggest a capacity for P6 peptide (FEHSLYKEMTHL) to be utilized as RON binding agent, and hence be used for various purposes, including design of drug delivery systems for transferring cytotoxic agents to RON-positive cancer cells, interfering with RON signaling, peptidomimetics design, and diagnostic imaging.     

Structural insight into interaction between C20 phenylalanyl derivative of tylosin and ribosomal tunnel

Macrolides are clinically important antibiotics that inhibit protein biosynthesis on ribosomes by binding to ribosomal tunnel. Tylosin belongs to the group of 16-membered macrolides. It is a potent inhibitor of translation whose activity is largely due to reversible covalent binding of its aldehyde group with the base of A2062 in 23S ribosomal RNA. It is known that the conversion of the aldehyde group of tylosin to methyl or carbinol groups dramatically reduces its inhibitory activity. However, earlier we obtained several derivatives of tylosin having comparable activity in spite of the fact that the aldehyde group of tylosin in these compounds was substituted with an amino acid or a peptide residue. Details of the interaction of these compounds with the ribosome that underlies their high inhibitory activity were not known. In the present work, the structure of the complex of tylosin derivative containing in position 20 the residue of ethyl ester of 2-imino(oxy)acetylphenylalanine with the tunnel of the E. coli ribosome was identified by means of molecular dynamics simulations, which could explain high biological activity of this compound.     

HIV-1跨膜蛋白gp41与N-取代吡咯衍生物的结合模式研究

采用分子对接,分子动力学(MD)模拟和分子力学/泊松-波尔兹曼溶剂可有面积方法与分子力学/广义伯恩溶剂可及面积方法(MM-PBSA/MM-GBSA),预测两种N-取代吡咯衍生物与HIV-1 跨膜蛋白gp41疏水口袋的结合模式与作用机理．分子对接采用多种受体构象,并从结果中选取几种可能的结合模式进行MD 模拟,然后通过MM-PBSA计算结合能的方法识别最优的结合模式. MM-PBSA计算结果表明,范德华相互作用是结合的主要驱动力,而极性相互作用决定了配体在结合过程中的取向．进一步的结合能分解显示,配体的羧基与gp41残基Arg579的静电相互作用对结合有重要贡献．上述工作为进一步优化N-取代吡咯衍生物类的HIV-1融合抑制剂建立了良好的理论基础．     

Molecular docking and molecular dynamics simulation approaches for identifying new lead compounds as potential AChE inhibitors

To provide hints for the design of novel acetylcholinesterase (AChE) inhibitors with higher potency and specificity, the binding modes of the (RS, S)-17b and (RS, R)-17b enantiomers on AChE were chosen to investigate by molecular docking and molecular dynamics simulation. The results show that the binding modes of (RS, S)-17b and (RS, R)-17b are clearly different from each other. In particular, the (RS, S)-17b and (RS, R)-17b enantiomers tend to be planar and bend conformations to interact with AChE, respectively. Furthermore, based on the binding mode on AChE and structure modification of (RS, S)-17b, two novel inhibitors (1 and 2) with higher inhibitory activity were designed. Our design strategy suggests that the number of N and O atoms should be increased, the 5, 6-dimethoxy should be transformed into ring and the indanone moiety should be ring-opening, which would result in generating potent and selective AChE inhibitors.     

Correlation between biological activity and binding energy in systems of integrin with cyclic RGD-containing binders: a QM/MM molecular dynamics study

We here report a combined quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) study on the binding interactions between the α(V)β(3) integrin and eight cyclic arginine-glycine-aspartate (RGD) containing peptides. The initial conformation of each peptide within the binding site of the integrin was determined by docking the ligand to the reactive site of the integrin crystal structure with the aid of docking software FRED. The subsequent QM/MM MD simulations of the complex structures show that these eight cyclic RGD-peptides have a generally similar interaction mode with the binding site of the integrin to the cyclo(RGDf-N[M]V) analog found in the crystal structure. Still, there are subtle differences in the interactions of peptide ligands with the integrin, which contribute to the different inhibition activities. The averaged QM/MM protein-ligand interaction energy (IE) is remarkably correlated to the biological activity of the ligand. The IE, as well as a three-variable model which is somewhat interpretable, thus can be used to predict the bioactivity of a new ligand quantitatively, at least within a family of analogs. The present study establishes a helpful protocol for advancing lead compounds to potent inhibitors.     

Tripeptide analogues of MG132 as protease inhibitors

The 26S proteasome and calpain are linked to a number of important human diseases. Here, we report a series of analogues of the prototypical tripeptide aldehyde inhibitor MG132 that show a unique combination of high activity and selectivity for calpains over proteasome. Tripeptide aldehydes (1–3) with an aromatic P3 substituent show enhanced activity and selectivity against ovine calpain 2 relative to chymotrypsin-like activity of proteasome. Docking studies reveal the key contacts between inhibitors and calpain to confirm the importance of the S3 pocket with respect to selectivity between calpains 1 and 2 and the proteasome.     

Structure-Based Mutational Studies of Substrate Inhibition of Betaine Aldehyde Dehydrogenase BetB from Staphylococcus aureus

Inhibition of enzyme activity by high concentrations of substrate and/or cofactor is a general phenomenon demonstrated in many enzymes, including aldehyde dehydrogenases. Here we show that the uncharacterized protein BetB (SA2613) from Staphylococcus aureus is a highly specific betaine aldehyde dehydrogenase, which exhibits substrate inhibition at concentrations of betaine aldehyde as low as 0.15 mM. In contrast, the aldehyde dehydrogenase YdcW from Escherichia coli, which is also active against betaine aldehyde, shows no inhibition by this substrate. Using the crystal structures of BetB and YdcW, we performed a structure-based mutational analysis of BetB and introduced the YdcW residues into the BetB active site. From a total of 32 mutations, those in five residues located in the substrate binding pocket (Val288, Ser290, His448, Tyr450, and Trp456) greatly reduced the substrate inhibition of BetB, whereas the double mutant protein H448F/Y450L demonstrated a complete loss of substrate inhibition. Substrate inhibition was also reduced by mutations of the semiconserved Gly234 (to Ser, Thr, or Ala) located in the BetB NAD⁺ binding site, suggesting some cooperativity between the cofactor and substrate binding sites. Substrate docking analysis of the BetB and YdcW active sites revealed that the wild-type BetB can bind betaine aldehyde in both productive and nonproductive conformations, whereas only the productive binding mode can be modeled in the active sites of YdcW and the BetB mutant proteins with reduced substrate inhibition. Thus, our results suggest that the molecular mechanism of substrate inhibition of BetB is associated with the nonproductive binding of betaine aldehyde.     

Caspase-3 binds diverse P4 residues in peptides as revealed by crystallography and structural modeling

Caspase-3 recognition of various P4 residues in its numerous protein substrates was investigated by crystallography, kinetics, and calculations on model complexes. Asp is the most frequent P4 residue in peptide substrates, although a wide variety of P4 residues are found in the cellular proteins cleaved by caspase-3. The binding of peptidic inhibitors with hydrophobic P4 residues, or no P4 residue, is illustrated by crystal structures of caspase-3 complexes with Ac-IEPD-Cho, Ac-WEHD-Cho, Ac-YVAD-Cho, and Boc-D(OMe)-Fmk at resolutions of 1.9–2.6 Å. The P4 residues formed favorable hydrophobic interactions in two separate hydrophobic regions of the binding site. The side chains of P4 Ile and Tyr form hydrophobic interactions with caspase-3 residues Trp206 and Trp214 within a non-polar pocket of the S4 subsite, while P4 Trp interacts with Phe250 and Phe252 that can also form the S5 subsite. These interactions of hydrophobic P4 residues are distinct from those for polar P4 Asp, which indicates the adaptability of caspase-3 for binding diverse P4 residues. The predicted trends in peptide binding from molecular models had high correlation with experimental values for peptide inhibitors. Analysis of structural models for the binding of 20 different amino acids at P4 in the aldehyde peptide Ac-XEVD-Cho suggested that the majority of hydrophilic P4 residues interact with Phe250, while hydrophobic residues interact with Trp206, Phe250, and Trp214. Overall, the S4 pocket of caspase-3 exhibits flexible adaptation for different residues and the new structures and models, especially for hydrophobic P4 residues, will be helpful for the design of caspase-3 based drugs.     

Conformer and pharmacophore based identification of peptidomimetic inhibitors of chikungunya virus nsP2 protease

Chikungunya virus nsP2 replication protein is a cysteine protease, which cleaves the nonstructural nsP1234 polyprotein into functional replication components. The cleavage and processing of nsP1234 by nsP2 protease is essential for the replication and proliferation of the virus. Thus, ChikV nsP2 protease is a promising target for antiviral drug discovery. In this study, the crystal structure of the C-terminal domain of ChikV nsP2 protease (PDB ID: 4ZTB) was used for structure based identification and rational designing of peptidomimetic inhibitors against nsP2 protease. The interactions of the junction residues of nsP3/4 polyprotein in the active site of nsP2 protease have been mimicked to identify and design potential inhibitory molecules. Molecular docking of the nsP3/4 junction peptide in the active site of ChikV nsP2 protease provided the structural insight of the probable binding mode of nsP3/4 peptide and pigeonholed the molecular interactions critical for the substrate binding. Further, the shape and pharmacophoric properties of the viral nsP3/4 substrate peptide were taken into consideration and the mimetic molecules were identified and designed. The designed mimetic compounds were then analyzed by docking and their binding affinity was assessed by molecular dynamics simulations.     

How to choose relevant multiple receptor conformations for virtual screening: a test case of Cdk2 and normal mode analysis

Better treatment of protein flexibility is essential in structure-based drug design projects such as virtual screening and protein-ligand docking. Diversity in ligand-binding mechanisms and receptor conformational changes makes it difficult to treat dynamic features of the receptor during the docking simulation. Thus, the use of pregenerated multiple receptor conformations is applied today in virtual screening studies. However, generation of a small relevant set of receptor conformations remains challenging. To address this problem, we propose a new protocol for the generation of multiple receptor conformations via normal mode analysis and for the selection of several receptor conformations suitable for docking/virtual screening. We validated this protocol on cyclin-dependent kinase 2, which possesses a binding site located at the interface between two subdomains and is known to undergo significant conformational changes in the active site region upon ligand binding. We believe that the suggested rules for the choice of suitable receptor conformations can be applied to other targets when dealing with in silico screening on flexible receptors.     

Molecular determinants of thapsigargin binding by SERCA Ca2+-ATPase: a computational docking study

Thapsigargin (TG) is a potent and commonly used inhibitor of the ion transport activity of sarco/endoplasmic reticulum Ca2+-ATPases (SERCA). Based on the recently published crystal structures of rabbit muscle SERCA1a in the Ca2+/E1 (E1) and TG/E2 (E2) conformations, we performed computational docking studies to characterize the molecular interactions that govern binding of TG and TG-analogs by the enzyme. Using the program GOLD (genetic optimization for ligand docking) in combination with the scoring function ChemScore, TG was docked into the binding site of the E1 and E2 conformations of SERCA1a. The docking results revealed a consensus ligand-binding mode consistent with the crystal structure and showed that hydrophobic interactions are the primary driving force of TG binding by SERCA. Moreover, it was shown that the conformational changes accompanying the E2 to E1 transition in the enzyme likely displace TG from its favored orientation in the binding site, thereby substantially reducing its binding affinity. This finding illustrates on the molecular level how TG may exert its inhibitory effect in binding tightly to the E2 form and preventing it from converting into its E1 form, a requirement for catalytic function. We also docked 9 TG analogs into the E2 conformation of the enzyme. Eight of the analogs adopted a binding mode very similar to that of TG, whereas one compound preferred a different orientation in the binding site. Analysis of the predicted binding affinities showed a good correlation with the experimentally observed inhibitory potencies of the analogs. Docking was also performed with several modeled mutants of SERCA1a, whose phenylalanine residue in position 256 (Phe256) had been modified. The experimentally observed declines in TG sensitivity in most of the Phe256 mutants was qualitatively accounted for and appears, at least in part, be due to a slightly altered TG-binding mode.     

Modeling of Arylamide Helix Mimetics in the p53 Peptide Binding Site of hDM2 Suggests Parallel and Anti-Parallel Conformations Are Both Stable

The design of novel α-helix mimetic inhibitors of protein-protein interactions is of interest to pharmaceuticals and chemical genetics researchers as these inhibitors provide a chemical scaffold presenting side chains in the same geometry as an α-helix. This conformational arrangement allows the design of high affinity inhibitors mimicking known peptide sequences binding specific protein substrates. We show that GAFF and AutoDock potentials do not properly capture the conformational preferences of α-helix mimetics based on arylamide oligomers and identify alternate parameters matching solution NMR data and suitable for molecular dynamics simulation of arylamide compounds. Results from both docking and molecular dynamics simulations are consistent with the arylamides binding in the p53 peptide binding pocket. Simulations of arylamides in the p53 binding pocket of hDM2 are consistent with binding, exhibiting similar structural dynamics in the pocket as simulations of known hDM2 binders Nutlin-2 and a benzodiazepinedione compound. Arylamide conformations converge towards the same region of the binding pocket on the 20 ns time scale, and most, though not all dihedrals in the binding pocket are well sampled on this timescale. We show that there are two putative classes of binding modes for arylamide compounds supported equally by the modeling evidence. In the first, the arylamide compound lies parallel to the observed p53 helix. In the second class, not previously identified or proposed, the arylamide compound lies anti-parallel to the p53 helix.