Unfolding and Aggregation Pathways of Variable Domains from Immunoglobulin Light Chains.

Light chain amyloidosis is the most common form of systemic amyloidosis. This disease is caused by the formation and deposition of amyloid fibers made from immunoglobulin light chains. Environmental conditions such as pH and temperature can affect protein structure and induce the development of these fibers. Several studies have shed light on the native state, stability, dynamics, and final amyloid state of these proteins; however, the initiation process and the fibril formation pathway remain poorly understood structurally and kinetically. To study this, we analyzed the unfolding and aggregation process of the 6aJL2 protein under acidic conditions, with temperature changes, and upon mutation, using biophysical and computational techniques. Our results suggest that the differences in amyloidogenicity displayed by 6aJL2 under these conditions are caused by traversing different aggregation pathways, including unfolded intermediates and the formation of oligomers.

Clinical Manifestations, Mutational Analysis, and Immunological Phenotype in Patients with RAG1/2 Mutations: First Cases Series from Mexico and Description of Two Novel Mutations.

Mutations in recombinase activating genes 1 and 2 (RAG1/2) result in human severe combined immunodeficiency (SCID). The products of these genes are essential for V(D)J rearrangement of the antigen receptors during lymphocyte development. Mutations resulting in null-recombination activity in RAG1 or RAG2 are associated with the most severe clinical and immunological phenotypes, whereas patients with hypomorphic mutations may develop leaky SCID, including Omenn syndrome (OS). A group of previously unrecognized clinical phenotypes associated with granulomata and/or autoimmunity have been described as a consequence of hypomorphic mutations. Here, we present six patients from unrelated families with missense variants in RAG1 or RAG2. Phenotypes observed in these patients ranged from OS to severe mycobacterial infections and granulomatous disease. Moreover, we report the first evidence of two variants that had not been associated with immunodeficiency. This study represents the first case series of RAG1- or RAG2-deficient patients from Mexico and Latin America.

Revealing the Structural Contributions to Thermal Adaptation of the TATA-Box Binding Protein: Molecular Dynamics and QSPR Analyses.

The TATA-box binding protein (TBP) is an important element of the transcription machinery in archaea and eukaryotic organisms. TBP is expressed in organisms adapted to different temperatures, indicating a robust structure, and experimental studies have shown that the mid-unfolding temperature (Tm) of TBP is directly correlated with the optimal growth temperature (OGT) of the organism. To understand which are the relevant structural requirements for its stability, we present the first structural and dynamic computational study of TBPs, combining molecular dynamics (MD) simulations and a quantitative structure-property relationship (QSPR) over a set of TBPs of organisms adapted to different temperatures. We found that the main structural properties of TBP used to adapt to high temperatures are an increase in the ease of desolvation of charged residues at the surface, an increase in the local resiliency, the presence of Leu clusters in the protein core, and an increase in the loss of hydrophobic packing in the N-terminal subdomain. In view of our results, we consider that TBP is a good model to study thermal adaptation, and our analysis opens the possibility of performing protein engineering on TBPs to study transcription at high or low temperatures.

Exploring the Role of Phenylalanine Residues in Modulating the Flexibility and Topography of the Active Site in the Peroxygenase Variant PaDa-I.

Unspecific peroxygenases (UPOs) are fungal heme-thiolate enzymes able to catalyze a wide range of oxidation reactions, such as peroxidase-like, catalase-like, haloperoxidase-like, and, most interestingly, cytochrome P450-like. One of the most outstanding properties of these enzymes is the ability to catalyze the oxidation a wide range of organic substrates (both aromatic and aliphatic) through cytochrome P450-like reactions (the so-called peroxygenase activity), which involves the insertion of an oxygen atom from hydrogen peroxide. To catalyze this reaction, the substrate must access a channel connecting the bulk solution to the heme group. The composition, shape, and flexibility of this channel surely modulate the catalytic ability of the enzymes in this family. In order to gain an understanding of the role of the residues comprising the channel, mutants derived from PaDa-I, a laboratory-evolved UPO variant from Agrocybe aegerita, were obtained. The two phenylalanine residues at the surface of the channel, which regulate the traffic towards the heme active site, were mutated by less bulky residues (alanine and leucine). The mutants were experimentally characterized, and computational studies (i.e., molecular dynamics (MD)) were performed. The results suggest that these residues are necessary to reduce the flexibility of the region and maintain the topography of the channel.

In silico analysis of missense mutations in exons 1-5 of the F9 gene that cause hemophilia B.

