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.

Identification of a Family of Glycoside Derivatives Biologically Active against Acinetobacter baumannii and Other MDR Bacteria Using a QSPR Model.

As the rate of discovery of new antibacterial compounds for multidrug-resistant bacteria is declining, there is an urge for the search for molecules that could revert this tendency. Acinetobacter baumannii has emerged as a highly virulent Gram-negative bacterium that has acquired multiple resistance mechanisms against antibiotics and is considered of critical priority. In this work, we developed a quantitative structure-property relationship (QSPR) model with 592 compounds for the identification of structural parameters related to their property as antibacterial agents against A. baumannii. QSPR mathematical validation (R2 = 70.27, RN = -0.008, a(R2) = 0.014, and δK = 0.021) and its prediction ability (Q2LMO= 67.89, Q2EXT = 67.75, a(Q2) = -0.068, δQ = 0.0, rm2¯ = 0.229, and Δrm2 = 0.522) were obtained with different statistical parameters; additional validation was done using three sets of external molecules (R2 = 72.89, 71.64 and 71.56). We used the QSPR model to perform a virtual screening on the BIOFACQUIM natural product database. From this screening, our model showed that molecules 32 to 35 and 54 to 68, isolated from different extracts of plants of the Ipomoea sp., are potential antibacterials against A. baumannii. Furthermore, biological assays showed that molecules 56 and 60 to 64 have a wide antibacterial activity against clinically isolated strains of A. baumannii, as well as other multidrug-resistant bacteria, including Staphylococcus aureus, Escherichia coli, Klebsiella pneumonia, and Pseudomonas aeruginosa. Finally, we propose 60 as a potential lead compound due to its broad-spectrum activity and its structural simplicity. Therefore, our QSPR model can be used as a tool for the investigation and search for new antibacterial compounds against A. baumannii.

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.

Colombian Scorpion Centruroides margaritatus: Purification and Characterization of a Gamma Potassium Toxin with Full-Block Activity on the hERG1 Channel.

The Colombian scorpion Centruroides margaritatus produces a venom considered of low toxicity. Nevertheless, there are known cases of envenomation resulting in cardiovascular disorders, probably due to venom components that target ion channels. Among them, the humanether-à-go-go-Related gene (hERG1) potassium channels are critical for cardiac action potential repolarization and alteration in its functionality are associated with cardiac disorders. This work describes the purification and electrophysiological characterization of a Centruroides margaritatus venom component acting on hERG1 channels, the CmERG1 toxin. This novel peptide is composed of 42 amino acids with a MW of 4792.88 Da, folded by four disulfide bonds and it is classified as member number 10 of the γ-KTx1 toxin family. CmERG1 inhibits hERG1 currents with an IC50 of 3.4 ± 0.2 nM. Despite its 90.5% identity with toxin É£-KTx1.1, isolated from Centruroides noxius, CmERG1 completely blocks hERG1 current, suggesting a more stable plug of the hERG channel, compared to that formed by other É£-KTx.

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.

Identification of an α-MoRF in the Intrinsically Disordered Region of the Escargot Transcription Factor.

Molecular recognition features (MoRFs) are common in intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). MoRFs are in constant order-disorder structural transitions and adopt well-defined structures once they are bound to their targets. Here, we study Escargot (Esg), a transcription factor in Drosophila melanogaster that regulates multiple cellular functions, and consists of a disordered N-terminal domain and a group of zinc fingers at its C-terminal domain. We analyzed the N-terminal domain of Esg with disorder predictors and identified a region of 45 amino acids with high probability to form ordered structures, which we named S2. Through 54 µs of molecular dynamics (MD) simulations using CHARMM36 and implicit solvent (generalized Born/surface area (GBSA)), we characterized the conformational landscape of S2 and found an α-MoRF of ∼16 amino acids stabilized by key contacts within the helix. To test the importance of these contacts in the stability of the α-MoRF, we evaluated the effect of point mutations that would impair these interactions, running 24 µs of MD for each mutation. The mutations had mild effects on the MoRF, and in some cases, led to gain of residual structure through long-range contacts of the α-MoRF and the rest of the S2 region. As this could be an effect of the force field and solvent model we used, we benchmarked our simulation protocol by carrying out 32 µs of MD for the (AAQAA)3 peptide. The results of the benchmark indicate that the global amount of helix in shorter peptides like (AAQAA)3 is reasonably predicted. Careful analysis of the runs of S2 and its mutants suggests that the mutation to hydrophobic residues may have nucleated long-range hydrophobic and aromatic interactions that stabilize the MoRF. Finally, we have identified a set of residues that stabilize an α-MoRF in a region still without functional annotations in Esg.

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.

