SUMMARY The target profiles of many drugs are established early in their development and are not systematically revisited at the time of FDA approval. Thus, it is often unclear whether therapeutics with the same nominal targets but different chemical structures are functionally equivalent. In this paper we use five different phenotypic and biochemical assays to compare approved inhibitors of cyclin-dependent kinases 4/6 – collectively regarded as breakthroughs in the treatment of hormone receptor-positive breast cancer. We find that transcriptional, proteomic, and phenotypic changes induced by palbociclib, ribociclib, and abemaciclib differ significantly; abemaciclib in particular has advantageous activities partially overlapping those of alvocidib, an older polyselective CDK inhibitor. In cells and mice, abemaciclib inhibits kinases other than CDK4/6 including CDK2/Cyclin A/E – implicated in resistance to CDK4/6 inhibition – and CDK1/Cyclin B. The multi-faceted experimental and computational approaches described here therefore uncover under-appreciated differences in CDK4/6 inhibitor activities with potential importance in treating human patients.

Word models (natural language descriptions of molecular mechanisms) are a common currency in spoken and written communication in biomedicine but are of limited use in predicting the behavior of complex biological networks. We present an approach to building computational models directly from natural language using automated assembly. Molecular mechanisms described in simple English are read by natural language processing algorithms, converted into an intermediate representation, and assembled into executable or network models. We have implemented this approach in the Integrated Network and Dynamical Reasoning Assembler (INDRA), which draws on existing natural language processing systems as well as pathway information in Pathway Commons and other online resources. We demonstrate the use of INDRA and natural language to model three biological processes of increasing scope: (i) p53 dynamics in response to DNA damage, (ii) adaptive drug resistance in BRAF‐V600E‐mutant melanomas, and (iii) the RAS signaling pathway. The use of natural language makes the task of developing a model more efficient and it increases model transparency, thereby promoting collaboration with the broader biology community.

Graphical Abstract Highlights d Implementing FAIR data standards requires identification of experimental confounders d Five labs performed the same experiment on mammalian cells and compared results d Several factors affecting reproducibility were explored d Biological context had an unexpected impact on the robustness of cell-based assays