CNA web server: rigidity theory-based thermal unfolding simulations of proteins for linking structure, (thermo-)stability, and function

Exploring Bacillus species xylanases for industrial applications: screening via thermostability and reaction modelling

CONTEXT

Xylanases derived from Bacillus species hold significant importance in various large-scale production sectors, with increasing demand driven by biofuel production. However, despite their potential, the extreme environmental conditions often encountered in production settings have led to their underutilisation. To address this issue and enhance their efficacy under adverse conditions, we conducted a theoretical investigation on a group of five Bacillus species xylanases belonging to the glycoside hydrolase GH11 family. Bacillus sp. NCL 87-6-10 (sp_NCL 87-6-10) emerged as a potent candidate among the selected biocatalysts; this Bacillus strain exhibited high thermal stability and achieved a transition state with minimal energy requirements, thereby accelerating the biocatalytic reaction process. Our approach aims to provide support for experimentalists in the industrial sector, encouraging them to employ structural-based reaction modelling scrutinisation to predict the ability of targeted xylanases.

METHODS

Utilising crystal structure data available in the Carbohydrate-Active enzymes database, we aimed to analyse their structural capabilities in terms of thermal-stability and activity. Our investigation into identifying the most prominent Bacillus species xylanases unfolds with the help of the semi-empirical quantum mechanics MOPAC method integrated with the DRIVER program is used in calculations of reaction pathways to understand the activation energy. Additionally, we scrutinised the selected xylanases using various analyses, including constrained network analyses, intermolecular interactions of the enzyme-substrate complex and molecular orbital assessments calculated using the AM1 method with the MO-G model (MO-G AM1) to validate their reactivity.

Conformational Stability of the N-Terminal Region of MDM2

MDM2 is an E3 ubiquitin ligase which is crucial for the degradation and inhibition of the key tumor-suppressor protein p53. In this work, we explored the stability and the conformational features of the N-terminal region of MDM2 (N-MDM2), through which it binds to the p53 protein as well as other protein partners. The isolated domain possessed a native-like conformational stability in a narrow pH range (7.0 to 10.0), as shown by intrinsic and 8-anilinonapthalene-1-sulfonic acid (ANS) fluorescence, far-UV circular dichroism (CD), and size exclusion chromatography (SEC). Guanidinium chloride (GdmCl) denaturation followed by intrinsic and ANS fluorescence, far-UV CD and SEC at physiological pH, and differential scanning calorimetry (DSC) and thermo-fluorescence experiments showed that (i) the conformational stability of isolated N-MDM2 was very low; and (ii) unfolding occurred through the presence of several intermediates. The presence of a hierarchy in the unfolding intermediates was also evidenced through DSC and by simulating the unfolding process with the help of computational techniques based on constraint network analysis (CNA). We propose that the low stability of this protein is related to its inherent flexibility and its ability to interact with several molecular partners through different routes.

Analyzing Current Trends and Possible Strategies to Improve Sucrose Isomerases' Thermostability

Due to their ability to produce isomaltulose, sucrose isomerases are enzymes that have caught the attention of researchers and entrepreneurs since the 1950s. However, their low activity and stability at temperatures above 40 °C have been a bottleneck for their industrial application. Specifically, the instability of these enzymes has been a challenge when it comes to their use for the synthesis and manufacturing of chemicals on a practical scale. This is because industrial processes often require biocatalysts that can withstand harsh reaction conditions, like high temperatures. Since the 1980s, there have been significant advancements in the thermal stabilization engineering of enzymes. Based on the literature from the past few decades and the latest achievements in protein engineering, this article systematically describes the strategies used to enhance the thermal stability of sucrose isomerases. Additionally, from a theoretical perspective, we discuss other potential mechanisms that could be used for this purpose.

Protein engineering and iterative multimodule optimization for vitamin B6 production in Escherichia coli

Vitamin B6 is an essential nutrient with extensive applications in the medicine, food, animal feed, and cosmetics industries. Pyridoxine (PN), the most common commercial form of vitamin B6, is currently chemically synthesized using expensive and toxic chemicals. However, the low catalytic efficiencies of natural enzymes and the tight regulation of the metabolic pathway have hindered PN production by the microbial fermentation process. Here, we report an engineered Escherichia coli strain for PN production. Parallel pathway engineering is performed to decouple PN production and cell growth. Further, protein engineering is rationally designed including the inefficient enzymes PdxA, PdxJ, and the initial enzymes Epd and Dxs. By the iterative multimodule optimization strategy, the final strain produces 1.4 g/L of PN with productivity of 29.16 mg/L/h by fed-batch fermentation. The strategies reported here will be useful for developing microbial strains for the production of vitamins and other bioproducts having inherently low metabolic fluxes.

Single substitution in α-helix of active center enhanced thermostability of Aspergillus awamori exo-inulinase

Exo-inulinases are applied in inulin hydrolysis and production of feed additives and need to be stable at temperatures of 60-95 °C. Aspergillus awamori exo-inulinase Inu1 is considerably thermostable, with a T_m of 73.2 °C. However, the thermostability of the enzyme should be improved. A single substitution G338A in α-helix in the active center of the enzyme provided a 3.5 °C improvement in T_m. The time of half-life at 70 °C and 80 °C was increased in 5.7- and 2.7-times, respectively, compared to wild-type. Molecular dynamics simulations demonstrated that the substitution G338A caused a decrease in RMSF not only for the α-helix 337-YAANI-341, but also for the catalytically active residues D41 and E241 and the amino acid residues forming the cleft of the active center. Calculations with Constraint Network Analysis for the variant G338A showed the increase in the stability of intramolecular clusters.

Enzyme adaptation to habitat thermal legacy shapes the thermal plasticity of marine microbiomes

Microbial communities respond to temperature with physiological adaptation and compositional turnover. Whether thermal selection of enzymes explains marine microbiome plasticity in response to temperature remains unresolved. By quantifying the thermal behaviour of seven functionally-independent enzyme classes (esterase, extradiol dioxygenase, phosphatase, beta-galactosidase, nuclease, transaminase, and aldo-keto reductase) in native proteomes of marine sediment microbiomes from the Irish Sea to the southern Red Sea, we record a significant effect of the mean annual temperature (MAT) on enzyme response in all cases. Activity and stability profiles of 228 esterases and 5 extradiol dioxygenases from sediment and seawater across 70 locations worldwide validate this thermal pattern. Modelling the esterase phase transition temperature as a measure of structural flexibility confirms the observed relationship with MAT. Furthermore, when considering temperature variability in sites with non-significantly different MATs, the broadest range of enzyme thermal behaviour and the highest growth plasticity of the enriched heterotrophic bacteria occur in samples with the widest annual thermal variability. These results indicate that temperature-driven enzyme selection shapes microbiome thermal plasticity and that thermal variability finely tunes such processes and should be considered alongside MAT in forecasting microbial community thermal response.