Site-Specific Interactions with Copper Promote Amyloid Fibril Formation for λ6aJL2-R24G.

Light-chain amyloidosis (AL) is one of the most common systemic amyloidoses, and it is characterized by the deposition of immunoglobulin light chain (LC) variable domains as insoluble amyloid fibers in vital organs and tissues. The recombinant protein 6aJL2-R24G contains λ6a and JL2 germline genes and also contains the Arg24 by Gly substitution. This mutation is present in 25% of all amyloid-associated λ6 LC cases, reduces protein stability, and increases the propensity to form amyloid fibers. In this study, it was found that the interaction of 6aJL2-R24G with Cu(II) decreases the thermal stability of the protein and accelerates the amyloid fibril formation, as observed by fluorescence spectroscopy. Isothermal calorimetry titration showed that Cu(II) binds to the protein with micromolar affinity. His99 may be one of the main Cu(II) interaction sites, as observed by nuclear magnetic resonance spectroscopy. The binding of Cu(II) to His99 induces larger fluctuations of the CDR1 and loop C″, as shown by molecular dynamics simulations. Thus, Cu(II) binding may be inducing the loss of interactions between CDR3 and CDR1, making the protein less stable and more prone to form amyloid fibers. This study provides insights into the mechanism of metal-induced aggregation of the 6aJL2-R24G protein and sheds light on the bio-inorganic understanding of AL disease.

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.

The TATA-binding Protein DNA-binding domain of eukaryotic parasites is a potentially druggable target.

The TATA-binding protein (TBP) is a central transcription factor in eukaryotes that interacts with a large number of different transcription factors; thus, affecting these interactions will be lethal for any living being. In this work, we present the first structural and dynamic computational study of the surface properties of the TBP DNA-binding domain for a set of parasites involved in diseases of worldwide interest. The sequence and structural differences of these TBPs, as compared with human TBP, were proposed to select representative ensembles generated from molecular dynamics simulations and to evaluate their druggability by molecular ensemble-based docking of drug-like molecules. We found that potential druggable sites correspond to the NC2-binding site, N-terminal tail, H2 helix, and the interdomain region, with good selectivity for Plasmodium falciparum, Necator americanus, Entamoeba histolytica, Candida albicans, and Taenia solium TBPs. The best hit compounds share structural similarity among themselves and have predicted dissociation constants ranging from nM to µM. These can be proposed as initial scaffolds for experimental testing and further optimization. In light of the obtained results, we propose TBP as an attractive therapeutic target for treatment of parasitic diseases.

A Review on Viral Metagenomics in Extreme Environments.

Viruses are the most abundant biological entities in the biosphere, and have the ability to infect Bacteria, Archaea, and Eukaryotes. The virome is estimated to be at least ten times more abundant than the microbiome with 107 viruses per milliliter and 109 viral particles per gram in marine waters and sediments or soils, respectively. Viruses represent a largely unexplored genetic diversity, having an important role in the genomic plasticity of their hosts. Moreover, they also play a significant role in the dynamics of microbial populations. In recent years, metagenomic approaches have gained increasing popularity in the study of environmental viromes, offering the possibility of extending our knowledge related to both virus diversity and their functional characterization. Extreme environments represent an interesting source of both microbiota and their virome due to their particular physicochemical conditions, such as very high or very low temperatures and >1 atm hydrostatic pressures, among others. Despite the fact that some progress has been made in our understanding of the ecology of the microbiota in these habitats, few metagenomic studies have described the viromes present in extreme ecosystems. Thus, limited advances have been made in our understanding of the virus community structure in extremophilic ecosystems, as well as in their biotechnological potential. In this review, we critically analyze recent progress in metagenomic based approaches to explore the viromes in extreme environments and we discuss the potential for new discoveries, as well as methodological challenges and perspectives.

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.

