Getting to know each other: PPIMem, a novel approach for predicting transmembrane protein-protein complexes.

Because of their considerable number and diversity, membrane proteins and their macromolecular complexes represent the functional units of cells. Their quaternary structure may be stabilized by interactions between the α-helices of different proteins in the hydrophobic region of the cell membrane. Membrane proteins equally represent potential pharmacological targets par excellence for various diseases. Unfortunately, their experimental 3D structure and that of their complexes with other intramembrane protein partners are scarce due to technical difficulties. To overcome this key problem, we devised PPIMem, a computational approach for the specific prediction of higher-order structures of α-helical transmembrane proteins. The novel approach involves proper identification of the amino acid residues at the interface of molecular complexes with a 3D structure. The identified residues compose then nonlinear interaction motifs that are conveniently expressed as mathematical regular expressions. These are efficiently implemented for motif search in amino acid sequence databases, and for the accurate prediction of intramembrane protein-protein complexes. Our template interface-based approach predicted 21,544 binary complexes between 1,504 eukaryotic plasma membrane proteins across 39 species. We compare our predictions to experimental datasets of protein-protein interactions as a first validation method. The online database that results from the PPIMem algorithm with the annotated predicted interactions are implemented as a web server and can be accessed directly at https://transint.univ-evry.fr.

Lin28, a major translation reprogramming factor, gains access to YB-1-packaged mRNA through its cold-shock domain.

The RNA-binding protein Lin28 (Lin28a) is an important pluripotency factor that reprograms translation and promotes cancer progression. Although Lin28 blocks let-7 microRNA maturation, Lin28 also binds to a large set of cytoplasmic mRNAs directly. However, how Lin28 regulates the processing of many mRNAs to reprogram global translation remains unknown. We show here, using a structural and cellular approach, a mixing of Lin28 with YB-1 (YBX1) in the presence of mRNA owing to their cold-shock domain, a conserved ß-barrel structure that binds to ssRNA cooperatively. In contrast, the other RNA binding-proteins without cold-shock domains tested, HuR, G3BP-1, FUS and LARP-6, did not mix with YB-1. Given that YB-1 is the core component of dormant mRNPs, a model in which Lin28 gains access to mRNPs through its co-association with YB-1 to mRNA may provide a means for Lin28 to reprogram translation. We anticipate that the translational plasticity provided by mRNPs may contribute to Lin28 functions in development and adaptation of cancer cells to an adverse environment.

Molecular dynamics of the histamine H3 membrane receptor reveals different mechanisms of GPCR signal transduction.

In this work, we studied the mechanisms of classical activation and inactivation of signal transduction by the histamine H3 receptor, a 7-helix transmembrane bundle G-Protein Coupled Receptor through long-time-scale atomistic molecular dynamics simulations of the receptor embedded in a hydrated double layer of dipalmitoyl phosphatidyl choline, a zwitterionic polysaturated ordered lipid. Three systems were prepared: the apo receptor, representing the constitutively active receptor; and two holo-receptors-the receptor coupled to the antagonist/inverse agonist ciproxifan, representing the inactive state of the receptor, and the receptor coupled to the endogenous agonist histamine and representing the active state of the receptor. An extensive analysis of the simulation showed that the three states of H3R present significant structural and dynamical differences as well as a complex behavior given that the measured properties interact in multiple and interdependent ways. In addition, the simulations described an unexpected escape of histamine from the orthosteric binding site, in agreement with the experimental modest affinities and rapid off-rates of agonists.

YB-1, an abundant core mRNA-binding protein, has the capacity to form an RNA nucleoprotein filament: a structural analysis.

