Selective recognition and discrimination of single isomeric changes in peptide strands with a host : guest sensing array.

An indirect competitive binding mechanism can be exploited to allow a combination of cationic fluorophores and water-soluble synthetic receptors to selectively recognize and discriminate peptide strands containing a single isomeric residue in the backbone. Peptide isomerization occurs in long-lived proteins and has been linked with diseases such as Alzheimer's, cataracts and cancer, so isomers are valuable yet underexplored targets for selective recognition. Planar cationic fluorophores can selectively bind hydrophobic, Trp-containing peptide strands in solution, and when paired with receptors that provide a competitive host for the fluorophore, can form a differential sensing array that enables selective discrimination of peptide isomers. Residue variations such as D- and L-Asp, D- and L-isoAsp, D-Ser and D-Glu can all be recognized, simply by their effects on the folded structure of the flexible peptide. Molecular dynamics simulations were applied to determine the most favorable conformation of the peptide : fluorophore conjugate, indicating that favorable π-stacking with internal tryptophan residues in a folded binding pocket enables micromolar binding affinity.

MYC-Targeting Inhibitors Generated from a Stereodiversified Bicyclic Peptide Library.

Here, we present the second generation of our bicyclic peptide library (NTB), featuring a stereodiversified structure and a simplified construction strategy. We utilized a tandem ring-opening metathesis and ring-closing metathesis reaction (ROM-RCM) to cyclize the linear peptide library in a single step, representing the first reported instance of this reaction being applied to the preparation of macrocyclic peptides. Moreover, the resulting bicyclic peptide can be easily linearized for MS/MS sequencing with a one-step deallylation process. We employed this library to screen against the E363-R378 epitope of MYC and identified several MYC-targeting bicyclic peptides. Subsequent in vitro cell studies demonstrated that one candidate, NT-B2R, effectively suppressed MYC transcription activities and cell proliferation.

An Integrated Pharmacological, Structural, and Genetic Analysis of Extracellular Versus Intracellular ROS Production in Neutrophils.

The neutrophil NADPH oxidase produces both intracellular and extracellular reactive oxygen species (ROS). Although oxidase activity is essential for microbial killing, and ROS can act as signaling molecules in the inflammatory process, excessive extracellular ROS directly contributes to inflammatory tissue damage, as well as to cancer progression and immune dysregulation in the tumor microenvironment. How specific signaling pathways contribute to ROS localization is unclear. Here we used a systems pharmacology approach to identify the specific Class I PI3-K isoform p110ß, and PLD1, but not PLD2, as critical regulators of extracellular, but not intracellular ROS production in primary neutrophils. Combined crystallographic and molecular dynamics analysis of the PX domain of the oxidase component p47phox, which binds the lipid products of PI 3-K and PLD, was used to clarify the membrane-binding mechanism and guide the design of mutant mice whose p47phox is unable to bind 3-phosphorylated inositol phospholipids. Neutrophils from these K43A mutant animals were specifically deficient in extracellular, but not intracellular, ROS production, and showed increased dependency on signaling through the remaining PLD1 arm. These findings identify the PX domain of p47phox as a critical integrator of PLD1 and p110ß signaling for extracellular ROS production, and as a potential therapeutic target for modulating tissue damage and extracellular signaling during inflammation.

Discovery of antimicrobial agent targeting tryptophan synthase.

Antibiotic resistance is a continually growing challenge in the treatment of various bacterial infections worldwide. New drugs and new drug targets are necessary to curb the threat of infectious diseases caused by multidrug-resistant pathogens. The tryptophan biosynthesis pathway is essential for bacterial growth but is absent in higher animals and humans. Drugs that can inhibit the bacterial biosynthesis of tryptophan offer a new class of antibiotics. In this work, we combined a structure-based strategy using in silico docking screening and molecular dynamics (MD) simulations to identify compounds targeting the α subunit of tryptophan synthase with experimental methods involving the whole-cell minimum inhibitory concentration (MIC) test, solution state NMR, and crystallography to confirm the inhibition of L-tryptophan biosynthesis. Screening 1,800 compounds from the National Cancer Institute Diversity Set I against α subunit revealed 28 compounds for experimental validation; four of the 28 hit compounds showed promising activity in MIC testing. We performed solution state NMR experiments to demonstrate that a one successful inhibitor, 3-amino-3-imino-2-phenyldiazenylpropanamide (Compound 1) binds to the α subunit. We also report a crystal structure of Salmonella enterica serotype Typhimurium tryptophan synthase in complex with Compound 1 which revealed a binding site at the αß interface of the dimeric enzyme. MD simulations were carried out to examine two binding sites for the compound. Our results show that this small molecule inhibitor could be a promising lead for future drug development.