Rational Design of a Thermostable 2'-Deoxyribosyltransferase for Nelarabine Production by Prediction of Disulfide Bond Engineering Sites

One of the major drawbacks of the industrial implementation of enzymatic processes is the low operational stability of the enzymes under tough industrial conditions. In this respect, the use of thermostable enzymes in the industry is gaining ground during the last decades. Herein, we report a structure-guided approach for the development of novel and thermostable 2′-deoxyribosyltransferases (NDTs) based on the computational design of disulfide bonds on hot spot positions. To this end, a small library of NDT variants from Lactobacillus delbrueckii (LdNDT) with introduced cysteine pairs was created. Among them, LdNDTS104C (100% retained activity) was chosen as the most thermostable variant, displaying a six- and two-fold enhanced long-term stability when stored at 55 °C (t1/255 °C ≈ 24 h) and 60 °C (t1/260 °C ≈ 4 h), respectively. Moreover, the biochemical characterization revealed that LdNDTS104C showed >60% relative activity across a broad range of temperature (30−90 °C) and pH (5−7). Finally, to study the potential application of LdNDTS104C as an industrial catalyst, the enzymatic synthesis of nelarabine was successfully carried out under different substrate conditions (1:1 and 3:1) at different reaction times. Under these experimental conditions, the production of nelarabine was increased up to 2.8-fold (72% conversion) compared with wild-type LdNDT.

Structure and Function of Redox-Sensitive Superfolder Green Fluorescent Protein Variant

Aims: Genetically encoded green fluorescent protein (GFP)-based redox biosensors are widely used to monitor specific and dynamic redox processes in living cells. Over the last few years, various biosensors for a variety of applications were engineered and enhanced to match the organism and cellular environments, which should be investigated. In this context, the unicellular intraerythrocytic parasite Plasmodium, the causative agent of malaria, represents a challenge, as the small size of the organism results in weak fluorescence signals that complicate precise measurements, especially for cell compartment-specific observations. To address this, we have functionally and structurally characterized an enhanced redox biosensor superfolder roGFP2 (sfroGFP2). Results: SfroGFP2 retains roGFP2-like behavior, yet with improved fluorescence intensity (FI) in cellulo. SfroGFP2-based redox biosensors are pH insensitive in a physiological pH range and show midpoint potentials comparable with roGFP2-based redox biosensors. Using crystallography and rigidity theory, we identified the superfolding mutations as being responsible for improved structural stability of the biosensor in a redox-sensitive environment, thus explaining the improved FI in cellulo. Innovation: This work provides insight into the structure and function of GFP-based redox biosensors. It describes an improved redox biosensor (sfroGFP2) suitable for measuring oxidizing effects within small cells where applicability of other redox sensor variants is limited. Conclusion: Improved structural stability of sfroGFP2 gives rise to increased FI in cellulo. Fusion to hGrx1 (human glutaredoxin-1) provides the hitherto most suitable biosensor for measuring oxidizing effects in Plasmodium. This sensor is of major interest for studying glutathione redox changes in small cells, as well as subcellular compartments in general. Antioxid. Redox Signal. 37, 1-18.

F/G Region Rigidity is Inversely Correlated to Substrate Promiscuity of Human CYP Isoforms Involved in Metabolism

Of 57 human cytochrome P450 (CYP) enzymes, 12 metabolize 90% of xenobiotics. To our knowledge, no study has addressed the relation between enzyme dynamics and substrate promiscuity for more than three CYPs. Here, we show by constraint dilution simulations with the Constraint Network Analysis for the 12 isoforms that structural rigidity of the F/G region is significantly inversely correlated to the enzymes' substrate promiscuity. This highlights the functional importance of structural dynamics of the substrate tunnel.

Aqueous ionic liquids redistribute local enzyme stability via long-range perturbation pathways

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for biocatalysis due to their unique properties. On the other hand, the incubation of enzymes in IL or aIL often reduces enzyme activity. Recent studies proposed various aIL-induced effects to explain the reduction, classified as direct effects, e.g., local dehydration or competitive inhibition, and indirect effects, e.g., structural perturbations or disturbed catalytic site integrity. However, the molecular origin of indirect effects has largely remained elusive. Here we show by multi-μs long molecular dynamics simulations, free energy computations, and rigidity analyses that aIL favorably interact with specific residues of Bacillus subtilis Lipase A (BsLipA) and modify the local structural stability of this model enzyme by inducing long-range perturbations of noncovalent interactions. The perturbations percolate over neighboring residues and eventually affect the catalytic site and the buried protein core. Validation against a complete experimental site saturation mutagenesis library of BsLipA (3620 variants) reveals that the residues of the perturbation pathways are distinguished sequence positions where substitutions highly likely yield significantly improved residual activity. Our results demonstrate that identifying these perturbation pathways and specific IL ion-residue interactions there effectively predicts focused variant libraries with improved aIL tolerance.

An engineered cellobiohydrolase I for sustainable degradation of lignocellulosic biomass

This study provides computational-assisted engineering of the cellobiohydrolase I (CBH-I) from Penicillium verruculosum with simultaneous enhanced thermostability and tolerance in ionic liquids, deep eutectic solvent, and concentrated seawater without affecting its wild-type activity. Engineered triple variant CBH-I R1 (A65R-G415R-S181F) showed 2.48-fold higher thermostability in terms of relative activity at 65°C after 1 h of incubation when compared with CBH-I wild type. CBH-I R1 exhibited 1.87-fold, 1.36-fold, and 1.57-fold higher specific activities compared with CBH-I wild type in [Bmim]Cl (50 g/L), [Ch]Cl (50 g/L), and two-fold concentrated seawater, respectively. In the multicellulases mixture, CBH-I R1 showed higher hydrolytic efficiency to hydrolyze aspen wood compared with CBH-I wild type in the buffer, [Bmim]Cl (50 g/L), and two-fold concentrated seawater, respectively. Structural analysis revealed a molecular basis for the higher stability of the CBH-I structure in which A65R and G415R substitutions form salt bridges (D64 … R65, E411 … R415) and S181F forms π-π interaction (Y155 … F181), leading to stabilize surface-exposed flexible α-helixes and loop in the multidomain β-jelly roll fold structure, respectively. In conclusion, the variant CBH-I R1 could enable efficient lignocellulosic biomass degradation as a cost-effective alternative for the sustainable production of biofuels and value-added chemicals.

