Search

Abstract Accurate molecular representation of compounds is a fundamental challenge for prediction of drug targets and molecular properties. In this study, we present a molecular video-based foundation model, named VideoMol, pretrained on 120 million frames of 2 million unlabeled drug-like and bioactive molecules. VideoMol renders each molecule as a video with 60-frame and designs three self-supervised learning strategies on molecular videos to capture molecular representation. We show high performance of VideoMol in predicting molecular targets and properties across 43 drug discovery benchmark datasets. VideoMol achieves high accuracy in identifying antiviral molecules against common diverse disease-specific drug targets (i.e., BACE1 and EP4). Drugs screened by VideoMol show better binding affinity than molecular docking, revealing the effectiveness in understanding the three-dimensional structure of molecules. We further illustrate interpretability of VideoMol using key chemical substructures.

Abstract Development of potent and broad-spectrum antimicrobial peptides (AMPs) could help overcome the antimicrobial resistance crisis. We develop a peptide language-based deep generative framework (deepAMP) for identifying potent, broad-spectrum AMPs. Using deepAMP to reduce antimicrobial resistance and enhance the membrane-disrupting abilities of AMPs, we identify, synthesize, and experimentally test 18 T1-AMP (Tier 1) and 11 T2-AMP (Tier 2) candidates in a two-round design and by employing cross-optimization-validation. More than 90% of the designed AMPs show a better inhibition than penetratin in both Gram-positive (i.e., S. aureus) and Gram-negative bacteria (i.e., K. pneumoniae and P. aeruginosa). T2-9 shows the strongest antibacterial activity, comparable to FDA-approved antibiotics. We show that three AMPs (T1-2, T1-5 and T2-10) significantly reduce resistance to S. aureus compared to ciprofloxacin and are effective against skin wound infection in a female wound mouse model infected with P. aeruginosa. In summary, deepAMP expedites discovery of effective, broad-spectrum AMPs against drug-resistant bacteria.

Abstract RHOA mutations are found at diverse residues in various cancer types, implying mutation- and cell-specific mechanisms of tumorigenesis. Here, we focus on the underlying mechanisms of two gain-of-function RHOA mutations, A161P and A161V, identified in adult T-cell leukemia/lymphoma. We find that RHOAA161P and RHOAA161V are both fast-cycling mutants with increased guanine nucleotide dissociation/association rates compared with RHOAWT and show reduced GTP-hydrolysis activity. Crystal structures reveal an altered nucleotide association in RHOAA161P and an open nucleotide pocket in RHOAA161V. Both mutations perturb the dynamic properties of RHOA switch regions and shift the conformational landscape important for RHOA activity, as shown by 31P NMR and molecular dynamics simulations. Interestingly, RHOAA161P and RHOAA161V can interact with effectors in the GDP-bound state. 1H-15N HSQC NMR spectra support the existence of an active population in RHOAA161V-GDP. The distinct interaction mechanisms resulting from the mutations likely favor an RHOAWT-like “ON” conformation, endowing GDP-bound state effector binding activity.

Abstract Nuclear factor erythroid 2-related factor 2 (NRF2) hyperactivation has been established as an oncogenic driver in a variety of human cancers, including non-small cell lung cancer (NSCLC). However, despite massive efforts, no specific therapy is currently available to target NRF2 hyperactivation. Here, we identify peptidylprolyl isomerase A (PPIA) is required for NRF2 protein stability. Ablation of PPIA promotes NRF2 protein degradation and blocks NRF2-driven growth in NSCLC cells. Mechanistically, PPIA physically binds to NRF2 and blocks the access of ubiquitin/Kelch Like ECH Associated Protein 1 (KEAP1) to NRF2, thus preventing ubiquitin-mediated degradation. Our X-ray co-crystal structure reveals that PPIA directly interacts with a NRF2 interdomain linker via a trans-proline 174-harboring hydrophobic sequence. We further demonstrate that an FDA-approved drug, cyclosporin A (CsA), impairs the interaction of NRF2 with PPIA, inducing NRF2 ubiquitination and degradation. Interestingly, CsA interrupts glutamine metabolism mediated by the NRF2/KLF5/SLC1A5 pathway, consequently suppressing the growth of NRF2-hyperactivated NSCLC cells. CsA and a glutaminase inhibitor combination therapy significantly retard tumor progression in NSCLC patient-derived xenograft (PDX) models with NRF2 hyperactivation. Our study demonstrates that targeting NRF2 protein stability is an actionable therapeutic approach to treat NRF2-hyperactivated NSCLC.

Pediatric low grade brain tumors and neurodevelopmental disorders share proteins, signaling pathways, and networks. They also share germline mutations and an impaired prenatal differentiation origin. They may differ in the timing of the events and proliferation. We suggest that their pivotal distinct, albeit partially overlapping, outcomes relate to the cell states, which depend on their spatial location, and timing of gene expression during brain development. These attributes are crucial as the brain develops sequentially, and single-cell spatial organization influences cell state, thus function. Our underlying premise is that the root cause in neurodevelopmental disorders and pediatric tumors is impaired prenatal differentiation. Data related to pediatric brain tumors, neurodevelopmental disorders, brain cell (sub)types, locations, and timing of expression in the developing brain are scant. However, emerging single cell technologies, including transcriptomic, spatial biology, spatial high-resolution imaging performed over the brain developmental time, could be transformational in deciphering brain pathologies thereby pharmacology.