Proteome-Wide Characterizations of N6-Methyl-Adenosine Triphosphate- and N6-Furfuryl-Adenosine Triphosphate-Binding Capabilities of Kinases.

Kinases catalyze the transfer of the γ-phosphate group from adenosine triphosphate (ATP) to their protein and small-molecule substrates, and this phosphorylation is a crucial element of multiple cell signaling pathways. Herein, we employed isotope-coded ATP acyl-phosphate probes, in conjunction with a multiple-reaction monitoring (MRM)-based targeted proteomic method for proteome-wide identifications of endogenous kinases that can bind to two N6-modified ATP derivatives, N6-methyl-ATP (N6-Me-ATP), and N6-furfuryl-ATP (a.k.a. kinetin triphosphate, KTP). We found that, among the ∼300 quantified kinases, 27 and 18 are candidate kinases that can bind to KTP and N6-Me-ATP, respectively. Additionally, GSK3α and GSK3ß are among the kinases that can bind to both ATP analogues. Moreover, the in vitro biochemical assay showed that GSK3ß could employ N6-Me-ATP but not KTP as the phosphate group donor to phosphorylate its substrate peptide. Molecular modeling studies provided insights into the differences between N6-Me-ATP and KTP in enabling the GSK3ß-mediated phosphorylation. Together, our chemoproteomic approach led to the identification of endogenous kinases that can potentially be targeted by the two ATP analogues. The approach should be generally applicable for assessing endogenous kinases targeted by other ATP and purine analogues.

Structure and noncanonical Cdk8 activation mechanism within an Argonaute-containing Mediator kinase module.

The Cdk8 kinase module (CKM) in Mediator, comprising Med13, Med12, CycC, and Cdk8, regulates RNA polymerase II transcription through kinase-dependent and -independent functions. Numerous pathogenic mutations causative for neurodevelopmental disorders and cancer congregate in CKM subunits. However, the structure of the intact CKM and the mechanism by which Cdk8 is non-canonically activated and functionally affected by oncogenic CKM alterations are poorly understood. Here, we report a cryo-electron microscopy structure of Saccharomyces cerevisiae CKM that redefines prior CKM structural models and explains the mechanism of Med12-dependent Cdk8 activation. Med12 interacts extensively with CycC and activates Cdk8 by stabilizing its activation (T-)loop through conserved Med12 residues recurrently mutated in human tumors. Unexpectedly, Med13 has a characteristic Argonaute-like bi-lobal architecture. These findings not only provide a structural basis for understanding CKM function and pathological dysfunction, but also further impute a previously unknown regulatory mechanism of Mediator in transcriptional modulation through its Med13 Argonaute-like features.

Inhibiting Matrix Metalloproteinase-2 Activation by Perturbing Protein-Protein Interactions Using a Cyclic Peptide.

We report on a cyclic peptide that inhibits matrix metalloproteinase-2 (MMP2) activation with a low-nM-level potency. This inhibitor specifically binds to the D570-A583 epitope on proMMP2 and interferes with the protein-protein interaction (PPI) between proMMP2 and tissue inhibitor of metalloproteinases-2 (TIMP2), thereby preventing the TIMP2-assisted proMMP2 activation process. We developed this cyclic peptide inhibitor through an epitope-targeted library screening process and validated its binding to proMMP2. Using a human melanoma cell line, we demonstrated the cyclic peptide's ability to modulate cellular MMP2 activities and inhibit cell migration. These results provide the first successful example of targeting the PPI between proMMP2 and TIMP2, confirming the feasibility of an MMP2 inhibition strategy that has been sought after for 2 decades.

