Molecular Evolution of RAMOSA1 (RA1) in Land Plants.

RAMOSA1 (RA1) is a Cys2-His2-type (C2H2) zinc finger transcription factor that controls plant meristem fate and identity and has played an important role in maize domestication. Despite its importance, the origin of RA1 is unknown, and the evolution in plants is only partially understood. In this paper, we present a well-resolved phylogeny based on 73 amino acid sequences from 48 embryophyte species. The recovered tree topology indicates that, during grass evolution, RA1 arose from two consecutive SUPERMAN duplications, resulting in three distinct grass sequence lineages: RA1-like A, RA1-like B, and RA1; however, most of these copies have unknown functions. Our findings indicate that RA1 and RA1-like play roles in the nucleus despite lacking a traditional nuclear localization signal. Here, we report that copies diversified their coding region and, with it, their protein structure, suggesting different patterns of DNA binding and protein-protein interaction. In addition, each of the retained copies diversified regulatory elements along their promoter regions, indicating differences in their upstream regulation. Taken together, the evidence indicates that the RA1 and RA1-like gene families in grasses underwent subfunctionalization and neofunctionalization enabled by gene duplication.

Characteristics of Mycobacterium tuberculosis PtpA interaction and activity on the alpha subunit of human mitochondrial trifunctional protein, a key enzyme of lipid metabolism.

During Mycobacterium tuberculosis (Mtb) infection, the virulence factor PtpA belonging to the protein tyrosine phosphatase family is delivered into the cytosol of the macrophage. PtpA interacts with numerous eukaryotic proteins modulating phagosome maturation, innate immune response, apoptosis, and potentially host-lipid metabolism, as previously reported by our group. In vitro, the human trifunctional protein enzyme (hTFP) is a bona fide PtpA substrate, a key enzyme of mitochondrial ß-oxidation of long-chain fatty acids, containing two alpha and two beta subunits arranged in a tetramer structure. Interestingly, it has been described that the alpha subunit of hTFP (ECHA, hTFPα) is no longer detected in mitochondria during macrophage infection with the virulent Mtb H37Rv. To better understand if PtpA could be the bacterial factor responsible for this effect, in the present work, we studied in-depth the PtpA activity and interaction with hTFPα. With this aim, we performed docking and in vitro dephosphorylation assays defining the P-Tyr-271 as the potential target of mycobacterial PtpA, a residue located in the helix-10 of hTFPα, previously described as relevant for its mitochondrial membrane localization and activity. Phylogenetic analysis showed that Tyr-271 is absent in TFPα of bacteria and is present in more complex eukaryotic organisms. These results suggest that this residue is a specific PtpA target, and its phosphorylation state is a way of regulating its subcellular localization. We also showed that phosphorylation of Tyr-271 can be catalyzed by Jak kinase. In addition, we found by molecular dynamics that PtpA and hTFPα form a stable protein complex through the PtpA active site, and we determined the dissociation equilibrium constant. Finally, a detailed study of PtpA interaction with ubiquitin, a reported PtpA activator, showed that additional factors are required to explain a ubiquitin-mediated activation of PtpA. Altogether, our results provide further evidence supporting that PtpA could be the bacterial factor that dephosphorylates hTFPα during infection, potentially affecting its mitochondrial localization or ß-oxidation activity.

The high identity of the Trypanosoma cruzi Group-I of trans-sialidases points them as promising vaccine immunogens.