The connection and causality between cancer and neurodevelopmental disorders have been puzzling. How can the same cellular pathways, proteins, and mutations lead to pathologies with vastly different clinical presentations? And why do individuals with neurodevelopmental disorders, such as autism and schizophrenia, face higher chances of cancer emerging throughout their lifetime? Our broad review emphasizes the multi-scale aspect of this type of reasoning. As these examples demonstrate, rather than focusing on a specific organ system or disease, we aim at the new understanding that can be gained. Within this framework, our review calls attention to computational strategies which can be powerful in discovering connections, causalities, predicting clinical outcomes, and are vital for drug discovery. Thus, rather than centering on the clinical features, we draw on the rapidly increasing data on the molecular level, including mutations, isoforms, three-dimensional structures, and expression levels of the respective disease-associated genes. Their integrated analysis, together with chromatin states, can delineate how, despite being connected, neurodevelopmental disorders and cancer differ, and how the same mutations can lead to different clinical symptoms. Here, we seek to uncover the emerging connection between cancer, including pediatric tumors, and neurodevelopmental disorders, and the tantalizing questions that this connection raises.

Abstract Epidemiological studies suggest that individuals with neurodevelopmental disorders (NDDs) are more prone to develop certain types of cancer. Notably, however, the case statistics can be impacted by late discovery of cancer in individuals afflicted with NDDs, such as intellectual disorders, autism, and schizophrenia, which may bias the numbers. As to NDD-associated mutations, in most cases, they are germline while cancer mutations are sporadic, emerging during life. However, somatic mosaicism can spur NDDs, and cancer-related mutations can be germline. NDDs and cancer share proteins, pathways, and mutations. Here we ask (i) exactly which features they share, and (ii) how, despite their commonalities, they differ in clinical outcomes. To tackle these questions, we employed a statistical framework followed by network analysis. Our thorough exploration of the mutations, reconstructed disease-specific networks, pathways, and transcriptome levels and profiles of autism spectrum disorder (ASD) and cancers, point to signaling strength as the key factor: strong signaling promotes cell proliferation in cancer, and weaker (moderate) signaling impacts differentiation in ASD. Thus, we suggest that signaling strength, not activating mutations, can decide clinical outcome.

Abstract Proper organization of intracellular assemblies is fundamental for efficient promotion of biochemical processes and optimal assembly functionality. Although advances in imaging technologies have shed light on how the centrosome is organized, how its constituent proteins are coherently architected to elicit downstream events remains poorly understood. Using multidisciplinary approaches, we showed that two long coiled-coil proteins, Cep63 and Cep152, form a heterotetrameric building block that undergoes a stepwise formation into higher molecular weight complexes, ultimately generating a cylindrical architecture around a centriole. Mutants defective in Cep63•Cep152 heterotetramer formation displayed crippled pericentriolar Cep152 organization, polo-like kinase 4 (Plk4) relocalization to the procentriole assembly site, and Plk4-mediated centriole duplication. Given that the organization of pericentriolar materials (PCM) is evolutionarily conserved, this work could serve as a model for investigating the structure and function of PCM in other species, while offering a new direction in probing the organizational defects of PCM-related human diseases.

Abstract Blood flow produces shear stress exerted on the endothelial layer of the vessels. Spatial characterization of the endothelial proteome is required to uncover the mechanisms of endothelial activation by shear stress, as blood flow varies in the vasculature. An integrative ubiquitinome and proteome analysis of shear‐stressed endothelial cells demonstrated that the non‐degradative ubiquitination of several GTPases is regulated by mechano‐signaling. Spatial analysis reveals increased ubiquitination of the small GTPase RAP1 in the descending aorta, a region exposed to laminar shear stress. The ubiquitin ligase WWP2 is identified as a novel regulator of RAP1 ubiquitination during shear stress response. Non‐degradative ubiquitination fine‐tunes the function of GTPases by modifying their interacting network. Specifically, WWP2‐mediated RAP1 ubiquitination at lysine 31 switches the balance from the RAP1/ Talin 1 (TLN1) toward RAP1/ Afadin (AFDN) or RAP1/ RAS Interacting Protein 1 (RASIP1) complex formation, which is essential to suppress shear stress‐induced reactive oxygen species (ROS) production and maintain endothelial barrier integrity. Increased ROS production in endothelial cells in the descending aorta of endothelial‐specific Wwp2‐knockout mice leads to increased levels of oxidized lipids and inflammation. These results highlight the importance of the spatially regulated non‐degradative ubiquitination of GTPases in endothelial mechano‐activation.

A statistical methodology is used to identify potential ‘latent driver’ mutations from public cancer bulk sequencing data.

RNA conformational heterogeneity is important to diverse functions. Here, the authors use AFM to directly visualize individual RNA molecules that are in various conformational states under near physiological solution conditions for the first time.

