Comparative molecular dynamics simulation studies for determining factors contributing to the thermostability of chemotaxis protein “CheY”

Postharvest ripening of newly harvested corn: Weakened interactions between starch and protein

The effects of postharvest ripening of corn on the mechanisms of starch and protein interactions were investigated using molecular dynamics and several chemical substances. Sodium dodecyl sulfate (SDS) treatment all significantly affected the starch content, molecular weight of proteins, relative crystallinity, pasting characteristics and dynamic viscoelasticity in samples before and after postharvest ripening. In the corn that had not undergone postharvest ripening, there were also significant electrostatic interactions and hydrogen bonds between starch and protein. In addition, molecular dynamics had demonstrated that the forces between starch and protein in corn were mainly hydrophobic interactions, electrostatic interaction, and hydrogen bonds. Compared with zein, corn glutelin was more tightly bound to starch. The binding energy of starch to both proteins was reduced in after postharvest-ripened corn. This study laid a rationale for investigating the change mechanism of corn postharvest ripening quality and improving processing property and edible quality of corn.

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

Understanding the basis of thermostability for enzyme "Nanoluc" towards designing industry-competent engineered variants

As a leading contender in the study of luminescence, nanoluciferase has recently attracted attention and proven effective in a wide variety of research areas. Although numerous attempts have been made to improve activity, there has yet to be a thorough exploration of further possibilities to improve thermostability. In this study, protein engineering in tandem with molecular dynamics simulation at various temperatures (300 K, 400 K, 450 K and 500 K) was used to improve our understanding of nanoluciferase dynamics and identification of factors that could significantly enhance the thermostability. Based on these, three novel mutations have been narrowed down, which were hypothesised to improve thermostability. Root mean square deviation and root mean square fluctuation studies confirmed higher stability of mutant at high temperature. Solvent-accessible surface area and protein unfolding studies revealed a decreased tendency of mutant to unfold at higher temperatures. Further free energy landscape and principal component analysis was adapted to get deeper insights into the thermodynamic and structural behavior of these proteins at elevated temperature. Thus, this study provides a deeper insight into the dynamic factors for thermostability and introduces a novel, enhanced nanoluciferase candidate with potential use in industry.Communicated by Ramaswamy H. Sarma.

Simulation and practical investigation of carbonic anhydrase stability in an industrial solvent system of methyl diethanolamine for carbon dioxide capture

Carbonic anhydrase owing to its potential as an industrial biocatalyst for carbon dioxide sequestration from flue gas has attracted considerable attention in solving global warming problems. A large body of research has been conducted to increase the thermal stability of carbonic anhydrase from different sources against the harsh operational conditions of CO2 capture systems. In contrast to cost-intensive protein engineering methods, solvation with aqueous-organic binary mixtures offers a convenient and economical alternative strategy for retention of protein structure and stability. This study aimed to examine the stabilizing effect of methyl diethanolamine (MDEA) as a component of an aqueous-organic solvent mixture on human carbonic anhydrase II (HCA II) at extreme temperatures. Computational and also spectroscopic examinations were employed for tracking conformational changes and stability evaluation of HCA II in 50:50 (vol %) water: MDEA binary mixture at high temperature. Molecular dynamic (MD) simulation studies predicted the high thermal stability of HCA II in the presence of MDEA. UV absorbance spectra confirmed the thermo-stabilizing effect of the binary solvent mixture on HCA II. While the enzymatic activity of HCA II at 25 °C in the presence of 10, 25, and 50 (vol%) of MDEA was substantially increased, no obvious effect on retention of HCA II activity in the water-MDEA binary solvent mixture at 85 °C was seen. It is shown that the solvation of HCA II in the presence of MDEA could result in the prevention of aggregate formation in high temperatures but not functional stability.Communicated by Ramaswamy H. Sarma.

Potentiality of bioactive compounds as inhibitor of M protein and F protein function of human respiratory syncytial virus

