Glycoprotein 63 (gp63) genes show gene conversion and reveal the evolution of Old World Leishmania

Species of the subgenus Leishmania (Leishmania) cause the debilitating disease leishmaniasis on four continents. Species grouped within the Leishmania donovani complex cause visceral leishmaniasis, a life-threatening disease, often associated with poverty, and affecting some 0.5 million people each year. The Leishmania glycoprotein GP63, or major surface protease, is a metalloprotease involved in parasite survival, infectivity and virulence. Here, we show that evolution of the gp63 multigene family is influenced by mosaic or fragmental gene conversion. This is a major evolutionary force for both homogenisation and for generating diversity, even in the absence of sexual reproduction. We propose here that the high GC content at the third codon position in the gp63 family of Old World Leishmania may be higher in multicopy regions, under the biased gene conversion model, because increased copy numbers may lead to increased rates of recombination. We confirm that one class of gp63 genes with an extended 3'end signal, gp63(EXT), reveals genetic groups within the complex and gives insights into evolution and host associations. Gp63(EXT) genes can also provide the basis for rapid and reliable genotyping of strains in the L. donovani complex. Our results confirmed that a more stringent definition of Leishmania infantum is required and that the species Leishmania archibaldi should be suppressed.     

Processing and trafficking of Leishmania mexicana GP63. Analysis using GP18 mutants deficient in glycosylphosphatidylinositol protein anchoring

GPI8 is a clan CD, family C13 cysteine protease and the catalytic core of the GPI-protein transamidase complex. In Leishmania mexicana, GPI8 is nonessential, and Deltagpi8 mutants lack the GPI-anchored metalloprotease GP63, which is the major surface protein of promastigotes. We have identified the active site histidine and cysteine residues of leishmanial GPI8 and generated Deltagpi8 lines expressing modified GPI8 proteins. This has allowed us to study the processing and trafficking of GP63 in wild type and Deltagpi8 mutants. We show using pulse-chase labeling that in Deltagpi8 non-GPI-anchored GP63 was glycosylated and secreted without further processing from the cell with a t(12) of 120 min. This secretion was prevented by growth of cells in the presence of tunicamycin, indicating that glycosylation is necessary for secretion of non-GPI-anchored proteins. In contrast, in wild type cells the majority of GP63 was rapidly glycosylated, GPI-anchored, and trafficked to the surface with defined processing intermediate forms. Tunicamycin inhibited glycosylation but did not prevent GPI anchor addition or trafficking. These results show that GPI-anchored and unanchored GP63 are trafficked via different pathways. In addition, the balance between GPI anchor addition and secretion of GP63 in Leishmania can vary depending on the activity of the GPI-protein transamidase, which has implications for the host-parasite interaction.     

Molecular dynamics of the full-length p53 monomer

The p53 protein is frequently mutated in a very large proportion of human tumors, where it seems to acquire gain-of-function activity that facilitates tumor onset and progression. A possible mechanism is the ability of mutant p53 proteins to physically interact with other proteins, including members of the same family, namely p63 and p73, inactivating their function. Assuming that this interaction might occurs at the level of the monomer, to investigate the molecular basis for this interaction, here, we sample the structural flexibility of the wild-type p53 monomeric protein. The results show a strong stability up to 850 ns in the DNA binding domain, with major flexibility in the N-terminal transactivations domains (TAD1 and TAD2) as well as in the C-terminal region (tetramerization domain). Several stable hydrogen bonds have been detected between N-terminal or C-terminal and DNA binding domain, and also between N-terminal and C-terminal. Essential dynamics analysis highlights strongly correlated movements involving TAD1 and the proline-rich region in the N-terminal domain, the tetramerization region in the C-terminal domain; Lys120 in the DNA binding region. The herein presented model is a starting point for further investigation of the whole protein tetramer as well as of its mutants.     

Differential stability of proteolytically active and inactive recombinant metalloproteinase in Chinese hamster ovary cells

