Discovery of a hidden transient state in all bromodomain families

ARIP: A Tool for Precise Interatomic Contact Area and Volume Calculation in Proteins

The interplay patterns of amino acid residues are pivotal in determining the tertiary structure and flexibility of proteins, which in turn are intricately linked to their functionality and interactions with other molecules. Here, we introduce ARIP, a novel tool designed to identify contact residues within proteins. ARIP employs a modified version of the dr_sasa algorithm and an atomic overlap weighted algorithm to directly calculate the contact area and volume between atoms based on their van der Waals radius. It also allows for the selection of solvent radii, recognizing that not every atom in proteins can interact with water molecules. The solvent parameters were derived from the analysis of approximately 5000 protein and nucleic acid structures with water molecules determined using X-ray crystallography. One advantage of the modified algorithm is its capability to analyze multiple models within a single PDB file, making it suitable for molecular dynamic capture. The contact volume is symmetrically distributed between the interacting atoms, providing more informative results than contact area for the analysis of intra- and intermolecular interactions and the development of scoring functions. Furthermore, ARIP has been applied to four distinct cases: capturing key residue-residue contacts in NMR structures of P4HB, protein-drug binding of CYP17A1, protein-DNA binding of SPI1, and molecular dynamic simulations of BRD4.

Toward physics-based precision medicine: Exploiting protein dynamics to design new therapeutics and interpret variants

The goal of precision medicine is to utilize our knowledge of the molecular causes of disease to better diagnose and treat patients. However, there is a substantial mismatch between the small number of food and drug administration (FDA)-approved drugs and annotated coding variants compared to the needs of precision medicine. This review introduces the concept of physics-based precision medicine, a scalable framework that promises to improve our understanding of sequence-function relationships and accelerate drug discovery. We show that accounting for the ensemble of structures a protein adopts in solution with computer simulations overcomes many of the limitations imposed by assuming a single protein structure. We highlight studies of protein dynamics and recent methods for the analysis of structural ensembles. These studies demonstrate that differences in conformational distributions predict functional differences within protein families and between variants. Thanks to new computational tools that are providing unprecedented access to protein structural ensembles, this insight may enable accurate predictions of variant pathogenicity for entire libraries of variants. We further show that explicitly accounting for protein ensembles, with methods like alchemical free energy calculations or docking to Markov state models, can uncover novel lead compounds. To conclude, we demonstrate that cryptic pockets, or cavities absent in experimental structures, provide an avenue to target proteins that are currently considered undruggable. Taken together, our review provides a roadmap for the field of protein science to accelerate precision medicine.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Biochemical and cellular insights into the Baz2B protein, a non-catalytic subunit of the chromatin remodeling complex

Baz2B is a regulatory subunit of the ATP-dependent chromatin remodeling complexes BRF1 and BRF5, which control access to DNA during DNA-templated processes. Baz2B has been implicated in several diseases and also in unhealthy ageing, however limited information is available on the domains and cellular roles of Baz2B. To gain more insight into the Baz2B function, we biochemically characterized the TAM (Tip5/ARBP/MBD) domain with the auxiliary AT-hook motifs and the bromodomain (BRD). We observed alterations in histone code recognition in bromodomains carrying cancer-associated point mutations, suggesting their potential involvement in disease. Furthermore, the depletion of Baz2B in the Hap1 cell line resulted in altered cell morphology, reduced colony formation and perturbed transcriptional profiles. Despite that, super-resolution microscopy images revealed no changes in the overall chromatin structure in the absence of Baz2B. These findings provide insights into the biological function of Baz2B.

Selectivity Mechanism of Pyrrolopyridone Analogues Targeting Bromodomain 2 of Bromodomain-Containing Protein 4 from Molecular Dynamics Simulations

