Androgen receptor (AR) is a ligand-activated transcription factor that plays a pivotal role in the development and progression of many severe diseases such as prostate cancer, muscle atrophy, and osteoporosis. Binding of ligands to AR triggers the conformational changes in AR that may affect the recruitment of coactivators and downstream response of AR signaling pathway. Therefore, AR ligands have great potential to treat these diseases. In this study, we searched for novel AR ligands by performing a docking-based virtual screening (VS) on the basis of the crystal structure of the AR ligand binding domain (LBD) in complex with its agonist. A total of 58 structurally diverse compounds were selected and subjected to LBD affinity assay, with five of them (HBP1-3, HBP1-17, HBP1-38, HBP1-51, and HBP1-58) exhibiting strong binding to AR-LBD. The IC50 values of HBP1-51 and HBP1-58 are 3.96 µM and 4.92 µM, respectively, which are even lower than that of enzalutamide (Enz, IC50 = 13.87 µM), a marketed second-generation AR antagonist. Further bioactivity assays suggest that HBP1-51 is an AR agonist, whereas HBP1-58 is an AR antagonist. In addition, molecular dynamics (MD) simulations and principal components analysis (PCA) were carried out to reveal the binding principle of the newly-identified AR ligands toward AR. Our modeling results indicate that the conformational changes of helix 12 induced by the bindings of antagonist and agonist are visibly different. In summary, the current study provides a highly efficient way to discover novel AR ligands, which could serve as the starting point for development of new therapeutics for AR-related diseases.

Phosphorylation is one of the most common post-translational modifications. The phosphorylation of the kinase-inducible domain (KID), which is an intrinsically disordered protein (IDP), promotes the folding of KID and binding with the KID-interacting domain (KIX). However, the regulation mechanism of the phosphorylation on KID is still elusive. In this study, the structural ensembles and binding process of pKID and KIX are studied by all-atom enhanced sampling technologies. The results show that more hydrophobic interactions are formed in pKID, which promote the formation of the special hydrophobic residue cluster (HRC). The pre-formed HRC promotes binding to the correct sites of KIX and further lead the folding of pKID. Consequently, a flexible conformational selection model is proposed to describe the binding and folding process of intrinsically disordered proteins. The binding mechanism revealed in this work provides new insights into the dynamic interactions and phosphorylation regulation of proteins.

The aggregation of human islet amyloid polypeptide (hIAPP) can damage the membrane of the β-cells in the pancreatic islets and induce type 2 diabetes (T2D). Growing evidences indicated that the major toxic species are small oligomers of IAPP. Due to the fast aggregation nature, it is hard to characterize the structures of IAPP oligomers by experiments, especially in the complex membrane environment. On the other side, molecular dynamics simulation can provide atomic details of the structure and dynamics of the aggregation of IAPP. In this study, all-atom bias-exchange metadynamics (BE-Meta) and unbiased molecular dynamics simulations were employed to study the structural properties of IAPP dimer in the membranes environments. A number of intermediates, including α-helical states, β-sheet states, and fully disordered states, are identified. The formation of N-terminal β-sheet structure is prior to the C-terminal β-sheet structure towards the final fibril-like structures. The α-helical intermediates have lower propensity in the dimeric hIAPP and are off-pathway intermediates. The simulations also demonstrate that the β-sheet intermediates induce more perturbation on the membrane than the α-helical and disordered states and thus pose higher disruption ability.