Trans-sialidase (TS) superfamily of proteins comprises eight subgroups, being the proteins of Group-I (TS-GI) promising immunogens in vaccine approaches against Trypanosoma cruzi. Strikingly, TS-GI antigenic variability among parasite lineages and their influence on vaccine development has not been previously analyzed. Here, a search in GenBank detects 49 TS-GI indexed sequences, whereas the main infecting human different parasite discrete typing units (DTU) are represented. In silico comparison among these sequences indicate that they share an identity above 92%. Moreover, the antigenic regions (T-cell and B-cell epitopes) are conserved in most sequences or present amino acid substitutions that scarcely may alter the antigenicity. Additionally, since the generic term TS is usually used to refer to different immunogens of this broad family, a further in silico analysis of the TS-GI-derived fragments tested in preclinical vaccines was done to determine the coverage and identity among them, showing that overall amino acid identity of vaccine immunogens is high, but the segment coverage varies widely. Accordingly, strong H-2K, H-2I, and B-cell epitopes are dissimilarly represented among vaccine TS-derived fragments depending on the extension of the TG-GI sequence used. Moreover, bioinformatic analysis detected a set of 150 T-cell strong epitopes among the DTU-indexed sequences that strongly bind human HLA-I supertypes. In all currently reported experimental vaccines based on TS-GI fragments, mapping these 150 epitopes showed that they are moderately represented. However, despite vaccine epitopes do not present all the substitutions observed in the DTUs, these regions of the proteins are equally recognized by the same HLAs.  Interestingly, the predictions regarding global and South American population coverage estimated in these 150 epitopes are similar to the estimations in experimental vaccines when the complete sequence of TS-GI is used as an antigen. In silico prediction also shows that a number of these MHC-I restricted T-cell strong epitopes could be also cross-recognized by HLA-I supertypes and H-2Kb or H-2Kd backgrounds, indicating that these mice may be used to improve and facilitate the development of new TS-based vaccines and suggesting an immunogenic and protective potential in humans. Further molecular docking analyses were performed to strengthen these results. Taken together, different strategies that would cover more or eventually fully of these T-cell and also B-cell epitopes to reach a high level of coverage are considered.

Elucidating paramylon and other carbohydrate metabolism in Euglena gracilis: Kinetic characterization, structure and cellular localization of UDP-glucose pyrophosphorylase.

Many oligo and polysaccharides (including paramylon) are critical in the Euglena gracilis life-cycle and they are synthesized by glycosyl transferases using UDP-glucose as a substrate. Herein, we report the molecular cloning of a gene putatively coding for a UDP-glucose pyrophosphorylase (EgrUDP-GlcPPase) in E. gracilis. After heterologous expression of the gene in Escherichia coli, the recombinant enzyme was characterized structural and functionally. Highly purified EgrUDP-GlcPPase exhibited a monomeric structure, able to catalyze synthesis of UDP-glucose with a Vmax of 3350 U.mg-1. Glucose-1P and UTP were the preferred substrates, although the enzyme also used (with lower catalytic efficiency) TTP, galactose-1P and mannose-1P. Oxidation by hydrogen peroxide inactivated the enzyme, an effect reversed by reduction with dithiothreitol or thioredoxin. The redox process would involve sulfenic acid formation, since no pair of the 7 cysteine residues is close enough in the 3D structure of the protein to form a disulfide bridge. Electrophoresis studies suggest that, after oxidation, the enzyme arranges in many enzymatically inactive structural conformations; which were also detected in vivo. Finally, confocal fluorescence microscopy provided evidence for a cytosolic (mainly in the flagellum) localization of the enzyme.

Human telomerase protein: Understanding how the catalytic activity is suppressed under single substitutions of some conserved residues. A computational study.