The human respiratory syncytial virus (RSV) creates a pandemic every year in several countries in the world. Lack of target therapeutics and absence of vaccines have prompted scientists to create novel vaccines or small chemical treatments against RSV's numerous targets. The matrix (M) protein and fusion (F) glycoprotein of RSV are well characterized and attractive drug targets. Five bioactive compounds from Alnus japonica (Thunb.) Steud. were taken into consideration as lead compounds. Drug-likeness characters of them showed the drugs are non-toxic and non-mutagenic and mostly lipophobic. Molecular docking reveals that all bioactive compounds have better binding and better inhibitory effect than ribavirin which is currently used against RSV. Praecoxin A appeared as the best lead compound between them. It creates 7 different types of bonds with amino acids of M protein and 5 different types of bonds with amino acids of F protein. Van der Waals interactions highly influenced the binding energies. Molecular dynamic simulations represent the non-deviated and less fluctuating nature of praecoxin A. Principal Component Analysis showed praecoxin A complex with RSV matrix protein is more stable than ribavirin complex. This study will help to develop a new drug to inhibit RSV. All ligands were minimized through semi-empirical PM3 process with MOPAC. Toxicity was tested by ProTox-II server. Molecular docking studies were carried out using AutoDock 4.2. Molecular dynamics simulations for 100 ns were carried out through GROMACS 5.12 MD and GROMOS96 43a1 force field. The graphs were produced by GROMACS's XMGrace program.

Graphical abstract

Genomic analysis of Paenibacillus sp. MDMC362 from the Merzouga desert leads to the identification of a potentially thermostable catalase

Microorganisms in hot deserts face heat and other environmental conditions, such as desiccation, UV radiation, or low nutrient availability. Therefore, this hostile environment harbour microorganisms with acquired characteristics related to survival in their habitat, which can be exploited in biotechnology. In this work, the genome of Paenibacillus sp. MDMC362 isolated from the Merzouga desert in Morocco was sequenced to understand its survival strategy's genetic basis; and to evaluate the thermostability of a catalase extracted from genomic annotation files using molecular dynamics. Paenibacillus sp. MDMC362 genome was rich in genetic elements involved in the fight against different stresses, notably temperature stress, UV radiations, osmotic stress, carbon starvation, and oxidative stress. Indeed, we could identify genes of the operons groES-groEL and hrcA-grpE-dnaK and those involved in the different stages of sporulation, which can help the bacteria to survive the high temperatures imposed by a desertic environment. We also observed the genetic components of the UvrABC system and additional mechanisms involved in DNA repair, which help overcome UV radiation damage. Other genes have been identified in the genome, like those coding for ectoine and proline, that aids fight osmotic stress and desiccation. Catalase thermostability investigation using molecular dynamics showed that the protein reached stability and conserved its compactness at temperatures up to 373.15 K. These results suggest a potential thermostability of the enzyme. Since the studied protein is a core protein, thermostability could be conserved among Paenibacillus sp. MDMC362 closely related strains; however, bacteria from harsh environments may have a slight advantage regarding protein stability.

In silico modeling revealed new insights into the mechanism of action of enzyme 2'-5'-oligoadenylate synthetase in cattle

The innate immune system has an important role in developing the initial resistance to virus infection, and the ability of oligoadenylate synthetase to overcome viral evasion and enhance innate immunity is already established in humans. In the present study, we have tried to explore the molecular and structural variations present in Sahiwal (indigenous) and crossbred (Frieswal) cattle to identify the molecular mechanism of action of OAS1 gene in activation of innate immune response. The significant changes in structural alignment in terms of orientation of loops, shortening of β-sheets and formation of 3-10 α-helix was noticed in Sahiwal and Frieswal cattle. Further, it has been observed that OAS1 from Sahiwal had better binding with APC and DTP ligand than Frieswal OAS1. A remarkable change was seen in orientation at the nucleoside base region of both the ligands, which are bound with OAS1 protein from Frieswal and Sahiwal cattle. The Molecular Dynamic study of apo and ligand complex structures was provided more insight towards the stability of OAS1 from both cattle. This analysis displayed that the Sahiwal cattle protein has more steady nature throughout the simulation and has better binding towards Frieswal in terms of APC and DTP binding. Thus, OAS1 protein is the potential target for explaining the innate immune response in Sahiwal than Frieswal.Communicated by Ramaswamy H. Sarma.

In silico study of cox protein from P2 type enteric bacteriophages based on sequence, structure and dynamics to understand its functional integrity

Cox protein plays a critical role in deciding the lytic-lysogenic switch of P2 enteric phages. This phenomenon makes Cox protein one of the most important candidates in developing novel phage-based therapeutics against antibacterial resistant pathogens. The principle focus concerning protein and its decision making is a DNA binding event, which helps to regulate differential promoter expression. In the current study, we have attempted to understand the sequence, structural and dynamic features associated with Cox protein and its DNA binding. Unavailability of information was a big burden in further proceedings. We have done an extensive literature search to develop a database of Cox with relevant information. That information coupled with the methods of Sequence-based phylogenetic and conservation studies, Homology Modelling, Atomic-level Docking and Molecular Dynamics (MD) Simulation (50 ns each for 10 systems, i.e. total of 500 ns) were performed in the current study. Analysis of those extensive studies has provided us the required sequence to structure to dynamics to functional understanding. Our present study would indeed be very helpful in understanding the biochemical mechanism of Cox activation as well as designing potential phage therapeutics.

