Association study of leptin receptor polymorphisms in women with obesity and their impact on protein domains: a case-control study and in silico analyses.

Leptin receptor (LEPR) is a member of the class I cytokine receptor family that receives and transmits leptin signals. It is primarily involved in the regulation of energy expenditure and food intake. This study aimed to evaluate the association of LEPR gene polymorphisms, Lys109Arg, Gln223Arg and Lys656Asn, with obesity in Moroccan women and to explore the structural and functional consequences of these SNPs. The variants were genotyped using the Sanger sequencing method. The three-dimensional structures of LEPR extracellular domains were determined using a template-based tertiary structure modeling web server and the protein variants were generated using in silico mutagenesis. The amino acids conservation analysis in the variants region was performed based on a protein's evolutionary profile. The molecular dynamics simulations of the wild-types and variants N-terminal, cytokine receptor homology I and fibronectin type III domains of LEPR protein were performed to investigate their impact on the domain structures. We identified that only Lys656Asn polymorphism is associated with obesity in Moroccan women (P = 0.024). In silico analyses revealed that Lys109, Gln223 and Lys656 are exposed residues and their substitution leads to changes in protein structure through loss or gain of hydrogen bonds and hydrophobic interactions. Lys656Asn increases the stability and decreased flexibility of the fibronectin type III domain. Lys109Arg highly decreases the stability and increases flexibility and the overall dimension of N-terminal and cytokine receptor homology I domains. Gln223Arg increases the stability and the compaction level of these domains. These results provide insight into the involvement of LEPR variants in obesity development.Communicated by Ramaswamy H. Sarma.

Influenza A Virus NS1 Protein Structural Flexibility Analysis According to Its Structural Polymorphism Using Computational Approaches.

Influenza A viruses are highly contagious RNA viruses that cause respiratory tract infections in humans and animals. Their non-structural protein NS1, a homodimer of two 230-residue chains, is the main viral factor in counteracting the antiviral defenses of the host cell. Its RNA-binding domain is an obligate dimer that is connected to each of the two effector domains by a highly flexible unstructured linker region of ten amino acids. The flexibility of NS1 is a key property that allows its effector domains and its RNA binding domain to interact with several protein partners or RNAs. The three-dimensional structures of full-length NS1 dimers revealed that the effector domains could adopt three distinct conformations as regards their mutual interactions and their orientation relative to the RNA binding domain (closed, semi-open and open). The origin of this structural polymorphism is currently being investigated and several hypotheses are proposed, among which one posits that it is a strain-specific property. In the present study, we explored through computational molecular modeling the dynamic and flexibility properties of NS1 from three important influenza virus A strains belonging to three distinct subtypes (H1N1, H6N6, H5N1), for which at least one conformation is available in the Protein Data Bank. In order to verify whether NS1 is stable in three forms for the three strains, we constructed homology models if the corresponding forms were not available in the Protein Data Bank. Molecular dynamics simulations were performed in order to predict the stability over time of the three distinct sequence variants of NS1, in each of their three distinct conformations. Our results favor the co-existence of three stable structural forms, regardless of the strain, but also suggest that the length of the linker, along with the presence of specific amino acids, modulate the dynamic properties and the flexibility of NS1.

Druggable Pockets at the RNA Interface Region of Influenza A Virus NS1 Protein Are Conserved across Sequence Variants from Distinct Subtypes.

Influenza A viruses still represent a major health issue, for both humans and animals. One of the main viral proteins of interest to target is the NS1 protein, which counters the host immune response and promotes viral replication. NS1 is a homodimer composed of a dimeric RNA-binding domain (RBD), which is structurally stable and conserved in sequence, and two effector domains that are tethered to the RBD by linker regions. This linker flexibility leads to NS1 polymorphism and can therefore exhibit different forms. Previously, we identified a putative drug-binding site, located in the RBD interface in a crystal structure of NS1. This pocket could be targeted to block RNA binding and inhibit NS1 activities. The objective of the present study is to confirm the presence of this druggable site, whatever the sequence variants, in order to develop a universal therapeutic compound that is insensitive to sequence variations and structural flexibility. Using a set of four NS1 full-length structures, we combined different bioinformatics approaches such as pocket tracking along molecular dynamics simulations, druggability prediction and classification. This protocol successfully confirmed a frequent large binding-site that is highly druggable and shared by different NS1 forms, which is promising for developing a robust NS1-targeted therapy.

Impacts of drug resistance mutations on the structural asymmetry of the HIV-2 protease.