Discovery of CDK8/CycC Ligands with a New Virtual Screening Tool.

Selective inhibition of cyclin-dependent kinase 8 and cyclin C (CDK8/CycC) has been suggested as a promising strategy for decreasing mitogenic signals in cancer cells with reduced toxicity toward normal cells. We developed a novel virtual screening protocol for drug development and applied it to the discovery of new CDK8/CycC type II ligands, which is likely to achieve long residence time and specificity. We first analyzed the binding thermodynamics of 11 published pyrazolourea ligands using molecular dynamics simulations and a free-energy calculation method, VM2, and extracted the key binding information to assist virtual screening. The urea moiety was found to be the critical structural contributor of the reference ligands. Starting with the urea moiety, we conducted substructure-based searches with our newly developed superposition and single-point energy evaluation method, followed by free-energy calculations, and singled out three purchasable compounds for bioassay testing. The ranking from the experimental results is completely consistent with the predicted rankings. A potent drug-like compound was found to have a Kd value of 42.5 nm, which is similar to those of the most potent reference ligands; this provided a good starting point for further improvement. This study shows that our novel virtual screening protocol is an accurate and efficient tool for drug development.

Lecithin-Stabilized Polymeric Micelles (LsbPMs) for Delivering Quercetin: Pharmacokinetic Studies and Therapeutic Effects of Quercetin Alone and in Combination with Doxorubicin.

In this study, lecithin-stabilized polymeric micelles (LsbPMs) were prepared to load quercetin (QUE) in order to improve its bioavailability and increase its antitumor activity. Its combination with doxorubicin (DOX) to minimize DOX-mediated cardiac toxicity and increase the antitumor activity of QUE-loaded LsbPMs was also examined. LsbPMs were prepared following a previously reported procedure. Results demonstrated that optimal QUE-loaded LsbPMs contained quercetin, D-α-tocopheryl polyethylene glycol succinate, and lecithin at a weight ratio of 6:40:80. Drug-release studies showed that QUE released from LsbPMs followed a controlled release pattern. A cytotoxicity assay revealed that QUE-loaded LsbPMs had significant anticancer activities against MCF-7, SKBR-3, and MDA-MB-231 human breast cancer cells and CT26 mouse colon cancer cells. In animal studies, intravenous administration of QUE-loaded LsbPMs resulted in efficient growth inhibition of CT26 colon cancer cells in a Balb/c mice model. In a pharmacokinetics study compared to free QUE, intravenous and oral administration of QUE-loaded LsbPMs was found to have significantly increased the relative bioavailability to 158% and 360%, respectively, and the absolute bioavailability to 5.13%. The effect of QUE-loaded LsbPMs in combination with DOX resulted in efficient growth inhibition of CT26 colon cancer cells and reduced cardiac toxicity in the Balb/c mice model.

Construction of a Sequenceable Protein Mimetic Peptide Library with a True 3D Diversifiable Chemical Space.

We present here a library of protein mimetic bicyclic peptides. These nanosized structures exhibit rigid backbones and spatially diversifiable side chains. They present modular amino acids on all three linkages, providing access to a true 3D diversifiable chemical space. These peptides are synthesized through a Cu-catalyzed click reaction and a Ru-catalyzed ring-closing metathesis reaction. Their bicyclic topology can be reduced to a linear one, using Edman degradation and Pd-catalyzed deallylation reactions. The linearization approaches allow de novo sequencing through mass spectrometry methods. We demonstrate the function of a particular peptide that was identified through a high throughput screening against the E363-R378 epitope on the intrinsically disordered c-Myc oncoprotein. Intracellular delivery of this peptide could interfere with the c-Myc-mediated transcription and inhibit proliferation in a human glioblastoma cell line.

A molecular dynamics investigation of CDK8/CycC and ligand binding: conformational flexibility and implication in drug discovery.