BACKGROUND: Missense mutations in the first five exons of F9, which encodes factor FIX, represent 40% of all mutations that cause hemophilia B. To address the ongoing debate regarding in silico identification of disease-causing mutations at these exons, we analyzed 215 missense mutations from www.factorix.org using six in silico prediction tools, which are the most common used programs for analysis prediction of impact of mutations on the protein structure and function, with further advantage of using similar approaches. We developed different algorithms to integrate multiple predictions from such tools. In order to approach a structural analysis on FIX we performed a modeling of five selected pathogenic mutations. RESULTS: SIFT, PolyPhen-2 HumDiv, SNAP2, and MutationAssessor were the most successful in identifying true non-causative and causative mutations. A proposed function integrating these algorithms (wgP4) was the most sensitive (90.1%), specific (22.6%), and accurate (87%) than similar functions, and identified 187 variants as deleterious. Clinical phenotype was significantly associated with predicted causative mutations at all five exons. However, PolyPhen-2 HumDiv was more successful in linking clinical severity to specific exons, while functions that integrate 4-6 predictions were more successful in linking phenotype to genotypes at the light chain (exons 3-5). The most important value of integrating multiple predictions is the inclusion of scores derived from different approaches. Modeling of protein structure showed the effects of pathogenic nsSNPs on structure and function of FIX. CONCLUSIONS: A simple function that integrates information from different in silico programs yields the best prediction of mutated phenotypes. However, the specificity, sensitivity, and accuracy of genotype-phenotype predictions depend on specific characteristics of the protein domain and the disease of interest as we validated by the structural analysis of selected pathogenic F9 mutations. The proposed function integrating algorithm (wgP4) might be useful for the analysis of nsSNPs impact on other genes.

Different Dynamics in 6aJL2 Proteins Associated with AL Amyloidosis, a Conformational Disease.

Light-chain amyloidosis (AL) is the most common systemic amyloidosis and is caused by the deposition of mainly insoluble immunoglobulin light chain amyloid fibrils in multiple organs, causing organ failure and eventually death. The germ-line λ6a has been implicated in AL, where a single point mutant at amino acid 24 (6aJL2-R24G) has been observed in around 25% of patient samples. Structural analysis has shown only subtle differences between both proteins; nevertheless, 6aJL2-R24G is more prone to form amyloid fibrils. To improve our understanding of the role of protein flexibility in amyloid fibril formation, we have used a combination of solution nuclear magnetic resonance spectroscopy and molecular dynamics simulations to complement the structural insight with dynamic knowledge. Fast timescale dynamics (ps-ns) were equivalent for both proteins, but suggested exchange events for some residues. Even though most of the intermediate dynamics (µs-ms) occurred at a similar region for both proteins, the specific characteristics are very different. A minor population detected in the dispersion experiments could be associated with the formation of an off-pathway intermediate that protects from fiber formation more efficiently in the germ-line protein. Moreover, we found that the hydrogen bond patterns for both proteins are similar, but the lifetime for the mutant is significantly reduced; as a consequence, there is a decrease in the stability of the tertiary structure that extends throughout the protein and leads to an increase in the propensity to form amyloid fibers.

Localized conformational changes trigger the pH-induced fibrillogenesis of an amyloidogenic λ light chain protein.

Solvent conditions modulate the expression of the amyloidogenic potential of proteins. In this work the effect of pH on the fibrillogenic behavior and the conformational properties of 6aJL2, a model protein of the highly amyloidogenic variable light chain λ6a gene segment, was examined. Ordered aggregates showing the ultrastructural and spectroscopic properties observed in amyloid fibrils were formed in the 2.0-8.0 pH range. At pH <3.0 a drastic decrease in lag time and an increase in fibril formation rate were found. In the 4.0-8.0 pH range there was no spectroscopic evidence for significant conformational changes in the native state. Likewise, heat capacity measurements showed no evidence for residual structure in the unfolded state. However, at pH <3.0 stability is severely decreased and the protein suffers conformational changes as detected by circular dichroism, tryptophan and ANS fluorescence, as well as by NMR spectroscopy. Molecular dynamics simulations indicate that acid-induced conformational changes involve the exposure of the loop connecting strands E and F. These results are compatible with pH-induced changes in the NMR spectra. Overall, the results indicate that the mechanism involved in the acid-induced increase in the fibrillogenic potential of 6aJL2 is profoundly different to that observed in κ light chains, and is promoted by localized conformational changes in a region of the protein that was previously not known to be involved in acid-induced light chain fibril formation. The identification of this region opens the potential for the design of specific inhibitors.

QSAR and Molecular Docking Studies of the Inhibitory Activity of Novel Heterocyclic GABA Analogues over GABA-AT.