The reverse transcriptase domain in telomerase proteins contains the essential conserved residues to catalyze the addition of a single nucleotide to the ends of DNA strands of most eukaryotic cells. In human telomerase protein, mutations in the conserved residues K902, R631, K626, D712, D868, and D869 are known to suppress catalytic activity. To understand these results, a computational model was constructed to simulate a ternary complex consisting of a model of the protein reverse transcriptase domain, a DNA/RNA double helix, an incoming dNTP, and two Mg2+ ions. Three independent Molecular Dynamics Simulations were performed for the wild type and the mutated K902N, R631Q, D712A, D868A, and D869A complexes. Binding Free Energies and alanine-scanning studies were also performed. Using the two-metal-ion mechanism for the nucleotide addition, deviations from the wild type which stop the activity of the human protein, were identified in each mutated complex. The K902N and R631Q mutations might stop the catalytic activity preventing the exit of the pyrophosphate from the catalytic pocket. Additionally, evidence that the same mechanism probably applies to the K626A, R631A, and K902A mutations is presented. For D712A mutation, the pentavalent intermediate state might not form; therefore, the catalytic reaction might not even begin. For the D868A mutation, the O3'-hydroxyl might lose coordination with the Mg ion and the reaction might not either start. Finally, from the limited sampling carried out in this work, we did not obtain any evidence to identify the mechanism by which the D869A mutation cancels the activity of telomerase.

An unusual cysteine VL87 affects the antibody fragment conformations without interfering with the disulfide bond formation.

The Cys residues are almost perfectly conserved in all antibodies. They contribute significantly to the antibody fragment stability. The relevance of two natural contiguous Cys residues of an anti-recombinant human-follicle stimulation hormone (rhFSH) in a format of single-chain variable fragment (scFv) was studied. This scFv contains 5 Cys residues: VH22 and VH92 in the variable heavy chain (VH) and VL23, VL87 and VL88 in the variable light chain (VL). The influence of two unusual contiguous Cys at positions VL87 and VL88 was studied by considering the wild type fragment and mutant variants: VL-C88S, VL-C87S, VL-C87Y. The analysis was carried out using antigen-binding ability measurement by indirect specific ELISA and a detailed molecular modeling that comprises homology methods, long molecular dynamics simulations and docking. We found that VL-C87 affected the antibody fragment stability without interfering with the disulfide bond formation. The effect of mutating the VL-C87 by a usual residue at this position like Tyr caused distant structural changes at the VH region that confers a higher mobility to the VH-CDR2 and VH-CDR3 loops improving the scFv binding to the antigen.

Singular Interaction between an Antimetastatic Agent and the Lipid Bilayer: The Ohmline Case.

SK3 channels are abnormaly expressed in metastatic cells, and Ohmline (OHM), an ether lipid, has been shown to reduce the activity of SK3 channels and the migration capacity of cancer cells. OHM incorporation into the plasma membrane is proposed to dissociate the protein complex formed between SK3 and Orai1, a potassium and a calcium channel, respectively, and would lead to a modification in the lipid environment of both the proteins. Here, we report the synthesis of deuterated OHM that affords the determination, through solid-state NMR, of its entire partitioning into membranes mimicking the SK3 environment. Use of deuterated lipids affords the demonstration of an OHM-induced membrane disordering, which is dose-dependent and increases with increasing amounts of cholesterol (CHOL). Molecular dynamics simulations comfort the disordering action and show that OHM interacts with the carbonyl and phosphate groups of stearoylphosphatidylcholine and sphingomyelin and to a minor extent with CHOL. OHM is thus proposed to remove the CHOL OH moieties away from their main binding sites, forcing a new rearrangement with other lipid groups. Such an interaction takes its origin at the lipid-water interface, but it propagates toward the entire lipid molecules and leads to a cooperative destabilization of the lipid acyl chains, that is, membrane disordering. The consequences of this reorganization of the lipid phases are discussed in the context of the OHM-induced inhibition of SK3 channels.

Structural properties of CHAPS micelles, studied by molecular dynamics simulations.