Abstract Blocked cellular differentiation is a critical pathologic hallmark of acute myeloid leukemia (AML). Here, we showed that genetic activation of the orphan GPCR GPR132 significantly induced cell differentiation of AML both in vitro and in vivo, indicating that GPR132 is a potential trigger of myeloid differentiation. To explore the therapeutic potential of GPR132 signaling, we screened and validated a natural product 8-gingerol (8GL) as a GPR132 agonist. Notably, GPR132 activation by 8GL promoted differentiation and reduced colony formation in human AML cell lines with diverse genetic profiles. Mechanistic studies revealed that 8GL treatment inhibits the activation of the mammalian target of rapamycin (mTOR), a regulator of AML cell differentiation blockade, via activating GPR132-Gs-PKA pathway. We further showed that the combination of 8GL and an mTOR inhibitor synergistically elicited AML cell differentiation in vitro. Importantly, 8GL alone or in combination with an mTOR inhibitor remarkably impaired tumor growth and extended mouse survival in an AML xenograft model accompanied by enhanced cell differentiation. Notably, genetic or pharmacological activation of GPR132 triggered the differentiation of human primary AML cells. In summary, this study demonstrated that activation of orphan GPR132 represents a potential strategy for inducing myeloid differentiation in AML patients.

Recognition of pathogen-derived nucleic acids by host cells is an effective host strategy to detect pathogenic invasion and trigger immune responses. In the context of pathogen-specific pharmacology, there is a growing interest in mapping the interactions between pathogen-derived nucleic acids and host proteins. Insight into the principles of the structural and immunological mechanisms underlying such interactions and their roles in host defense is necessary to guide therapeutic intervention. Here, we discuss the newest advances in studies of molecular interactions involving pathogen nucleic acids and host factors, including their drug design, molecular structure and specific patterns. We observed that two groups of nucleic acid recognizing molecules, Toll-like receptors (TLRs) and the cytoplasmic retinoic acid-inducible gene (RIG)-I-like receptors (RLRs) form the backbone of host responses to pathogen nucleic acids, with additional support provided by absent in melanoma 2 (AIM2) and DNA-dependent activator of Interferons (IFNs)-regulatory factors (DAI) like cytosolic activity. We review the structural, immunological, and other biological aspects of these representative groups of molecules, especially in terms of their target specificity and affinity and challenges in leveraging host-pathogen protein-nucleic acid interactions (HP-PNI) in drug discovery.

c-Abl, a non-receptor tyrosine kinase, regulates cell growth and survival in healthy cells and causes chronic myeloid leukemia (CML) when fused by Bcr. Its activity is blocked in the assembled inactive state, where the SH3 and SH2 domains dock into the kinase domain, reducing its conformational flexibility, resulting in the autoinhibited state. It is active in an extended ‘open’ conformation. Allostery governs the transitions between the autoinhibited and active states. Even though experiments revealed the structural hallmarks of the two states, a detailed grasp of the determinants of c-Abl autoinhibition and activation at the atomic level, which may help innovative drug discovery, is still lacking. Here, using extensive molecular dynamics simulations, we decipher exactly how these determinants regulate it. Our simulations confirm and extend experimental data that the myristoyl group serves as the switch for c-Abl inhibition/activation. Its dissociation from the kinase domain promotes the SH2-SH3 release, initiating c-Abl activation. We show that the precise SH2/N-lobe interaction is required for full activation of c-Abl. It stabilizes a catalysis-favored conformation, priming it for catalytic action. Bcr-Abl allosteric drugs elegantly mimic the endogenous myristoyl-mediated autoinhibition state of c-Abl 1b. Allosteric activating mutations shift the ensemble to the active state, blocking ATP-competitive drugs. Allosteric drugs alter the active-site conformation, shifting the ensemble to re-favor ATP-competitive drugs. Our work provides a complete mechanism of c-Abl activation and insights into critical parameters controlling at the atomic level c-Abl inactivation, leading us to propose possible strategies to counter reemergence of drug resistance.

Host-pathogen interactions (HPIs) affect and involve multiple mechanisms in both the pathogen and the host. Pathogen interactions disrupt homeostasis in host cells, with their toxins interfering with host mechanisms, resulting in infections, diseases, and disorders, extending from AIDS and COVID-19, to cancer. Studies of the three-dimensional (3D) structures of host-pathogen complexes aim to understand how pathogens interact with their hosts. They also aim to contribute to the development of rational therapeutics, as well as preventive measures. However, structural studies are fraught with challenges toward these aims. This review describes the state-of-the-art in protein-protein interactions (PPIs) between the host and pathogens from the structural standpoint. It discusses computational aspects of predicting these PPIs, including machine learning (ML) and artificial intelligence (AI)-driven, and overviews available computational methods and their challenges. It concludes with examples of how theoretical computational approaches can result in a therapeutic agent with a potential of being used in the clinics, as well as future directions.