Stability of heterologous protein expression during production is critical for regulatory approval of vaccine and therapeutic products. Leishmania GP63, a zinc metalloproteinase that is a potential vaccine candidate, has been expressed on the surface of Chinese hamster ovary (CHO) cells. Flow cytometry was used to follow the stability of GP63 expression. Expression of proteolytically active GP63 (GP63WT) was unstable whether or not methotrexate (MTX) selection was maintained. In contrast, expression of an active site mutant (GP63E265D) was stable under MTX selection. In the absence of selection, the decline in GP63E265D expression was more gradual than the loss of GP63WT expression. Different molecular mechanisms accounted for these losses and resulted in higher growth rate nonproducer populations. A dynamic population model was used to calculate the conversion rates of GP63WT producers to nonproducers. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 53: 594-600, 1997.     

Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations. Proteins 31:116–127, 1998. © 1998 Wiley-Liss, Inc.     

Leishmania parasites act as a Trojan horse that paralyzes the translation system of host macrophages

GP63 is an abundant GPI-anchored surface metalloprotease of Leishmania. Jaramillo et al. (2011) show that GP63 manipulates the translation system of host macrophages by cleaving mTOR, which leads to 4E-BP1 dephosphorylation. This study pioneers the observation that Leishmania parasites metabolically paralyze their host cells using an elegant translation shutoff mechanism.     

Allostery in Orai1 binding to calmodulin revealed from conformational thermodynamics

Here, we study microscopic mechanism of complex formation between Ca²⁺-bound calmodulin (holoCaM) and Orai1 that regulates Ca²⁺-dependent inactivation process in eukaryotic cells. We compute conformational thermodynamic changes in holoCaM with respect to complex of Orai1 bound to C-terminal domain of holoCaM using histograms of dihedral angles of the proteins over trajectories from molecular dynamics simulations. Our analysis shows that the N-terminal domain residues L4, T5, Q41, N42, T44 and E67 of holoCaM get destabilized and disordered due to Orai1 binding to C-terminal domain of calmodulin affect the N-terminal domain residues. Among these residues, polar T44, having maximum destabilization and disorder via backbone fluctuations, shows the largest change in solvent exposure. This suggests that N-terminal domain is allosterically regulated via T44 by the binding of Orai1 to the C-terminal domain.     

Structure of G57W mutant of human γS-crystallin and its involvement in cataract formation

A recently identified mutant of human γS-crystallin, G57W is associated with dominant congenital cataracts, the familial determinate of childhood blindness worldwide. To investigate the structural and functional changes that mediate the effect of this cataract-related mutant to compromise eye lens transparency and cause lens opacification in children, we recently reported complete sequence-specific resonance assignments of γS-G57W using a suite of heteronuclear NMR experiments. As a follow up, we have determined the 3D structure of γS-G57W and studied its conformational dynamics by solution NMR spectroscopy. Our structural dynamics results reveal greater flexibility of the N-terminal domain, which undergoes site-specific structural changes to accommodate W57, than its C-terminal counterpart. Our structural inferences that the unusual solvent exposure of W57 is associated with rearrangement of the N-terminal domain suggest an efficient pathway for increased aggregation in γS-G57W and illuminates the molecular dynamics underlying cataractogenic aggregation of lens crystallins in particular and aggregation of proteins in general.     

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.     

Specific recognition and cleavage of galectin-3 by Leishmania major through species-specific polygalactose epitope

Lipophosphoglycan is a major surface molecule of Leishmania, protozoa parasites, which are the causative agents of leishmaniasis, a disease that annually afflicts millions of people worldwide. The oligosaccharide structures of lipophosphoglycan varies among species, and epitopes of these species-specific oligosaccharides are suggested to be implicated in the interaction of Leishmania with macrophages as well as species-specific tissue tropism observed in leishmaniasis. The recognition of the species-specific variation of oligosaccharides is likely to be mediated by host carbohydrate-binding proteins, lectins, but the identities of the lectins remain elusive. Galectin-3 is a mammalian soluble beta-galactoside-binding lectin and is expressed in macrophages, dendritic cells, and keratinocytes, as well as fibroblasts, all of which are present in the site of Leishmania infection. In this paper, we found that galectin-3 binds to lipophosphoglycan of Leishmania major but not to those of Leishmania donovani through L. major-specific polygalactose epitopes. Association of galectin-3 with L. major led to the cleavage of galectin-3, resulting in truncated galectin-3 containing the C-terminal lectin domain but lacking the N-terminal domain implicated in lectin oligomerization. This cleavage was inhibited by the galectin-3 antagonist lactose, as well as 1,10-ortho-phenanthroline, suggesting that galectin-3 is cleaved by zinc metalloproteases after its binding to lipophosphoglycans. The modulation of various innate immunity reactions by galectin-3 is affected by its oligomerization; therefore, we propose the L. major-specific truncation of galectin-3 may contribute to the species-specific immune responses induced by Leishmania.     