Detergents are essential tools to study biological membranes, and they are frequently used to solubilize lipids and integral membrane proteins. Particularly the nondenaturing zwitterionic detergent usually named CHAPS was designed for membrane biochemistry and integrates the characteristics of the sulfobetaine-type detergents and bile salts. Despite the available experimental data little is known about the molecular structure of its micelles. In this work, molecular dynamics simulations were performed to study the aggregation in micelles of several numbers of CHAPS (≤ 18) starting from a homogeneous water dilution. The force field parameters to describe the interactions of the molecule were developed and validated. After 50 ns of simulation almost all the systems result in the formation of stable micelles. The molecular shape (gyration radii, volume, surface) and the molecular structure (RDF, salt bridges, H-bonds, SAS) of the micelles were characterized. It was found that the main interactions that lead to the stability of the micelles are the electrostatic ones among the polar groups of the tails and the OH's from the ring moiety. Unlike micelles of other compounds, CHAPS show a grainlike heterogeneity with hydrophobic micropockets. The results are in complete agreement with the available experimental information from NMR, TEM, and SAXS studies, allowing the modeling of the molecular structure of CHAPS micelles. Finally, we hope that the new force field parameters for this detergent will be a significant contribution to the knowledge of such an interesting molecule.

Molecular dynamics study of the interaction of arginine with phosphatidylcholine and phosphatidylethanolamine bilayers.

In this work, the differential interaction of zwitterionic arginines with fully hydrated dimyristoylphosphatidylcholine (DMPC) and dimyristoylphosphatidylethanolamine (DMPE) bilayers was analyzed by molecular dynamics simulations. In both systems, arginine binds to lipids with the carboxylate moiety oriented toward the aqueous phase, in agreement with previous experimental determinations of ζ potential of DMPC and DMPE liposomes. The guanidinium groups are found at different depths within the bilayers indicating that some arginines are buried, especially in DMPE. We observe, in the DMPE system, that the strongest interaction occurs between the guanidinium group and the carbonyl oxygen of the lipid. In the case of DMPC membranes, the strongest interaction is found between the guanidinium groups of the arginines and the phosphate groups of the lipids. Unexpectedly, arginine zwitterions are stabilized through the creation of hydrogen bonds (HB), either with water or with polar groups of the lipids. The mechanisms of interaction seem to be different in both membranes. In DMPE bilayers, arginines insert by breaking the inner HB network of the polar head groups, consequently increasing the occupied area per lipid molecule. In the DMPC bilayers the arginines insert by replacing the already present water molecules within the membrane, without significant effects on the area per lipid.

Similar structures but different mechanisms: Prediction of FABPs-membrane interaction by electrostatic calculation.

The role of fatty acid binding proteins as intracellular fatty acid transporters may require their direct interaction with membranes. In this way different mechanisms have been previously characterized through experimental studies suggesting different models for FABPs-membrane association, although the process in which the molecule adsorbs to the membrane remains to be elucidated. To estimate the importance of the electrostatic energy in the FABP-membrane interaction, we computationally modeled the interaction of different FABPs with both anionic and neutral membranes. Free Electrostatic Energy of Binding (dE), was computed using Finite Difference Poisson Boltzmann Equation (FDPB) method as implemented in APBS (Adaptive Poisson Boltzmann Solver). Based on the computational analysis, it is found that recruitment to membranes is facilitated by non-specific electrostatic interactions. Also energetic analysis can quantitatively differentiate among the mechanisms of membrane association proposed and determinate the most energetically favorable configuration for the membrane-associated states of different FABPs. This type of calculations could provide a starting point for further computational or experimental analysis.

Structure and dynamics of nano-sized raft-like domains on the plasma membrane.

Cell membranes are constitutively composed of thousands of different lipidic species, whose specific organization leads to functional heterogeneities. In particular, sphingolipids, cholesterol and some proteins associate among them to form stable nanoscale domains involved in recognition, signaling, membrane trafficking, etc. Atomic-detail information in the nanometer/second scale is still elusive to experimental techniques. In this context, molecular simulations on membrane systems have provided useful insights contributing to bridge this gap. Here we present the results of a series of simulations of biomembranes representing non-raft and raft-like nano-sized domains in order to analyze the particular structural and dynamical properties of these domains. Our results indicate that the smallest (5 nm) raft domains are able to preserve their distinctive structural and dynamical features, such as an increased thickness, higher ordering, lower lateral diffusion, and specific lipid-ion interactions. The insertion of a transmembrane protein helix into non-raft, extended raft-like, and raft-like nanodomain environments result in markedly different protein orientations, highlighting the interplay between the lipid-lipid and lipid-protein interactions.