The small GTPase Ran is the main regulator of the nucleo-cytoplasmic import and export through the nuclear pore complex. It functions as a molecular switch cycling between the GDP-bound inactive and GTP-bound active state. It consists of a globular (G) domain and a C-terminal region, which is bound to the G-domain in the inactive, GDP-bound states. Crystal structures of the GTP-bound active form complexed with Ran binding proteins (RanBP) show that the C-terminus undergoes a large conformational change, embracing Ran binding domains (RanBD). Whereas in the crystal structures of macromolecular complexes not containing RanBDs the structure of the C-terminal segment remains unresolved, indicating its large conformational flexibility. This movement could not have been followed either by experimental or simulation methods. Here, starting from the crystal structure of Ran in both GDP- and GTP-bound forms we show how rigid the C-terminal region in the inactive structure is during molecular dynamics (MD) simulations. Furthermore, we show how MD simulations of the active form are incapable of mapping the open conformations of the C-terminus. By using the MDeNM (Molecular Dynamics with excited Normal Modes) method, we were able to widely map the conformational surface of the C-terminus of Ran in the active GTP-bound form, which allows us to envisage how it can embrace RanBDs.

Summary: Immunity could be viewed as the common factor in neurodevelopmental disorders and cancer. The immune and nervous systems coevolve as the embryo develops. Immunity can release cytokines that activate MAPK signaling in neural cells. In specific embryonic brain cell types, dysregulated signaling that results from germline or embryonic mutations can promote changes in chromatin organization and gene accessibility, and thus expression levels of essential genes in neurodevelopment. In cancer, dysregulated signaling can emerge from sporadic somatic mutations during human life. Neurodevelopmental disorders and cancer share similarities. In neurodevelopmental disorders, immunity, and cancer, there appears an almost invariable involvement of small GTPases (e.g., Ras, RhoA, and Rac) and their pathways. TLRs, IL-1, GIT1, and FGFR signaling pathways, all can be dysregulated in neurodevelopmental disorders and cancer. Although there are signaling similarities, decisive differentiating factors are timing windows, and cell type specific perturbation levels, pointing to chromatin reorganization. Finally, we discuss drug discovery.

Abstract Which signaling pathway and protein to select to mitigate the patient’s expected drug resistance? The number of possibilities facing the physician is massive, and the drug combination should fit the patient status. Here, we briefly review current approaches and data and map an innovative patient-specific strategy to forecast drug resistance targets that centers on parallel (or redundant) proliferation pathways in specialized cells. It considers the availability of each protein in each pathway in the specific cell, its activating mutations, and the chromatin accessibility of its encoding gene. The construction of the resulting Proliferation Pathway Network Atlas will harness the emerging exascale computing and advanced artificial intelligence (AI) methods for therapeutic development. Merging the resulting set of targets, pathways, and proteins, with current strategies will augment the choice for the attending physicians to thwart resistance.

Abstract Massive genome sequencing data have inspired new challenges in personalized treatments and facilitated oncological drug discovery. We present a comprehensive database, My Personal Mutanome (MPM), for accelerating the development of precision cancer medicine protocols. MPM contains 490,245 mutations from over 10,800 tumor exomes across 33 cancer types in The Cancer Genome Atlas mapped to 94,563 structure-resolved/predicted protein-protein interaction interfaces (“edgetic”) and 311,022 functional sites (“nodetic”), including ligand-protein binding sites and 8 types of protein posttranslational modifications. In total, 8884 survival results and 1,271,132 drug responses are obtained for these mapped interactions. MPM is available at https://mutanome.lerner.ccf.org .

Oncogenic mutations in the serine/threonine kinase B-Raf, particularly the V600E mutation, are frequent in cancer, making it a major drug target. Although much is known about B-Raf’s active and inactive states, questions remain about the mechanism by which the protein changes between these two states. Here, we utilize molecular dynamics to investigate both wild-type and V600E B-Raf to gain mechanistic insights into the impact of the Val to Glu mutation. The results show that the wild-type and mutant follow similar activation pathways involving an extension of the activation loop and an inward motion of the αC-helix. The V600E mutation, however, destabilizes the inactive state by disrupting hydrophobic interactions present in the wild-type structure while the active state is stabilized through the formation of a salt bridge between Glu600 and Lys507. Additionally, when the activation loop is extended, the αC-helix is able to move between an inward and outward orientation as long as the DFG motif adopts a specific orientation. In that orientation Phe595 rotates away from the αC-helix, allowing the formation of a salt bridge between Lys483 and Glu501. These mechanistic insights have implications for the development of new Raf inhibitors.