Abnormal activity of cyclin-dependent kinase 8 (CDK8) along with its partner protein cyclin C (CycC) is a common feature of many diseases including colorectal cancer. Using molecular dynamics (MD) simulations, this study determined the dynamics of the CDK8-CycC system and we obtained detailed breakdowns of binding energy contributions for four type-I and five type-II CDK8 inhibitors. We revealed system motions and conformational changes that will affect ligand binding, confirmed the essentialness of CycC for inclusion in future computational studies, and provide guidance in development of CDK8 binders. We employed unbiased all-atom MD simulations for 500 ns on twelve CDK8-CycC systems, including apoproteins and protein-ligand complexes, then performed principal component analysis (PCA) and measured the RMSF of key regions to identify protein dynamics. Binding pocket volume analysis identified conformational changes that accompany ligand binding. Next, H-bond analysis, residue-wise interaction calculations, and MM/PBSA were performed to characterize protein-ligand interactions and find the binding energy. We discovered that CycC is vital for maintaining a proper conformation of CDK8 to facilitate ligand binding and that the system exhibits motion that should be carefully considered in future computational work. Surprisingly, we found that motion of the activation loop did not affect ligand binding. Type-I and type-II ligand binding is driven by van der Waals interactions, but electrostatic energy and entropic penalties affect type-II binding as well. Binding of both ligand types affects protein flexibility. Based on this we provide suggestions for development of tighter-binding CDK8 inhibitors and offer insight that can aid future computational studies.

Novel application of pluronic lecithin organogels (PLOs) for local delivery of synergistic combination of docetaxel and cisplatin to improve therapeutic efficacy against ovarian cancer.

The synergistic combination of docetaxel (DTX) and cisplatin (CIS) by local drug delivery with a pluronic lecithin organogel (PLO) to facilitate high drug concentrations at tumor sites and less nonspecific distribution to normal organs is thought to be beneficial in chemotherapy. In this study, using Capryol-90 (C90) with the addition of lecithin as the oil phase was developed to carry DTX, which was then incorporated into a PLO-containing CIS to formulate a dual-drug injectable PLO for local delivery. An optimal PLO composite, P13L0.15O1.5, composed of PF127:lecithin:C90 at a 13:0.15:1.5 weight ratio was obtained. The sol-gel transition temperature of P13L0.15O1.5 was found to be 33 °C. Tumor inhibition studies illustrated that DTX/CIS-loaded P13L0.15O1.5 could efficiently suppress tumor growth by both intratumoral and peritumoral injections in SKOV-3 xenograft mouse model. Pharmacokinetic studies showed that subcutaneous administration of P13L0.15O1.5 was able to sustain the release of DTX and CIS leading to their slow absorption into the systemic circulation resulting in lower area under the plasma concentration curve at 0-72 h (AUC0-72) and maximum concentration (Cmax) values but longer half-life (T1/2) and mean residence time (MRT) values. An in vivo biodistribution study showed lower DTX and CIS concentrations in organs compared to other treatment groups after IT administration of the dual drug-loaded P13L0.15O1.5. It was concluded that the local co-delivery of DTX and CIS by PLOs may be a promising and effective platform for local anticancer drug delivery with minimal systemic toxicities.

Immunoproteasome inhibition and bioactivity of thiasyrbactins.

A family of macrodilactam natural products, the syrbactins, are known proteasome inhibitors. A small group of syrbactin analogs was prepared with a sulfur-for-carbon substitution to enhance synthetic accessibility and facilitate modulation of their solubility. Two of these compounds surprisingly proved to be inhibitors of the trypsin-like catalytic site, including of the immunoproteasome. Their bound and free conformations suggest special properties of the thiasyrbactin ring are responsible for this unusual preference, which may be exploited to develop drug-like immunoproteasome inhibitors. These compounds show greater selectivity than earlier compounds used to infer phenotypes of immunoproteasome inhibition, like ONX-0914.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design.

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Codelivery of doxorubicin-containing thermosensitive hydrogels incorporated with docetaxel-loaded mixed micelles enhances local cancer therapy.