BACKGROUND: Drug resistance is a severe problem in HIV treatment. HIV protease is a common target for the design of new drugs for treating HIV infection. Previous studies have shown that the crystallographic structures of the HIV-2 protease (PR2) in bound and unbound forms exhibit structural asymmetry that is important for ligand recognition and binding. Here, we investigated the effects of resistance mutations on the structural asymmetry of PR2. Due to the lack of structural data on PR2 mutants, the 3D structures of 30 PR2 mutants of interest have been modeled using an in silico protocol. Structural asymmetry analysis was carried out with an in-house structural-alphabet-based approach. RESULTS: The systematic comparison of the asymmetry of the wild-type structure and a large number of mutants highlighted crucial residues for PR2 structure and function. In addition, our results revealed structural changes induced by PR2 flexibility or resistance mutations. The analysis of the highlighted structural changes showed that some mutations alter protein stability or inhibitor binding. CONCLUSIONS: This work consists of a structural analysis of the impact of a large number of PR2 resistant mutants based on modeled structures. It suggests three possible resistance mechanisms of PR2, in which structural changes induced by resistance mutations lead to modifications in the dimerization interface, ligand recognition or inhibitor binding.

DHH pathogenic variants involved in 46,XY disorders of sex development differentially impact protein self-cleavage and structural conformation.

In humans, pathogenic variants in the DHH gene underlie cases of 46,XY gonadal dysgenesis. DHH is part of the Hedgehog family of proteins, which require extensive processing, including self-cleavage of the precursor for efficient signalling. In our work, we have assessed the effect of several human DHH pathogenic variants involved in recessive complete or partial gonadal dysgenesis, on protein processing and sub-cellular localization. We found that a subset of variants was unable to perform self-cleavage, which correlated albeit not perfectly with an altered subcellular localization of the resulting proteins. For the processing-proficient variants, we used structural modelling tools and molecular dynamic (MD) simulations to predict the potential impact of the variants on protein conformation and/or interaction with partners. Our study contributes to a better understanding of the molecular mechanisms involved in DHH dysfunction leading to 46,XY disorders of sex development.

Molecular Dynamics Simulations of Influenza A Virus NS1 Reveal a Remarkably Stable RNA-Binding Domain Harboring Promising Druggable Pockets.

The non-structural protein NS1 of influenza A viruses is considered to be the major antagonist of the interferon system and antiviral defenses of the cell. It could therefore represent a suitable target for novel antiviral strategies. As a first step towards the identification of small compounds targeting NS1, we here investigated the druggable potential of its RNA-binding domain since this domain is essential to the biological activities of NS1. We explored the flexibility of the full-length protein by running molecular dynamics simulations on one of its published crystal structures. While the RNA-binding domain structure was remarkably stable along the simulations, we identified a flexible site at the two extremities of the "groove" that is delimited by the antiparallel α-helices that make up its RNA-binding interface. This groove region is able to form potential binding pockets, which, in 60% of the conformations, meet the druggability criteria. We characterized these pockets and identified the residues that contribute to their druggability. All the residues involved in the druggable pockets are essential at the same time to the stability of the RNA-binding domain and to the biological activities of NS1. They are also strictly conserved across the large sequence diversity of NS1, emphasizing the robustness of this search towards the identification of broadly active NS1-targeting compounds.

Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures.

HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.Communicated by Ramaswamy H. Sarma.

Novel STAG3 mutations in a Caucasian family with primary ovarian insufficiency.

Primary ovarian insufficiency (POI) affects ~ 1-3, 7% of women under forty and is a public health problem. Most causes are unknown, but an increasing number of genetic causes have been identified recently. The identification of such causes is essential for genetic and therapeutic counseling in patients and their families. We performed whole exome sequencing in two Caucasian sisters displaying non syndromic POI and their unaffected mother. We identified two novel pathogenic variants in STAG3 encoding a meiosis-specific subunit of the cohesin ring, which ensures correct sister chromatid cohesion: a c.3052delC truncating mutation in exon 28 yielding p.Arg1018Aspfs*14, and a c.659T > G substitution in exon seven yielding p.Leu220Arg. Leu220, highly conserved throughout species, belongs to the STAG domain conserved with other mitotic subunits of the cohesion complex STAG1 and 2. In silico analysis reveals that this substitution markedly impacts the structure of this domain. The truncation removes the last 206 C-terminal residues, not conserved in STAG1 and 2, supporting an important specific role in STAG3, especially meiosis. This is the first occurrence of STAG3 mutations in a Caucasian family. Very little is known about the function of STAG proteins domains. The "knock out-like" phenotype described here supports the crucial role of a single residue in the STAG domain and of the C-terminal region in STAG3 function. In conclusion, this observation shows the necessity to perform the genetic study of POI worldwide including STAG3. This could lead to appropriate genetic counseling and long term follow-up since these patients may develop ovarian tumors.