Abstract Immune checkpoint blockade (ICB) has a limited effect on colorectal cancer, underlining the requirement of co‐targeting the complementary mechanisms. Here, we identified prostaglandin E2 (PGE2) receptor 4 (EP4) as the master regulator of immunosuppressive myeloid cells (IMCs), which are the major driver of resistance to ICB therapy. PGE2‐bound EP4 promotes the differentiation of immunosuppressive M2 macrophages and myeloid‐derived suppressor cells (MDSCs) and reduces the expansion of immunostimulated M1 macrophages. To explore the immunotherapeutic role of EP4 signaling, we developed a novel and selective EP4 antagonist TP‐16. TP‐16 effectively blocked the function of IMCs and enhanced cytotoxic T‐cell‐mediated tumor elimination in vivo. Cell co‐culture experiments revealed that TP‐16 promoted T‐cell proliferation, which was impaired by tumor‐derived CD11b+ myeloid cells. Notably, TP‐16 and anti‐PD‐1 combination therapy significantly impeded tumor progression and prolonged mice survival. We further demonstrated that TP‐16 increased responsiveness to anti‐PD‐1 therapy in an IMC‐related spontaneous colorectal cancer mouse model. In summary, this study demonstrates that inhibition of EP4‐expressing IMCs may offer a potential strategy for enhancing the efficacy of immunotherapy for colorectal cancer.

Proximally located in the membrane, oncogenic Ras dimers (or nanoclusters) can recruit and promote Raf dimerization and MAPK (Raf/MEK/ERK) signaling. Among Ras isoforms, KRas4B is the most frequently mutated. Recent data on the binary KRas4B–Raf-1 complex suggested that Raf-1 CRD not only executes membrane anchorage, but also supports the high-affinity interaction of Raf-1 RBD with KRas4B catalytic domain. For a detailed mechanistic picture of Raf activation at the membrane, we employ explicit MD simulations of the quaternary KRas4B–Raf-1 complex. The complex contains two active GTP-bound KRas4B proteins forming a dimer through the allosteric lobe interface and two tandem RBD-CRD segments of Raf-1 interacting with the effector lobes at both ends of the KRas4B dimer. We show that Raf-1 RBD-CRD supports stable KRas4B dimer at preferred interface and orientation at the membrane, thereby cooperatively enhancing the affinity of the KRas4B–Raf-1 interaction. We propose that a Ras dimer at the membrane can increase the population of proximal Raf kinase domains, promoting kinase domain dimerization in the cytoplasm. Collectively, the dynamic Ras–Raf assembly promotes Raf activation not by allostery; instead, Ras activates Raf by shifting its ensemble toward kinase domain-accessible states through enhanced affinity at the membrane.

Over the past decade, Markov State Models (MSM) have emerged as powerful methodologies to build discrete models of dynamics over structures obtained from Molecular Dynamics trajectories. The identification of macrostates for the MSM is a central decision that impacts the quality of the MSM but depends on both the selected representation of a structure and the clustering algorithm utilized over the featurized structures. Motivated by a large molecular system in its free and bound state, this paper investigates two directions of research, further reducing the representation dimensionality in a non-parametric, data-driven manner and including more structures in the computation. Rigorous evaluation of the quality of obtained MSMs via various statistical tests in a comparative setting firmly shows that fewer dimensions and more structures result in a better MSM. Many interesting findings emerge from the best MSM, advancing our understanding of the relationship between antibody dynamics and antibody–antigen recognition.

Summary: Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP3 produced by PI3Ks. Excess PIP3 promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP2 and PIP3. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP3. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP3 allosterically unfolding the N-terminal PIP2 binding domain, positioning it favorably for the polybasic motif interaction with PIP2. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

Summary: Drug repurposing or repositioning is a technique whereby existing drugs are used to treat emerging and challenging diseases, including COVID-19. Drug repurposing has become a promising approach because of the opportunity for reduced development timelines and overall costs. In the big data era, artificial intelligence (AI) and network medicine offer cutting-edge application of information science to defining disease, medicine, therapeutics, and identifying targets with the least error. In this Review, we introduce guidelines on how to use AI for accelerating drug repurposing or repositioning, for which AI approaches are not just formidable but are also necessary. We discuss how to use AI models in precision medicine, and as an example, how AI models can accelerate COVID-19 drug repurposing. Rapidly developing, powerful, and innovative AI and network medicine technologies can expedite therapeutic development. This Review provides a strong rationale for using AI-based assistive tools for drug repurposing medications for human disease, including during the COVID-19 pandemic.

K-Ras is one of the most frequently mutated oncogenes in human tumor cells. It consists of a well-conserved globular catalytic domain and a flexible tail-like hypervariable region (HVR) at its C-terminal end. It plays a key role in signaling networks in proliferation, differentiation, and survival, undergoing a conformational switch between the active and inactive states. It is regulated through the GDP-GTP cycle of the inactive GDP-bound and active GTP-bound states. Here, without imposing any prior constraints, we mapped the interaction pattern between the catalytic domain and the HVR using Molecular Dynamics with excited Normal Modes (MDeNM) starting from an initially extended HVR conformation for both states. Our sampling captured similar interaction patterns in both GDP- and GTP-bound states with shifted populations depending on the bound nucleotide. In the GDP-bound state, the conformations where the HVR interacts with the effector lobe are more populated than in the GTP-bound state, forming a buried thus autoinhibited catalytic site; in the GTP-bound state conformations where the HVR interacts with the allosteric lobe are more populated, overlapping the α3/α4 dimerization interface. The interaction of the GTP with Switch I and Switch II is stronger than that of the GDP in line with a decrease in the fluctuation upon GTP binding.