Bromodomain and extra-terminal domain (BET) proteins play an important role in epigenetic regulation and are linked to several diseases; therefore, they are interesting targets. BET has two bromodomains: bromodomain 1 (BD1) and BD2. Selective targeting of BD1 or BD2 may produce different activities and greater effects than pan-BD inhibitors. However, the selective mechanism of the specific core must be studied at the atomic level. This study determined the effectiveness of pyrrolopyridone analogues to selectively inhibit BD2 using a pan-BD inhibitor (ABBV-075) and a selective-BD2 inhibitor (ABBV-744). Molecular dynamics simulations and calculations of binding free energies were used to systematically study the selectivity of BD2 inhibition by the pyrrolopyridone analogues. Overall, the pyrrolopyridone analogue inhibitors targeting BD2 interacted mainly with the following amino acid pairs between bromodomain-containing protein 4 (BRD4)-BD1 and BRD4-BD2 complexes: I146/V439, N140/N433, D144/H437, P82/P375, V87/V380, D88/D381, and Y139/Y432. The pyrrolopyridone analogues targeting BRD4-BD2 were divided into five regions based on selectivity mechanism. These results suggest that the R3 and R5 regions of pyrrolopyridone analogues can be modified to improve the selectivity between BRD4-BD1 and BRD4-BD2. The selectivity of BD2 inhibition by pyrrolopyridone analogues can be used to design novel BD2 inhibitors based on a pyrrolopyridone core.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long-time scale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We do not, however, find evidence for the existence of canonical conformational selection or induced fit binding pathways. We observe four kinetically separated native-like bound states that interconvert on time scales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Predicting locations of cryptic pockets from single protein structures using the PocketMiner graph neural network

Cryptic pockets expand the scope of drug discovery by enabling targeting of proteins currently considered undruggable because they lack pockets in their ground state structures. However, identifying cryptic pockets is labor-intensive and slow. The ability to accurately and rapidly predict if and where cryptic pockets are likely to form from a structure would greatly accelerate the search for druggable pockets. Here, we present PocketMiner, a graph neural network trained to predict where pockets are likely to open in molecular dynamics simulations. Applying PocketMiner to single structures from a newly curated dataset of 39 experimentally confirmed cryptic pockets demonstrates that it accurately identifies cryptic pockets (ROC-AUC: 0.87) >1,000-fold faster than existing methods. We apply PocketMiner across the human proteome and show that predicted pockets open in simulations, suggesting that over half of proteins thought to lack pockets based on available structures likely contain cryptic pockets, vastly expanding the potentially druggable proteome.

Targeting protein conformations with small molecules to control protein complexes

Dynamic protein complexes function in all cellular processes, from signaling to transcription, using distinct conformations that regulate their activity. Conformational switching of proteins can turn on or off their activity through protein-protein interactions, catalytic function, cellular localization, or membrane interaction. Recent advances in structural, computational, and chemical methodologies have enabled the discovery of small-molecule activators and inhibitors of conformationally dynamic proteins by using a more rational design than a serendipitous screening approach. Here, we discuss such recent examples, focusing on the mechanism of protein conformational switching and its regulation by small molecules. We emphasize the rational approaches to control protein oligomerization with small molecules that offer exciting opportunities for investigation of novel biological mechanisms and drug discovery.

Bromodomain Interactions with Acetylated Histone 4 Peptides in the BRD4 Tandem Domain: Effects on Domain Dynamics and Internal Flexibility

The bromodomain and extra-terminal (BET) protein BRD4 regulates gene expression via recruitment of transcriptional regulatory complexes to acetylated chromatin. Like other BET proteins, BRD4 contains two bromodomains, BD1 and BD2, that can interact cooperatively with target proteins and designed ligands, with important implications for drug discovery. Here, we used nuclear magnetic resonance (NMR) spectroscopy to study the dynamics and interactions of the isolated bromodomains, as well as the tandem construct including both domains and the intervening linker, and investigated the effects of binding a tetra-acetylated peptide corresponding to the tail of histone 4. The peptide affinity is lower for both domains in the tandem construct than for the isolated domains. Using ¹⁵N spin relaxation, we determined the global rotational correlation times and residue-specific order parameters for BD1 and BD2. Isolated BD1 is monomeric in the apo state but apparently dimerizes upon binding the tetra-acetylated peptide. Isolated BD2 partially dimerizes in both the apo and peptide-bound states. The backbone order parameters reveal marked differences between BD1 and BD2, primarily in the acetyl-lysine binding site where the ZA loop is more flexible in BD2. Peptide binding reduces the order parameters of the ZA loop in BD1 and the ZA and BC loops in BD2. The AB loop, located distally from the binding site, shows variable dynamics that reflect the different dimerization propensities of the domains. These results provide a basis for understanding target recognition by BRD4.