SAFlex: A structural alphabet extension to integrate protein structural flexibility and missing data information.

In this paper, we describe SAFlex (Structural Alphabet Flexibility), an extension of an existing structural alphabet (HMM-SA), to better explore increasing protein three dimensional structure information by encoding conformations of proteins in case of missing residues or uncertainties. An SA aims to reduce three dimensional conformations of proteins as well as their analysis and comparison complexity by simplifying any conformation in a series of structural letters. Our methodology presents several novelties. Firstly, it can account for the encoding uncertainty by providing a wide range of encoding options: the maximum a posteriori, the marginal posterior distribution, and the effective number of letters at each given position. Secondly, our new algorithm deals with the missing data in the protein structure files (concerning more than 75% of the proteins from the Protein Data Bank) in a rigorous probabilistic framework. Thirdly, SAFlex is able to encode and to build a consensus encoding from different replicates of a single protein such as several homomer chains. This allows localizing structural differences between different chains and detecting structural variability, which is essential for protein flexibility identification. These improvements are illustrated on different proteins, such as the crystal structure of an eukaryotic small heat shock protein. They are promising to explore increasing protein redundancy data and obtain useful quantification of their flexibility.

Exploration of the effect of sequence variations located inside the binding pocket of HIV-1 and HIV-2 proteases.

HIV-2 protease (PR2) is naturally resistant to most FDA (Food and Drug Administration)-approved HIV-1 protease inhibitors (PIs), a major antiretroviral class. In this study, we compared the PR1 and PR2 binding pockets extracted from structures complexed with 12 ligands. The comparison of PR1 and PR2 pocket properties showed that bound PR2 pockets were more hydrophobic with more oxygen atoms and fewer nitrogen atoms than PR1 pockets. The structural comparison of PR1 and PR2 pockets highlighted structural changes induced by their sequence variations and that were consistent with these property changes. Specifically, substitutions at residues 31, 46, and 82 induced structural changes in their main-chain atoms that could affect PI binding in PR2. In addition, the modelling of PR1 mutant structures containing V32I and L76M substitutions revealed a cooperative mechanism leading to structural deformation of flap-residue 45 that could modify PR2 flexibility. Our results suggest that substitutions in the PR1 and PR2 pockets can modify PI binding and flap flexibility, which could underlie PR2 resistance against PIs. These results provide new insights concerning the structural changes induced by PR1 and PR2 pocket variation changes, improving the understanding of the atomic mechanism of PR2 resistance to PIs.

Analysis of the HIV-2 protease's adaptation to various ligands: characterization of backbone asymmetry using a structural alphabet.

The HIV-2 protease (PR2) is a homodimer of 99 residues with asymmetric assembly and binding various ligands. We propose an exhaustive study of the local structural asymmetry between the two monomers of all available PR2 structures complexed with various inhibitors using a structural alphabet approach. On average, PR2 exhibits asymmetry in 31% of its positions-i.e., exhibiting different backbone local conformations in the two monomers. This asymmetry was observed all along its structure, particularly in the elbow and flap regions. We first differentiated structural asymmetry conserved in most PR2 structures from the one specific to some PR2. Then, we explored the origin of the detected asymmetry in PR2. We localized asymmetry that could be induced by PR2's flexibility, allowing transition from the semi-open to closed conformations and the asymmetry potentially induced by ligand binding. This latter could be important for the PR2's adaptation to diverse ligands. Our results highlighted some differences between asymmetry of PR2 bound to darunavir and amprenavir that could explain their differences of affinity. This knowledge is critical for a better description of PR2's recognition and adaptation to various ligands and for a better understanding of the resistance of PR2 to most PR2 inhibitors, a major antiretroviral class.

Exploring the potential of a structural alphabet-based tool for mining multiple target conformations and target flexibility insight.

Protein flexibility is often implied in binding with different partners and is essential for protein function. The growing number of macromolecular structures in the Protein Data Bank entries and their redundancy has become a major source of structural knowledge of the protein universe. The analysis of structural variability through available redundant structures of a target, called multiple target conformations (MTC), obtained using experimental or modeling methods and under different biological conditions or different sources is one way to explore protein flexibility. This analysis is essential to improve the understanding of various mechanisms associated with protein target function and flexibility. In this study, we explored structural variability of three biological targets by analyzing different MTC sets associated with these targets. To facilitate the study of these MTC sets, we have developed an efficient tool, SA-conf, dedicated to capturing and linking the amino acid and local structure variability and analyzing the target structural variability space. The advantage of SA-conf is that it could be applied to divers sets composed of MTCs available in the PDB obtained using NMR and crystallography or homology models. This tool could also be applied to analyze MTC sets obtained by dynamics approaches. Our results showed that SA-conf tool is effective to quantify the structural variability of a MTC set and to localize the structural variable positions and regions of the target. By selecting adapted MTC subsets and comparing their variability detected by SA-conf, we highlighted different sources of target flexibility such as induced by binding partner, by mutation and intrinsic flexibility. Our results support the interest to mine available structures associated with a target using to offer valuable insight into target flexibility and interaction mechanisms. The SA-conf executable script, with a set of pre-compiled binaries are available at http://www.mti.univ-paris-diderot.fr/recherche/plateformes/logiciels.