Prefusion spike protein conformational changes are slower in SARS-CoV-2 than in SARS-CoV-1

Within the last 2 decades, severe acute respiratory syndrome coronaviruses 1 and 2 (SARS-CoV-1 and SARS-CoV-2) have caused two major outbreaks; yet, for reasons not fully understood, the coronavirus disease 2019 pandemic caused by SARS-CoV-2 has been significantly more widespread than the 2003 SARS epidemic caused by SARS-CoV-1, despite striking similarities between these two viruses. The SARS-CoV-1 and SARS-CoV-2 spike proteins, both of which bind to host cell angiotensin-converting enzyme 2, have been implied to be a potential source of their differential transmissibility. However, the mechanistic details of prefusion spike protein binding to angiotensin-converting enzyme 2 remain elusive at the molecular level. Here, we performed an extensive set of equilibrium and nonequilibrium microsecond-level all-atom molecular dynamics simulations of SARS-CoV-1 and SARS-CoV-2 prefusion spike proteins to determine their differential dynamic behavior. Our results indicate that the active form of the SARS-CoV-2 spike protein is more stable than that of SARS-CoV-1 and the energy barrier associated with the activation is higher in SARS-CoV-2. These results suggest that not only the receptor-binding domain but also other domains such as the N-terminal domain could play a crucial role in the differential binding behavior of SARS-CoV-1 and SARS-CoV-2 spike proteins.     

MARCKS-related protein (MRP) is a substrate for the Leishmania major surface protease leishmanolysin (gp63).

Myristoylated alanine-rich C kinase substrate (MARCKS) and MARCKS-related protein (MRP; MacMARCKS) are protein kinase C substrates in diverse cell types. Activation of murine macrophages by cytokines increases MRP expression, but infection with Leishmania promastigotes during activation results in MRP depletion. We therefore examined the effect of Leishmania major LV39 on recombinant MRP. Both live promastigotes and a soluble fraction of LV39 lysates degraded MRP to yield lower molecular weight fragments. Degradation was independent of MRP myristoylation and was inhibited by protein kinase C-dependent phosphorylation of MRP. MRP was similarly degraded by purified leishmanolysin (gp63), a Leishmania surface metalloprotease. Degradation was evident at low enzyme/substrate ratios, over a broad pH range, and was inhibited by 1,10-phenanthroline and by a hydroxamate dipeptide inhibitor of leishmanolysin. Using mass spectrometric analysis, cleavage was shown to occur within the effector domain of MRP between Ser(92) and Phe(93), in accordance with the substrate specificity of leishmanolysin. Moreover, an MRP construct in which the effector domain had been deleted was resistant to cleavage. Thus, Leishmania infection may result in leishmanolysin-dependent hydrolysis of MRP, a major protein kinase C substrate in macrophages.     

Leishmania repression of host translation through mTOR cleavage is required for parasite survival and infection

The protozoan parasite Leishmania alters the activity of its host cell, the macrophage. However, little is known about the effect of Leishmania infection on host protein synthesis. Here, we show that the Leishmania protease GP63 cleaves the mammalian/mechanistic target of rapamycin (mTOR), a serine/threonine kinase that regulates the translational repressor 4E-BP1. mTOR cleavage results in the inhibition of mTOR complex 1 (mTORC1) and concomitant activation of 4E-BP1 to promote Leishmania proliferation. Consistent with these results, pharmacological activation of 4E-BPs with rapamycin, results in a dramatic increase in parasite replication. In contrast, genetic deletion of 4E-BP1/2 reduces parasite load in macrophages ex vivo and decreases susceptibility to cutaneous leishmaniasis in vivo. The parasite resistant phenotype of 4E-BP1/2 double-knockout mice involves an enhanced type I IFN response. This study demonstrates that Leishmania evolved a survival mechanism by activating 4E-BPs, which serve as major targets for host translational control.     

