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

Discovery and Validation of Lmj_04_BRCT Domain, a Novel Therapeutic Target: Identification of Candidate Drugs for Leishmaniasis

Since many of the currently available antileishmanial treatments exhibit toxicity, low effectiveness, and resistance, search and validation of new therapeutic targets allowing the development of innovative drugs have become a worldwide priority. This work presents a structure-based drug discovery strategy to validate the Lmj_04_BRCT domain as a novel therapeutic target in Leishmania spp. The structure of this domain was explored using homology modeling, virtual screening, and molecular dynamics studies. Candidate compounds were validated in vitro using promastigotes of Leishmania major, L. amazonensis, and L. infantum, as well as primary mouse macrophages infected with L. major. The novel inhibitor CPE2 emerged as the most active of a group of compounds against Leishmania, being able to significantly reduce the viability of promastigotes. CPE2 was also active against the intracellular forms of the parasites and significantly reduced parasite burden in murine macrophages without exhibiting toxicity in host cells. Furthermore, L. major promastigotes treated with CPE2 showed significant lower expression levels of several genes (α-tubulin, Cyclin CYCA, and Yip1) related to proliferation and treatment resistance. Our in silico and in vitro studies suggest that the Lmj_04_BRCT domain and its here disclosed inhibitors are new potential therapeutic options against leishmaniasis.

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.

Promiscuous proteins are commonly observed in biological systems, such as modular domains that recognize phosphopeptides during signal transduction. The use of phosphopeptides and compounds with phosphate groups as inhibitors to protein–protein interactions have attracted increasing interest for years. By using atomistic molecular dynamics simulations, we are able to perform detailed analyses of the dihedral space to explore protein fluctuation upon ligand binding to better understand promiscuous molecular recognition. Free energy calculation can further provide insights into the mechanism of binding, including both enthalpic and entropic contributions for molecular recognition, which assist in inhibitor design. Our calculation results show that pre-rigidifying a ligand is not always advantageous, suggesting the challenge in retaining optimized intermolecular interactions in pre-rigidified ligand. Instead, certain flexible ligands with multiple binding conformations can reduce entropic penalty, and therefore improves binding affinity. According to our computations, we can introduce new intermolecular interactions to flexible ligand to strengthen attractions while maintaining small entropic penalty by retaining its plasticity in the bound conformation. The study might cast light on a new general strategy for designing inhibitors targeting promiscuous modular domains and protein–protein interactions.

A Thermodynamic Model for Genome Packaging in Hepatitis B Virus

Understanding the fundamentals of genome packaging in viral capsids is important for finding effective antiviral strategies and for utilizing benign viral particles for gene therapy. While the structure of encapsidated genomic materials has been routinely characterized with experimental techniques such as cryo-electron microscopy and x-ray diffraction, much less is known about the molecular driving forces underlying genome assembly in an intracellular environment and its in vivo interactions with the capsid proteins. Here we study the thermodynamic basis of the pregenomic RNA encapsidation in human Hepatitis B virus in vivo using a coarse-grained molecular model that captures the essential components of nonspecific intermolecular interactions. The thermodynamic model is used to examine how the electrostatic interaction between the packaged RNA and the highly charged C-terminal domains (CTD) of capsid proteins regulate the nucleocapsid formation. The theoretical model predicts optimal RNA content in Hepatitis B virus nucleocapsids with different CTD lengths in good agreement with mutagenesis measurements, confirming the predominant role of electrostatic interactions and molecular excluded-volume effects in genome packaging. We find that the amount of encapsidated RNA is not linearly correlated with the net charge of CTD tails as suggested by earlier theoretical studies. Our thermodynamic analysis of the nucleocapsid structure and stability indicates that ∼10% of the CTD residues are free from complexation with RNA, resulting in partially exposed CTD tails. The thermodynamic model also predicts the free energy of complex formation between macromolecules, which corroborates experimental results for the impact of CTD truncation on the nucleocapsid stability.

Achieving peptide binding specificity and promiscuity by loops: case of the forkhead-associated domain

The regulation of a series of cellular events requires specific protein–protein interactions, which are usually mediated by modular domains to precisely select a particular sequence from diverse partners. However, most signaling domains can bind to more than one peptide sequence. How do proteins create promiscuity from precision? Moreover, these complex interactions typically occur at the interface of a well-defined secondary structure, α helix and β sheet. However, the molecular recognition primarily controlled by loop architecture is not fully understood. To gain a deep understanding of binding selectivity and promiscuity by the conformation of loops, we chose the forkhead-associated (FHA) domain as our model system. The domain can bind to diverse peptides via various loops but only interact with sequences containing phosphothreonine (pThr). We applied molecular dynamics (MD) simulations for multiple free and bound FHA domains to study the changes in conformations and dynamics. Generally, FHA domains share a similar folding structure whereby the backbone holds the overall geometry and the variety of sidechain atoms of multiple loops creates a binding surface to target a specific partner. FHA domains determine the specificity of pThr by well-organized binding loops, which are rigid to define a phospho recognition site. The broad range of peptide recognition can be attributed to different arrangements of the loop interaction network. The moderate flexibility of the loop conformation can help access or exclude binding partners. Our work provides insights into molecular recognition in terms of binding specificity and promiscuity and helpful clues for further peptide design.