Tumor-specific genomic alterations allow systematic identification of genetic interactions that promote tumorigenesis and tumor vulnerabilities, offering novel strategies for development of targeted therapies for individual patients. We develop an Individualized Network-based Co-Mutation (INCM) methodology by inspecting over 2.5 million nonsynonymous somatic mutations derived from 6,789 tumor exomes across 14 cancer types from The Cancer Genome Atlas. Our INCM analysis reveals a higher genetic interaction burden on the significantly mutated genes, experimentally validated cancer genes, chromosome regulatory factors, and DNA damage repair genes, as compared to human pan-cancer essential genes identified by CRISPR-Cas9 screenings on 324 cancer cell lines. We find that genes involved in the cancer type-specific genetic subnetworks identified by INCM are significantly enriched in established cancer pathways, and the INCM-inferred putative genetic interactions are correlated with patient survival. By analyzing drug pharmacogenomics profiles from the Genomics of Drug Sensitivity in Cancer database, we show that the network-predicted putative genetic interactions (e.g., BRCA2-TP53) are significantly correlated with sensitivity/resistance of multiple therapeutic agents. We experimentally validated that afatinib has the strongest cytotoxic activity on BT474 (IC50 = 55.5 nM, BRCA2 and TP53 co-mutant) compared to MCF7 (IC50 = 7.7 μM, both BRCA2 and TP53 wild type) and MDA-MB-231 (IC50 = 7.9 μM, BRCA2 wild type but TP53 mutant). Finally, drug-target network analysis reveals several potential druggable genetic interactions by targeting tumor vulnerabilities. This study offers a powerful network-based methodology for identification of candidate therapeutic pathways that target tumor vulnerabilities and prioritization of potential pharmacogenomics biomarkers for development of personalized cancer medicine.

Abstract How cytochrome C is released from the mitochondria to the cytosol via Bax oligomeric pores, a process which is required for apoptosis, is still a mystery. Based on experimentally measured residue-residue distances, we recently solved the first atomic model for Bax oligomeric pores at the membranes using computational approaches. Here, we investigate the mechanism at the microsecond time- and nanometer space- scale using MD simulations. Our free energy landscape depicts a low barrier for the permeation of cytochrome C into the Bax C-terminal mouth, with the pathway proceeding to the inner cavity and exiting via the N-terminal mouth. Release is guided by organized charged/hydrophilic surfaces. The hydrophilicity and negative charge of the pore surface gradually increase along the release pathway from the pore entry to the exit opening. Rather than inert passing of the cytochrome C through a rigid pore, the flexible pore may selectively aid the cytochrome C passage. Once the Bax pore is formed in the membrane, with a low energy barrier, the release of cytochrome C may be readily achieved through energy fluctuations. Collectively, our work provides mechanistic insight in atomic detail into the release of cytochrome C through Bax oligomeric pores.

Abstract DNA shape is emerging as an important determinant of transcription factor binding beyond just the DNA sequence. The only tool for large scale DNA shape estimates, DNAshape was derived from Monte-Carlo simulations and predicts four broad and static DNA shape features, Propeller twist, Helical twist, Minor groove width and Roll. The contributions of other shape features e.g. Shift, Slide and Opening cannot be evaluated using DNAshape. Here, we report a novel method DynaSeq, which predicts molecular dynamics-derived ensembles of a more exhaustive set of DNA shape features. We compared the DNAshape and DynaSeq predictions for the common features and applied both to predict the genome-wide binding sites of 1312 TFs available from protein interaction quantification (PIQ) data. The results indicate a good agreement between the two methods for the common shape features and point to advantages in using DynaSeq. Predictive models employing ensembles from individual conformational parameters revealed that base-pair opening - known to be important in strand separation - was the best predictor of transcription factor-binding sites (TFBS) followed by features employed by DNAshape. Of note, TFBS could be predicted not only from the features at the target motif sites, but also from those as far as 200 nucleotides away from the motif.

The mechanism through which oncogenic Ras activates its effectors is vastly important to resolve. If allostery is at play, then targeting allosteric pathways could help in quelling activation of MAPK (Raf/MEK/ERK) and PI3K (PI3K/Akt/mTOR) cell proliferation pathways. On the face of it, allosteric activation is reasonable: Ras binding perturbs the conformational ensembles of its effectors. Here, however, we suggest that at least for Raf, PI3K, and NORE1A (RASSF5), that is unlikely. Raf's long disordered linker dampens effective allosteric activation. Instead, we suggest that the high-affinity Ras–Raf binding relieves Raf's autoinhibition, shifting Raf's ensemble from the inactive to the nanocluster-mediated dimerized active state, as Ras also does for NORE1A. PI3K is recruited and allosterically activated by RTK (e.g., EGFR) at the membrane. Ras restrains PI3K's distribution and active site orientation. It stabilizes and facilitates PIP2 binding at the active site and increases the PI3K residence time at the membrane. Thus, RTKs allosterically activate PI3Kα; however, merging their action with Ras accomplishes full activation. Here we review their activation mechanisms in this light and draw attention to implications for their pharmacology.

