Increased Understanding of the Biochemistry and Biosynthesis of MUC2 and Other Gel-Forming Mucins Through the Recombinant Expression of Their Protein Domains

The gel-forming mucins are large and heavily O-glycosylated proteins which build up mucus gels. The recombinant production of full-length gel-forming mucins has not been possible to date. In order to study mucin biosynthesis and biochemistry, we and others have taken the alternative approach of constructing different recombinant proteins consisting of one or several domains of these large proteins and expressing them separately in different cell lines. Using this approach, we have determined that MUC2, the intestinal gel-forming mucin, dimerizes via its C-terminal cysteine-knot domain and also trimerizes via one of the N-terminal von Willebrand D domains. Both of these interactions are disulfide bond mediated. Via this assembly, a molecular network is built by which the mucus gel is formed. Here we discuss not only the functional understanding obtained from studies of the recombinant proteins, but also highlight the difficulties encountered when these proteins were produced recombinantly. We often found an accumulation of the proteins in the ER and consequently no secretion. This was especially apparent when the cysteine-rich domains of the N- and C-terminal parts of the mucins were expressed. Other proteins that we constructed were either not secreted or not expressed at all. Despite these problems, the knowledge of mucin biosynthesis and assembly has advanced considerably through the studies of these recombinant proteins.     

Conformational stability and three-dimensional model of the δ-opioid pharmacophore for the extended antiparallel dimer structure of Met-enkephalin in water

The conformational stability of the extended antiparallel dimer structure of Met-enkephalin in water was analyzed by examining the hydration structure of enkephalin using molecular dynamics simulations. The result shows that, despite of the hydrophicility of the terminal atoms in the pentapeptide, the main contributor for the stability of the dimer in water is the four intermolecular hydrogen bonds between the Gly² and Phe⁴ groups. The three-dimensional model of the δ-opioid pharmacophore for this dimer structure was also established. Such a model was demonstrated to match the δ-opioid pharmacophore query derived from the non-peptides SIOM, TAN-67, and OMI perfectly. This result thus strongly supports the assumption that the dimer structure of Met-enkephalin is a possible δ-receptor binding conformation. Figure Schematic model of the extended antiparallel dimer structure of Met-enkephalin     

MUC16 is produced in tracheal surface epithelium and submucosal glands and is present in secretions from normal human airway and cultured bronchial epithelial cells

The gel-forming MUC5AC and MUC5B mucins have been identified as major components of human airway mucus but it is not known whether additional mucin species, possibly with other functions, are also present. MUC16 mucin is a well-known serum marker for ovarian cancer, but the molecule has also been found on the ocular surface and in cervical secretions suggesting that it may play a role on the normal mucosal surface. In this investigation, the LUM16-2 antiserum (raised against a sequence in the N-terminal repeat domain) recognized MUC16 in goblet and submucosal gland mucous cells as well as on the epithelial surface of human tracheal tissue suggesting that the mucin originates from secretory cells. MUC16 mucin was present in 'normal' respiratory tract mucus as well as in secretions from normal human bronchial epithelial (NHBE) cells. MUC16 from NHBE cells was a high-molecular-mass, monomeric mucin which gave rise to large glycopeptides after proteolysis. N- and C-terminal fragments of the molecule were separated on gel electrophoresis showing that the MUC16 apoprotein undergoes a cleavage between these domains, possibly in the SEA domain as demonstrated for other transmembrane mucins; MUC1 and MUC3. After metabolic labeling of NHBE cells, most of the secreted monomeric, high-molecular-mass [(35)S]sulphate-labelled molecules were immunoprecipitated with the OC125 antibody indicating that MUC16 is the major [(35)S]sulphate-labelled mucin in NHBE cell secretions.     

Proton transfer in some periodic molecular systems