GNE-064: A Potent, Selective, and Orally Bioavailable Chemical Probe for the Bromodomains of SMARCA2 and SMARCA4 and the Fifth Bromodomain of PBRM1

Bromodomains are acetyllysine recognition domains present in a variety of human proteins. Bromodomains also bind small molecules that compete with acetyllysine, and therefore bromodomains have been targets for drug discovery efforts. Highly potent and selective ligands with good cellular permeability have been proposed as chemical probes for use in exploring the functions of many of the bromodomain proteins. We report here the discovery of a class of such inhibitors targeting the family VIII bromodomains of SMARCA2 (BRM) and SMARCA4 (BRG1), and PBRM1 (polybromo-1) bromodomain 5. We propose one example from this series, GNE-064, as a chemical probe for the bromodomains SMARCA2, SMARCA4, and PBRM1(5) with the potential for in vivo use.

Crystal structure of the BAZ2B TAM domain

BAZ2B is a regulatory subunit of the ISWI (Imitation Switch) remodeling complex and engages in nucleosome remodeling. Loss-of-function and haploinsufficiency of BAZ2B are associated with different diseases. BAZ2B is a large multidomain protein. In addition to the epigenetic reader domains plant homeodomain (PHD) and bromodomain (BRD), BAZ2B also has a Tip5/ARBP/MBD (TAM) domain. Sequence alignment revealed that the TAM domains of BAZ2A and BAZ2B share 53% sequence identity. How the BAZ2A TAM domain bound with DNA has been characterized recently, however, the DNA binding ability and methylation preference, as well as the structural basis of the BAZ2B TAM domain are not studied yet. In this study, we measured the DNA binding affinity of the TAM domain of BAZ2B, and also determined its apo crystal structure. We found that the TAM domains of BAZ2A and BAZ2B adopt almost the same fold, and like BAZ2A, the BAZ2B TAM domain also binds to dsDNA without methyl-cytosine preference, implying that the BAZ2B TAM domain might recognize DNA in a similar binding mode to that of the BAZ2A TAM domain. These results provide clues for the biological function study of BAZ2B in the future.

A cryptic pocket in Ebola VP35 allosterically controls RNA binding

Protein-protein and protein-nucleic acid interactions are often considered difficult drug targets because the surfaces involved lack obvious druggable pockets. Cryptic pockets could present opportunities for targeting these interactions, but identifying and exploiting these pockets remains challenging. Here, we apply a general pipeline for identifying cryptic pockets to the interferon inhibitory domain (IID) of Ebola virus viral protein 35 (VP35). VP35 plays multiple essential roles in Ebola's replication cycle but lacks pockets that present obvious utility for drug design. Using adaptive sampling simulations and machine learning algorithms, we predict VP35 harbors a cryptic pocket that is allosterically coupled to a key dsRNA-binding interface. Thiol labeling experiments corroborate the predicted pocket and mutating the predicted allosteric network supports our model of allostery. Finally, covalent modifications that mimic drug binding allosterically disrupt dsRNA binding that is essential for immune evasion. Based on these results, we expect this pipeline will be applicable to other proteins.

Dynamics of the Histone Acetyltransferase Lysine-Rich Loop in the Catalytic Core of the CREB-Binding Protein