Study of potentiality of dexamethasone and its derivatives against Covid-19

In December 2019, COVID-19 epidemic was reported in Wuhan, China, and subsequently the infection has spread all over the world and became pandemic. The death toll associated with the pandemic is increasing day by day in a high rate. Herein, an effort has been made to identify the potentiality of commercially available drugs and also their probable derivatives for creation of better opportunity to make more powerful drugs against coronavirus. This study involves the in-silico interactions of dexamethasone and its derivatives against the multiple proteins of SARS-CoV-2 with the help of various computational methods. Descriptor parameters revealed their non-toxic effect in the human body. Ultimately docking studies and molecular dynamic simulation on those target protein by dexamethasone and its derivatives showed a high binding energy. Dexamethasone showed -9.8 kcal/mol and its derivative D5 showed -12.1 kcal/mol binding energy. Those scores indicate that dexamethasone has more therapeutic effect on SARS CoV-2 than other currently used drugs. Derivatives give the clue for the synthesis of a novel drug to remove SARS CoV-2. Until then, dexamethasone will be used as a potential inhibitor of SARS CoV-2.Communicated by Ramaswamy H. Sarma.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Interdisciplinary in silico studies to understand in-depth molecular level mechanism of drug resistance involving NS3-4A protease of HCV

Hepatitis C virus (HCV) causes hepatitis, a life-threatening disease responsible for liver cirrhosis. Urgent measures have been taken to develop therapeutics against this deadly pathogen. NS3/4A protease is an extremely important target. A series of inhibitors have been developed against this viral protease including Faldaprevir. Unfortunately, the error-prone viral RNA polymerase causes the emergence of resistance, thereby causing reduced effectiveness of those peptidomimetic inhibitors. Among the drug resistant variants, three single amino acid residues (R155, A156 and D168) are notable for their presence in clinical isolates and also their effectivity against most of the known inhibitors in clinical development. Therefore, it is crucial to understand the mechanistic role of those drug resistant variants while designing potent novel inhibitors. In this communication, we have deeply analyzed through using in silico studies to understand the molecular mechanism of alteration of inhibitor binding between wild type and its R155K, A156V and D168V variants. Principal component analysis was carried to identify the backbone fluctuations of important residues in HCV NS3/4A responsible for the inhibitor binding and maintaining drug resistance. Free energy landscape as a function of the principal components has been used to identify the stability and conformation of the key residues that regulate inhibitor binding and their impact in developing drug resistance. Our findings are consistent with the trend of experimental results. The observations are also true in case of other Faldaprevir-like peptidomimetic inhibitors. Understanding this binding mechanism would be significant for the development of novel inhibitors with less susceptibility towards drug resistance.

Optimization of IL-1RA structure to achieve a smaller protein with a higher affinity to its receptor

Interleukine-1 family cytokines are key orchestrators of innate and adaptive immunity. In particular, up-regulation of IL-1R1 via its agonistic ligands consisting of IL-1β and IL-1α is implicated in a variety of human diseases, such as rheumatoid arthritis, psoriasis, type I diabetes, amyotrophic lateral sclerosis, and dry-eye disease. Until now, there are no small-molecule inhibitors of the IL-1R1 with increased antagonistic potency to be used for the treatment of peripheral inflammation. The objective of this study was to engineer a low-molecular-weight version of IL-1RA with increased affinity and enhanced antagonistic activity for potential therapeutic use. To develop a smaller protein–ligand with a better affinity to IL-1R, we used bioinformatics studies and in silico simulations to anticipate non-binding areas on IL-1RA. In this study, we have identified a 41aa (F57-F98) non-binding site of IL-1RA. Overall RMSF of the Truncated complex (1.5 nm) was lower than the Native complex (2 nm), which could prove higher stability of the Truncated complex. The free binding energy of the T-IL-1RA (− 1087.037 kJ/mol) was significantly lower than the IL-1RA (− 836.819 kJ/mol) which could demonstrate a higher binding affinity of the truncated ligand with its receptor as a result of new important interactions. These findings have demonstrated a higher binding affinity of the T-IL-1RA with its receptor than the native protein. These results should: have an impact on the development of new treatments that block IL-1 signaling, although more research is needed in vitro and in vivo.

Synthesis of Dihydrobenzofuro[3,2-b]chromenes as a potential 3CLpro inhibitors of SARS-CoV-2: A molecular docking and dynamics simulation study