Promiscuous Esterases Counterintuitively Are Less Flexible than Specific Ones

Understanding mechanisms of promiscuity is increasingly important from a fundamental and application point of view. As to enzyme structural dynamics, more promiscuous enzymes generally have been recognized to also be more flexible. However, examples for the opposite received much less attention. Here, we exploit comprehensive experimental information on the substrate promiscuity of 147 esterases tested against 96 esters together with computationally efficient rigidity analyses to understand the molecular origin of the observed promiscuity range. Unexpectedly, our data reveal that promiscuous esterases are significantly less flexible than specific ones, are significantly more thermostable, and have a significantly increased specific activity. These results may be reconciled with a model according to which structural flexibility in the case of specific esterases serves for conformational proofreading. Our results signify that an esterase sequence space can be screened by rigidity analyses for promiscuous esterases as starting points for further exploration in biotechnology and synthetic chemistry.

Thermostability Engineering of a Class II Pyruvate Aldolase from Escherichia coli by in Vivo Folding Interference

The use of enzymes in industrial processes is often limited by the unavailability of biocatalysts with prolonged stability. Thermostable enzymes allow increased process temperature and thus higher substrate and product solubility, reuse of expensive biocatalysts, resistance against organic solvents, and better "evolvability" of enzymes. In this work, we have used an activity-independent method for the selection of thermostable variants of any protein in Thermus thermophilus through folding interference at high temperature of a thermostable antibiotic reporter protein at the C-terminus of a fusion protein. To generate a monomeric folding reporter, we have increased the thermostability of the moderately thermostable Hph5 variant of the hygromycin B phosphotransferase from Escherichia coli to meet the method requirements. The final Hph17 variant showed 1.5 °C higher melting temperature (T _m) and 3-fold longer half-life at 65 °C compared to parental Hph5, with no changes in the steady-state kinetic parameters. Additionally, we demonstrate the validity of the reporter by stabilizing the 2-keto-3-deoxy-l-rhamnonate aldolase from E. coli (YfaU). The most thermostable multiple-mutated variants thus obtained, YfaU99 and YfaU103, showed increases of 2 and 2.9 °C in T _m compared to the wild-type enzyme but severely lower retro-aldol activities (150- and 120-fold, respectively). After segregation of the mutations, the most thermostable single variant, Q107R, showed a T _m 8.9 °C higher, a 16-fold improvement in half-life at 60 °C and higher operational stability than the wild-type, without substantial modification of the kinetic parameters.

Can constraint network analysis guide the identification phase of KnowVolution? A case study on improved thermostability of an endo-β-glucanase

Cellulases are industrially important enzymes, e.g., in the production of bioethanol, in pulp and paper industry, feedstock, and textile. Thermostability is often a prerequisite for high process stability and improving thermostability without affecting specific activities at lower temperatures is challenging and often time-consuming. Protein engineering strategies that combine experimental and computational are emerging in order to reduce experimental screening efforts and speed up enzyme engineering campaigns. Constraint Network Analysis (CNA) is a promising computational method that identifies beneficial positions in enzymes to improve thermostability. In this study, we compare CNA and directed evolution in the identification of beneficial positions in order to evaluate the potential of CNA in protein engineering campaigns (e.g., in the identification phase of KnowVolution). We engineered the industrially relevant endoglucanase EGLII from Penicillium verruculosum towards increased thermostability. From the CNA approach, six variants were obtained with an up to 2-fold improvement in thermostability. The overall experimental burden was reduced to 40% utilizing the CNA method in comparison to directed evolution. On a variant level, the success rate was similar for both strategies, with 0.27% and 0.18% improved variants in the epPCR and CNA-guided library, respectively. In essence, CNA is an effective method for identification of positions that improve thermostability.

Computational Design of Nitrile Hydratase from Pseudonocardia thermophila JCM3095 for Improved Thermostability

High thermostability and catalytic activity are key properties for nitrile hydratase (NHase, EC 4.2.1.84) as a well-industrialized catalyst. In this study, rational design was applied to tailor the thermostability of NHase from Pseudonocardia thermophila JCM3095 (PtNHase) by combining FireProt server prediction and molecular dynamics (MD) simulation. Site-directed mutagenesis of non-catalytic residues provided by the rational design was subsequentially performed. The positive multiple-point mutant, namely, M10 (αI5P/αT18Y/αQ31L/αD92H/βA20P/βP38L/βF118W/βS130Y/βC189N/βC218V), was obtained and further analyzed. The Melting temperature (T_m) of the M10 mutant showed an increase by 3.2 °C and a substantial increase in residual activity of the enzyme at elevated temperatures was also observed. Moreover, the M10 mutant also showed a 2.1-fold increase in catalytic activity compared with the wild-type PtNHase. Molecular docking and MD simulations demonstrated better substrate affinity and improved thermostability for the mutant.

Systematically Scrutinizing the Impact of Substitution Sites on Thermostability and Detergent Tolerance for Bacillus subtilis Lipase A

Improving an enzyme's (thermo-)stability or tolerance against solvents and detergents is highly relevant in protein engineering and biotechnology. Recent developments have tended toward data-driven approaches, where available knowledge about the protein is used to identify substitution sites with high potential to yield protein variants with improved stability, and subsequently, substitutions are engineered by site-directed or site-saturation (SSM) mutagenesis. However, the development and validation of algorithms for data-driven approaches have been hampered by the lack of availability of large-scale data measured in a uniform way and being unbiased with respect to substitution types and locations. Here, we extend our knowledge on guidelines for protein engineering following a data-driven approach by scrutinizing the impact of substitution sites on thermostability or/and detergent tolerance for Bacillus subtilis lipase A (BsLipA) at very large scale. We systematically analyze a complete experimental SSM library of BsLipA containing all 3439 possible single variants, which was evaluated as to thermostability and tolerances against four detergents under respectively uniform conditions. Our results provide systematic and unbiased reference data at unprecedented scale for a biotechnologically important protein, identify consistently defined hot spot types for evaluating the performance of data-driven protein-engineering approaches, and show that the rigidity theory and ensemble-based approach Constraint Network Analysis yields hot spot predictions with an up to ninefold gain in precision over random classification.

Data-driven computational analysis of allosteric proteins by exploring protein dynamics, residue coevolution and residue interaction networks

BACKGROUND

Computational studies of allosteric interactions have witnessed a recent renaissance fueled by the growing interest in modeling of the complex molecular assemblies and biological networks. Allosteric interactions in protein structures allow for molecular communication in signal transduction networks.

METHODS

