A major problem in virtual screening concerns the accuracy of the binding free energy between a target protein and a putative ligand. Here we report an example supporting the outperformance of the AutoDock scoring function in virtual screening in comparison to the other popular docking programs. The original AutoDock program is in itself inefficient to be used in virtual screening because the grids of interaction energy have to be calculated for each putative ligand in chemical database. However, the automation of the AutoDock program with the potential grids defined in common for all putative ligands leads to more than twofold increase in the speed of virtual database screening. The utility of the automated AutoDock in virtual screening is further demonstrated by identifying the actual inhibitors of various target enzymes in chemical databases with accuracy higher than the other docking tools including DOCK and FlexX. These results exemplify the usefulness of the automated AutoDock as a new promising tool in structure-based virtual screening.

Memapsin 2 (BACE) is an aspartyl protease known as beta-secretase that acts on the production of the beta-amyloid peptide in the human brain, a key event in the pathogenesis of Alzheimer's disease. Although it is expected that the net charge of the catalytic Asp diad would be -1 as in other kinds of aspartyl proteases, the exact protonation states of Asp32 and Asp228 have not been known without ambiguity. Two independent molecular dynamics (MD) simulations of BACE in complex with the potent inhibitor OM99-2 are carried out to determine the preferred protonation state of the Asp diad in the context that is consistent with the previous X-ray crystal structure. The results show that a strong hydrogen bond between the inhibitor hydroxyl group and Asp228 can be maintained only when Asp32 is neutral and Asp228 is ionized. The preference of this protonation state is further supported from the energetic and structural features found in the docking experiment of a novel potent inhibitor with the BACE active site. Thus, both MD and docking studies suggest that the role of hydrogen bond acceptor for the hydroxyl and piperazine groups of the inhibitors should be played by Asp228 instead of Asp32. This may be a key piece of information for the structure-based design/discovery of new inhibitor drugs.

No abstract

A peptide-cleaving catalyst selective for peptide deformylase (PDF) was obtained from a library containing about 15 000 catalyst candidates. The catalyst cleaved the polypeptide backbone of PDF at Gln(152)-Arg(153). Docking simulations suggested multiple modes of interactions in the complex formed between the catalyst and PDF.

Cytochrome P450 (CYP) 3A4 is responsible for the oxidative degradation of more than 50% of clinically used drugs. By means of molecular dynamics simulations with the newly developed force field parameters for the heme-thiolate group and its dioxygen adduct, we examine the differences in structural and dynamic properties between CYP3A4 in the resting form and its complexes with the substrate progesterone and the inhibitor metyrapone. The results indicate that the broad substrate specificity of CYP3A4 stems from the malleability of a loop (residues 211-218) that resides in the vicinity of the channel connecting the active site and bulk solvent. However, the high-amplitude motion of the flexible loop is found to be damped out upon binding of the inhibitor or the substrate in the active site. In the resting form of CYP3A4, a structural water molecule is bound to the sixth coordination position of the heme iron, stabilizing the octahedral coordination geometry. In addition to the direct coordination of metyrapone to the heme iron, the hydrogen bond interaction between the inhibitor carbonyl group and the side chain of Ser119 also contributes significantly to stabilizing the CYP3A4-metyrapone complex. On the other hand, progesterone is stabilized in the active site by the formation of two hydrogen bonds with Ser119 and Arg106, as well as by the van der Waals interactions with the heme and hydrophobic residues. The structural and dynamic features of the CYP3A4-progesterone complex indicate that the oxidative degradation of progesterone occurs through hydroxylation at the C16 position by the reactive oxygen coordinated to the heme iron.