The structural rearrangements accompanying mRNA during translation in mammalian cells remain poorly understood. Here, we discovered that YB-1 (YBX1), a major partner of mRNAs in the cytoplasm, forms a linear nucleoprotein filament with mRNA, when part of the YB-1 unstructured C-terminus has been truncated. YB-1 possesses a cold-shock domain (CSD), a remnant of bacterial cold shock proteins that have the ability to stimulate translation under the low temperatures through an RNA chaperone activity. The structure of the nucleoprotein filament indicates that the CSD of YB-1 preserved its chaperone activity also in eukaryotes and shows that mRNA is channeled between consecutive CSDs. The energy benefit needed for the formation of stable nucleoprotein filament relies on an electrostatic zipper mediated by positively charged amino acid residues in the YB-1 C-terminus. Thus, YB-1 displays a structural plasticity to unfold structured mRNAs into extended linear filaments. We anticipate that our findings will shed the light on the scanning of mRNAs by ribosomes during the initiation and elongation steps of mRNA translation.

Multiplex epithelium dysfunction due to CLDN10 mutation: the HELIX syndrome.

PurposeWe aimed to identify the genetic cause to a clinical syndrome encompassing hypohidrosis, electrolyte imbalance, lacrimal gland dysfunction, ichthyosis, and xerostomia (HELIX syndrome), and to comprehensively delineate the phenotype.MethodsWe performed homozygosity mapping, whole-genome sequencing, gene sequencing, expression studies, functional tests, protein bioinformatics, and histological characterization in two unrelated families with HELIX syndrome.ResultsWe identified biallelic missense mutations (c.386C>T, p.S131L and c.2T>C, p.M1T) in CLDN10B in six patients from two unrelated families. CLDN10B encodes Claudin-10b, an integral tight junction (TJ) membrane-spanning protein expressed in the kidney, skin, and salivary glands. All patients had hypohidrosis, renal loss of NaCl with secondary hyperaldosteronism and hypokalemia, as well as hypolacrymia, ichthyosis, xerostomia, and severe enamel wear. Functional testing revealed that patients had a decreased NaCl absorption in the thick ascending limb of the loop of Henle and a severely decreased secretion of saliva. Both mutations resulted in reduced or absent Claudin-10 at the plasma membrane of epithelial cells.ConclusionCLDN10 mutations cause a dysfunction in TJs in several tissues and, subsequently, abnormalities in renal ion transport, ectodermal gland homeostasis, and epidermal integrity.

Deciphering the GPER/GPR30-agonist and antagonists interactions using molecular modeling studies, molecular dynamics, and docking simulations.

The G-protein coupled estrogen receptor 1 GPER/GPR30 is a transmembrane seven-helix (7TM) receptor involved in the growth and proliferation of breast cancer. Due to the absence of a crystal structure of GPER/GPR30, in this work, molecular modeling studies have been carried out to build a three-dimensional structure, which was subsequently refined by molecular dynamics (MD) simulations (up to 120 ns). Furthermore, we explored GPER/GPR30's molecular recognition properties by using reported agonist ligands (G1, estradiol (E2), tamoxifen, and fulvestrant) and the antagonist ligands (G15 and G36) in subsequent docking studies. Our results identified the E2 binding site on GPER/GPR30, as well as other receptor cavities for accepting large volume ligands, through GPER/GPR30 π-π, hydrophobic, and hydrogen bond interactions. Snapshots of the MD trajectory at 14 and 70 ns showed almost identical binding motifs for G1 and G15. It was also observed that C107 interacts with the acetyl oxygen of G1 (at 14 ns) and that at 70 ns the residue E275 interacts with the acetyl group and with the oxygen from the other agonist whereas the isopropyl group of G36 is oriented toward Met141, suggesting that both C107 and E275 could be involved in the protein activation. This contribution suggest that GPER1 has great structural changes which explain its great capacity to accept diverse ligands, and also, the same ligand could be recognized in different binding pose according to GPER structural conformations.

Mutations in zinc finger 407 [ZNF407] cause a unique autosomal recessive cognitive impairment syndrome.