In this work, we performed a large scale comprehensive and multi-faceted analysis of >300 diverse allosteric proteins and complexes with allosteric modulators. By modeling and exploring coarse-grained dynamics, residue coevolution, and residue interaction networks for allosteric proteins, we have determined unifying molecular signatures shared by allosteric systems.

RESULTS

The results of this study have suggested that allosteric inhibitors and allosteric activators may differentially affect global dynamics and network organization of protein systems, leading to diverse allosteric mechanisms. By using structural and functional data on protein kinases, we present a detailed case study that that included atomic-level analysis of coevolutionary networks in kinases bound with allosteric inhibitors and activators.

CONCLUSIONS

We have found that coevolutionary networks can form direct communication pathways connecting functional regions and can recapitulate key regulatory sites and interactions responsible for allosteric signaling in the studied protein systems. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of known regulatory hotspots in protein kinases.

GENERAL SIGNIFICANCE

This study has shown that allosteric inhibitors and allosteric activators can have a different effect on residue interaction networks and can exploit distinct regulatory mechanisms, which could open up opportunities for probing allostery and new drug combinations with broad range of activities.

Integration of Random Forest Classifiers and Deep Convolutional Neural Networks for Classification and Biomolecular Modeling of Cancer Driver Mutations

Development of machine learning solutions for prediction of functional and clinical significance of cancer driver genes and mutations are paramount in modern biomedical research and have gained a significant momentum in a recent decade. In this work, we integrate different machine learning approaches, including tree based methods, random forest and gradient boosted tree (GBT) classifiers along with deep convolutional neural networks (CNN) for prediction of cancer driver mutations in the genomic datasets. The feasibility of CNN in using raw nucleotide sequences for classification of cancer driver mutations was initially explored by employing label encoding, one hot encoding, and embedding to preprocess the DNA information. These classifiers were benchmarked against their tree-based alternatives in order to evaluate the performance on a relative scale. We then integrated DNA-based scores generated by CNN with various categories of conservational, evolutionary and functional features into a generalized random forest classifier. The results of this study have demonstrated that CNN can learn high level features from genomic information that are complementary to the ensemble-based predictors often employed for classification of cancer mutations. By combining deep learning-generated score with only two main ensemble-based functional features, we can achieve a superior performance of various machine learning classifiers. Our findings have also suggested that synergy of nucleotide-based deep learning scores and integrated metrics derived from protein sequence conservation scores can allow for robust classification of cancer driver mutations with a limited number of highly informative features. Machine learning predictions are leveraged in molecular simulations, protein stability, and network-based analysis of cancer mutations in the protein kinase genes to obtain insights about molecular signatures of driver mutations and enhance the interpretability of cancer-specific classification models.

Three dimensional structure prediction of panomycocin, a novel Exo-β-1,3-glucanase isolated from Wickerhamomyces anomalus NCYC 434 and the computational site-directed mutagenesis studies to enhance its thermal stability for therapeutic applications

Panomycocin is a naturally produced potent antimycotic/antifungal protein secreted by the yeast Wickerhamomyces anomalus NCYC 434 with an exo-β-1,3-glucanase activity. In this study the three dimensional structure of panomycocin was predicted and the computational site-directed mutagenesis was performed to enhance its thermal stability in liquid formulations over the body temperature for topical therapeutic applications. Homology modeling was performed with MODELLER and I-TASSER. Among the generated models, the model with the lowest energy and DOPE score was selected for further loop modeling. The loop model was optimized and the reliability of the model was confirmed with ERRAT, Verify 3D and Ramachandran plot values. Enhancement of the thermal stability of the model was done using contemporary servers and programs such as SPDBViewer, CNA, I-Mutant2.0, Eris, AUTO-MUTE and MUpro. In the region outside the binding site of the model Leu52 Arg, Phe223Arg and Gly254Arg were found to be the best thermostabilizing mutations with 6.26 K, 6.26 K and 8.27 K increases, respectively. In the binding site Glu186Arg was found to be the best thermostabilizer mutation with a 9.58 K temperature increase. The results obtained in this study led us to design a mutant panomycocin that can be used as a novel antimycotic/antifungal drug in a liquid formulation for topical applications over the normal body temperature.

Atomistic Modeling of the ABL Kinase Regulation by Allosteric Modulators Using Structural Perturbation Analysis and Community-Based Network Reconstruction of Allosteric Communications

In this work, we have examined the molecular mechanisms of allosteric regulation of the ABL tyrosine kinase at the atomic level. Atomistic modeling of the ABL complexes with a panel of allosteric modulators has been performed using a combination of molecular dynamics simulations, structural residue perturbation scanning, and a novel community analysis of the residue interaction networks. Our results have indicated that allosteric inhibitors and activators may exert a differential control on allosteric signaling between the kinase binding sites and functional regions. While the inhibitor binding can strengthen the closed ABL state and induce allosteric communications directed from the allosteric pocket to the ATP binding site, the DPH activator may induce a more dynamic open form and activate allosteric couplings between the ATP and substrate binding sites. By leveraging a network-centric theoretical framework, we have introduced a novel community analysis method and global topological parameters that have unveiled the hierarchical modularity and the intercommunity bridging sites in the residue interaction network. We have found that allosteric functional hotspots responsible for the kinase regulation may serve the intermodular bridges in the global interaction network. The central conclusion from this analysis is that the regulatory switch centers play a fundamental role in the modular network organization of ABL as the unique intercommunity bridges that connect the SH2 and SH3 domains with the catalytic core into a functional kinase assembly. The hierarchy of network organization in the ABL regulatory complexes may allow for the synergistic action of dense intercommunity links required for the robust signal transfer in the catalytic core and sparse network bridges acting as the regulatory control points that orchestrate allosteric transitions between the inhibited and active kinase forms.

Disulfide Bond Engineering of an Endoglucanase from Penicillium verruculosum to Improve Its Thermostability

Endoglucanases (EGLs) are important components of multienzyme cocktails used in the production of a wide variety of fine and bulk chemicals from lignocellulosic feedstocks. However, a low thermostability and the loss of catalytic performance of EGLs at industrially required temperatures limit their commercial applications. A structure-based disulfide bond (DSB) engineering was carried out in order to improve the thermostability of EGLII from Penicillium verruculosum. Based on in silico prediction, two improved enzyme variants, S127C-A165C (DSB2) and Y171C-L201C (DSB3), were obtained. Both engineered enzymes displayed a 15⁻21% increase in specific activity against carboxymethylcellulose and β-glucan compared to the wild-type EGLII (EGLII-wt). After incubation at 70 °C for 2 h, they retained 52⁻58% of their activity, while EGLII-wt retained only 38% of its activity. At 80 °C, the enzyme-engineered forms retained 15⁻22% of their activity after 2 h, whereas EGLII-wt was completely inactivated after the same incubation time. Molecular dynamics simulations revealed that the introduced DSB rigidified a global structure of DSB2 and DSB3 variants, thus enhancing their thermostability. In conclusion, this work provides an insight into DSB protein engineering as a potential rational design strategy that might be applicable for improving the stability of other enzymes for industrial applications.