We have previously reported the synthesis, in vitro and in silico activities of new GABA analogues as inhibitors of the GABA-AT enzyme from Pseudomonas fluorescens, where the nitrogen atom at the γ-position is embedded in heterocyclic scaffolds. With the goal of finding more potent inhibitors, we now report the synthesis of a new set of GABA analogues with a broader variation of heterocyclic scaffolds at the γ-position such as thiazolidines, methyl-substituted piperidines, morpholine and thiomorpholine and determined their inhibitory potential over the GABA-AT enzyme from Pseudomonas fluorescens. These structural modifications led to compound 9b which showed a 73% inhibition against this enzyme. In vivo studies with PTZ-induced seizures on male CD1 mice show that compound 9b has a neuroprotective effect at a 0.50 mmole/kg dose. A QSAR study was carried out to find the molecular descriptors associated with the structural changes in the GABA scaffold to explain their inhibitory activity against GABA-AT. Employing 3D molecular descriptors allowed us to propose the GABA analogues enantiomeric active form. To evaluate the interaction with Pseudomonas fluorescens and human GABA-AT by molecular docking, the constructions of homology models was carried out. From these calculations, 9b showed a strong interaction with both GABA-AT enzymes in agreement with experimental results and the QSAR model, which indicates that bulky ligands tend to be the better inhibitors especially those with a sulfur atom on their structure.

Novel-Substituted Heterocyclic GABA Analogues. Enzymatic Activity against the GABA-AT Enzyme from Pseudomonas fluorescens and In Silico Molecular Modeling.

γ-Aminobutyric acid (GABA) is the most important inhibitory neurotransmitter in the central nervous system, and a deficiency of GABA is associated with serious neurological disorders. Due to its low lipophilicity, there has been an intensive search for new molecules with increased lipophilicity to cross the blood-brain barrier to raise GABA concentrations. We have designed and evaluated in vitro and in silico some new analogues of GABA, where the nitrogen atom at the γ-position is embedded in heterocyclic scaffolds and determined their inhibitory potential over the GABA-AT enzyme from Pseudomonas fluorescens. These modifications lead to compounds with inhibitory activity as it occurs with compounds 18a and 19a. The construction of Pseudomonas fluorescens and human GABA-AT models were carried out by homology modeling. Docking assays were done for these compounds over the GABA-AT enzyme models where 19a showed a strong interaction with both GABA-AT enzymes.

A family 13 thioesterase isolated from an activated sludge metagenome: Insights into aromatic compounds metabolism.

Activated sludge is produced during the treatment of sewage and industrial wastewaters. Its diverse chemical composition allows growth of a large collection of microbial phylotypes with very different physiologic and metabolic profiles. Thus, activated sludge is considered as an excellent environment to discover novel enzymes through functional metagenomics, especially activities related with degradation of environmental pollutants. Metagenomic DNA was isolated and purified from an activated sludge sample. Metagenomic libraries were subsequently constructed in Escherichia coli. Using tributyrin hydrolysis, a screening by functional analysis was conducted and a clone that showed esterase activity was isolated. Blastx analysis of the sequence of the cloned DNA revealed, among others, an ORF that encodes a putative thioesterase with 47-64% identity to GenBank CDS reported genes, similar to those in the hotdog fold thioesterase superfamily. On the basis of its amino acid similarity and its homology-modelled structure we deduced that this gene encodes an enzyme (ThYest_ar) that belongs to family TE13, with a preference for aryl-CoA substrates and a novel catalytic residue constellation. Plasmid retransformation in E. coli confirmed the clone's phenotype, and functional complementation of a paaI E. coli mutant showed preference for phenylacetate over chlorobenzene as a carbon source. This work suggests a role for TE13 family thioesterases in swimming and degradation approaches for phenyl acetic acid. Proteins 2017; 85:1222-1237. © 2017 Wiley Periodicals, Inc.

Dissecting the protein architecture of DNA-binding transcription factors in bacteria and archaea.