BACKGROUND: A consanguineous Arab family is affected by an apparently novel autosomal recessive disorder characterized by cognitive impairment, failure-to-thrive, hypotonia and dysmorphic features including bilateral ptosis and epicanthic folds, synophrys, midface hypoplasia, downturned mouth corners, thin upper vermillion border and prominent ears, bilateral 5th finger camptodactyly, bilateral short 4th metatarsal bones, and limited knee mobility bilaterally. METHODS: The family was studied by homozygosity mapping, candidate gene mutation screening and whole Exome Next Generation Sequencing of a single affected member to identify the offending gene and mutation. The mutated gene product was studied by structural bioinformatics methods. RESULTS: A damaging c.C5054G mutation affecting an evolutionary highly conserved amino acid p.S1685W was identified in the ZNF407 gene at 18q23. The Serine to Tryptophane mutation affects two of the three ZNF407 isoforms and is located in the last third of the protein, in a linker peptide adjoining two zinc-finger domains. Structural analyses of this mutation shows disruption of an H-bond that locks the relative spatial position of the two fingers, leading to a higher flexibility of the linker and thus to a decreased probability of binding to the target DNA sequence essentially eliminating the functionality of downstream domains and interfering with the expression of various genes under ZNF407 control during fetal brain development. CONCLUSIONS: ZNF407 is a transcription factor with an essential role in brain development. When specific and limited in number homozygosity intervals exist that harbor the offending gene in consanguineous families, Whole Exome Sequencing of a single affected individual is an efficient approach to gene mapping and mutation identification.

Homology modeling of the structure of acyl coA:isopenicillin N-acyltransferase (IAT) from Penicillium chrysogenum. IAT interaction studies with isopenicillin-N, combining molecular dynamics simulations and docking.

In the last step of penicillin biosynthesis, acyl-CoA:isopenicillin N acyltransferase (IAT) (E.C. 2.3.1.164) catalyzes the conversion of isopenicillin N (IPN) to penicillin G. IAT substitutes the α-aminoadipic acid side chain of IPN by a phenylacetic acid phenolate group (from phenylacetyl-CoA). Having a three-dimensional (3D) structure of IAT helps to determine the steps involved in side chain exchange by identifying the atomic details of substrate recognition. We predicted the IAT 3-D structure (α- and ß-subunits), as well as the manner of IPN and phenylacetyl-CoA bind to the mature enzyme (ß-subunit). The 3D IAT prediction was achieved by homology modeling and molecular docking in different snapshots, and refined by molecular dynamic simulations. Our model can reasonably interpret the results of a number of experiments, where key residues for IAT processing as well as strictly conserved residues most probably involved with enzymatic activity were mutated. Based on the results of docking studies, energies associated with the complexes, and binding constants calculated, we identified a site located in the region generated by ß1, ß2 and ß5 strands, which forms part of the central structure of ß-subunit, as the potential binding site of IPN. The site comprises the amino acid residues Cys103, Asp121, Phe122, Phe123, Ala168, Leu169, His170, Gln172, Phe212, Arg241, Leu262, Asp264, Arg302, Ser309, and Arg310. Through hydrogen bonds, the IPN binding site establishes interactions with Cys103, Leu169, Gln172, Asp264 and Arg310. Our model is also validated by a recently revealed crystal structure of the mature enzyme.

Bioinformatics and functional analysis of an Entamoeba histolytica mannosyltransferase necessary for parasite complement resistance and hepatical infection.

The glycosylphosphatidylinositol (GPI) moiety is one of the ways by which many cell surface proteins, such as Gal/GalNAc lectin and proteophosphoglycans (PPGs) attach to the surface of Entamoeba histolytica, the agent of human amoebiasis. It is believed that these GPI-anchored molecules are involved in parasite adhesion to cells, mucus and the extracellular matrix. We identified an E. histolytica homolog of PIG-M, which is a mannosyltransferase required for synthesis of GPI. The sequence and structural analysis led to the conclusion that EhPIG-M1 is composed of one signal peptide and 11 transmembrane domains with two large intra luminal loops, one of which contains the DXD motif, involved in the enzymatic catalysis and conserved in most glycosyltransferases. Expressing a fragment of the EhPIG-M1 encoding gene in antisense orientation generated parasite lines diminished in EhPIG-M1 levels; these lines displayed reduced GPI production, were highly sensitive to complement and were dramatically inhibited for amoebic abscess formation. The data suggest a role for GPI surface anchored molecules in the survival of E. histolytica during pathogenesis.