The recent emergence of pandemic of coronavirus (COVID‐19) caused by SARS‐CoV‐2 has raised significant global health concerns. More importantly, there is no specific therapeutics currently available to combat against this deadly infection. The enzyme 3‐chymotrypsin‐like cysteine protease (3CLpro) is known to be essential for viral life cycle as it controls the coronavirus replication. 3CLpro could be a potential drug target as established before in the case of severe acute respiratory syndrome coronavirus (SARS‐CoV) and Middle East respiratory syndrome coronavirus (MERS‐CoV). In the current study, we wanted to explore the potential of fused flavonoids as 3CLpro inhibitors. Fused flavonoids (5a,10a‐dihydro‐11H‐benzofuro[3,2‐b]chromene) are unexplored for their potential bioactivities due to their low natural occurrences. Their synthetic congeners are also rare due to unavailability of general synthetic methodology. Here we designed a simple strategy to synthesize 5a,10a‐dihydro‐11H‐benzofuro[3,2‐b]chromene skeleton and it's four novel derivatives. Our structural bioinformatics study clearly shows excellent potential of the synthesized compounds in comparison to experimentally validated inhibitor N3. Moreover, in‐silico ADMET study displays excellent druggability and extremely low level of toxicity of the synthesized molecules. Further, for better understanding, the molecular dynamic approach was implemented to study the change in dynamicity after the compounds bind to the protein. A detailed investigation through clustering analysis and distance calculation gave us sound comprehensive data about their molecular interaction. In summary, we anticipate that the currently synthesized molecules could not only be a potential set of inhibitors against 3CLpro but also the insights acquired from the current study would be instrumental in further developing novel natural flavonoid based anti‐COVID therapeutic spectrums.

Docking-guided rational engineering of a macrolide glycosyltransferase glycodiversifies epothilone B

Glycosyltransferases typically display acceptor substrate flexibility but more stringent donor specificity. BsGT-1 is a highly effective glycosyltransferase to glycosylate macrolides, including epothilones, promising antitumor compounds. Here, we show that BsGT-1 has three major regions significantly influencing the glycodiversification of epothilone B based on structural molecular docking, "hot spots" alanine scanning, and site saturation mutagenesis. Mutations in the PSPG-like motif region and the C2 loop region are more likely to expand donor preference; mutations of the flexible N3 loop region located at the mouth of the substrate-binding cavity produce novel epothilone oligosaccharides. These "hot spots" also functioned in homologues of BsGT-1. The glycosides showed significantly enhanced water solubility and decreased cytotoxicity, although the glycosyl appendages of epothilone B also reduced drug permeability and attenuated antitumor efficacy. This study laid a foundation for the rational engineering of other GTs to synthesize valuable small molecules.

N-terminal truncation of VC0395_0300 protein from Vibrio cholerae does not lead to loss of diguanylate cyclase activity

The bacterial secondary messenger bis-(3',5')-cyclic-dimeric-guanosine monophosphate (c-di-GMP) has been implicated in the pathogenesis of Vibrio cholerae, due to its significant role in regulating the virulence, biofilm formation and motility of the host organism. The VC0395_0300 protein from V. cholerae, possessing a GGEEF sequence has been established as a diguanylate cyclase (DGC) capable of catalyzing the conversion of two GTP molecules to form cyclic-di-GMP. This in turn, plays a crucial role in allowing the organism to adopt a dual lifestyle, thriving both in human and aquatic systems. The difficulty in procuring sufficient amounts of homogenous soluble protein for structural assessment of the GGDEF domain in VC0395_0300 and the lack of soluble protein yield, prompted the truncation into smaller constructs (Sebox31 and Sebox32) carrying the GGDEF domain. The truncates retained their diguanylate cyclase activity comparable to the wild type, and were able to form biofilms as well. Fluorescence and circular dichroism spectroscopy measurements revealed that the basic structural elements do not show significant changes in the truncated proteins as compared to the full-length. This has also been confirmed using homology modeling and molecular docking of the wild type and truncates. This led us to conclude that the truncated constructs retain their activity in spite of the deletions in the N terminal region. This is supportive of the fact that DGC activity in GGDEF proteins is predominantly dependent on the presence of the conserved GGD(/E)EF domain and its interaction with GTP.

Comparative modelling studies of fruit bromelain using molecular dynamics simulation