Salt induced asymmetry in membrane simulations by partial restriction of ionic motion.

The specific ionic composition differs considerably at both sides of biological membranes and specific lipid/electrolyte interactions may be essential for their structure, stability and function. Hence, explicit consideration of the ionic asymmetry is important to achieve an accurate description of lipid bilayers. Molecular dynamics simulations have proven to be a reliable tool to study biomembranes at atomic detail. Nevertheless, the use of periodic boundary conditions allows ions to diffuse rapidly and reach both sides of the bilayer. Therefore, ad hoc simulation schemes have to be applied to take into account ionic asymmetry. In this work we present a simple implementation to overcome this problem. A more realistic description of the biomembranes can be achieved by partially restricting the ionic motion in the direction normal to the membrane within a region of the space near to only one of the leaflets. This creates two different situations: one leaflet is highly exposed to ions while the second one can be completely or partially depleted of them. Comparison between this new method and control simulations performed using a previously proposed approach consisting of a double-membrane setup yielded an excellent agreement with a speed-up of nearly 60%. The performance of the method with different ionic species is explored and remaining limitations are examined.

Inhibition of alpha-synuclein fibrillization by dopamine is mediated by interactions with five C-terminal residues and with E83 in the NAC region.

The interplay between dopamine and alpha-synuclein (AS) plays a central role in Parkinson's disease (PD). PD results primarily from a severe and selective devastation of dopaminergic neurons in substantia nigra pars compacta. The neuropathological hallmark of the disease is the presence of intraneuronal proteinaceous inclusions known as Lewy bodies within the surviving neurons, enriched in filamentous AS. In vitro, dopamine inhibits AS fibril formation, but the molecular determinants of this inhibition remain obscure. Here we use molecular dynamic (MD) simulations to investigate the binding of dopamine and several of its derivatives onto conformers representative of an NMR ensemble of AS structures in aqueous solution. Within the limitations inherent to MD simulations of unstructured proteins, our calculations suggest that the ligands bind to the (125)YEMPS(129) region, consistent with experimental findings. The ligands are further stabilized by long-range electrostatic interactions with glutamate 83 (E83) in the NAC region. These results suggest that by forming these interactions with AS, dopamine may affect AS aggregation and fibrillization properties. To test this hypothesis, we investigated in vitro the effects of dopamine on the aggregation of mutants designed to alter or abolish these interactions. We found that point mutations in the (125)YEMPS(129) region do not affect AS aggregation, which is consistent with the fact that dopamine interacts non-specifically with this region. In contrast, and consistent with our modeling studies, the replacement of glutamate by alanine at position 83 (E83A) abolishes the ability of dopamine to inhibit AS fibrillization.

On the oligomeric state of DJ-1 protein and its mutants associated with Parkinson Disease. A combined computational and in vitro study.

Mutations in the DJ-1 protein are present in patients suffering from familial Parkinson disease. Here we use computational methods and biological assays to investigate the relationship between DJ-1 missense mutations and the protein oligomeric state. Molecular dynamics calculations suggest that: (i) the structure of DJ-1 wild type (WT) in aqueous solution, in both oxidized and reduced forms, is similar to the crystal structure of the reduced form; (ii) the Parkinson disease-causing M26I variant is structurally similar to the WT, consistent with the experimental evidence showing the protein is a dimer as WT; (iii) R98Q is structurally similar to the WT, consistent with the fact that this is a physiological variant; and (iv) the L166P monomer rapidly evolves toward a conformation significantly different from WT, suggesting a change in its ability to oligomerize. Our combined computational and experimental approach is next used to identify a mutant (R28A) that, in contrast to L166P, destabilizes the dimer subunit-subunit interface without significantly changing secondary structure elements.