Statistical Profiling of One Promiscuous Protein Binding Site: Illustrated by Urokinase Catalytic Domain.

While recent literature focuses on drug promiscuity, the characterization of promiscuous binding sites (ability to bind several ligands) remains to be explored. Here, we present a proteochemometric modeling approach to analyze diverse ligands and corresponding multiple binding sub-pockets associated with one promiscuous binding site to characterize protein-ligand recognition. We analyze both geometrical and physicochemical profile correspondences. This approach was applied to examine the well-studied druggable urokinase catalytic domain inhibitor binding site, which results in a large number of complex structures bound to various ligands. This approach emphasizes the importance of jointly characterizing pocket and ligand spaces to explore the impact of ligand diversity on sub-pocket properties and to establish their main profile correspondences. This work supports an interest in mining available 3D holo structures associated with a promiscuous binding site to explore its main protein-ligand recognition tendency.

A novel variant of DHH in a familial case of 46,XY disorder of sex development: Insights from molecular dynamics simulations.

OBJECTIVE: Disorders of sex development (DSD) are a heterogeneous group of conditions affecting the differentiation and development of the internal and external genitalia. Here, we aimed at identifying the genetic cause of DSD in two 46,XY sisters from a consanguineous family. DESIGN: We performed a whole-exome sequencing of two 46,XY female individuals. Sanger sequencing was used to validate the most likely candidate variant, affecting the desert hedgehog (DHH) gene. Molecular dynamics simulations were performed to get insights into the impact of the variant on protein structure and on its interaction with the protein partner BOC (brother of CDO/cell adhesion molecule, downregulated by oncogenes). PATIENTS: The index patient presented with a female phenotype, primary amenorrhoea (low oestradiol and testosterone and high FSH and LH). She also had an apparent absence of intra-abdominal gonads and uterus, facial dysmorphy, psychomotor retardation and neuropathy. Her sister displayed a similar gonadal and endocrinological picture, without dysmorphy or psychomotor retardation. RESULTS: Whole-exome sequencing revealed a homozygous variant in DHH leading to the p.Trp173Cys substitution. The relevant Trp residue is conserved, and its alteration was predicted to be deleterious. Molecular dynamics simulations showed that the mutation increases the conformational flexibility of the protein and potentially alters its interaction with BOC, a positive regulator of Hedgehog signalling. We do not exclude an interference of the mutation with DHH-intein-mediated auto-processing. CONCLUSIONS: This report increases the number of described homozygous DHH variants and highlights the importance of advanced bioinformatic tools to better understand the pathogenicity of human variants.

Global vision of druggability issues: applications and perspectives.

During the preliminary stage of a drug discovery project, the lack of druggability information and poor target selection are the main causes of frequent failures. Elaborating on accurate computational druggability prediction methods is a requirement for prioritizing target selection, designing new drugs and avoiding side effects. In this review, we describe a survey of recently reported druggability prediction methods mainly based on networks, statistical pocket druggability predictions and virtual screening. An application for a frequent mutation of p53 tumor suppressor is presented, illustrating the complementarity of druggability prediction approaches, the remaining challenges and potential new drug development perspectives.

Screening for protein-protein interactions using Förster resonance energy transfer (FRET) and fluorescence lifetime imaging microscopy (FLIM).

We present a high content multiwell plate cell-based assay approach to quantify protein interactions directly in cells using Förster resonance energy transfer (FRET) read out by automated fluorescence lifetime imaging (FLIM). Automated FLIM is implemented using wide-field time-gated detection, typically requiring only 10 s per field of view (FOV). Averaging over biological, thermal and shot noise with 100's to 1000's of FOV enables unbiased quantitative analysis with high statistical power. Plotting average donor lifetime vs. acceptor/donor intensity ratio clearly identifies protein interactions and fitting to double exponential donor decay models provides estimates of interacting population fractions that, with calibrated donor and acceptor fluorescence intensities, can yield dissociation constants. We demonstrate the application to identify binding partners of MST1 kinase and estimate interaction strength among the members of the RASSF protein family, which have important roles in apoptosis via the Hippo signalling pathway. KD values broadly agree with published biochemical measurements.