Fruit bromelain is a cysteine protease accumulated in pineapple fruits. This proteolytic enzyme has received high demand for industrial and therapeutic applications. In this study, fruit bromelain sequences QIM61759, QIM61760 and QIM61761 were retrieved from the National Center for Biotechnology Information (NCBI) Genbank Database. The tertiary structure of fruit bromelain QIM61759, QIM61760 and QIM61761 was generated by using MODELLER. The result revealed that the local stereochemical quality of the generated models was improved by using multiple templates during modelling process. Moreover, by comparing with the available papain model, structural analysis provides an insight on how pro-peptide functions as a scaffold in fruit bromelain folding and contributing to inactivation of mature protein. The structural analysis also disclosed the similarities and differences between these models. Lastly, thermal stability of fruit bromelain was studied. Molecular dynamics simulation of fruit bromelain structures at several selected temperatures demonstrated how fruit bromelain responds to elevation of temperature.

Understanding structure-based dynamic interactions of antihypertensive peptides extracted from food sources

Functional foods are emerging as essential healthy nutritional component due to their abundant wellbeing benefits. Especially the food-derived peptides are considered as key components for playing their biologically active roles. One such robust therapeutics that already exploited with food peptides that help treating high blood pressure via targeting Angiotensin-Converting Enzyme (ACE). This in silico study demonstrated the inhibitory potential of antihypertensive peptides derived from food sources. This study involves an intensive structure-based analysis of enzyme-peptide interactions using Molecular Dynamics (MD) simulations. Interestingly, this study will help us to get deeper understanding on how food peptides achieve successful inhibition of ACE. In this study, the peptide-enzyme complexes revealed two binding pockets, A and B, on either side of the active site Zn atom. Pocket B has a smaller binding site volume than pocket A, comprised of β-sheets and the active site opening cleft. The interface of the binding sites showed that the enzyme structure was negative to neutral charge, and the peptide structure was positive to neutral charge. The dynamics of complex structures of seven highly potential peptides were performed for 20 ns each at 300 K. Comparative analysis of RMSD, RMSF and binding energies show the enzyme-peptide complexes and the overall stability of apo-enzyme. Importantly, two peptides AFKAWAVAR and IWHHTF showed the highest variation in their RMSD as compared to the apo-enzyme. This study will further be useful for the assessment of the characteristics to predict novel inhibitory peptides that can be generated from food proteins.Communicated by Ramaswamy H. Sarma.

PI3Kα-regulated gelsolin activity is a critical determinant of cardiac cytoskeletal remodeling and heart disease

Biomechanical stress and cytoskeletal remodeling are key determinants of cellular homeostasis and tissue responses to mechanical stimuli and injury. Here we document the increased activity of gelsolin, an actin filament severing and capping protein, in failing human hearts. Deletion of gelsolin prevents biomechanical stress-induced adverse cytoskeletal remodeling and heart failure in mice. We show that phosphatidylinositol (3,4,5)-triphosphate (PIP3) lipid suppresses gelsolin actin-severing and capping activities. Accordingly, loss of PI3Kα, the key PIP3-producing enzyme in the heart, increases gelsolin-mediated actin-severing activities in the myocardium in vivo, resulting in dilated cardiomyopathy in response to pressure-overload. Mechanical stretching of adult PI3Kα-deficient cardiomyocytes disrupts the actin cytoskeleton, which is prevented by reconstituting cells with PIP3. The actin severing and capping activities of recombinant gelsolin are effectively suppressed by PIP3. Our data identify the role of gelsolin-driven cytoskeletal remodeling in heart failure in which PI3Kα/PIP3 act as negative regulators of gelsolin activity.

Developing Hispolon-based novel anticancer therapeutics against human (NF-κβ) using in silico approach of modelling, docking and protein dynamics

Hispolon is a polyphenolic compound derived from black hoof mushroom (Phellinus linteus) or shaggy bracket mushroom (Inonotus hispidus) which induces the inhibition of cancer-promoting nuclear factor-kappa beta (NF-κβ) complex. To develop more potent lead molecules with enhanced anticancer efficiency, the mechanism of hispolon-mediated nuclear factor-κβ inhibition has been investigated by molecular modelling and docking. Ten derivatives of hispolon (DRG1-10) have been developed by pharmacophore-based design with a view to enhance the anticancer efficacy. Hispolon and its derivatives were further screened for different pharmacological parameters like binding free energy, drug likeliness, absorption-digestion-metabolism-excretion (ADME), permeability, mutagenicity, toxicity and inhibitory concentration 50 (IC50) to find a potent lead molecule. Based on pharmacological validation, comparative molecular dynamics (MD) simulations have been performed for three lead molecules: Hispolon, DRG2 and DRG7 complexed with human NF-κβ up to 50 ns. By analysing different factors like root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA) and principal component analysis (PCA), Gibb's free energy plots DRG2 have more binding efficiency compared to hispolon and DRG7. In RMSD plot, hispolon-bound NF-κβ has the most deviation within a range between 0.125 and 0.45 nm, and DRG2-bound complex showed the range between 0.125 and 0.25 nm. The residues of NF-κβ responsible for hydrophobic interactions with ligand, e.g. Met469, Leu522 and Cys533, have the lowest fluctuation values in DRG2-bound complex. The average Rg fluctuation for DRG2-bound NF-κβ has been recorded under 2.025 nm for most of the simulation time which is much less compared to hispolon and DRG7. Gibb's free energy plots also define the highest stability of DRG2-bound NF-κβ. Communicated by Ramaswamy H. Sarma.