Doxorubicin (DOX) thermosensitive hydrogels (TSHs) incorporated with docetaxel (DOC)-loaded mixed micelles were developed to co-deliver these two drugs through a TSH system, DH700kMF-13.5/M-DocLF, to improve local cancer therapy and reduce side effects. First, Pluronics-based DOC-loaded mixed micelles were developed and optimized. The optimal formulation designated as M-DocLF was composed of 1mg/g docetaxel, 15mg/g Pluronic F127 (PF127), and 45mg/g Pluronic L121 (PL121). Rheological tests showed that DH700kMF-13.5/M-DocLF was an injectable flowing solution, which formed a nonflowing gel at body temperature. After intratumoral (IT) or peritumoral (PT) administration, DH700kMF-13.5/M-DocLF demonstrated efficient growth inhibition of CT-26 tumors in a Balb/c mice model. The tumor inhibitory rate after IT administration of DH700kMF-13.5/M-DocLF was 92.4%, followed by 85.8%, 75.6%, 62.9%, 50.6%, and 49.5% for DH700kMF-15, free DOX, F-13.5/M-DocLF, Tynen (DOC solution), and M-DocLF, respectively. Furthermore, PT administration of DH700kMF-13.5/M-DocLF resulted in similar efficacies. Pharmacokinetic and biodistribution studies showed that after subcutaneous (SC) and IT administration of the designated formulations, smaller amounts of DOX and DOC were absorbed from the local SC or tumor sites into systemic circulation, probably reducing their systemic toxicity. Tumor retention of DOX and DOC in biodistribution studies further revealed that co-delivery of these two drugs in DH700KMF-13.5/M-DocLF potentially enhanced the efficacy of tumor inhibition. In conclusion, our in situ injectable DOX and DOC TSH is a potential dual drug delivery system, which can enhance the efficacy of cancer chemotherapy with minimal side effects and reduced chemoresistance.

Syrbactin Structural Analog TIR-199 Blocks Proteasome Activity and Induces Tumor Cell Death.

Multiple myeloma is an aggressive hematopoietic cancer of plasma cells. The recent emergence of three effective FDA-approved proteasome-inhibiting drugs, bortezomib (Velcade®), carfilzomib (Kyprolis®), and ixazomib (Ninlaro®), confirms that proteasome inhibitors are therapeutically useful against neoplastic disease, in particular refractory multiple myeloma and mantle cell lymphoma. This study describes the synthesis, computational affinity assessment, and preclinical evaluation of TIR-199, a natural product-derived syrbactin structural analog. Molecular modeling and simulation suggested that TIR-199 covalently binds each of the three catalytic subunits (ß1, ß2, and ß5) and revealed key interaction sites. In vitro and cell culture-based proteasome activity measurements confirmed that TIR-199 inhibits the proteasome in a dose-dependent manner and induces tumor cell death in multiple myeloma and neuroblastoma cells as well as other cancer types in the NCI-60 cell panel. It is particularly effective against kidney tumor cell lines, with >250-fold higher anti-tumor activities than observed with the natural product syringolin A. In vivo studies in mice revealed a maximum tolerated dose of TIR-199 at 25 mg/kg. The anti-tumor activity of TIR-199 was confirmed in hollow fiber assays in mice. Adverse drug reaction screens in a kidney panel revealed no off-targets of concern. This is the first study to examine the efficacy of a syrbactin in animals. Taken together, the results suggest that TIR-199 is a potent new proteasome inhibitor with promise for further development into a clinical drug for the treatment of multiple myeloma and other forms of cancer.

Endosidin2 targets conserved exocyst complex subunit EXO70 to inhibit exocytosis.

The exocyst complex regulates the last steps of exocytosis, which is essential to organisms across kingdoms. In humans, its dysfunction is correlated with several significant diseases, such as diabetes and cancer progression. Investigation of the dynamic regulation of the evolutionarily conserved exocyst-related processes using mutants in genetically tractable organisms such as Arabidopsis thaliana is limited by the lethality or the severity of phenotypes. We discovered that the small molecule Endosidin2 (ES2) binds to the EXO70 (exocyst component of 70 kDa) subunit of the exocyst complex, resulting in inhibition of exocytosis and endosomal recycling in both plant and human cells and enhancement of plant vacuolar trafficking. An EXO70 protein with a C-terminal truncation results in dominant ES2 resistance, uncovering possible distinct regulatory roles for the N terminus of the protein. This study not only provides a valuable tool in studying exocytosis regulation but also offers a potentially new target for drugs aimed at addressing human disease.