The electronic structure of representative hydrogen bonded systems: hydrogen cyanide, imidazole and malonic acid have been studied at the non-empirical level. The role of the dimensionality on the potential barrier for the proton transfer has been examined. It was shown that it depends on the crystal structure and only in some cases like hydrogen cyanide or imidazole the relevant crystals may be considered as one-dimensional. However, for more complicated crystallographic structures, e.g. malonic acid, the evaluated barrier is strongly dependent on the dimensionality taken into account in our calculations. Figure Density of states for the imidazole dimer: MD-EQ- standard proton position (black line), MD-TR- proton position in the middle of the hydrogen bond (red line)     

The human MUC2 intestinal mucin has cysteine-rich subdomains located both upstream and downstream of its central repetitive region.

The human MUC2 mucin is a large secretory glycoconjugate that coats the epithelia of the intestines, airways, and other mucus membrane-containing organs. Previous work has shown that this mucin contains an extended tandem repeat-containing domain rich in Thr and Pro. In the present work we describe two additional regions of this mucin located both upstream and downstream of the tandem repeat array. The carboxyl-terminal domain contains 984 residues and can be divided into mucin-like (139 residues) and cysteine-rich (845 residues) subdomains. This latter subdomain exhibits varying degrees of sequence similarity to a wide range of mucins and mucin-like proteins including those isolated from rats, pigs, cows, and frogs. We also report here the sequence of 1270 residues lying immediately upstream of the tandem repeats. This region contains a repetitive, mucin-like subdomain and a second cysteine-rich stretch of more than 700 residues. Both cysteine-rich subdomains of this mucin have sequence similarity with von Willebrand factor, a serum protein that exists as a disulfide-linked polymer. This suggests that these cysteine-rich subdomains are important in the catenation of mucin monomers into oligomers, the structures that confer viscoelasticity upon mucus.     

Identification of a nonmucin glycoprotein (gp-340) from a purified respiratory mucin preparation: evidence for an association involving the MUC5B mucin.

Rate-zonal centrifugation of a reduced and alkylated respiratory mucin preparation identified a protein-rich fraction. This was subjected to trypsin treatment and one of the many liberated peptides was purified and its N-terminal sequence determined. The peptide was identical to a 14 amino acid sequence from the scavenger receptor cysteine-rich domain containing glycoprotein gp-340. A polyclonal antiserum, raised against the peptide, stained the serous cells in the submucosal glands of human tracheal tissue. The glycoprotein was purified from respiratory mucus by density-gradient centrifugation, gel chromatography, and anion exchange chromatography. The molecule exhibited a heterogeneous distribution of buoyant density (1.28-1.46 g/ml) that overlapped with the gel-forming mucins, was included on Sepharose CL-2B and was quite highly anionic. SDS-PAGE indicated a mass greater than 208 kDa and measurements performed across the molecular size distribution indicated an average M(r) of 5 x 10(5) with a range of M(r) from 2 x 10(5) to 1 x 10(6). Gel chromatography of respiratory mucus extracts ("associative" and "dissociative") indicated that this glycoprotein forms complexes that may involve the large gel-forming mucins MUC5AC and MUC5B. Rate zonal centrifugation suggested such complexes are more likely to involve MUC5B rather than MUC5AC mucins.     

Cellular interaction through LewisX cluster: theoretical studies

It is well known that cell surface carbohydrates play a role in cell–cell adhesion and communication. LewisX glycosphingolipids form microdomains on cell surfaces. Homotypic and calcium-mediated LewisX–LewisX (LeX-LeX) interactions were proposed to be responsible for the initial steps of cell adhesion, and to mediate embryogenesis and metastasis. Various techniques have been used to investigate such interactions, but little information is available on the geometry and the mechanism of dimerisation. To better understand these interactions, a new molecular model was developed to simulate homotypic interactions in explicit solvent with and without calcium ions. Accurate analysis of both trajectories yielded valuable information about the energetics of LeX-LeX dimerisation. Detailed interpretation of the hydrogen bond network and the presence of calcium ions along the trajectory provide valuable insights into the role of calcium ions in this carbohydrate–carbohydrate interaction. Figure Calcium population density around the LewisX carbohydrate (after the trajectory has been fitted to the primary unit cell). All central dimer coordinates are fitted along the time axis, whereas calcium ion positions are recorded and represented as points. The clouds of points indicate that the ions are not randomly placed around the dimer but take up preferred positions     