FRET-based analysis and molecular modeling of the human GPN-loop GTPases 1 and 3 heterodimer unveils a dominant-negative protein complex.

GPN-loop GTPases 1 and 3 are required for RNA polymerase II (RNAPII) nuclear import. Gpn1 and Gpn3 display some sequence similarity, physically associate, and their protein expression levels are mutually dependent in human cells. We performed here Fluorescence Resonance Energy Transfer (FRET), molecular modeling, and cell biology experiments to understand, and eventually disrupt, human Gpn1-Gpn3 interaction in live HEK293-AD cells. Transiently expressed EYFP-Gpn1 and Gpn3-CFP generated a strong FRET signal, indicative of a very close proximity, in the cytoplasm of HEK293-AD cells. Molecular modeling of the human Gpn1-Gpn3 heterodimer based on the crystallographic structure of Npa3, the Saccharomyces cerevisiae Gpn1 ortholog, revealed that human Gpn1 and Gpn3 associate through a large interaction surface formed by internal α-helix 7, insertion 2, and the GPN-loop from each protein. In site-directed mutagenesis experiments of interface residues, we identified the W132D and M227D EYFP-Gpn1 mutants as defective to produce a FRET signal when coexpressed with Gpn3-CFP. Simultaneous but not individual expression of Gpn1 and Gpn3, with either or both proteins fused to EYFP, retained RNAPII in the cytoplasm and markedly inhibited global transcription in HEK293-AD cells. Interestingly, the W132D and M227D Gpn1 mutants that showed an impaired ability to interact with Gpn3 by FRET were also unable to delocalize RNAPII in this assay, indicating that an intact Gpn1-Gpn3 interaction is required to display the dominant-negative effect on endogenous Gpn1/Gpn3 function we described here. Altogether, our results suggest that a Gpn1-Gpn3 strong interaction is critical for their cellular function.

In silico structure-based design of GABAB receptor agonists using a combination of docking and QSAR.

The study of γ-aminobutyric acid B receptor (GABAB ) activation is of great interest for several brain disorders. The search of new GABAB receptor agonists has been carried out by many research groups. As a result, Baclofen has become the prototypical GABAB receptor agonist. However, several attempts have been made to modify its structure to generate derivatives with improved activity. In this work, we carried out a theoretical and computational study for a wide range of GABAB receptor agonists reported in the literature. Molecular docking and QSAR techniques were combined by using the interaction energies of the agonists with the key residues of GABAB receptor, as molecular descriptors for the QSAR construction. The resulting mathematical model suggests that the activity of GABAB receptor agonists is influenced by three factors: their shape and molecular size (PW5 and PJI2), their constitutional features (ELUMO and T(N…O)) and the energy interaction with GABAB receptor (ETRP278 ). This model was validated by the QUIK, REDUNDANCY and OVERFITTING rules, and its predicted ability was tasted by the QLOO , QASYM , R02 and rm2 rules. Finally, six new compounds are proposed (35-40) with high potential to be used as GABAB receptor agonists.

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.

The human adenovirus type 5 E1B 55kDa protein interacts with RNA promoting timely DNA replication and viral late mRNA metabolism.

The E1B 55kDa produced by human adenovirus type 5 is a multifunctional protein that participates in the regulation of several steps during the viral replication cycle. Previous studies suggest this protein plays an important role in postranscriptional regulation of viral and cellular gene expression, as it is required for the selective accumulation of maximal levels of viral late mRNA in the cytoplasm of the infected cell; however the molecular mechanisms that are altered or regulated by this protein have not been elucidated. A ribonucleoprotein motif that could implicate the direct interaction of the protein with RNA was initially predicted and tested in vitro, but the interaction with RNA could not be detected in infected cells, suggesting the interaction may be weak or transient. Here it was determined that the E1B 55kDa interacts with RNA in the context of the viral infection in non-transformed human cells, and its contribution to the adenovirus replication cycle was evaluated. Using recombinant adenoviruses with amino acid substitutions or a deletion in the ribonucleoprotein motif the interaction of E1B 55kDa with RNA was found to correlate with timely and efficient viral DNA replication and viral late mRNA accumulation and splicing.

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.