A Hot-spot of In-frame Duplications Activates the Oncoprotein AKT1 in Juvenile Granulosa Cell Tumors.

BACKGROUND: Ovarian granulosa cell tumors are the most common sex-cord stromal tumors and have juvenile (JGCTs) and adult forms. In a previous study we reported the occurrence of activating somatic mutations of Gαs, which transduces mitogenic signals, in 30% of the analyzed JGCTs. METHODS: We have searched for alterations in other proteins involved in ovarian mitogenic signaling. We focused on the PI3K-AKT axis. As we found mutations in AKT1, we analyzed the subcellular localization of the mutated proteins and performed functional explorations using Western-blot and luciferase assays. FINDINGS: We detected in-frame duplications affecting the pleckstrin-homology domain of AKT1 in more than 60% of the tumors occurring in girls under 15 years of age. The somatic status of the mutations was confirmed when peritumoral DNA was available. The JGCTs without duplications carried point mutations affecting highly conserved residues. Several of these substitutions were somatic lesions. The mutated proteins carrying the duplications had a non-wild-type subcellular distribution, with a marked enrichment at the plasma membrane. This led to a striking degree of AKT1 activation demonstrated by a strong phosphorylation level and by reporter assays. INTERPRETATION: Our study incriminates somatic mutations of AKT1 as a major event in the pathogenesis of JGCTs. The existence of AKT inhibitors currently tested in clinical trials opens new perspectives for targeted therapies for these tumors, which are currently treated with standard non-specific chemotherapy protocols.

Comparative analysis of interactions of RASSF1-10.

Members of the RASSF family (RASSF1-10) have been identified as candidate tumour suppressors that are frequently downregulated by promoter hypermethylation in cancers. These proteins carry a common Ras-association (RA) and SARAH domain (RASSF1-6) that can potentially bind Ras oncoproteins and mediate protein-protein interactions with other SARAH domain proteins. However, there is a notable lack of comparative characterisation of the RASSF family, as well as molecular and structural information that facilitate their tumour suppressive functions. As part of our comparative analysis, we modelled the RA and SARAH domains of the RASSF members based on existing structures and predicted their potential interactions. These in silico predictions were compared to in vitro interaction studies with Ras and MST kinase (a SARAH domain-containing protein). Our data shows a diversity of interaction within the RASSF family RA domain, whereas the SARAH domain-mediated interactions for RASSF1-6 are consistent with the predictions. This suggests that different members, despite shared general architecture, could have distinct functional properties. Additionally, we identify a new interacting partner for MST kinase in the form of RASSF7. Current data supports an interaction model where RASSF serves as an adaptor for the assembly of multiple protein complexes and further functional interactions, involving MST kinases and other SARAH domain proteins, which could be regulated by Ras.

Detection and architecture of small heat shock protein monomers.

BACKGROUND: Small Heat Shock Proteins (sHSPs) are chaperone-like proteins involved in the prevention of the irreversible aggregation of misfolded proteins. Although many studies have already been conducted on sHSPs, the molecular mechanisms and structural properties of these proteins remain unclear. Here, we propose a better understanding of the architecture, organization and properties of the sHSP family through structural and functional annotations. We focused on the Alpha Crystallin Domain (ACD), a sandwich fold that is the hallmark of the sHSP family. METHODOLOGY/PRINCIPAL FINDINGS: We developed a new approach for detecting sHSPs and delineating ACDs based on an iterative Hidden Markov Model algorithm using a multiple alignment profile generated from structural data on ACD. Using this procedure on the UniProt databank, we found 4478 sequences identified as sHSPs, showing a very good coverage with the corresponding PROSITE and Pfam profiles. ACD was then delimited and structurally annotated. We showed that taxonomic-based groups of sHSPs (animals, plants, bacteria) have unique features regarding the length of their ACD and, more specifically, the length of a large loop within ACD. We detailed highly conserved residues and patterns specific to the whole family or to some groups of sHSPs. For 96% of studied sHSPs, we identified in the C-terminal region a conserved I/V/L-X-I/V/L motif that acts as an anchor in the oligomerization process. The fragment defined from the end of ACD to the end of this motif has a mean length of 14 residues and was named the C-terminal Anchoring Module (CAM). CONCLUSIONS/SIGNIFICANCE: This work annotates structural components of ACD and quantifies properties of several thousand sHSPs. It gives a more accurate overview of the architecture of sHSP monomers.