Alpha particle chemistry. On the formation of stable complexes between He2+ and other simple species: implications for atmospheric and interstellar chemistry

The possibility that stable complexes may be formed between alpha particles (He²⁺) and small molecules is investigated using QCISD quantum mechanical calculations. Implications for their presence in the terrestrial atmosphere and/or in interstellar space are discussed. Figure Optimized structure of a stable H₂OHe²⁺ complex     

Three-dimensional structure of the catalytic domain of the yeast β-(1,3)-glucan transferase Gas1: a molecular modeling investigation

The three-dimensional (3D) structure of the catalytic domain of Gas1p, a protein belonging to the only family of β-(1,3)-glucan transferases so far identified in yeasts and some pathogenic fungi (family GH-72), has been predicted by combining results derived from threading methods, multiple sequence alignments and secondary-structure predictions. The 3D model has allowed the identification of several residues that are predicted to play a crucial role in structural integrity, substrate recognition and catalysis. In particular, the model of the catalytic domain can be useful for designing site-directed mutagenesis experiments and for developing inhibitors of Gas1p enzymatic activity. Figure Three-dimensional models of the Gas1p catalytic domain as predicted using as template 7A3H (PDB code) protein Electronic Supplementary Material Supplementary material is available for this article at     

Neighboring group stabilization by σ-holes

We have used density-functional theory to investigate the neighboring-group stabilization of iodine, arsenic, and phosphorus-centered oxyanion moieties in species such as deprotonated 2-iodoxybenzoic acid (IBX) and its analogs. The magnitudes of different stabilizing effects and further candidates for analogous stabilization are analyzed.      

Coarse-grained modeling of mucus barrier properties

We designed a simple coarse-grained model of the glycocalyx layer, or adhesive mucus layer (AML), covered by mucus gel (luminal mucus layer) using a polymer lattice model and stochastic sampling (replica exchange Monte Carlo) for canonical ensemble simulations. We assumed that mucin MUC16 is responsible for the structural properties of the AML. Other mucins that are much smaller in size and less relevant for layer structure formation were not included. We further assumed that the system was in quasi-equilibrium. For systems with surface coverage and concentrations of model mucins mimicking physiological conditions, we determined the equilibrium distribution of inert nanoparticles within the mucus layers using an efficient replica exchange Monte Carlo sampling procedure. The results show that the two mucus layers penetrate each other only marginally, and the bilayer imposes a strong barrier for nanoparticles, with the AML layer playing a crucial role in the mucus barrier.     

Large scale identification of proteins, mucins, and their O-glycosylation in the endocervical mucus during the menstrual cycle

The mucus filling the human cervical opening blocks the entry to the uterus, but this has to be relative and allow for the sperm to penetrate at ovulation. We studied this mucus, its content of proteins and mucins, and the mucin O-glycosylation in cervical secretions before, during, and after ovulation. Cervical mucosal secretions from 12 subjects were collected, reduced-alkylated, separated with polyacrylamide or agarose/polyacrylamide gel electrophoresis, and stained with silver, Alcian blue, or Coomassie Blue stain. Protein and mucin bands from before and during ovulation were digested and subsequently analyzed by nano-LC-FT-ICR MS and MS/MS. We identified 194 proteins after searches against the NCBI non-redundant protein database and an in-house mucin database. Three gel-forming (MUC5B, MUC5AC, and MUC6) and two transmembrane mucins (MUC16 and MUC1) were identified. For the analysis of mucin O-glycosylation, separated mucins from six individuals were blotted to PVDF membranes, and the O-glycans were released by reductive beta-elimination and analyzed with capillary HPLC-MS and -MS/MS. At least 50 neutral, sialic acid-, and sulfate-containing oligosaccharides were found. An increase of GlcNAc-6GalNAcol Core 2 structures and a relative decrease of NeuAc residues are typical for ovulation, and NeuAc-6GalNAcol and NeuAc-3Gal- epitopes are typical for the non-ovulatory phases. The cervical mucus at ovulation is thus characterized by a relative increase in neutral fucosylated oligosaccharides. This comprehensive characterization of the mucus during the menstrual cycle suggests mucin glycosylation as the major alteration at ovulation, but the relation to the altered physicochemical properties and sperm penetrability is still not understood.     