Comparison of SARS and NL63 papain-like protease binding sites and binding site dynamics: inhibitor design implications

The human severe acute respiratory syndrome coronavirus (SARS-CoV) and the NL63 coronaviruses are human respiratory pathogens for which no effective antiviral treatment exists. The papain-like cysteine proteases encoded by the coronavirus (SARS-CoV: PLpro; NL63: PLP1 and PLP2) represent potential targets for antiviral drug development. Three recent inhibitor-bound PLpro structures highlight the role of an extremely flexible six-residue loop in inhibitor binding. The high binding site plasticity is a major challenge in computational drug discovery/design efforts. From conventional molecular dynamics and accelerated molecular dynamics (aMD) simulations, we find that with conventional molecular dynamics simulation, PLpro translationally samples the open and closed conformation of BL2 loop on a picosecond-nanosecond timescale but does not reproduce the peptide bond inversion between loop residues Tyr269 and Gln270 that is observed on inhibitor GRL0617 binding. Only aMD simulation, starting from the closed loop conformation, reproduced the 180° ?-ψ dihedral rotation back to the open loop state. The Tyr-Gln peptide bond inversion appears to involve a progressive conformational change of the full loop, starting at one side, and progressing to the other. We used the SARS-CoV apo X-ray structure to develop a model of the NL63-PLP2 catalytic site. Superimposition of the PLP2 model on the PLpro X-ray structure identifies binding site residues in PLP2 that contribute to the distinct substrate cleavage site specificities between the two proteases. The topological and electrostatic differences between the two protease binding sites also help explain the selectivity of non-covalent PLpro inhibitors.     

Human microtubule affinity-regulating kinase 4 is stable at extremes of pH

MAP/microtubule affinity-regulating kinase 4 (MARK4) is a member of adenosine monophosphate-activated protein kinases, directly associated with cancer and neurodegenerative diseases. Here, we have cloned, expressed, and purified two variants of MARK4 [the kinase domain (MARK4-F2), and kinase domain along with 59 N-terminal residues (MARK4-F1)] and compared their stability at varying pH range. Structural and functional changes were observed by incubating both forms of MARK4 in buffers of different pH. We measured the secondary structure of MARK4 using circular dichroism and tertiary structure by measuring intrinsic fluorescence and absorbance properties along with the size of proteins by dynamic light scattering. We observed that at extremes of pH (below pH 3.5 and above pH 9.0), MARK4 is quite stable. However, a remarkable aggregate formation was observed at intermediate pH (between pH 3.5 and 9.0). To further validate this result, we have modeled both forms of MARK4 and performed molecular dynamics simulation for 15 ns. The spectroscopic observations are in excellent agreement with the findings of molecular dynamics simulation. We also performed ATPase activity at varying pH and found a significant correlation of structure of MARK4 with its enzyme activity. It is interesting to note that both forms of MARK4 are showing a similar pattern of structure changes with reference to pH.     

The comparative fine structure and surface glycoconjugate expression of three life stages of Leishmania major

The cellular ultrastructure and surface glycoconjugate expression of three life stages of Leishmania major were compared. Noninfective logarithmic phase promastigotes (LP) are immature cells bearing a thin cell coat, short flagellum, small and empty flagellar pocket, and a loose cytoplasm filled with profiles of ER and large Golgi complex. LP also contain subpopulations of maturing cells containing less ER and Golgi and synthesizing cytoplasmic granules of different size, number, and electron-density. Infective or metacyclic promastigotes (MP) are fully differentiated nondividing forms with a thickened, prominent cell coat, long flagellum, distended flagellar pocket filled with secretory material, and few cytoplasmic organelles other than abundant electron-dense granules. Tissue amastigotes also contain electron-dense cytoplasmic granules, their flagellar pockets are also enlarged and contain secretory material, but they lack a discernable cell coat. Immunogold labeling of GP63 on the cell surface was extensive only on amastigotes. Promastigote GP63 appeared to be masked by the presence of a densely packed lipophosphoglycan (LPG) coat which was extensively labeled on the entire surface of MP and LP. An elongated, developmentally modified form of LPG was abundantly labeled only on MP. LPG was poorly labeled on amastigotes, arguing that the promastigote cell coat is a stage-specific structure which is lost during intracellular transformation.     