B genome specific polymorphism in the TdDRF1 gene is in relationship with grain yield

A and B genome copies of DRF1 gene in durum wheat were isolated and sequenced using gene variability. B genome specific polymorphism resulted, in a RIL population, in relationship with grain yield mainly in drought condition. Drought tolerance is one of the main components of yield potential and stability, and its improvement is a major challenge to breeders. Transcription factors are considered among the best candidate genes for developing functional markers, since they are components of the signal transduction pathways that coordinate the expression of several downstream genes. Polymorphisms of the Triticum durum dehydration responsive factor 1 (TdDRF1) gene that belongs to DREB2 transcription factor family were identified and specifically assigned to the A or B genome. A panel of primers was derived to selectively isolate the corresponding gene copies. These molecular information were also used to develop a new molecular marker: an allele-specific PCR assay discriminating two genotypes (Mohawk and Cocorit) was developed and used for screening a durum wheat recombinant inbred line population (RIL-pop) derived from the above genotypes. Phenotypic data from the RIL-pop grown during two seasons, under different environmental conditions, adopting an α-lattice design with two repetitions, were collected, analyzed and correlated with molecular data from the PCR assay. A significant association between a specific polymorphism in the B genome copy of the TdDRF1 gene and the grain yield in drought conditions were observed.

Molecular modeling and in silico characterization of GmABCC5: a phytate transporter and potential target for low-phytate crops

Designing low-phytate crops without affecting the developmental process in plants had led to the identification of ABCC5 gene in soybean. The GmABCC5 gene was identified and a partial gene sequence was cloned from popular Indian soybean genotype Pusa16. Conserved domains and motifs unique to ABC transporters were identified in the 30 homologous sequences retrieved by BLASTP analysis. The homologs were analyzed for their evolutionary relationship and physiochemical properties. Conserved domains, transmembrane architecture and secondary structure of GmABCC5 were predicted with the aid of computational tools. Analysis identified 53 alpha helices and 31 beta strands, predicting 60% residues in alpha conformation. A three-dimensional (3D) model for GmABCC5 was developed based on 5twv.1.B (Homo sapiens) template homology to gain better insight into its molecular mechanism of transport and sequestration. Spatio-temporal real-time PCR analysis identified mid-to-late seed developmental stages as the time window for the maximum GmABCC5 gene expression, a potential target stage for phytate reduction. Results of this study provide valuable insights into the structural and functional characteristics of GmABCC5, which may be further utilized for the development of nutritionally enriched low-phytate soybean with improved mineral bioavailability.

Atomistic simulations and network-based modeling of the Hsp90-Cdc37 chaperone binding with Cdk4 client protein: A mechanism of chaperoning kinase clients by exploiting weak spots of intrinsically dynamic kinase domains

A fundamental role of the Hsp90 and Cdc37 chaperones in mediating conformational development and activation of diverse protein kinase clients is essential in signal transduction. There has been increasing evidence that the Hsp90-Cdc37 system executes its chaperoning duties by recognizing conformational instability of kinase clients and modulating their folding landscapes. The recent cryo-electron microscopy structure of the Hsp90-Cdc37-Cdk4 kinase complex has provided a framework for dissecting regulatory principles underlying differentiation and recruitment of protein kinase clients to the chaperone machinery. In this work, we have combined atomistic simulations with protein stability and network-based rigidity decomposition analyses to characterize dynamic factors underlying allosteric mechanism of the chaperone-kinase cycle and identify regulatory hotspots that control client recognition. Through comprehensive characterization of conformational dynamics and systematic identification of stabilization centers in the unbound and client- bound Hsp90 forms, we have simulated key stages of the allosteric mechanism, in which Hsp90 binding can induce instability and partial unfolding of Cdk4 client. Conformational landscapes of the Hsp90 and Cdk4 structures suggested that client binding can trigger coordinated dynamic changes and induce global rigidification of the Hsp90 inter-domain regions that is coupled with a concomitant increase in conformational flexibility of the kinase client. This process is allosteric in nature and can involve reciprocal dynamic exchanges that exert global effect on stability of the Hsp90 dimer, while promoting client instability. The network-based rigidity analysis and emulation of thermal unfolding of the Cdk4-cyclin D complex and Hsp90-Cdc37-Cdk4 complex revealed weak spots of kinase instability that are present in the native Cdk4 structure and are targeted by the chaperone during client recruitment. Our findings suggested that this mechanism may be exploited by the Hsp90-Cdc37 chaperone to recruit and protect intrinsically dynamic kinase clients from degradation. The results of this investigation are discussed and interpreted in the context of diverse experimental data, offering new insights into mechanisms of chaperone regulation and binding.

Ensemble-based modeling and rigidity decomposition of allosteric interaction networks and communication pathways in cyclin-dependent kinases: Differentiating kinase clients of the Hsp90-Cdc37 chaperone

The overarching goal of delineating molecular principles underlying differentiation of protein kinase clients and chaperone-based modulation of kinase activity is fundamental to understanding activity of many oncogenic kinases that require chaperoning of Hsp70 and Hsp90 systems to attain a functionally competent active form. Despite structural similarities and common activation mechanisms shared by cyclin-dependent kinase (CDK) proteins, members of this family can exhibit vastly different chaperone preferences. The molecular determinants underlying chaperone dependencies of protein kinases are not fully understood as structurally similar kinases may often elicit distinct regulatory responses to the chaperone. The regulatory divergences observed for members of CDK family are of particular interest as functional diversification among these kinases may be related to variations in chaperone dependencies and can be exploited in drug discovery of personalized therapeutic agents. In this work, we report the results of a computational investigation of several members of CDK family (CDK5, CDK6, CDK9) that represented a broad repertoire of chaperone dependencies—from nonclient CDK5, to weak client CDK6, and strong client CDK9. By using molecular simulations of multiple crystal structures we characterized conformational ensembles and collective dynamics of CDK proteins. We found that the elevated dynamics of CDK9 can trigger imbalances in cooperative collective motions and reduce stability of the active fold, thus creating a cascade of favorable conditions for chaperone intervention. The ensemble-based modeling of residue interaction networks and community analysis determined how differences in modularity of allosteric networks and topography of communication pathways can be linked with the client status of CDK proteins. This analysis unveiled depleted modularity of the allosteric network in CDK9 that alters distribution of communication pathways and leads to impaired signaling in the client kinase. According to our results, these network features may uniquely define chaperone dependencies of CDK clients. The perturbation response scanning and rigidity decomposition approaches identified regulatory hotspots that mediate differences in stability and cooperativity of allosteric interaction networks in the CDK structures. By combining these synergistic approaches, our study revealed dynamic and network signatures that can differentiate kinase clients and rationalize subtle divergences in the activation mechanisms of CDK family members. The therapeutic implications of these results are illustrated by identifying structural hotspots of pathogenic mutations that preferentially target regions of the increased flexibility to enable modulation of activation changes. Our study offers a network-based perspective on dynamic kinase mechanisms and drug design by unravelling relationships between protein kinase dynamics, allosteric communications and chaperone dependencies.