Molecular dynamics study of the initial stages of catalyzed single-wall carbon nanotubes growth: force field development

Effective force fields for Ni-C interactions developed by Yamaguchi and Maruyama for the formation of metallofullerenes are modified to simulate the catalyzed growth of single-wall carbon nanotubes on Ni_n clusters with n >20, and the reactive empirical bond order Brenner potential for C-C interactions is also revised to include the effect of the metal atoms on such interactions. Figure Force field parameters for carbon-metal interactions obtained from DFT calculations in small clusters.     

CHIH-DFT determination of the molecular structure and IR and UV spectra of solanidine

Solanidine is the steroidal aglycon of some potato glycoalkaloids and a very important precursor for the synthesis of hormones and some pharmacologically active compounds. In this work, we make use of a new chemistry model within Density Functional Theory, called CHIH-DFT, to calculate the molecular structure of solanidine, as well to predict its infrared and ultraviolet spectra. The calculated values are compared with the experimental data available for this molecule as a means of validation of our proposed chemistry model. Figure Molecular structure of solanidine calculated with the CHIH(small) model chemistry     

Intestinal Gel-Forming Mucins Polymerize by Disulfide-Mediated Dimerization of D3 Domains

The mucin 2 glycoprotein assembles into a complex hydrogel that protects intestinal epithelia and houses the gut microbiome. A major step in mucin 2 assembly is further multimerization of preformed mucin dimers, thought to produce a honeycomb-like arrangement upon hydrogel expansion. Important open questions are how multiple mucin 2 dimers become covalently linked to one another and how mucin 2 multimerization compares with analogous processes in related polymers such as respiratory tract mucins and the hemostasis protein von Willebrand factor. Here we report the x-ray crystal structure of the mucin 2 multimerization module, found to form a dimer linked by two intersubunit disulfide bonds. The dimer structure calls into question the current model for intestinal mucin assembly, which proposes disulfide-mediated trimerization of the same module. Key residues making interactions across the dimer interface are highly conserved in intestinal mucin orthologs, supporting the physiological relevance of the observed quaternary structure. With knowledge of the interface residues, it can be demonstrated that many of these amino acids are also present in other mucins and in von Willebrand factor, further indicating that the stable dimer arrangement reported herein is likely to be shared across this functionally broad protein family. The mucin 2 module structure thus reveals the manner by which both mucins and von Willebrand factor polymerize, drawing deep structural parallels between macromolecular assemblies critical to mucosal epithelia and the vasculature.     

MUC genes: mucin or not mucin? That is the question

Mucins are macromolecules lying the cells in contact with external environment and protect the epithelium against constant attacks such as digestive fluids, microorganisms, pollutants, and toxins. Mucins are the main components of mucus and are synthesized and secreted by specialized cells of the epithelium (goblet cells, cells of mucous glands) or non mucin-secreting cells. Human mucin genes show common features: large size of their mRNAs, large nucleotide tandem repeat domains, complex expression both at tissular and cellular level. Since 1987, 21 MUC symbols have been used to designate genes encoding O-glycoproteins containing tandem repeat domains rich in serine, threonine and proline. Some of these genes encode true mucins while others encode non mucin adhesion O-glycoproteins. In this paper, we propose a classification based on sequence similarities and expression areas. Two main families can be distinguished: secreted mucins or gel-forming mucins (MUC2, MUC5AC, MUC5B, MUC6), and membrane-bound mucins (MUC1, MUC3, MUC4, MUC12, MUC17). Muc-deficient mice will provide important models in the study of functional relationships between these two mucin families.     

Quantitative MUC5AC and MUC6 mucin estimations in gastric mucus by a least-squares minimization method