The tight control of transcriptional coactivators is a fundamental aspect of gene expression in cells. The regulation of the CREB-binding protein (CBP) and p300 coactivators, two paralog multidomain proteins, involves an autoinhibitory loop (AIL) of the histone acetyltransferase (HAT) domain. There is experimental evidence for the AIL engaging with the HAT binding site, thus interrupting the acetylation of histone tails or other proteins. Both CBP and p300 contain a domain of about 110 residues (called the bromodomain) that recognizes histone tails with one or more acetylated lysine side chains. Here, we investigate by molecular dynamics simulations whether the AIL of CBP (residues 1556-1618) acetylated at the side chain of Lys1595 can bind to the bromodomain. The structural instability and fast unbinding kinetics of the AIL from the bromodomain pocket suggest that the AIL is not a ligand of the bromodomain on the same protein chain. This is further supported by the absence of strong and persistent contacts at the binding interface. Furthermore, the simulations of unbinding show an initial fast detachment of the acetylated lysine and a slower phase necessary for complete AIL dissociation. We provide further evidence for the instability of the AIL intramolecular binding by comparison with a natural ligand, the histone peptide H3K56ac, which shows higher stability in the pocket.

Divide and Conquer. Pocket-Opening Mixed-Solvent Simulations in the Perspective of Docking Virtual Screening Applications for Drug Discovery

The existence of a druggable binding pocket is a prerequisite for computational drug-target interaction studies including virtual screening. Retrospective studies have shown that extended sampling methods like Markov State Modeling and mixed-solvent simulations can identify cryptic pockets relevant for drug discovery. Here, we apply a combination of mixed-solvent molecular dynamics (MD) and time-structure independent component analysis (TICA) to four retrospective case studies: NPC2, the CECR2 bromodomain, TEM-1, and MCL-1. We compare previous experimental and computational findings to our results. It is shown that the successful identification of cryptic pockets depends on the system and the cosolvent probes. We used alternative TICA internal features such as the unbiased backbone coordinates or backbone dihedrals versus biased interatomic distances. We found that in the case of NPC2, TEM-1, and MCL-1, the use of unbiased features is able to identify cryptic pockets, although in the case of the CECR2 bromodomain, more specific features are required to properly capture a pocket opening. In the perspective of virtual screening applications, it is shown how docking studies with the parent ligands depend critically on the conformational state of the targets.

Differential BET Bromodomain Inhibition by Dihydropteridinone and Pyrimidodiazepinone Kinase Inhibitors

BRD4 and other members of the bromodomain and extraterminal (BET) family of proteins are promising epigenetic targets for the development of novel therapeutics. Among the reported BRD4 inhibitors are dihydropteridinones and benzopyrimidodiazepinones originally designed to target the kinases PLK1, ERK5, and LRRK2. While these kinase inhibitors were identified as BRD4 inhibitors, little is known about their binding potential and structural details of interaction with the other BET bromodomains. We comprehensively characterized a series of known and newly identified dual BRD4-kinase inhibitors against all eight individual BET bromodomains. A detailed analysis of 23 novel cocrystal structures of BET-kinase inhibitor complexes in combination with direct binding assays and cell signaling studies revealed significant differences in molecular shape complementarity and inhibitory potential. Collectively, the data offer new insights into the action of kinase inhibitors across BET bromodomains, which may aid the development of drugs to inhibit certain BET proteins and kinases differentially.

A distal regulatory strategy of enzymes: from local to global conformational dynamics