1H, 13C and 15N NMR chemical shift assignments of A. thaliana RCD1 RST

The A. thaliana RCD1 (radical-induced cell death1) protein is a cellular signaling hub protein which interacts with numerous plant transcription factors from different families. It consists of three conserved domains and intervening unstructured regions, the C-terminal RST domain being responsible for the interactions with the transcription factors. It has been shown that many partner proteins interact with RCD1 RST via their intrinsically disordered regions, and that the domain is able to house partners with divergent folds. We aim to structurally characterize the RCD1 RST domain and its complexes [complex with DREB2A]. Here we report the ¹H, ¹⁵N and ¹³C chemical shift assignments of the backbone and sidechain atoms for RCD1 (468-589) containing the RST (510-567) domain.

Cryo-EM structure of the spinach chloroplast ribosome reveals the location of plastid-specific ribosomal proteins and extensions

Ribosomes are the protein synthesizing machines of the cell. Recent advances in cryo-EM have led to the determination of structures from a variety of species, including bacterial 70S and eukaryotic 80S ribosomes as well as mitoribosomes from eukaryotic mitochondria, however, to date high resolution structures of plastid 70S ribosomes have been lacking. Here we present a cryo-EM structure of the spinach chloroplast 70S ribosome, with an average resolution of 5.4 Å for the small 30S subunit and 3.6 Å for the large 50S ribosomal subunit. The structure reveals the location of the plastid-specific ribosomal proteins (RPs) PSRP1, PSRP4, PSRP5 and PSRP6 as well as the numerous plastid-specific extensions of the RPs. We discover many features by which the plastid-specific extensions stabilize the ribosome via establishing additional interactions with surrounding ribosomal RNA and RPs. Moreover, we identify a large conglomerate of plastid-specific protein mass adjacent to the tunnel exit site that could facilitate interaction of the chloroplast ribosome with the thylakoid membrane and the protein-targeting machinery. Comparing the Escherichia coli 70S ribosome with that of the spinach chloroplast ribosome provides detailed insight into the co-evolution of RP and rRNA.

Role of Cys73 in the thermostability of farnesyl diphosphate synthase from Geobacillus stearothermophilus

Farnesyl diphosphate synthase (FPPase) is an enzyme that catalyzes the condensation between one molecule of dimethylallyl diphosphate (DMAPP) and two molecules of isopentenyl diphosphate (IPP) to produce farnesyl diphosphate (FPP). FPP is an important precursor in the isoprenoid synthesis pathway. In this study, the crystal structure of FPPase from Geobacillus stearothermophilus (GsFPPase) was determined at 2.31 Å resolution. The structure of GsFPPase shows a three-layered all α-helical fold and conserved functional domains similar to other prenyltransferases. We have analyzed the structural features of GsFPPase related to thermostability and compared it with those of human and avian mesophilic FPPases. “Semi-conserved” regions which appear to be possible features contributing to the thermostability of FPPase were found.

Genome-wide identification, in silico characterization and expression analysis of ZIP-like genes from Trichomonas vaginalis in response to Zinc and Iron

Trace elements such as Zinc and Iron are essential components of metalloproteins and serve as cofactors or structural elements for enzymes involved in several important biological processes in almost all organisms. Because either excess or insufficient levels of Zn and Fe can be harmful for the cells, the homeostatic levels of these trace minerals must be tightly regulated. The Zinc regulated transporter, Iron regulated transporter-like Proteins (ZIP) comprise a diverse family, with several paralogues in diverse organisms and are considered essential for the Zn and Fe uptake and homeostasis. Zn and Fe has been shown to regulate expression of proteins involved in metabolism and pathogenicity mechanisms in the protozoan pathogen Trichomonas vaginalis, in contrast high concentrations of these elements were also found to be toxic for T. vaginalis trophozoites. Nevertheless, Zn and Fe uptake and homeostasis mechanisms is not yet clear in this parasite. We performed a genome-wide analysis and localized the 8 members of the ZIP gene family in T. vaginalis (TvZIP1-8). The bioinformatic programs predicted that the TvZIP proteins are highly conserved and show similar properties to the reported in other ZIP orthologues. The expression patterns of TvZIP1, 3, 5 and 7 were diminished in presence of Zinc, while the rest of the TvZIP genes showed an unchanged profile in this condition. In addition, TvZIP2 and TvZIP4 showed a differential expression pattern in trophozoites growth under different Iron conditions. These results suggest that TvZIP genes encode membrane transporters that may be responsible for the Zn and Fe acquisition in T. vaginalis.

H255Y and K508R missense mutations in tumour suppressor folliculin (FLCN) promote kidney cell proliferation

Germline H255Y and K508R missense mutations in the folliculin (FLCN) gene have been identified in patients with bilateral multifocal (BMF) kidney tumours and clinical manifestations of Birt-Hogg-Dubé (BHD) syndrome, or with BMF kidney tumours as the only manifestation; however, their impact on FLCN function remains to be determined. In order to determine if FLCN H255Y and K508R missense mutations promote aberrant kidney cell proliferation leading to pathogenicity, we generated mouse models expressing these mutants using BAC recombineering technology and investigated their ability to rescue the multi-cystic phenotype of Flcn-deficient mouse kidneys. Flcn H255Y mutant transgene expression in kidney-targeted Flcn knockout mice did not rescue the multi-cystic kidney phenotype. However, expression of the Flcn K508R mutant transgene partially, but not completely, abrogated the phenotype. Notably, expression of the Flcn K508R mutant transgene in heterozygous Flcn knockout mice resulted in development of multi-cystic kidneys and cardiac hypertrophy in some mice. These results demonstrate that both FLCN H255Y and K508R missense mutations promote aberrant kidney cell proliferation, but to different degrees. Based on the phenotypes of our preclinical models, the FLCN H255Y mutant protein has lost it tumour suppressive function leading to the clinical manifestations of BHD, whereas the FLCN K508R mutant protein may have a dominant negative effect on the function of wild-type FLCN in regulating kidney cell proliferation and, therefore, act as an oncoprotein. These findings may provide mechanistic insight into the role of FLCN in regulating kidney cell proliferation and facilitate the development of novel therapeutics for FLCN-deficient kidney cancer.