Oncoviruses rewire host pathways to subvert host immunity and promote their survival and proliferation. However, exactly how is challenging to understand. Here, by employing the first and to date only interface-based host-microbe interaction (HMI) prediction method, we explore a pivotal strategy oncoviruses use to drive cancer: mimicking binding surfaces—interfaces—of human proteins. We show that oncoviruses can target key human network proteins and transform cells by acquisition of cancer hallmarks. Experimental large-scale mapping of HMIs is difficult and individual HMIs do not permit in-depth grasp of tumorigenic virulence mechanisms. Our computational approach is tractable and 3D structural HMI models can help elucidate pathogenesis mechanisms and facilitate drug design. We observe that many host proteins are unique targets for certain oncoviruses, whereas others are common to several, suggesting similar infectious strategies. A rough estimation of our false discovery rate based on the tissue expression of oncovirus-targeted human proteins is 25%.

Classically, phenotype is what is observed, and genotype is the genetic makeup. Statistical studies aim to project phenotypic likelihoods of genotypic patterns. The traditional genotype-to-phenotype theory embraces the view that the encoded protein shape together with gene expression level largely determines the resulting phenotypic trait. Here, we point out that the molecular biology revolution at the turn of the century explained that the gene encodes not one but ensembles of conformations, which in turn spell all possible gene-associated phenotypes. The significance of a dynamic ensemble view is in understanding the linkage between genetic change and the gained observable physical or biochemical characteristics. Thus, despite the transformative shift in our understanding of the basis of protein structure and function, the literature still commonly relates to the classical genotype-phenotype paradigm. This is important because an ensemble view clarifies how even seemingly small genetic alterations can lead to pleiotropic traits in adaptive evolution and in disease, why cellular pathways can be modified in monogenic and polygenic traits, and how the environment may tweak protein function.

[This corrects the article DOI: 10.1371/journal.pcbi.1006658.].

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor's genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses-all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

Recent advances in next-generation sequencing and computational technologies have enabled routine analysis of large-scale single-cell ribonucleic acid sequencing (scRNA-seq) data. However, scRNA-seq technologies have suffered from several technical challenges, including low mean expression levels in most genes and higher frequencies of missing data than bulk population sequencing technologies. Identifying functional gene sets and their regulatory networks that link specific cell types to human diseases and therapeutics from scRNA-seq profiles are daunting tasks. In this study, we developed a Component Overlapping Attribute Clustering (COAC) algorithm to perform the localized (cell subpopulation) gene co-expression network analysis from large-scale scRNA-seq profiles. Gene subnetworks that represent specific gene co-expression patterns are inferred from the components of a decomposed matrix of scRNA-seq profiles. We showed that single-cell gene subnetworks identified by COAC from multiple time points within cell phases can be used for cell type identification with high accuracy (83%). In addition, COAC-inferred subnetworks from melanoma patients' scRNA-seq profiles are highly correlated with survival rate from The Cancer Genome Atlas (TCGA). Moreover, the localized gene subnetworks identified by COAC from individual patients' scRNA-seq data can be used as pharmacogenomics biomarkers to predict drug responses (The area under the receiver operating characteristic curves ranges from 0.728 to 0.783) in cancer cell lines from the Genomics of Drug Sensitivity in Cancer (GDSC) database. In summary, COAC offers a powerful tool to identify potential network-based diagnostic and pharmacogenomics biomarkers from large-scale scRNA-seq profiles. COAC is freely available at https://github.com/ChengF-Lab/COAC.

Acute myeloid leukemia is a disorder characterized by abnormal differentiation of myeloid cells and a clonal proliferation derived from primitive hematopoietic stem cells. Interventions that overcome myeloid differentiation have been shown to be a promising therapeutic strategy for acute myeloid leukemia. In this study, we demonstrate that CRISPR/Cas9-mediated knockout of dihydroorotate dehydrogenase leads to apoptosis and normal differentiation of acute myeloid leukemia cells, indicating that dihydroorotate dehydrogenase is a potential differentiation regulator and a therapeutic target in acute myeloid leukemia. By screening a library of natural products, we identified a novel dihydroorotate dehydrogenase inhibitor, isobavachalcone, derived from the traditional Chinese medicine Psoralea corylifolia. Using enzymatic analysis, thermal shift assay, pull down, nuclear magnetic resonance, and isothermal titration calorimetry experiments, we demonstrate that isobavachalcone inhibits human dihydroorotate dehydrogenase directly, and triggers apoptosis and differentiation of acute myeloid leukemia cells. Oral administration of isobavachalcone suppresses subcutaneous HL60 xenograft tumor growth without obvious toxicity. Importantly, our results suggest that a combination of isobavachalcone and adriamycin prolonged survival in an intravenous HL60 leukemia model. In summary, this study demonstrates that isobavachalcone triggers apoptosis and differentiation of acute myeloid leukemia cells via pharmacological inhibition of human dihydroorotate dehydrogenase, offering a potential therapeutic strategy for acute myeloid leukemia.