The paradox of conformational constraint in the design of Cbl(TKB)-binding peptides.

Solving the crystal structure of Cbl(TKB) in complex with a pentapeptide, pYTPEP, revealed that the PEP region adopted a poly-L-proline type II (PPII) helix. An unnatural amino acid termed a proline-templated glutamic acid (ptE) that constrained both the backbone and sidechain to the bound conformation was synthesized and incorporated into the pYTPXP peptide. We estimated imposing structural constraints onto the backbone and sidechain of the peptide and preorganize it to the bound conformation in solution will yield nearly an order of magnitude improvement in activity. NMR studies confirmed that the ptE-containing peptide adopts the PPII conformation, however, competitive binding studies showed an order of magnitude loss of activity. Given the emphasis that is placed on imposing structural constraints, we provide an example to support the contrary. These results point to conformational flexibility at the interface, which have implications in the design of potent Cbl(TKB)-binding peptides.

Ligand-specific homology modeling of human cannabinoid (CB1) receptor.

Cannabinoid (CB1) receptor is a therapeutic drug target, and its structure and conformational changes after ligand binding are of great interest. To study the protein conformations in ligand bound state and assist in drug discovery, CB1 receptor homology models are needed for computer-based ligand screening. The known CB1 ligands are highly diverse structurally, so CB1 receptor may undergo considerable conformational changes to accept different ligands, which is challenging for molecular docking methods. To account for the flexibility of CB1 receptor, we constructed four CB1 receptor models based on four structurally distinct ligands, HU-210, ACEA, WIN55212-2 and SR141716A, using the newest X-ray crystal structures of human ß2 adrenergic receptor and adenosine A(2A) receptor as templates. The conformations of these four CB1-ligand complexes were optimized by molecular dynamics (MD) simulations. The models revealed interactions between CB1 receptor and known binders suggested by experiments and could successfully discriminate known ligands and non-binders in our docking assays. MD simulations were used to study the most flexible ligand, ACEA, in its free and bound states to investigate structural mobility achieved by the rearrangement of the fatty acid chain. Our models may capture important conformational changes of CB1 receptor to help improve accuracy in future CB1 drug screening.

Mechanistic insights into phosphopeptide--BRCT domain association: preorganization, flexibility, and phosphate recognition.

Promiscuous proteins are commonly observed in biological systems, for example, in modular domains that recognize phosphopeptides during signal transduction. This promiscuous recognition is of fundamental interest in chemistry and biology but is challenging when designing phosphopeptides in silico for cell biology studies. To investigate promiscuous recognition and binding processes of phosphopeptides and the modular domain, we selected a domain essential in breast cancer-the breast-cancer-associated protein 1 (BRCA1) C-terminal (BRCT) repeats as our model system. We performed molecular dynamics simulations and detailed analyses of the dihedral space to study protein fluctuation and conformational changes with phosphopeptide binding. We also studied the association processes of phosphorylated and unphosphorylated peptides using Brownian dynamics with a coarse-grained model. We found that the BRCT domain is preorganized for phosphopeptide binding but has a moderate arrangement of side chains to form complexes with various types of phosphopeptides. Phosphopeptide binding restricts the system motion in general, while the nonpolar phosphopeptide becomes more flexible in the bound state. Our analysis found that the BRCT domain utilizes different mechanisms, usually termed lock and key, induced-fit, and population-shift/conformational-selection models, to recognize peptides with different features. Brownian dynamics simulations revealed that the charged phosphate group may not always accelerate peptide association processes, but it helps the phosphopeptide orient into binding pockets accurately and stabilizes the complex. This work provides insights into molecular recognition in the promiscuous protein system.