Network-based modelling and percolation analysis of conformational dynamics and activation in the CDK2 and CDK4 proteins: dynamic and energetic polarization of the kinase lobes may determine divergence of the regulatory mechanisms

Network modeling and percolation analysis of conformational dynamics and energetics of regulatory mechanisms in cyclin-dependent kinases.

Residue Geometry Networks: A Rigidity-Based Approach to the Amino Acid Network and Evolutionary Rate Analysis

Amino acid networks (AANs) abstract the protein structure by recording the amino acid contacts and can provide insight into protein function. Herein, we describe a novel AAN construction technique that employs the rigidity analysis tool, FIRST, to build the AAN, which we refer to as the residue geometry network (RGN). We show that this new construction can be combined with network theory methods to include the effects of allowed conformal motions and local chemical environments. Importantly, this is done without costly molecular dynamics simulations required by other AAN-related methods, which allows us to analyse large proteins and/or data sets. We have calculated the centrality of the residues belonging to 795 proteins. The results display a strong, negative correlation between residue centrality and the evolutionary rate. Furthermore, among residues with high closeness, those with low degree were particularly strongly conserved. Random walk simulations using the RGN were also successful in identifying allosteric residues in proteins involved in GPCR signalling. The dynamic function of these residues largely remain hidden in the traditional distance-cutoff construction technique. Despite being constructed from only the crystal structure, the results in this paper suggests that the RGN can identify residues that fulfil a dynamical function.

Requirement of novel amino acid fragments of orphan nuclear receptor TR3/Nur77 for its functions in angiogenesis

Pathological angiogenesis is a hallmark of many diseases. We demonstrated that TR3/Nur77 is an excellent target for pro-angiogenesis and anti-angiogenesis therapies. Here, we report that TR3 transcriptionally regulates endothelial cell migration, permeability and the formation of actin stress fibers that is independent of RhoA GTPase. 1) Amino acid residues 344-GRR-346 and de-phosphorylation of amino acid residue serine 351 in the DNA binding domain, and 2) phosphorylation of amino acid residues in the 41-61 amino acid fragment of the transactivation domain, of TR3 are required for its induction of the formation of actin stress fibers, cell proliferation, migration and permeability. The 41-61 amino acid fragment contains one of the three potential protein interaction motifs in the transactivation domain of TR3, predicted by computational modeling and analysis. These studies further our understanding of the molecular mechanism, by which TR3 regulates angiogenesis, identify novel therapeutic targeted sites of TR3, and set the foundation for the development of high-throughput screening assays to identify compounds targeting TR3/Nur77 for pro-angiogenesis and anti-angiogenesis therapies.

Rv2744c Is a PspA Ortholog That Regulates Lipid Droplet Homeostasis and Nonreplicating Persistence in Mycobacterium tuberculosis

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), remains a significant cause of morbidity and mortality worldwide, despite the availability of a live attenuated vaccine and anti-TB antibiotics. The vast majority of individuals infected with M. tuberculosis develop an asymptomatic latent infection in which the bacterium survives within host-generated granulomatous lesions in a physiologically altered metabolic state of nonreplicating persistence. The granuloma represents an adverse environment, as M. tuberculosis is exposed to various stressors capable of disrupting the essential constituents of the bacterium. In Gram-negative and Gram-positive bacteria, resistance to cell envelope stressors that perturb the plasma membrane is mediated in part by proteins comprising the phage shock protein (Psp) system. PspA is an important component of the Psp system; in the presence of envelope stress, PspA localizes to the inner face of the plasma membrane, homo-oligomerizes to form a large scaffold-like complex, and helps maintain plasma membrane integrity to prevent a loss of proton motive force. M. tuberculosis and other members of the Mycobacterium genus are thought to encode a minimal functional unit of the Psp system, including an ortholog of PspA. Here, we show that Rv2744c possesses structural and physical characteristics that are consistent with its designation as a PspA family member. However, although Rv2744c is upregulated under conditions of cell envelope stress, loss of Rv2744c does not alter resistance to cell envelope stressors. Furthermore, Rv2744c localizes to the surface of lipid droplets in Mycobacterium spp. and regulates lipid droplet number, size, and M. tuberculosis persistence during anaerobically induced dormancy. Collectively, our results indicate that Rv2744c is a bona fide ortholog of PspA that may function in a novel role to regulate lipid droplet homeostasis and nonreplicating persistence (NRP) in M. tuberculosis

IMPORTANCE

Mycobacterium tuberculosis is the causative agent of tuberculosis, a disease associated with significant morbidity and mortality worldwide. M. tuberculosis is capable of establishing lifelong asymptomatic infections in susceptible individuals and reactivating during periods of immune suppression to cause active disease. The determinants that are important for persistent infection of M. tuberculosis or for reactivation of this organism from latency are poorly understood. In this study, we describe our initial characterizations of Rv2744c, an ortholog of phage shock protein A (PspA) that regulates the homeostasis of lipid bodies and nonreplicating persistence in M. tuberculosis This function of PspA in M. tuberculosis is novel and suggests that PspA may represent a unique bacterial target upon which to base therapeutic interventions against this organism.

The linker region of breast cancer resistance protein ABCG2 is critical for coupling of ATP-dependent drug transport

The ATP-binding cassette (ABC) transporters of class G display a different domain organisation than P-glycoprotein/ABCB1 and bacterial homologues with a nucleotide-binding domain preceding the transmembrane domain. The linker region connecting these domains is unique and its function and structure cannot be predicted. Sequence analysis revealed that the human ABCG2 linker contains a LSGGE sequence, homologous to the canonical C-motif/ABC signature present in all ABC nucleotide-binding domains. Predictions of disorder and of secondary structures indicated that this C2-sequence was highly mobile and located between an α-helix and a loop similarly to the C-motif. Point mutations of the two first residues of the C2-sequence fully abolished the transport-coupled ATPase activity, and led to the complete loss of cell resistance to mitoxantrone. The interaction with potent, selective and non-competitive, ABCG2 inhibitors was also significantly altered upon mutation. These results suggest an important mechanistic role for the C2-sequence of the ABCG2 linker region in ATP binding and/or hydrolysis coupled to drug efflux.

Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis

Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability.