Identification of differentially expressed microRNAs in Sahiwal (Bos indicus) breed of cattle during thermal stress

microRNAs (miRNAs) are a class of small non-coding RNAs that play key roles in post transcriptional gene regulation that influence various fundamental cellular processes, including the cellular responses during environmental stresses. However, perusal of literatures revealed few reports on the differential expression of miRNA during thermal stress in Indian native (Bos indicus) cattle breeds. The present investigation aimed to identify differentially expressed miRNAs during thermal stress in Sahiwal (Bos indicus) dairy cattle breed of India, adapted with tropical climate over a long period of time. Stress responses of the animals were characterized by determining various physiological as well as biochemical parameters and differential expression profile of major heat shock protein genes. Ion Torrent deep sequencing and CLC-genomic analysis identified a set of differentially expressed miRNAs during summer and winter seasons. Most of the identified differentially expressed miRNAs were found to target heat shock responsive genes especially members of heat shock protein (HSP) family. Real-time quantification-based analysis of selected miRNAs revealed that bta-mir-1248, bta-mir-2332, bta-mir-2478, and bta-mir-1839 were significantly (p < 0.01) over expressed while bta-mir-16a, bta-let-7b, bta-mir-142, and bta-mir-425 were significantly (p < 0.01) under expressed during summer in comparison to winter. The present study enlists differentially expressed miRNAs at different environmental temperatures in Sahiwal (Bos indicus) that may be importance for further understanding the role of miRNAs on thermo-regulatory mechanisms.

A consensus-guided approach yields a heat-stable alkane-producing enzyme and identifies residues promoting thermostability

Aldehyde-deformylating oxygenase (ADO) is an essential enzyme for production of long-chain alkanes as drop-in biofuels, which are compatible with existing fuel systems. The most active ADOs are present in mesophilic cyanobacteria, especially Nostoc punctiforme Given the potential applications of thermostable enzymes in biorefineries, here we generated a thermostable (Cts)-ADO based on a consensus of ADO sequences from several thermophilic cyanobacterial strains. Using an in silico design pipeline and a metagenome library containing 41 hot-spring microbial communities, we created Cts-ADO. Cts-ADO displayed a 3.8-fold increase in pentadecane production on raising the temperature from 30 to 42 °C, whereas ADO from N. punctiforme (Np-ADO) exhibited a 1.7-fold decline. 3D structure modeling and molecular dynamics simulations of Cts- and Np-ADO at different temperatures revealed differences between the two enzymes in residues clustered on exposed loops of these variants, which affected the conformation of helices involved in forming the ADO catalytic core. In Cts-ADO, this conformational change promoted ligand binding to its preferred iron, Fe2, in the di-iron cluster at higher temperature, but the reverse was observed in Np-ADO. Detailed mapping of residues conferring Cts-ADO thermostability identified four amino acids, which we substituted individually and together in Np-ADO. Among these substitution variants, A161E was remarkably similar to Cts-ADO in terms of activity optima, kinetic parameters, and structure at higher temperature. A161E was located in loop L6, which connects helices H5 and H6, and supported ligand binding to Fe2 at higher temperatures, thereby promoting optimal activity at these temperatures and explaining the increased thermostability of Cts-ADO.

Histidine substitution in the most flexible fragments of firefly luciferase modifies its thermal stability

Molecular dynamics (MD) at two temperatures of 300 and 340 K identified two histidine residues, His461 and His489, in the most flexible regions of firefly luciferase, a light emitting enzyme. We therefore designed four protein mutants H461D, H489K, H489D and H489M to investigate their enzyme kinetic and thermodynamic stability changes. Substitution of His461 by aspartate (H461D) decreased ATP binding affinity, reduced the melting temperature of protein by around 25 °C and shifted its optimum temperature of activity to 10 °C. In line with the common feature of psychrophilic enzymes, the MD data showed that the overall flexibility of H461D was relatively high at low temperature, probably due to a decrease in the number of salt bridges around the mutation site. On the other hand, substitution of His489 by aspartate (H489D) introduced a new salt bridge between the C-terminal and N-terminal domains and increased protein rigidity but only slightly improved its thermal stability. Similar changes were observed for H489K and, to a lesser degree, H489M mutations. Based on our results we conclude that the MD simulation-based rational substitution of histidines by salt-bridge forming residues can modulate conformational dynamics in luciferase and shift its optimal temperature activity.