Characterization of a human coagulation factor Xa-binding site on Viperidae snake venom phospholipases A2 by affinity binding studies and molecular bioinformatics.

BACKGROUND: The snake venom group IIA secreted phospholipases A2 (SVPLA2), present in the Viperidae snake family exhibit a wide range of toxic and pharmacological effects. They exert their different functions by catalyzing the hydrolysis of phospholipids (PL) at the membrane/water interface and by highly specific direct binding to: (i) presynaptic membrane-bound or intracellular receptors; (ii) natural PLA2-inhibitors from snake serum; and (iii) coagulation factors present in human blood. RESULTS: Using surface plasmon resonance (SPR) protein-protein interaction measurements and an in vitro biological test of inhibition of prothrombinase activity, we identify a number of Viperidae venom SVPLA2s that inhibit blood coagulation through direct binding to human blood coagulation factor Xa (FXa) via a non-catalytic, PL-independent mechanism. We classify the SVPLA2s in four groups, depending on the strength of their binding. Molecular electrostatic potentials calculated at the surface of 3D homology-modeling models show a correlation with inhibition of prothrombinase activity. In addition, molecular docking simulations between SVPLA2 and FXa guided by the experimental data identify the potential FXa binding site on the SVPLA2s. This site is composed of the following regions: helices A and B, the Ca2+ loop, the helix C-beta-wing loop, and the C-terminal fragment. Some of the SVPLA2 binding site residues belong also to the interfacial binding site (IBS). The interface in FXa involves both, the light and heavy chains. CONCLUSION: We have experimentally identified several strong FXa-binding SVPLA2s that disrupt the function of the coagulation cascade by interacting with FXa by the non-catalytic PL-independent mechanism. By theoretical methods we mapped the interaction sites on both, the SVPLA2s and FXa. Our findings may lead to the design of novel, non-competitive FXa inhibitors.

Snake venom serine proteinases: sequence homology vs. substrate specificity, a paradox to be solved.

Snake venom glands synthesize a variety of serine proteinases capable of affecting the haemostatic system. They act on macromolecular substrates of the coagulation, fibrinolytic, and kallikrein-kinin systems, and on platelets to cause an imbalance of the haemostatic system of the prey. In this review we describe their biochemical/biophysical characteristics, biological activities as well as aspects of their evolution and structure-activity relationship.

Entamoeba histolytica: inhibition of cellular functions by overexpression of EhGEF1, a novel Rho/Rac guanine nucleotide exchange factor.

The molecular, biochemical, and cellular characterization of EhGEF1 protein is described. Complete cDNA sequence of 1890 bp revealed an open reading frame that encodes a protein of 69 kDa. EhGEF1 is constituted of Dbl homology domain, pleckstrin homology domain, and several putative regulation sites. Studies of guanine nucleotide exchange activity of EhGEF1 on several GTPases from Entamoeba histolytica and Homo sapiens showed preferential activation on EhRacG, suggesting that EhGEF1 protein could be involved in mechanisms related to actin cytoskeleton activation, cytokinesis, capping, and uroid formation in trophozoite. Confocal microscopy studies of pExEhNeo/HSV-tagged-EhGEF1-transfected cells showed that trophozoites stimulated with ConA, EhGEF1, and EhRacG were localized at plasma membrane. Cellular studies showed that F-actin content of pExEhNeo/HSV-tagged-EhGEF1-transfected trophozoites as well as cellular migration and cell damage capacity were significantly altered. The observations suggest that EhRacG was the principal target of EhGEF1 and that EhGEF1 may provide a link between F-actin dynamics and EhRacG signaling.