Protein thermostability is a crucial factor for biotechnological enzyme applications. However, performance studies of computational approaches for predicting effects of mutations on protein (thermo)stability have suggested that there is still room for improvement. We describe and validate a novel and unique strategy to predict optimal amino acid substitutions at structural weak spots. At variance with other rational approaches, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state is the underlying cause for thermostability. When applied prospectively on lipase LipA from Bacillus subtilis, a high success rate in predicting thermostabilized lipase variants and a remarkably large increase in their thermostability is achieved. This demonstrates the value of the novel strategy, which extends the existing portfolio of methods for rational protein design.

Prediction of Phenotypic Effects of Variants Observed in LOC_Os04g36720 of FRO1 Gene in Rice (Oryza sativa L.)

In rice, ferric-chelate reductase-1 (FRO1) (LOC_Os04g36720) gene was present on chromosome number 4 and its beginning and ending coordinates where coding sequence lies are 22182599 and 22186943, respectively. It plays a vital role in metal homeostasis and iron transportation in plants. Based on the alignment results, location of single-nucleotide variants is located in open reading frame and their effects of variants were predicted using SIFT sequence tool. The non-synonymous variants at position 342 and 436 lies in helical and coil parts of the protein, respectively, as predicted by Psi-pred server. PSI-Blast which resulted in significant hits and the most similar protein sequence (Accession ID: NP_001052896) with available sequence features displayed 100 % identity with query cover of 99 %. Results suggest the non-synonymous variant at position 436 (Accession ID: TBGI204002) lies in FAD-binding domain and nsSNV at position 342 (Accession ID: TBGI203998) lies in periphery of NADP. The SNPs were also analyzed for the deleterious effect by PANTHER subPSEC scores and I-mutant score, and it was postulated that SNPs would be hampering on biological as well as molecular function of FRO 1 gene of rice. A cutoff of -3 corresponds to a 50 % probability that a score is deleterious. From this, the probability that a given variant will cause a deleterious effect on protein function is estimated by P _deleterious, such that a subPSEC score of -4 corresponds to a P _deleterious of 0.79. Hence, to study the phenotypic consequences of variant TBGI204002, we performed comparative molecular docking studies of native modeled protein and protein with induced mutation as receptors and FAD as ligand to be utilized for binding. The docking process was performed by AutoDock 4.2 software with Lamarckian Genetic algorithm as computational algorithm. Results suggest binding energies are higher in case of mutation-induced protein which suggests presence of variant TBGI204002 enhances binding of FAD ligand at FAD-binding domain site. In case of TBGI203998, similar comparative docking procedure was performed with FAD as binding ligand, which suggests presence of variant does not impact FAD binding at the domain site. We revealed impact of SNPs on the protein structure and its function using sequence-based tools.

Bovine papillomavirus type 7 in Italy: complete genomes and sequence variants

Two novel bovine papillomavirus type 7 (BPV-7) variants have been identified in teat cutaneous papillomas affecting dairy cows in northern Italy. The entire genome sequences of two BPV-7 Italian variants showed major sequence differences in the long control region (LCR) and in the L2 gene compared to the Japanese reference strain. In order to define the stability of these genetic variants, the L2 and LCR sequences of seven further BPV-7 positive isolates were characterized. An insertion of six amino acids in the L2 structural protein has been detected in all samples while different genetic variants have been identified for the LCR. These findings provide new insights on intra-type variability of BPVs and represent a starting point for future studies aimed at establishing the biological role of the different BPV genomic regions and investigating the pathogenic potential of papillomavirus variants.

Molecular Cloning, Expression Analysis, and Functional Characterization of the H(+)-Pyrophosphatase from Jatropha curcas

H(+)-pyrophosphatase (H(+)-PPase) is a primary pyrophosphate (PPi)-energized proton pump to generate electrochemical H(+) gradient for ATP production and substance translocations across membranes. It plays an important role in stress adaptation that was intensively substantiated by numerous transgenic plants overexpressing H(+)-PPases yet devoid of any correlated studies pointing to the elite energy plant, Jatropha curcas. Herein, we cloned the full length of J. curcas H(+)-PPase (JcVP1) complementary DNA (cDNA) by reverse transcription PCR, based on the assembled sequence of its ESTs highly matched to Hevea brasiliensis H(+)-PPase. This gene encodes a polypeptide of 765 amino acids that was predicted as a K(+)-dependent H(+)-PPase evolutionarily closest to those of other Euphorbiaceae plants. Many cis-regulatory elements relevant to environmental stresses, molecular signals, or tissue-specificity were identified by promoter prediction within the 1.5-kb region upstream of JcVP1 coding sequence. Meanwhile, the responses of JcVP1 expression to several common abiotic stresses (salt, drought, heat, cold) were characterized with a considerable accordance with the inherent stress tolerance of J. curcas. Moreover, we found that the heterologous expression of JcVP1 could significantly improve the salt tolerance in both recombinant Escherichia coli and Saccharomyces cerevisiae, and this effect could be further fortified in yeast by N-terminal addition of a vacuole-targeting signal peptide from the H(+)-PPase of Trypanosoma cruzi.

The Interaction of FABP with Kapα

Gene-activating lipophilic compounds are carried into the nucleus when loaded on fatty-acid-binding proteins (FABP). Some of these proteins are recognized by the α-Karyopherin (Kapα) through its nuclear localization signal (NLS) consisting of three positive residues that are not in a continuous sequence. The Importin system can distinguish between FABP loaded with activating and non-activating compounds. In the present study, we introduced molecular dynamics as a tool for clarifying the mechanism by which FABP4, loaded with activating ligand (linoleate) is recognized by Kapα. In the first phase, we simulated the complex between Kapα^ΔIBB (termed “Armadillo”) that was crystallized with two NLS hepta-peptides. The trajectory revealed that the crystal-structure orientation of the peptides is rapidly lost and new interactions dominate. Though, the NLS sequence of FABP4 is cryptic, since the functional residues are not in direct sequence, implicating more than one possible conformation. Therefore, four possible docked conformations were generated, in which the NLS of FABP4 is interacting with either the major or the minor sites of Kapα, and the N → C vectors are parallel or anti-parallel. Out of these four basic starting positions, only the FABP4-minor site complex exhibited a large number of contact points. In this complex, the FABP interacts with the minor and the major sites, suppressing the self-inhibitory interaction of the Kapα, rendering it free to react with Kapβ. Finally, we propose that the transportable conformation generated an extended hydrophobic domain which expanded out of the boundary of the FABP4, allowing the loaded linoleate to partially migrate out of the FABP into a joint complex in which the Kapα contributes part of a combined binding pocket.