Putative protein VC0395_0300 from Vibrio cholerae is a diguanylate cyclase with a role in biofilm formation

The hallmark of the lifecycle of Vibrio cholerae is its ability to switch between two lifestyles - the sessile, non-pathogenic form and the motile, infectious form in human hosts. One of these changes is in the formation of surface biofilms, when in sessile aquatic habitats. The cell-cell interactions within a V. cholerae biofilm are stabilized by the production of an exopolysachharide (EPS) matrix, which in turn is regulated by the ubiquitous secondary messenger, cyclic di-GMP (c-di-GMP), synthesized by proteins containing GGD(/E)EF domains in all prokaryotic systems. Here, we report the functional role of the VC0395_0300 protein (Sebox3) encoded by the chromosome I of V. cholerae, with a GGEEF signature sequence, in the formation of surface biofilms. In our study, we have shown that Escherichia coli containing the full-length Sebox3 displays enhanced biofilm forming ability with cellulose production as quantified and visualized by multiple assays, most notably using FEG-SEM. This has also been corroborated with the lack of motility of host containing Sebox3 in semi-solid media. Searching for the reasons for this biofilm formation, we have demonstrated in vitro that Sebox3 can synthesize c-di-GMP from GTP. The homology derived model of Sebox3 displayed significant conservation of the GGD(/E)EF architecture as well. Hence, we propose that the putative protein VC0395_0300 from V. cholerae is a diguanylate cyclase which has an active role in biofilm formation.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Two strategies to engineer flexible loops for improved enzyme thermostability

Flexible sites are potential targets for engineering the stability of enzymes. Nevertheless, the success rate of the rigidifying flexible sites (RFS) strategy is still low due to a limited understanding of how to determine the best mutation candidates. In this study, two parallel strategies were applied to identify mutation candidates within the flexible loops of Escherichia coli transketolase (TK). The first was a “back to consensus mutations” approach, and the second was computational design based on ΔΔG calculations in Rosetta. Forty-nine single variants were generated and characterised experimentally. From these, three single-variants I189H, A282P, D143K were found to be more thermostable than wild-type TK. The combination of A282P with H192P, a variant constructed previously, resulted in the best all-round variant with a 3-fold improved half-life at 60 °C, 5-fold increased specific activity at 65 °C, 1.3-fold improved k_cat and a T_m increased by 5 °C above that of wild type. Based on a statistical analysis of the stability changes for all variants, the qualitative prediction accuracy of the Rosetta program reached 65.3%. Both of the two strategies investigated were useful in guiding mutation candidates to flexible loops, and had the potential to be used for other enzymes.

Sequence-, structure-, and dynamics-based comparisons of structurally homologous CheY-like proteins

We recently introduced a physically based approach to sequence comparison, the property factor method (PFM). In the present work, we apply the PFM approach to the study of a challenging set of sequences-the bacterial chemotaxis protein CheY, the N-terminal receiver domain of the nitrogen regulation protein NT-NtrC, and the sporulation response regulator Spo0F. These are all response regulators involved in signal transduction. Despite functional similarity and structural homology, they exhibit low sequence identity. PFM sequence comparison demonstrates a statistically significant qualitative difference between the sequence of CheY and those of the other two proteins that is not found using conventional alignment methods. This difference is shown to be consonant with structural characteristics, using distance matrix comparisons. We also demonstrate that residues participating strongly in native contacts during unfolding are distributed differently in CheY than in the other two proteins. The PFM result is also in accord with dynamic simulation results of several types. Molecular dynamics simulations of all three proteins were carried out at several temperatures, and it is shown that the dynamics of CheY are predicted to differ from those of NT-NtrC and Spo0F. The predicted dynamic properties of the three proteins are in good agreement with experimentally determined B factors and with fluctuations predicted by the Gaussian network model. We pinpoint the differences between the PFM and traditional sequence comparisons and discuss the informatic basis for the ability of the PFM approach to detect physical differences between these sequences that are not apparent from traditional alignment-based comparison.

In Silico Designing of an Industrially Sustainable Carbonic Anhydrase Using Molecular Dynamics Simulation