Gene regulation at the transcriptional level is a central process in all organisms where DNA-binding transcription factors play a fundamental role. This class of proteins binds specifically at DNA sequences, activating or repressing gene expression as a function of the cell's metabolic status, operator context and ligand-binding status, among other factors, through the DNA-binding domain (DBD). In addition, TFs may contain partner domains (PaDos), which are involved in ligand binding and protein-protein interactions. In this work, we systematically evaluated the distribution, abundance and domain organization of DNA-binding TFs in 799 non-redundant bacterial and archaeal genomes. We found that the distributions of the DBDs and their corresponding PaDos correlated with the size of the genome. We also identified specific combinations between the DBDs and their corresponding PaDos. Within each class of DBDs there are differences in the actual angle formed at the dimerization interface, responding to the presence/absence of ligands and/or crystallization conditions, setting the orientation of the resulting helices and wings facing the DNA. Our results highlight the importance of PaDos as central elements that enhance the diversity of regulatory functions in all bacterial and archaeal organisms, and our results also demonstrate the role of PaDos in sensing diverse signal compounds. The highly specific interactions between DBDs and PaDos observed in this work, together with our structural analysis highlighting the difficulty in predicting both inter-domain geometry and quaternary structure, suggest that these systems appeared once and evolved with diverse duplication events in all the analysed organisms.

Convergent mechanisms favor fast amyloid formation in two lambda 6a Ig light chain mutants.

Extracellular deposition as amyloids of immunoglobulin light chains causes light chain amyloidosis. Among the light chain families, lambda 6a is one of the most frequent in light chain amyloidosis patients. Its germline protein, 6aJL2, and point mutants, R24G and P7S, are good models to study fibrillogenesis, because their stability and fibril formation characteristics have been described. Both mutations make the germline protein unstable and speed up its ability to aggregate. To date, there is no molecular mechanism that explains how these differences in amyloidogenesis can arise from a single mutation. To look into the structural and dynamical differences in the native state of these proteins, we carried out molecular dynamics simulations at room temperature. Despite the structural similarity of the germline protein and the mutants, we found differences in their dynamical signatures that explain the mutants' increased tendency to form amyloids. The contact network alterations caused by the mutations, though different, converge in affecting two anti-aggregation motifs present in light chain variable domains, suggesting a different starting point for aggregation in lambda chains compared to kappa chains.

Inhibition of Light Chain 6aJL2-R24G Amyloid Fiber Formation Associated with Light Chain Amyloidosis.

Light chain amyloidosis (AL) is a deadly disease characterized by the deposition of monoclonal immunoglobulin light chains as insoluble amyloid fibrils in different organs and tissues. Germ line λ VI has been closely related to this condition; moreover, the R24G mutation is present in 25% of the proteins of this germ line in AL patients. In this work, five small molecules were tested as inhibitors of the formation of amyloid fibrils from the 6aJL2-R24G protein. We have found by thioflavin T fluorescence and transmission electron microscopy that EGCG inhibits 6aJL2-R24G fibrillogenesis. Furthermore, using nuclear magnetic resonance spectroscopy, dynamic light scattering, and isothermal titration calorimetry, we have determined that the inhibition is due to binding to the protein in its native state, interacting mainly with aromatic residues.

Electrostatic analysis of bacterial expansins.

Expansins are a family of proteins with plant cell wall remodeling-activity, which bind cell wall components through hydrophobic and electrostatic interactions. A shallow area on the surface of the protein serves as the polysaccharide binding site (PBS) and it is composed of conserved residues. However, electric charge differences on the opposite face of the PBS produce basic, neutral, or acidic proteins. An analysis of forty-four bacterial expansins, homologues of BsEXLX1, revealed two main groups defined by: (a) the presence or absence of disulfide bonds; and (b) by the proteins isoelectric point (pI). We determined the location of the residues responsible for the pI on the structure of representative expansins. Our results suggest that the electric charge at the opposite site of the PBS may help in substrate differentiation among expansins from different species; in addition, electrostatic polarization between the front and the back of the molecule could affect expansin activity on cellulose.

Identification of a novel carbohydrate esterase from Bjerkandera adusta: structural and function predictions through bioinformatics analysis and molecular modeling.

A new gene from Bjerkandera adusta strain UAMH 8258 encoding a carbohydrate esterase (designated as BacesI) was isolated and expressed in Pichia pastoris. The gene had an open reading frame of 1410 bp encoding a polypeptide of 470 amino acid residues, the first 18 serving as a secretion signal peptide. Homology and phylogenetic analyses showed that BaCesI belongs to carbohydrate esterases family 4. Three-dimensional modeling of the protein and normal mode analysis revealed a breathing mode of the active site that could be relevant for esterase activity. Furthermore, the overall negative electrostatic potential of this enzyme suggests that it degrades neutral substrates and will not act on negative substrates such as peptidoglycan or p-nitrophenol derivatives. The enzyme shows a specific activity of 1.118 U mg(-1) protein on 2-naphthyl acetate. No activity was detected on p-nitrophenol derivatives as proposed from the electrostatic potential data. The deacetylation activity of the recombinant BaCesI was confirmed by measuring the release of acetic acid from several substrates, including oat xylan, shrimp shell chitin, N-acetylglucosamine, and natural substrates such as sugar cane bagasse and grass. This makes the protein very interesting for the biofuels production industry from lignocellulosic materials and for the production of chitosan from chitin.