Hundreds of different species colonize multicellular organisms making them "metaorganisms". A growing body of data supports the role of microbiota in health and in disease. Grasping the principles of host-microbiota interactions (HMIs) at the molecular level is important since it may provide insights into the mechanisms of infections. The crosstalk between the host and the microbiota may help resolve puzzling questions such as how a microorganism can contribute to both health and disease. Integrated superorganism networks that consider host and microbiota as a whole-may uncover their code, clarifying perhaps the most fundamental question: how they modulate immune surveillance. Within this framework, structural HMI networks can uniquely identify potential microbial effectors that target distinct host nodes or interfere with endogenous host interactions, as well as how mutations on either host or microbial proteins affect the interaction. Furthermore, structural HMIs can help identify master host cell regulator nodes and modules whose tweaking by the microbes promote aberrant activity. Collectively, these data can delineate pathogenic mechanisms and thereby help maximize beneficial therapeutics. To date, challenges in experimental techniques limit large-scale characterization of HMIs. Here we highlight an area in its infancy which we believe will increasingly engage the computational community: predicting interactions across kingdoms, and mapping these on the host cellular networks to figure out how commensal and pathogenic microbiota modulate the host signaling and broadly cross-species consequences.

p53 protein has about thirty phosphorylation sites located at the N- and C-termini and in the core domain. The phosphorylation sites are relatively less mutated than other residues in p53. To understand why and how p53 phosphorylation sites are rarely mutated in human cancer, using a bioinformatics approaches, we examined the phosphorylation site and its nearby flanking residues, focusing on the consensus phosphorylation motif pattern, amino-acid correlations within the phosphorylation motifs, the propensity of structural disorder of the phosphorylation motifs, and cancer mutations observed within the phosphorylation motifs. Many p53 phosphorylation sites are targets for several kinases. The phosphorylation sites match 17 consensus sequence motifs out of the 29 classified. In addition to proline, which is common in kinase specificity-determining sites, we found high propensity of acidic residues to be adjacent to phosphorylation sites. Analysis of human cancer mutations in the phosphorylation motifs revealed that motifs with adjacent acidic residues generally have fewer mutations, in contrast to phosphorylation sites near proline residues. p53 phosphorylation motifs are mostly disordered. However, human cancer mutations within phosphorylation motifs tend to decrease the disorder propensity. Our results suggest that combination of acidic residues Asp and Glu with phosphorylation sites provide charge redundancy which may safe guard against loss-of-function mutations, and that the natively disordered nature of p53 phosphorylation motifs may help reduce mutational damage. Our results further suggest that engineering acidic amino acids adjacent to potential phosphorylation sites could be a p53 gene therapy strategy.

Computational biology has made powerful advances. Among these, trends in human health have been uncovered through heterogeneous ‘big data’ integration, and disease-associated genes were identified and classified. Along a different front, the dynamic organization of chromatin is being elucidated to gain insight into the fundamental question of genome regulation. Powerful conformational sampling methods have also been developed to yield a detailed molecular view of cellular processes. when combining these methods with the advancements in the modeling of supramolecular assemblies, including those at the membrane, we are finally able to get a glimpse into how cells’ actions are regulated. Perhaps most intriguingly, a major thrust is on to decipher the mystery of how the brain is coded. Here, we aim to provide a broad, yet concise, sketch of modern aspects of computational biology, with a special focus on computational structural biology. We attempt to forecast the areas that computational structural biology will embrace in the future and the challenges that it may face. We skirt details, highlight successes, note failures, and map directions.

The concept of allostery has evolved in the past century. In this Editorial, we briefly overview the history of allostery, from the pre-allostery nomenclature era starting with the Bohr effect (1904) to the birth of allostery by Monod and Jacob (1961). We describe the evolution of the allostery concept, from a conformational change in a two-state model (1965, 1966) to dynamic allostery in the ensemble model (1999); from multi-subunit (1965) proteins to all proteins (2004). We highlight the current available methods to study allostery and their applications in studies of conformational mechanisms, disease, and allosteric drug discovery. We outline the challenges and future directions that we foresee. Altogether, this Editorial narrates the history of this fundamental concept in the life sciences, its significance, methodologies to detect and predict it, and its application in a broad range of living systems.

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

Aβ is the toxic amyloid polypeptide responsible for Alzheimer’s disease (AD). Prevention and elimination of the Aβ misfolded aggregates are the promising therapeutic strategies for the AD treatments. Gammabody, the Aβ-Specific Single-domain (VH) antibody, recognizes Aβ aggregates with high affinity and specificity and reduces their toxicities. Employing the molecular dynamics simulations, we studied diverse gammabody-Aβ recognition complexes to get insights into their structural and dynamic properties and gammabody-Aβ recognitions. Among many heterogeneous binding modes, we focused on two gammabody-Aβ recognition scenarios: recognition through Aβ β-sheet backbone and on sidechain surface. We found that the gammabody primarily uses the complementarity-determining region 3 (CDR3) loop with the grafted Aβ sequence to interact with the Aβ fibril, while CDR1/CDR2 loops have very little contact. The gammabody-Aβ complexes with backbone binding mode are more stable, explaining the gammabody’s specificity towards the C-terminal Aβ sequence.