Carbonic anhydrase (CA) is a family of metalloenzymes that has the potential to sequestrate carbon dioxide (CO₂) from the environment and reduce pollution. The goal of this study is to apply protein engineering to develop a modified CA enzyme that has both higher stability and activity and hence could be used for industrial purposes. In the current study, we have developed an in silico method to understand the molecular basis behind the stability of CA. We have performed comparative molecular dynamics simulation of two homologous α-CA, one of thermophilic origin (Sulfurihydrogenibium sp.) and its mesophilic counterpart (Neisseria gonorrhoeae), for 100 ns each at 300, 350, 400, and 500 K. Comparing the trajectories of two proteins using different stability-determining factors, we have designed a highly thermostable version of mesophilic α-CA by introducing three mutations (S44R, S139E, and K168R). The designed mutant α-CA maintains conformational stability at high temperatures. This study shows the potential to develop industrially stable variants of enzymes while maintaining high activity.

Structural and dynamic evolution of the amphipathic N-terminus diversifies enzyme thermostability in the glycoside hydrolase family 12

The N-terminus diversifies enzyme thermostability in the GH12 family, which was investigated by MD simulations, and provides potential applications in protein engineering.

Lysine-Based Site-Directed Mutagenesis Increased Rigid β-Sheet Structure and Thermostability of Mesophilic 1,3-1,4-β-Glucanase

1,3-1,4-β-Glucanase is widely applied in the food industry, while its low thermostability often reduces its performance. In a previous study, chemical modification of surface lysine residues was proved to increase the thermostability of β-glucanase. To improve the thermostability, the mesophilic β-glucanase from Bacillus terquilensis was rationally engineered through site-directed mutagenesis of the 12 lysines into serines. The results showed that the K20S, K117S, and K165S mutants could both enhance the specific activities and thermostability of β-glucanase. The triple mutant (K20S/K117S/K165S) could increase the optimal temperature and T50 value by 15 and 14 °C, respectively. Five percent more structured residues were observed in the mutant, which formed new β-sheet structures in the concave side. Molecular dynamics simulation analysis showed that the flexibility in the mutation regions was decreased, which resulted in the overall rigidity of the β-glucanase. Therefore, the lysine-based site-directed mutagenesis is a simple and effective method for improving the thermostability of β-glucanase.

Structure and functional characterization of pyruvate decarboxylase from Gluconacetobacter diazotrophicus

BACKGROUND

Bacterial pyruvate decarboxylases (PDC) are rare. Their role in ethanol production and in bacterially mediated ethanologenic processes has, however, ensured a continued and growing interest. PDCs from Zymomonas mobilis (ZmPDC), Zymobacter palmae (ZpPDC) and Sarcina ventriculi (SvPDC) have been characterized and ZmPDC has been produced successfully in a range of heterologous hosts. PDCs from the Acetobacteraceae and their role in metabolism have not been characterized to the same extent. Examples include Gluconobacter oxydans (GoPDC), G. diazotrophicus (GdPDC) and Acetobacter pasteutrianus (ApPDC). All of these organisms are of commercial importance.

RESULTS

This study reports the kinetic characterization and the crystal structure of a PDC from Gluconacetobacter diazotrophicus (GdPDC). Enzyme kinetic analysis indicates a high affinity for pyruvate (K M 0.06 mM at pH 5), high catalytic efficiencies (1.3 • 10(6) M(-1) • s(-1) at pH 5), pHopt of 5.5 and Topt at 45°C. The enzyme is not thermostable (T½ of 18 minutes at 60°C) and the calculated number of bonds between monomers and dimers do not give clear indications for the relatively lower thermostability compared to other PDCs. The structure is highly similar to those described for Z. mobilis (ZmPDC) and A. pasteurianus PDC (ApPDC) with a rmsd value of 0.57 Å for Cα when comparing GdPDC to that of ApPDC. Indole-3-pyruvate does not serve as a substrate for the enzyme. Structural differences occur in two loci, involving the regions Thr341 to Thr352 and Asn499 to Asp503.

CONCLUSIONS

This is the first study of the PDC from G. diazotrophicus (PAL5) and lays the groundwork for future research into its role in this endosymbiont. The crystal structure of GdPDC indicates the enzyme to be evolutionarily closely related to homologues from Z. mobilis and A. pasteurianus and suggests strong selective pressure to keep the enzyme characteristics in a narrow range. The pH optimum together with reduced thermostability likely reflect the host organisms niche and conditions under which these properties have been naturally selected for. The lack of activity on indole-3-pyruvate excludes this decarboxylase as the enzyme responsible for indole acetic acid production in G. diazotrophicus.