Modulating the distribution of various states in protein ensembles through distal sites may be promising in the evolution of enzymes in desired directions. However, the prediction of distal mutation hotspots that stabilize the favoured states from a computational perspective remains challenging. Here, we presented a strategy based on molecular dynamics (MD) and Markov state models (MSM) to predict distal mutation sites. Extensive MD combined with MSM was applied to determine the principally distributed metastable states interconverting at a slow timescale. Then, molecular docking was used to classify these states into active states and inactive ones. Distal mutation hotspots were targeted based on comparing the conformational features between active and inactive states, where mutations destabilize the inactive states and show little influence on the active state. The proposed strategy was used to explore the highly dynamic MHETase, which shows a potential application in the biodegradation of poly(ethylene terephthalate) (PET). Seven principally populated interrelated metastable states were identified, and the atomistic picture of their conformational changes was unveiled. Several residues at distal positions were found to adopt more H-bond occupancies in inactive states than active states, making them potential mutation hotspots for stabilizing the favoured conformations. In addition, the detailed mechanism revealed the significance of calcium ions at a distance from the catalytic centre in reshaping the free energy landscape. This study deepens the understanding of the conformational dynamics of α/β hydrolases containing a lid domain and advances the study of enzymatic plastic degradation.

Computational strategies for protein conformational ensemble detection

Protein function is constrained by the three-dimensional structure but is delineated by its dynamics. This framework must satisfy specificity of function along with adaptability to changing environments and evolvability under external constraints. The accessibility of the available regions of the energy landscape for a set of conditions and shifts in the populations upon their modulation have effects propagating across scales, from biomolecular interactions, to organisms, to populations. Developing the ability to detect and juggle protein conformations supplemented by a physics-based understanding has implications for not only in vivo problems but also for resistance impeding drug discovery and bionano-sensor design.

Design and synthesis of a new orthogonally protected glutamic acid analog and its use in the preparation of high affinity polo-like kinase 1 polo-box domain - binding peptide macrocycles

Targeting protein - protein interactions (PPIs) has emerged as an important area of discovery for anticancer therapeutic development. In the case of phospho-dependent PPIs, such as the polo-like kinase 1 (Plk1) polo-box domain (PBD), a phosphorylated protein residue can provide high-affinity recognition and binding to target protein hot spots. Developing antagonists of the Plk1 PBD can be particularly challenging if one relies solely on interactions within and proximal to the phospho-binding pocket. Fortunately, the affinity of phospho-dependent PPI antagonists can be significantly enhanced by taking advantage of interactions in both the phospho-binding site and hidden "cryptic" pockets that may be revealed on ligand binding. In our current paper, we describe the design and synthesis of macrocyclic peptide mimetics directed against the Plk1 PBD, which are characterized by a new glutamic acid analog that simultaneously serves as a ring-closing junction that provides accesses to a cryptic binding pocket, while at the same time achieving proper orientation of a phosphothreonine (pT) residue for optimal interaction in the signature phospho-binding pocket. Macrocycles prepared with this new amino acid analog introduce additional hydrogen-bonding interactions not found in the open-chain linear parent peptide. It is noteworthy that this new glutamic acid-based amino acid analog represents the first example of extremely high affinity ligands where access to the cryptic pocket from the pT-2 position is made possible with a residue that is not based on histidine. The concepts employed in the design and synthesis of these new macrocyclic peptide mimetics should be useful for further studies directed against the Plk1 PBD and potentially for ligands directed against other PPI targets.

Selectively targeting individual bromodomain: Drug discovery and molecular mechanisms

Bromodomain-containing proteins include bromodomain and extra-terminal (BET) and non-BET families. Due to the conserved bromodomain (BD) module between BD-containing proteins, and especially BETs with each member having two BDs (BD1 and BD2), the high degree of structural similarity makes BD-selective inhibitors much difficult to be designed. However, increasing evidences emphasized that individual BDs had distinct functions and different cellular phenotypes after pharmacological inhibition, and selectively targeting one of the BDs could result in a different efficacy and tolerability profile. This review is to summarize the pioneering progress of BD-selective inhibitors targeting BET and non-BET proteins, focusing on their structural features, biological activity, therapeutic application and experimental/theoretical mechanisms. The present proteolysis targeting chimeras (PROTAC) degraders targeting BDs, and clinical status of BD-selective inhibitors were also analyzed, providing a new insight into future direction of bromodomain-selective drug discovery.