We have determined the molar proportions of the MUC5AC and MUC6 mucus glycoproteins (mucins) in mucus from the normal and pathological human gastric antrum using a least-squares minimization analysis applied to amino acid compositions. We noted that the content of MUC5AC mucin in mucus from individuals without gastroduodenal disease was very high, suggesting that the integrity and barrier properties of the adherent gastric mucus layer are normally maintained by building-block structures formed from this mucin alone. We observed that the molar content of MUC6 mucin doubled (without significance) in mucus from patients with duodenal ulcer, and increased five times (with high significance) in mucus from patients with gastric ulcer, when compared with that in mucus from individuals without gastroduodenal disease.     

Theoretical study of 1-(4-hexylcyclohexyl)-4-isothiocyanatobenzene: molecular properties and spectral characteristics

The mesogenic species 4-(4-hexylcyclohexyl) isothiocyanatobenzene (6CHBT) was studied with density functional theory and molecular mechanics in order to investigate the molecular properties, interactions between dimers and to interpret the IR spectrum. Two types of calculations were performed for model systems containing single and double molecules of 6CHBT. Calculations (involving conformation analysis) for isolated species indicated that the trans isomer, in the equatorial–equatorial conformation, is the most energetically stable. The 6CHBT molecule is polar, with a rather high (4.43 D) dipole moment with negatively charged isothiocyanato (NCS) ligand. The dimer–dimer interaction energies show that the head-to-head configuration (where van der Waals attraction forces play the major role) is the most energetically stable. Vibrational analysis provided detailed assignment of the experimental infra-red (IR) spectrum. Figure Most favorite 6CHBT head to head interaction - ESP mapped to electron density surface Dedication  This paper is dedicated to the memory of Dr. Wacław Witko, who introduced us to research on mesogenic systems.     

The oligomerization of a family of four genetically clustered human gastrointestinal mucins

Mucins are synthesized and secreted by many epithelia. They are complex glycoproteins that offer cytoprotection. In their functional configuration, mucins form oligomers by a biosynthetic process that is poorly understood. A family of four human gastrointestinal mucin genes (MUC2, MUC5AC, MUC5B, and MUC6) is clustered to chromosome 11p15.5. To study oligomerization of these related mucins, we performed metabolic labeling experiments with [35S]amino acids in LS174T cells, and isolated mucin precursors by specific immunoprecipitations that were analyzed on SDS-PAGE. Each of the precursors of MUC2, MUC5AC, MUC5B, and MUC6 formed a single species of disulfide-linked homo-oligomer within 1 h after pulse labeling. Based on apparent molecular masses, these oligomeric precursors were most likely dimers. Inhibition of vesicular RER-to-Golgi transport, with brefeldin A and CCCP, did not affect the dimerization of MUC2 precursors, localizing dimerization to the RER. O-Glycosylation of MUC2 followed dimerization. Inhibition of N- glycosylation by tunicamycin retarded, but did not inhibit, dimerization, indicating that N-glycans play a role in efficient dimerization of MUC2 precursors. Based on sequence homology, the ability of MUC2, MUC5AC, MUC5B and MUC6 to dimerize most likely resides in their C-terminal domains. Thus, the RER-localized dimerization of secretory mucins likely proceeds by similar mechanisms, which is an essential step in the formation of the human gastrointestinal mucus- gels.      

DFT investigation of the 1-octene metathesis reaction mechanism with the Phobcat precatalyst

The productive self-metathesis reaction of 1-octene in the presence of the Phobcat precatalyst [RuCl₂(Phoban-Cy)₂(=CHPh)] using density functional theory was investigated and compared to the Grubbs 1 precatalyst [RuCl₂(PCy₃)₂(=CHPh)]. At the GGA-PW91/DNP level, the geometry optimization of all the participating species and the PES scans of the various activation and catalytic cycles in the dissociative mechanism were performed. The formation of the catalytically active heptylidene species is kinetically and thermodynamically favored, while the formation of trans-tetradecene is thermodynamically favored.