Human papillomaviruses (HPVs) are causative agents of various diseases associated with cellular hyperproliferation, including cervical cancer, one of the most prevalent tumors in women. E7 is one of the two HPV-encoded oncoproteins and directs recruitment and subsequent degradation of tumor-suppressive proteins such as retinoblastoma protein (pRb) via its LxCxE motif. E7 also triggers tumorigenesis in a pRb-independent pathway through its C-terminal domain, which has yet been largely undetermined, with a lack of structural information in a complex form with a host protein. Herein, we present the crystal structure of the E7 C-terminal domain of HPV18 belonging to the high-risk HPV genotypes bound to the catalytic domain of human nonreceptor-type protein tyrosine phosphatase 14 (PTPN14). They interact directly and potently with each other, with a dissociation constant of 18.2 nM. Ensuing structural analysis revealed the molecular basis of the PTPN14-binding specificity of E7 over other protein tyrosine phosphatases and also led to the identification of PTPN21 as a direct interacting partner of E7. Disruption of HPV18 E7 binding to PTPN14 by structure-based mutagenesis impaired E7’s ability to promote keratinocyte proliferation and migration. Likewise, E7 binding-defective PTPN14 was resistant for degradation via proteasome, and it was much more effective than wild-type PTPN14 in attenuating the activity of downstream effectors of Hippo signaling and negatively regulating cell proliferation, migration, and invasion when examined in HPV18-positive HeLa cells. These results therefore demonstrated the significance and therapeutic potential of the intermolecular interaction between HPV E7 and host PTPN14 in HPV-mediated cell transformation and tumorigenesis.

The design and discovery of selective cyclin-dependent kinase 4 (CDK4) inhibitors have been actively pursued in order to develop therapeutic cancer treatments. By means of a consecutive computational protocol involving homology modeling, docking experiments, and molecular dynamics simulations, we examine the characteristic structural and dynamic properties that distinguish CDK4 from CDK2 in its complexation with selective inhibitors. The results for all three CDK4-selective inhibitors under investigation show that the large-amplitude motion of a disordered loop of CDK4 is damped out in the presence of the inhibitors whereas their binding in the CDK2 active site has little effect on the loop flexibility. It is also found that the binding preference of CDK4- selective inhibitors for CDK4 over CDK2 stems from the reduced solvent accessibility in the active site of the former due to the formation of a stable hydrogen-bond triad by the Asp99, Arg101, and Thr102 side chains at the top of the active-site gorge. Besides the differences in loop flexibility and solvent accessibility, the dynamic stabilities of the hydrogen bonds between the inhibitors and the side chain of the lysine residue at the bottom of the active site also correlate well with the relative binding affinities of the inhibitors for the two CDKs. These results highlight the usefulness of this computational approach in evaluating the selectivity of a CDK inhibitor, and demonstrate the necessity of considering protein flexibility and solvent effects in designing new selective CDK4-selective inhibitors.

On the basis of ab initio calculations at the MP2/6-31G**//RHF/6-31G** level, we study the model reactions for the catalytic action of aspartatic proteinases. We elucidate the mechanistic features of two competing catalytic mechanisms by determining the reaction paths. In contrast to the previous theoretical studies which neglected the electron correlation, the concerted and the stepwise pathways are predicted to be almost equally favored in the general acid/general base mechanism. On the other hand, we find that a concerted reaction pathway is preferred to the stepwise one in the nucleophilic mechanism. We also find that both nucleophilic and general-acid/general-base mechanisms may be operative in a peptide hydrolysis by aspartic proteinases. For the model reaction under consideration, the former is energetically more favored if one considers only the potential energy profile along the intrinsic reaction coordinate, as has been done in previous theoretical studies. However, when the entropic contribution due to nuclear motions as well as the zero-point vibrational energies is included, the latter is predicted to be the preferred one by 1.9 kcal/mol. The covalent intermediate, which is the end-point minimum energy complex on the concerted nucleophilic pathway, is found to be a unstable one that is 18.1 kcal/mol higher in free energy than the incipient model enzyme-substrate complex. This fact is in accordance with the previous experimental implications that decomposition of the covalent intermediate is much faster than its formation. The present work provides a theoretical support for the persistent argument that the possibility of the nucleophilic mechanism cannot be excluded in the catalytic action of aspartic proteinases although the required experimental detection of the covalent intermediate has been unsuccessful so far. It is demonstrated that for all reaction pathways under consideration, the protonation of the nitrogen atom belonging to the peptidic bond is an essential step in crossing the activation barrier for the rupture of a peptide substrate. The relevance of the mechanistic features observed for the model reactions to an enzymatic reaction is discussed.