Insights from the analysis of a predicted model of gp63 in Leishmania donovani

Leishmaniasis is a protozoal disease of human that occurs in most parts of the world. By considering the progress of bioinformatics in molecular modeling, major surface glycoprotein of Leishmania donovani (gp63) structure was modeled using homology modeling with high accuracy based on the X-ray crystal structure of the Leishmania major gp63 as a template, and then analyzed 3D structure of gp63 which can reveal exact facts about its structure and interaction. The objective of this study was to find folding and three dimensional structure of the gp63 as potent antigen for human. In this project, we applied the theory of evolution method, including comparative modeling and threading. This study presented a simple protocol for rapid and precise finding 3D structure of gp63 and investigation of its structural properties. The translated amino acid sequence showed that Leishmania donovani gp63 contains 590 amino acids precursor protein consisting of an NH₂-terminal signal peptide of 39 amino acids for membrane targeting, a pro region of 48 amino acids, the mature protein of 478 amino acids containing glycosylation and putative catalytic sites, and a COOH-terminal signal peptide of 25 amino acids for GPI attachment. Based on our model, the protein consists of three domains: the N-terminal, central and C-terminal domains. Additionally, these results could guide future structure-function analyses of gp63 protein.     

NMR structure of the p63 SAM domain and dynamical properties of G534V and T537P pathological mutants, identified in the AEC syndrome

The p63 protein is crucial for epidermal development, and its mutations cause the extrodactyly ectodermal dysplasia and cleft lip/palate syndrome. The three-dimensional solution structure of the p63 sterile α-motif (SAM) domain (residues 505–579), a region crucial to explaining the human genetic disease ankyloblepharon-ectodermal dysplasia-clefting syndrome (AEC), has been determined by nuclear magnetic resonance spectroscopy. The structure indicates that the domain is a monomer with the characteristic five-helix bundle topology observed in other SAM domains. It includes five tightly packed helices with an extended hydrophobic core to form a globular and compact structure. The dynamics of the backbone and the global correlation time of the molecule have also been investigated and compared with the dynamical properties obtained through molecular dynamics simulation. Attempts to purify the pathological G534V and T537P mutants, originally identified in AEC, were not successful because of the occurrence of unspecific proteolytic degradation of the mutated SAM domains. Analysis of the structural dynamic properties of the G534V and T537P mutants through molecular dynamics simulation and comparison with the wild type permits detection of differences in the degree of free-dom of individual residues and discussion of the possible causes for the pathology.     

Molecular cloning and characterization of a novel repeat-containing Leishmania major gene, ppg1, that encodes a membrane-associated form of proteophosphoglycan with a putative glycosylphosphatidylinositol anchor.

Leishmania parasites secrete a variety of proteins that are modified by phosphoglycan chains structurally similar to those of the cell surface glycolipid lipophosphoglycan. These proteins are collectively called proteophosphoglycans. We report here the cloning and sequencing of a novel Leishmania major proteophosphoglycan gene, ppg1. It encodes a large polypeptide of approximately 2300 amino acids. The N-terminal domain of approximately 70 kDa exhibits 11 imperfect amino acid repeats that show some homology to promastigote surface glycoproteins of the psa2/gp46 complex. The large central domain apparently consists exclusively of approximately 100 repetitive peptides of the sequence APSASSSSA(P/S)SSSSS(+/-S). Gene fusion experiments demonstrate that these peptide repeats are the targets of phosphoglycosylation in Leishmania and that they form extended filamentous structures reminiscent of mammalian mucins. The C-terminal domain contains a functional glycosylphosphatidylinositol anchor addition signal sequence, which confers cell surface localization to a normally secreted Leishmania acid phosphatase, when fused to its C terminus. Antibody binding studies show that the ppg1 gene product is phosphoglycosylated by phosphoglycan repeats and cap oligosaccharides. In contrast to previously characterized proteophosphoglycans, the ppg1 gene product is predominantly membrane-associated and it is expressed on the promastigote cell surface. Therefore this membrane-bound proteophosphoglycan may be important for direct host-parasite interactions.