A novel TctA citrate transporter from an activated sludge metagenome: structural and mechanistic predictions for the TTT family.

We isolated a putative citrate transporter of the tripartite tricarboxylate transporter (TTT) class from a metagenomic library of activated sludge from a sewage treatment plant. The transporter, dubbed TctA_ar, shares ∼50% sequence identity with TctA of Comamonas testosteroni (TctA_ct) and other ß-Proteobacteria, and contains two 20-amino acid repeat signature sequences, considered a hallmark of this particular transporter class. The structures for both TctA_ar and TctA_ct were modeled with I-TASSER and two possible structures for this transporter family were proposed. Docking assays with citrate resulted in the corresponding sets of proposed critical residues for function. These models suggest functions for the 20-amino acid repeats in the context of the two different architectures. This constitutes the first attempt at structure modeling of the TTT family, to the best of our knowledge, and could aid functional understanding of this little-studied family.

Cu(II) binding to the λ6aJL2-R24G antibody light chain protein associated with light chain amyloidosis disease: The role of histidines.

Light chain amyloidosis is a conformational disease caused by the abnormal proliferation and deposition of antibody light chains as amyloid fibers in organs and tissues. The effect of Cu(II) binding to the model recombinant protein 6aJL2-R24G was previously characterized in our group, and we found an acceleration of the aggregation kinetics of the protein. In this study, in order to confirm the Cu(II) binding sites, histidine variants of 6aJL2-R24G were prepared and the effects of their interaction with Cu(II) were analyzed by circular dichroism, fluorescence spectroscopy, isothermal calorimetry titrations, and molecular dynamics simulations. Confirming our earlier work, we found that His8 and His99 are the highest affinity Cu(II) binding sites, and that Cu(II) binding to both sites is a cooperative event.

Characterization of a catalase-peroxidase variant (L333V-KatG) identified in an INH-resistant Mycobacterium tuberculosis clinical isolate.

Mycobacterium tuberculosis catalase-peroxidase (Mt-KatG) is a bifunctional heme-dependent enzyme that has been shown to activate isoniazid (INH), the widely used antibiotic against tuberculosis (TB). The L333V-KatG variant has been associated with INH resistance in clinical M. tuberculosis isolates from Mexico. To understand better the mechanisms of INH activation, its catalytic properties (catalase, peroxidase, and IN-NAD formation) and crystal structure were compared with those of the wild-type enzyme (WT-KatG). The rate of IN-NAD formation mediated by WT-KatG was 23% greater than L333V-KatG when INH concentration is varied. In contrast to WT-KatG, the crystal structure of the L333V-KatG variant has a perhydroxy modification of the indole nitrogen of W107 from MYW adduct. L333V-KatG shows most of the active site residues in a similar position to WT-KatG; only R418 is in the R-conformation instead of the double R and Y conformation present in WT-KatG. L333V-KatG shows a small displacement respect to WT-KatG in the helix from R385 to L404 towards the mutation site, an increase in length of the coordination bond between H270 and heme Fe, and a longer H-bond between proximal D381 and W321, compared to WT-KatG; these small displacements could explain the altered redox potential of the heme, and result in a less active and stable enzyme.

The Role of a Loop in the Non-catalytic Domain B on the Hydrolysis/Transglycosylation Specificity of the 4-α-Glucanotransferase from Thermotoga maritima.

The mechanism by which glycoside hydrolases control the reaction specificity through hydrolysis or transglycosylation is a key element embedded in their chemical structures. The determinants of reaction specificity seem to be complex. We looked for structural differences in domain B between the 4-α-glucanotransferase from Thermotoga maritima (TmGTase) and the α-amylase from Thermotoga petrophila (TpAmylase) and found a longer loop in the former that extends towards the active site carrying a W residue at its tip. Based on these differences we constructed the variants W131G and the partial deletion of the loop at residues 120-124/128-131, which showed a 11.6 and 11.4-fold increased hydrolysis/transglycosylation (H/T) ratio relative to WT protein, respectively. These variants had a reduction in the maximum velocity of the transglycosylation reaction, while their affinity for maltose as the acceptor was not substantially affected. Molecular dynamics simulations allow us to rationalize the increase in H/T ratio in terms of the flexibility near the active site and the conformations of the catalytic acid residues and their associated pKas.