Transcription factor protein interactomes reveal genetic determinants in heart disease.

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica).

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here, we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (<20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (YD) (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

Integrating Protein Interaction Surface Prediction with a Fragment-Based Drug Design: Automatic Design of New Leads with Fragments on Energy Surfaces.

Protein-protein interactions (PPIs) have emerged in the past years as significant pharmacological targets in the development of new therapeutics due to their key roles in determining pathological pathways. Herein, we present fragments on energy surfaces, a simple and general design strategy that integrates the analysis of the dynamic and energetic signatures of proteins to unveil the substructures involved in PPIs, with docking, selection, and combination of drug-like fragments to generate new PPI inhibitor candidates. Specifically, structural representatives of the target protein are used as inputs for the blind physics-based prediction of potential protein interaction surfaces using the matrix of low coupling energy decomposition method. The predicted interaction surfaces are subdivided into overlapping windows that are used as templates to direct the docking and combination of fragments representative of moieties typically found in active drugs. This protocol is then applied and validated using structurally diverse, important PPI targets as test systems. We demonstrate that our approach facilitates the exploration of the molecular diversity space of potential ligands, with no requirement of prior information on the location and properties of interaction surfaces or on the structures of potential lead compounds. Importantly, the hit molecules that emerge from our ab initio design share high chemical similarity with experimentally tested active PPI inhibitors. We propose that the protocol we describe here represents a valuable means of generating initial leads against difficult targets for further development and refinement.

SARS-CoV-2 Spike Protein Mutations and Escape from Antibodies: A Computational Model of Epitope Loss in Variants of Concern.

The SARS-CoV-2 spike (S) protein is exposed on the viral surface and is the first point of contact between the virus and the host. For these reasons it represents the prime target for Covid-19 vaccines. In recent months, variants of this protein have started to emerge. Their ability to reduce or evade recognition by S-targeting antibodies poses a threat to immunological treatments and raises concerns for their consequences on vaccine efficacy. To develop a model able to predict the potential impact of S-protein mutations on antibody binding sites, we performed unbiased multi-microsecond molecular dynamics of several glycosylated S-protein variants and applied a straightforward structure-dynamics-energy based strategy to predict potential changes in immunogenic regions on each variant. We recover known epitopes on the reference D614G sequence. By comparing our results, obtained on isolated S-proteins in solution, to recently published data on antibody binding and reactivity in new S variants, we directly show that modifications in the S-protein consistently translate into the loss of potentially immunoreactive regions. Our findings can thus be qualitatively reconnected to the experimentally characterized decreased ability of some of the Abs elicited against the dominant S-sequence to recognize variants. While based on the study of SARS-CoV-2 spike variants, our computational epitope-prediction strategy is portable and could be applied to study immunoreactivity in mutants of proteins of interest whose structures have been characterized, helping the development/selection of vaccines and antibodies able to control emerging variants.

The calcium-binding type III repeats domain of thrombospondin-2 binds to fibroblast growth factor 2 (FGF2).

Thrombospondin (TSP)-1 and TSP-2 share similar structures and functions, including a remarkable antiangiogenic activity. We have previously demonstrated that a mechanism of the antiangiogenic activity of TSP-1 is the interaction of its type III repeats domain with fibroblast growth factor-2 (FGF2), affecting the growth factor bioavailability and angiogenic activity. Since the type III repeats domain is conserved in TSP-2, this study aimed at investigating whether also TSP-2 retained the ability to interact with FGF2. The FGF2 binding properties of TSP-1 and TSP-2 and their recombinant domains were analyzed by solid-phase binding and surface plasmon resonance assays. TSP-2 bound FGF2 with high affinity (Kd = 1.3 nM). TSP-2/FGF2 binding was inhibited by calcium and heparin. The FGF2-binding domain of TSP-2 was located in the type III repeats and the minimal interacting sequence was identified as the GVTDEKD peptide in repeat 3C, corresponding to KIPDDRD, the active sequence of TSP-1. A second putative FGF2 binding sequence was also identified in repeat 11C of both TSPs. Computational docking analysis predicted that both the TSP-2 and TSP-1-derived heptapeptides interacted with FGF2 with comparable binding properties. Accordingly, small molecules based on the TSP-1 active sequence blocked TSP-2/FGF2 interaction. Binding of TSP-2 to FGF2 impaired the growth factor ability to interact with its cellular receptors, since TSP-2-derived fragments prevented the binding of FGF2 to both heparin (used as a structural analog of heparan sulfate proteoglycans) and FGFR-1. These findings identify TSP-2 as a new FGF2 ligand that shares with TSP-1 the same molecular requirements for interaction with the growth factor and a comparable capacity to block FGF2 interaction with proangiogenic receptors. These features likely contribute to TSP-2 antiangiogenic and antineoplastic activity, providing the rationale for future therapeutic applications.

Impact of Mutations on NPAC Structural Dynamics: Mechanistic Insights from MD Simulations.

NPAC is a cytokine-like nuclear factor involved in chromatin modification and regulation of gene expression. In humans, the C-terminal domain of NPAC has the conserved structure of the ß-hydroxyacid dehydrogenases (ß-HAD) protein superfamily, which forms a stable tetrameric core scaffold for demethylase enzymes and organizes multiple sites for chromatin interactions. In spite of the close structural resemblance to other ß-HAD family members, the human NPAC dehydrogenase domain lacks a highly conserved catalytic lysine, substituted by a methionine. The reintroduction of the catalytic lysine by M437 K mutation results in a significant decrease of stability of the tetramer. Here, we have computationally investigated the molecular determinants of the functional differences between methionine and lysine-containing NPAC proteins. We find that the single mutation can determine strong consequences in terms of dynamics, stability, and ultimately ability to assemble in supramolecular complexes: the higher stability and lower flexibility of the methionine variant structurally preorganizes the monomer for tetramerization, whereas lysine increases flexibility and favors conformations that, while catalytically active, are not optimal for tetrameric assembly. We combine structure-dynamics analysis to an evolutionary study of NPAC sequences, showing that the methionine mutation occurs in a specifically flexible region of the lysine-containing protein, flanked by two domains that concentrate most of the stabilizing interactions. In our model, such separation of stability nuclei and flexible regions appears to favor the functional innovability of the protein.

Expedient Access to 2-Benzazepines by Palladium-Catalyzed C-H Activation: Identification of a Unique Hsp90 Inhibitor Scaffold.

Bioactive 2-benzazepines were accessed in an atom- and step-economical manner through a versatile palladium-catalyzed C-H activation strategy. The C-H arylation required low catalyst loading and a mild base, which was reflected by a broad scope and high functional-group tolerance. The benzotriazolodiazepinones were identified as new heat shock protein 90 (Hsp90) inhibiting lead compounds, with considerable potential for anti-cancer applications.

Design of Allosteric Stimulators of the Hsp90 ATPase as New Anticancer Leads.

Allosteric compounds that stimulate Hsp90 adenosine triphosphatase (ATPase) activity were rationally designed, showing anticancer potencies in the low micromolar to nanomolar range. In parallel, the mode of action of these compounds was clarified and a quantitative model that links the dynamic ligand-protein cross-talk to observed cellular and in vitro activities was developed. The results support the potential of using dynamics-based approaches to develop original mechanism-based cancer therapeutics.

Activation of Hsp90 Enzymatic Activity and Conformational Dynamics through Rationally Designed Allosteric Ligands.

Hsp90 is a molecular chaperone of pivotal importance for multiple cell pathways. ATP-regulated internal dynamics are critical for its function and current pharmacological approaches block the chaperone with ATP-competitive inhibitors. Herein, a general approach to perturb Hsp90 through design of new allosteric ligands aimed at modulating its functional dynamics is proposed. Based on the characterization of a first set of 2-phenylbenzofurans showing stimulatory effects on Hsp90 ATPase and conformational dynamics, new ligands were developed that activate Hsp90 by targeting an allosteric site, located 65 Šfrom the active site. Specifically, analysis of protein responses to first-generation activators was exploited to guide the design of novel derivatives with improved ability to stimulate ATP hydrolysis. The molecules' effects on Hsp90 enzymatic, conformational, co-chaperone and client-binding properties were characterized through biochemical, biophysical and cellular approaches. These designed probes act as allosteric activators of the chaperone and affect the viability of cancer cell lines for which proper functioning of Hsp90 is necessary.

Antitumor activity of a novel homodimeric SMAC mimetic in ovarian carcinoma.

Treatment of ovarian carcinoma often fails to be curative because of drug resistance, and many efforts are directed to overcome tumor cell resistance by increasing apoptosis induction. The potential of second mitochondria-derived activator of caspases (SMAC) mimetics (SMACm) has appeared in preclinical studies, but novel proapoptotic agents of this class with improved pharmacological profile are needed. To identify novel treatment options for ovarian carcinoma by interfering with antiapoptotic factors, in the present study a novel homodimeric SMACm (SM83) was employed in preclinical models both in vitro and in vivo. An investigation of the structural features of dimeric SM83 as compared to a closely related reference compound indicated slight differences, likely because of the interaction between one of the terminal phenyl groups and triazole rings of SM83 with the BIR2 domain. Although SM83 per se did not inhibit cell proliferation, it displayed a synergistic effect in combination with TNF-related apoptosis inducing ligand (TRAIL) in cell sensitivity assays. Because the tumor microenvironment is a reservoir of cytokines that may act in conjunction with SMACm to affect tumor growth, the activity of the novel compound was tested in vivo in ovarian carcinoma cells subcutaneously xenografted into immunodeficient mice. A significant tumor volume inhibition was observed together with activation of caspase 3 and apoptotic cell death. A biochemical analysis of tumor necrosis factor (TNF) and TRAIL content in specimens from xenografted mice indicated that SM83 downmodulated the levels of human TNF in plasma samples and tended to upmodulate human TRAIL levels in tumors. Thus, TRAIL appears to contribute to the antitumor activity of novel SMACm SM83 in subcutaneously grown ovarian carcinoma. Overall, our results indicate that SM83 is an attractive candidate for further development.

Exploiting conformational dynamics in drug discovery: design of C-terminal inhibitors of Hsp90 with improved activities.

The interaction that occurs between molecules is a dynamic process that impacts both structural and conformational properties of the ligand and the ligand binding site. Herein, we investigate the dynamic cross-talk between a protein and the ligand as a source for new opportunities in ligand design. Analysis of the formation/disappearance of protein pockets produced in response to a first-generation inhibitor assisted in the identification of functional groups that could be introduced onto scaffolds to facilitate optimal binding, which allowed for increased binding with previously uncharacterized regions. MD simulations were used to elucidate primary changes that occur in the Hsp90 C-terminal binding pocket in the presence of first-generation ligands. This data was then used to design ligands that adapt to these receptor conformations, which provides access to an energy landscape that is not visible in a static model. The newly synthesized compounds demonstrated antiproliferative activity at ∼150 nM concentration. The method identified herein may be used to design chemical probes that provide additional information on structural variations of Hsp90 C-terminal binding site.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation.

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

N-Glycosylation as a Modulator of Protein Conformation and Assembly in Disease.

Glycosylation, a prevalent post-translational modification, plays a pivotal role in regulating intricate cellular processes by covalently attaching glycans to macromolecules. Dysregulated glycosylation is linked to a spectrum of diseases, encompassing cancer, neurodegenerative disorders, congenital disorders, infections, and inflammation. This review delves into the intricate interplay between glycosylation and protein conformation, with a specific focus on the profound impact of N-glycans on the selection of distinct protein conformations characterized by distinct interactomes-namely, protein assemblies-under normal and pathological conditions across various diseases. We begin by examining the spike protein of the SARS virus, illustrating how N-glycans regulate the infectivity of pathogenic agents. Subsequently, we utilize the prion protein and the chaperone glucose-regulated protein 94 as examples, exploring instances where N-glycosylation transforms physiological protein structures into disease-associated forms. Unraveling these connections provides valuable insights into potential therapeutic avenues and a deeper comprehension of the molecular intricacies that underlie disease conditions. This exploration of glycosylation's influence on protein conformation effectively bridges the gap between the glycome and disease, offering a comprehensive perspective on the therapeutic implications of targeting conformational mutants and their pathologic assemblies in various diseases. The goal is to unravel the nuances of these post-translational modifications, shedding light on how they contribute to the intricate interplay between protein conformation, assembly, and disease.

Structural determinants of cold activity and glucose tolerance of a family 1 glycoside hydrolase (GH1) from Antarctic Marinomonas sp. ef1.

Cold-active enzymes support life at low temperatures due to their ability to maintain high activity in the cold and can be useful in several biotechnological applications. Although information on the mechanisms of enzyme cold adaptation is still too limited to devise general rules, it appears that very diverse structural and functional changes are exploited in different protein families and within the same family. In this context, we studied the cold adaptation mechanism and the functional properties of a member of the glycoside hydrolase family 1 (GH1) from the Antarctic bacterium Marinomonas sp. ef1. This enzyme exhibits all typical functional hallmarks of cold adaptation, including high catalytic activity at 5 °C, broad substrate specificity, low thermal stability, and higher lability of the active site compared to the overall structure. Analysis of the here-reported crystal structure (1.8 Å resolution) and molecular dynamics simulations suggest that cold activity and thermolability may be due to a flexible region around the active site (residues 298-331), whereas the dynamic behavior of loops flanking the active site (residues 47-61 and 407-413) may favor enzyme-substrate interactions at the optimal temperature of catalysis (Topt) by tethering together protein regions lining the active site. Stapling of the N-terminus onto the surface of the ß-barrel is suggested to partly counterbalance protein flexibility, thus providing a stabilizing effect. The tolerance of the enzyme to glucose and galactose is accounted for by the presence of a "gatekeeping" hydrophobic residue (Leu178), located at the entrance of the active site.

Structures, dynamics, complexes, and functions: From classic computation to artificial intelligence.

Computational approaches can provide highly detailed insight into the molecular recognition processes that underlie drug binding, the assembly of protein complexes, and the regulation of biological functional processes. Classical simulation methods can bridge a wide range of length- and time-scales typically involved in such processes. Lately, automated learning and artificial intelligence methods have shown the potential to expand the reach of physics-based approaches, ushering in the possibility to model and even design complex protein architectures. The synergy between atomistic simulations and AI methods is an emerging frontier with a huge potential for advances in structural biology. Herein, we explore various examples and frameworks for these approaches, providing select instances and applications that illustrate their impact on fundamental biomolecular problems.

How aberrant N-glycosylation can alter protein functionality and ligand binding: An atomistic view.

Protein-assembly defects due to an enrichment of aberrant conformational protein variants are emerging as a new frontier in therapeutics design. Understanding the structural elements that rewire the conformational dynamics of proteins and pathologically perturb functionally oriented ensembles is important for inhibitor development. Chaperones are hub proteins for the assembly of multiprotein complexes and an enrichment of aberrant conformers can affect the cellular proteome, and in turn, phenotypes. Here, we integrate computational and experimental tools to investigte how N-glycosylation of specific residues in glucose-regulated protein 94 (GRP94) modulates internal dynamics and alters the conformational fitness of regions fundamental for the interaction with ATP and synthetic ligands and impacts substructures important for the recognition of interacting proteins. N-glycosylation plays an active role in modulating the energy landscape of GRP94, and we provide support for leveraging the knowledge on distinct glycosylation variants to design molecules targeting GRP94 disease-associated conformational states and assemblies.

Conformational Behavior of SARS-Cov-2 Spike Protein Variants: Evolutionary Jumps in Sequence Reverberate in Structural Dynamic Differences.

SARS-CoV-2 has evolved rapidly in the first 3 years of pandemic diffusion. The initial evolution of the virus appeared to proceed through big jumps in sequence changes rather than through the stepwise accumulation of point mutations on already established variants. Here, we examine whether this nonlinear mutational process reverberates in variations of the conformational dynamics of the SARS-CoV-2 Spike protein (S-protein), the first point of contact between the virus and the human host. We run extensive microsecond-scale molecular dynamics simulations of seven distinct variants of the protein in their fully glycosylated state and set out to elucidate possible links between the mutational spectrum of the S-protein and the structural dynamics of the respective variant, at global and local levels. The results reveal that mutation-dependent structural and dynamic modulations mostly consist of increased coordinated motions in variants that acquire stability and in an increased internal flexibility in variants that are less stable. Importantly, a limited number of functionally important substructures (the receptor binding domain, in particular) share the same time of movements in all variants, indicating efficient preorganization for functional regions dedicated to host interactions. Our results support a model in which the internal dynamics of the S-proteins from different strains varies in a way that reflects the observed random and non-stepwise jumps in sequence evolution, while conserving the functionally oriented traits of conformational dynamics necessary to support productive interactions with host receptors.

A NMR and computational study of Smac mimics targeting both the BIR2 and BIR3 domains in XIAP protein.

In this paper we report an extensive NMR analysis of small ligands (Smac mimics) complexed with different constructs of XIAP. The mimics-binding site of XIAP is known as the BIR3 domain - primary, and the linker BIR2 region - secondary site. Interactions between the BIR3 domain and Smac mimics have been extensively studied by X-ray but, as of today, there are scarce data about the interaction between BIR2, or the whole linker-BIR2-BIR3 construct, and Smac mimics. In order to characterize our Smac mimics, we performed a STD NMR study between our 4-substituted, 1-aza-2-oxobicyclo[5.3.0]decane scaffold-based molecules and three different XIAP fragments: single BIR2 and BIR3 domains, and bifunctional linker-BIR2-BIR3. The results were integrated with docking calculations and molecular dynamics simulations. NMR data, which are consistent with biological tests, indicated that the two BIR subunits interact differently with our Smac mimics and suggest that the ligands enter into more intimate contact with the linker-BIR2-BIR3. In conclusion, we observe that the SMAC mimics showed with the construct linker-BIR2-BIR3 a series of NOE contacts that were not observed in the mono-domain ligand:BIR2 or :BIR3 complexes. So, in agreement with the computational models we believe that the linker moieties of the binding site play a key role in the stability of the protein complex.

Dimeric Smac mimetics/IAP inhibitors as in vivo-active pro-apoptotic agents. Part II: Structural and biological characterization.

Novel pro-apoptotic, homodimeric and heterodimeric Smac mimetics/IAPs inhibitors connected through head-head (8), tail-tail (9) or head-tail linkers (10), were biologically and structurally characterized. In vitro characterization (binding to BIR3 and linker-BIR2-BIR3 domains from XIAP and cIAP1, cytotoxicity assays) identified early leads from each dimer family. Computational models and structural studies (crystallography, NMR, gel filtration) partially rationalized the observed properties for each dimer class. Tail-tail dimer 9a was shown to be active in a breast and in an ovary tumor model, highlighting the potential of dimeric Smac mimetics/IAP inhibitors based on the N-AVPI-like 4-substituted 1-aza-2-oxobicyclo[5.3.0]decane scaffold as potential antineoplastic agents.

Studying the Dynamics of a Complex G-Quadruplex System: Insights into the Comparison of MD and NMR Data.

Molecular dynamics (MD) simulations are coming of age in the study of nucleic acids, including specific tertiary structures such as G-quadruplexes. While being precious for providing structural and dynamic information inaccessible to experiments at the atomistic level of resolution, MD simulations in this field may still be limited by several factors. These include the force fields used, different models for ion parameters, ionic strengths, and water models. We address various aspects of this problem by analyzing and comparing microsecond-long atomistic simulations of the G-quadruplex structure formed by the human immunodeficiency virus long terminal repeat (HIV LTR)-III sequence for which nuclear magnetic resonance (NMR) structures are available. The system is studied in different conditions, systematically varying the ionic strengths, ion numbers, and water models. We comparatively analyze the dynamic behavior of the G-quadruplex motif in various conditions and assess the ability of each simulation to satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints and structural parameters. The conditions taking into account K+-ions to neutralize the system charge, mimicking the intracellular ionic strength, and using the four-atom water model are found to be the best in reproducing the experimental NMR constraints and data. Our analysis also reveals that in all of the simulated environments residues belonging to the duplex moiety of HIV LTR-III exhibit the highest flexibility.