Molecular dynamics simulations as a tool for improving protein stability

Robust ω-Transaminases by Computational Stabilization of the Subunit Interface

Transaminases are attractive catalysts for the production of enantiopure amines. However, the poor stability of these enzymes often limits their application in biocatalysis. Here, we used a framework for enzyme stability engineering by computational library design (FRESCO) to stabilize the homodimeric PLP fold type I ω-transaminase from Pseudomonas jessenii. A large number of surface-located point mutations and mutations predicted to stabilize the subunit interface were examined. Experimental screening revealed that 10 surface mutations out of 172 tested were indeed stabilizing (6% success), whereas testing 34 interface mutations gave 19 hits (56% success). Both the extent of stabilization and the spatial distribution of stabilizing mutations showed that the subunit interface was critical for stability. After mutations were combined, 2 very stable variants with 4 and 6 mutations were obtained, which in comparison to wild type (T _m ^app = 62 °C) displayed T _m ^app values of 80 and 85 °C, respectively. These two variants were also 5-fold more active at their optimum temperatures and tolerated high concentrations of isopropylamine and cosolvents. This allowed conversion of 100 mM acetophenone to (S)-1-phenylethylamine (>99% enantiomeric excess) with high yield (92%, in comparison to 24% with the wild-type transaminase). Crystal structures mostly confirmed the expected structural changes and revealed that the most stabilizing mutation, I154V, featured a rarely described stabilization mechanism: namely, removal of steric strain. The results show that computational interface redesign can be a rapid and powerful strategy for transaminase stabilization.

Learning the changes of barnase mutants thermostability from structural fluctuations obtained using anisotropic network modeling

In biotechnology applications, rational design of new proteins with improved physico-chemical properties includes a number of important tasks. One of the greatest practical and fundamental challenges is the design of highly thermostable protein enzymes that maintain catalytic activity at high temperatures. This problem may be solved by introducing mutations into the wild-type enzyme protein. In this work, to predict the impact of such mutations in barnase protein we applied the anisotropic network modeling approach, revealing atomic fluctuations in structural regions that are changed in mutants compared to the wild-type protein. A regression model was constructed based on these structural features that can allow one to predict the thermal stability of new barnase mutants. Moreover, the analysis of regression model provides a mechanistic explanation of how the structural features can contribute to the thermal stability of barnase mutants.

Parallel Strategy Increases the Thermostability and Activity of Glutamate Decarboxylase

Glutamate decarboxylase (GAD; EC 4.1.1.15) is a unique pyridoxal 5-phosphate (PLP)-dependent enzyme that specifically catalyzes the decarboxylation of L-glutamic acid to produce γ-aminobutyric acid (GABA), which exhibits several well-known physiological functions. However, glutamate decarboxylase from different sources has the common problem of poor thermostability that affects its application in industry. In this study, a parallel strategy comprising sequential analysis and free energy calculation was applied to identify critical amino acid sites affecting thermostability of GAD and select proper mutation contributing to improve structure rigidity of the enzyme. Two mutant enzymes, D203E and S325A, with higher thermostability were obtained, and their semi-inactivation temperature (T₅₀¹⁵) values were 2.3 °C and 1.4 °C higher than the corresponding value of the wild-type enzyme (WT), respectively. Moreover, the mutant, S325A, exhibited enhanced activity compared to the wild type, with a 1.67-fold increase. The parallel strategy presented in this work proved to be an efficient tool for the reinforcement of protein thermostability.

Molecular dynamics simulations suggest stabilizing mutations in a de novo designed α/β protein

Designing functional proteins that can withstand extreme heat is beneficial for industrial and protein therapeutic applications. Thus, elucidating the atomic-level determinants of thermostability is a major interest for rational protein design. To that end, we compared the structure and dynamics of a set of previously designed, thermostable proteins based on the activation domain of human procarboxypeptidase A2 (AYEwt). The mutations in these designed proteins were intended to increase hydrophobic core packing and inter-secondary-structure interactions. To evaluate whether these design strategies were successfully deployed, we performed all-atom, explicit-solvent molecular dynamics (MD) simulations of AYEwt and three designed variants at both 25 and 100°C. Our MD simulations agreed with the relative experimental stabilities of the designs based on their secondary structure content, Cα root-mean-square deviation/fluctuation, and buried-residue solvent accessible surface area. Using a contact analysis, we found that the designs stabilize inter-secondary structure interactions and buried hydrophobic surface area, as intended. Based on our analysis, we designed three additional variants to test the role of helix stabilization, core packing, and a Phe → Met mutation on thermostability. We performed the additional MD simulations and analysis on these variants, and these data supported our predictions.

From in vitro to in silico: Modeling and recombinant production of DT-Diaphorase enzyme

DT-Diaphorase (DTD) belonging to the oxidoreductase family, is among the most important enzymes and is of great significance in present-day biotechnology. Also, it has potential applications in glucose and pyruvate biosensors. Another important role of the DTD enzyme is in the detection of Phenylketonuria disease. According to the above demands, at first, we tried to study molecular cloning and production of recombinant DTD in E. coli BL21 strain. We have successfully cloned, expressed, and purified functionally active diaphorase. The amount of enzyme was increased in 10-h using IPTG induction, and the recombinant protein was purified by Ni-NTA agarose affinity chromatography. After that, the kinetic and thermodynamic parameters of the enzyme, optimum temperature and pH were also investigated to find more in-depth information. In the end, to represent the connections between the structures and function of this enzyme, the molecular dynamics simulations have been considered at two temperatures in which DTD had maximum and minimum activity (310 and 293 K, respectively). The results of MD simulations indicated that the interaction between NADH with phenylalanine 232 residue at 310 K is more severe than other residues. So, to investigate the interaction details of NADH/PHE 232 the DFT calculations were done.

Structure based designing of benzimidazole/benzoxazole derivatives as anti-leishmanial agents

Folates are essential biomolecules required to carry out many crucial processes in leishmania parasite. Dihydrofolate reductase-thymidylate synthase (DHFR-TS) and pteridine reductase 1 (PTR1) involved in folate biosynthesis in leishmania have been established as suitable targets for development of chemotherapy against leishmaniasis. In the present study, various computational tools such as homology modelling, pharmacophore modelling, docking, molecular dynamics and molecular mechanics have been employed to design dual DHFR-TS and PTR1 inhibitors. Two designed molecules, i.e. 2-(4-((4-nitrobenzyl)oxy)phenyl)-1H-benzo[d]imidazole and 2-(4-((2,4-dichlorobenzyl)oxy)phenyl)-1H-benzo[d]oxazolemolecules were synthesized. MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay was performed to evaluate in vitro activity of molecules against promastigote form of Leishmania donovani using Miltefosine as standard. 2-(4-((4-nitrobenzyl)oxy)phenyl)-1H-benzo[d]imidazole and 2-(4-((2,4-dichlorobenzyl)oxy)phenyl)-1H-benzo[d]oxazolemolecules were found to be moderately active with showed IC₅₀ = 68 ± 2.8 µM and 57 ± 4.2 µM, respectively.

Comparative antiradical activity and molecular Docking/Dynamics analysis of octopamine and norepinephrine: the role of OH groups

Octopamine is a neurotransmitter in invertebrates and a phenol analog of norepinephrine. The crystallographic and spectral (UV-visUV, and NMR) characteristics of octopamine were investigated experimentally and theoretically by applying appropriate level of theory, B3LYP-D3BJ/6-311++G(d,p), which reproduced well the experimental bond lengths and angles. The intramolecular interactions governing the stability of conformers were described by NBO and QTAIM analyses. The antiradical potencies of octopamine and norepinephrine towards DPPH and ABTS⁺ were examined with special emphasis on the preferred mechanism and effect of catechol moiety. Several techniques were used to distinguish Hydrogen Atom Transfer (HAT) and Proton Coupled Electron Transfer (PCET) mechanisms for reaction with DPPH. The calculated rate constants of the reactions with both radicals showed that Sequential Proton Loss Electron Transfer (SPLET) mechanism was dominant both thermodynamically and kinetically, with values of thermodynamic functions and rate constants clearly proving the importance of the second hydroxyl group in structure. The Molecular Docking and afterward Molecular Dynamics calculations of formed complexes between octopamine/norepinephrine with β1- and β2- adrenergic receptors examined in details the interactions that lead to the formation of stable complexes. The number of strong interactions of amino acids with norepinephrine was higher, but the absence of hydroxyl group in octopamine did not lead to a significant change in the type of interactions and stability. The formed complexes showed higher flexibility of amino acids, similar compactness of structure as proteins and increased interatomic distances of the backbone when compared to pure proteins.

Designing self-assembled peptide nanovaccine against Streptococcus pneumoniae: An in silico strategy

Streptococcus pneumoniae is the main cause of diseases such as meningitis, pneumoniae and sepsis, especially in children and old people. Due to costly antibiotic treatment, and increasing resistance of pneumococcus, developing high-efficient protective vaccine against this pathogen is an urgent need. Although the pneumoniae polysaccharide vaccine (PPV) and pneumonia conjugate vaccines (PCV) are the efficient pneumococcal vaccine in children and adult groups, but the serotype replacement of S. pneumoniae strains causes the reduction in efficacy of such vaccines. For overcoming the aforesaid drawbacks epitope-based vaccines are introduced as the relevant alternative. In our previous research, the epitope vaccine was designed based on immunodominant epitopes from PspA, CbpA antigens as cellular stimulants and PhtD, PiuA as humoral stimulants. Because the low immunogenicity is the main disadvantage of epitope vaccine, in the current study, we applied coiled-coil self-assembled structures for developing our vaccine. Recently, self-assembled peptide nanoparticles (SAPNs) have gained much attention in the field of vaccine development due to their multivalency, self-adjuvanticity, biocompatibility, and size similarity to pathogen. In this regard, the final designed vaccine is comprised of cytotoxic T lymphocytes (CTL) epitopes from PspA and CbpA, helper T lymphocytes (HTL) epitopes from PhtD and PiuA, the pentamer and trimmer oligomeric domains form 5-stranded and 3-stranded coiled-coils as self-assembled scaffold, Diphtheria toxoids (DTD) as a universal T-helper, which fused to each other with appropriate linkers. The four different arrangements based on the order of above-mentioned compartments were constructed, and each of them were modeled, and validated to find the 3D structure. The structural, physicochemical, and immunoinformatics analyses of final vaccine construct represented that our vaccine could stimulate potent immune response against S. pneumoniae; however, the potency of that should be approved via various in vivo and in vitro immunological tests.

•

Stimulating cellular and especially humoral immunities are essential for protection against Streptococcus Pneumoniae.

•

Immunodominant epitopes were selected from highly protective antigens of S. pneumoniae: PspA, CbpA, PiuA, PhtD.

•

In order to bypass the low immunogenicity of epitope-based peptide vaccine the self-assembled motifs, coiled-coil structure, was applied as the vaccine scaffold.

•

The structural, physicochemical, and immunoinformatics results indicate that the designed vaccine can incite strong immune response against S. pneumoniae.

Structure based virtual screening of natural products to disrupt the structural integrity of TRAF6 C-terminal domain homotrimer

Tumor necrosis factor receptor-associated factor 6 (TRAF6) is an E3 ligase which takes part in different cellular pathways. TRAF6 is seen to be highly expressed in various cancers and most importantly is known to drive cancer metastasis. This makes TRAF6 a potential therapeutic target. In our previous studies, we observed that the C-terminal domain of TRAF6 forms a mushroom shaped trimer structure. Lys340 and Glu345 were identified to be the most critical residues in the trimer interface. In this current work, we screened for more than 14000 small molecules derived from various natural sources and they were screened against TRAF6 C-terminal trimer interaction interface to prevent the formation of the interface. All the obtained molecules were tested for their drug-likeliness properties. The ligands which qualified the filter were considered for protein-ligand docking or structure based virtual screening in GOLD 5.2. Pose selection was carried out on the basis of GoldScore and ChemScore function of GOLD 5.2. Top 20 molecules binding the TRAF6 trimeric interface were tested for their ADME properties. From the top 20 molecules, top 3 ligands were chosen based on their abilities to pass the maximum numbers of ADME filters.

Insight into the binding of a synthetic nitro-flavone derivative with human poly (ADP-ribose) polymerase 1

Flavones are important bioactive compounds, many of which are effective in cancer therapy for their ability to target enzymes related to DNA repair and cell proliferation. In this report, the interaction of a synthetic nitroflavone, 2,4-nitrophenylchromen-4-one (4NCO) with human poly (ADP-ribose) polymerase 1 (hPARP1) was investigated to explore its inhibitory action. Its interaction with hPARP1 was compared with that of other inhibitors through molecular docking studies. Further insight into the 4NCO-hPARP1 interaction was obtained from competitive docking and molecular dynamic simulation studies. In silico mutagenesis studies and per-residue interaction energy calculations were carried out. Quantitative Structure Activity Relationship analysis was also performed to calculate its predictive percent inhibitory activity. Our results indicated that 4NCO exhibited competitive mode of binding to hPARP1. It formed a stable interaction with the protein thereby hindering any further molecular interaction to render it inactive with a predictive inhibition of 96%. It also had good ADMET properties and showed best Autodock binding free energy values compared to other known inhibitors. 4NCO showed good hPARP1 inhibitory properties with higher bioavailability and lower probability of getting effluxed. Development of inhibitors against hPARP1 is important for cell proliferative disorders, where 4NCO can be predicted as a potential new drug.

In silico analysis of PFN1 related to amyotrophic lateral sclerosis

Profilin 1 (PFN1) protein plays key roles in neuronal growth and differentiation, membrane trafficking, and regulation of the actin cytoskeleton. Four natural variants of PFN1 were described as related to ALS, the most common adult-onset motor neuron disorder. However, the pathological mechanism of PFN1 in ALS is not yet completely understood. The goal of this work is to thoroughly analyze the effects of the ALS-related mutations on PFN1 structure and function using computational simulations. Here, PhD-SNP, PMUT, PolyPhen-2, SIFT, SNAP, SNPS&GO, SAAP, nsSNPAnalyzer, SNPeffect4.0 and I-Mutant2.0 were used to predict the functional and stability effects of PFN1 mutations. ConSurf was used for the evolutionary conservation analysis, and GROMACS was used to perform the MD simulations. The mutations C71G, M114T, and G118V, but not E117G, were predicted as deleterious by most of the functional prediction algorithms that were used. The stability prediction indicated that the ALS-related mutations could destabilize PFN1. The ConSurf analysis indicated that the mutation C71G, M114T, E117G, and G118V occur in highly conserved positions. The MD results indicated that the studied mutations could affect the PFN1 flexibility at the actin and PLP-binding domains, and consequently, their intermolecular interactions. It may be therefore related to the functional impairment of PFN1 upon C71G, M114T, E117G and G118V mutations, and their involvement in ALS development. We also developed a database, SNPMOL (http://www.snpmol.org/), containing the results presented on this paper for biologists and clinicians to exploit PFN1 and its natural variants.

Changes of Thermostability, Organic Solvent, and pH Stability in Geobacillus zalihae HT1 and Its Mutant by Calcium Ion

Thermostable T1 lipase from Geobacillus zalihae has been crystallized using counter-diffusion method under space and Earth conditions. The comparison of the three-dimensional structures from both crystallized proteins show differences in the formation of hydrogen bond and ion interactions. Hydrogen bond and ion interaction are important in the stabilization of protein structure towards extreme temperature and organic solvents. In this study, the differences of hydrogen bond interactions at position Asp43, Thr118, Glu250, and Asn304 and ion interaction at position Glu226 was chosen to imitate space-grown crystal structure, and the impact of these combined interactions in T1 lipase-mutated structure was studied. Using space-grown T1 lipase structure as a reference, subsequent simultaneous mutation D43E, T118N, E226D, E250L, and N304E was performed on recombinant wild-type T1 lipase (wt-HT1) to generate a quintuple mutant term as 5M mutant lipase. This mutant lipase shared similar characteristics to its wild-type in terms of optimal pH and temperature. The stability of mutant 5M lipase improved significantly in acidic and alkaline pH as compared to wt-HT1. 5M lipase was highly stable in organic solvents such as dimethyl sulfoxide (DMSO), methanol, and n-hexane compared to wt-HT1. Both wild-type and mutant lipases were found highly activated in calcium as compared to other metal ions due to the presence of calcium-binding site for thermostability. The presence of calcium prolonged the half-life of mutant 5M and wt-HT1, and at the same time increased their melting temperature (T_m). The melting temperature of 5M and wt-HT1 lipases increased at 8.4 and 12.1 °C, respectively, in the presence of calcium as compared to those without. Calcium enhanced the stability of mutant 5M in 25% (v/v) DMSO, n-hexane, and n-heptane. The lipase activity of wt-HT1 also increased in 25% (v/v) ethanol, methanol, acetonitrile, n-hexane, and n-heptane in the presence of calcium. The current study showed that the accumulation of amino acid substitutions D43E, T118N, E226D, E250L, and N304E produced highly stable T1 mutant when hydrolyzing oil in selected organic solvents such as DMSO, n-hexane, and n-heptane. It is also believed that calcium ion plays important role in regulating lipase thermostability.

In Silico insights on enhancing thermostability and activity of a plant Fructosyltransferase from Pachysandra terminalis via introduction of disulfide bridges

Drawbacks of industrially-used fructosyltransferases (FTs) such as low optimum temperature and low fructooligosaccharides (FOS) yield necessitates the search for engineered FTs that are highly thermostable and active. With the availability of the first plant FT crystal structure from Pachysandra terminalis (PDB ID: 3UGH), computer-aided protein engineering of plant FT is now feasible. To obtain insights on the effect of specific mutations i.e. disulfide bridge introduction, wild-type and mutant FTs were subjected to a 15 μs Martini Coarse-grained Molecular Dynamics (CGMD) simulations at 303 K and 334 K. We report here the five mutants, M31C-Q49C, L144C-S193C, P34C-W300C, S219C-L226C and V470C-S498C with enhanced thermostabilities and/or activities relative to the wild type. Interestingly, M31C-Q49C, which is located within the catalytic-carrying blade of the catalytic domain, has an activity enhancement at both temperatures. At 334 K, three mutations, L144C-S193C, P34C-W300C and V470C-S498C, achieved thermostability relative to the wild type. Intriguingly, both activity and stability enhancement exhibited only at 334 K can be achieved provided that the mutation is located either on the catalytic-carrying residue blade of the catalytic domain or on the non-catalytic domain. Our results suggest that V470C-S498C and L144C-S193C are promising mutants and that domain-specific approach may be exploited to customize enzyme properties.

Construction of stabilized (R)-selective amine transaminase from Aspergillus terreus by consensus mutagenesis

Amine transaminases are a class of efficient and industrially-desired biocatalysts for the production of chiral amines. In this study, stabilized variants of the (R)-selective amine transaminase from Aspergillus terreus (AT-ATA) were constructed by consensus mutagenesis. Using Consensus Finder (http://cbs-kazlab.oit.umn.edu/), six positions with the most prevalent amino acid (over 60% threshold) among the homologous family members were identified. Subsequently, these six residues were individually mutated to match the consensus sequence (I77 L, Q97E, H210N, N245D, G292D, and I295 V) using site-directed mutagenesis. Compared to that of the wild-type, the thermostability of all six single variants was improved. The H210N variant displayed the largest shift in thermostability, with a 3.3-fold increase in half-life (t_1/2) at 40 °C, and a 4.6 °C increase in T₅₀¹⁰ among the single variants. In addition, the double mutant H210N/I77L displayed an even larger shift with 6.1-fold improvement of t_1/2 at 40 °C, and a 6.6 °C increase in T₅₀¹⁰. Furtherly, the H210N/I77L mutation was introduced into the previously engineered thermostable AT-ATA by the introduction of disulfide bonds, employing B-factor and folding free energy (ΔΔG^fold) calculations. Our results showed that the combined variant H210N/I77L/M150C-M280C had the largest shift in thermostability, with a 16.6-fold improvement of t_1/2 and a 11.8 °C higher T₅₀¹⁰.

Cysteine Disulfide Traps Reveal Distinct Conformational Ensembles in Dengue Virus NS2B-NS3 Protease

The dengue virus protease (NS2B-NS3pro) plays a critical role in the dengue viral life cycle, making it an attractive drug target for dengue-related pathologies, including dengue hemorrhagic fever. A number of studies indicate that NS2B-NS3pro undergoes a transition between two widely different conformational states: an "open" (inactive) conformation and a "closed" (active) conformation. For the past several years, the equilibrium between these states and the resting conformation of NS2B-NS3pro have been debated, although a strong consensus is emerging. To investigate the importance of such conformational states, we developed versions of NS2B-NS3pro that allow us to trap the enzyme in various distinct conformations. Our data from these variants suggest that the enzymatic activity appears to be dependent on the movement of NS2B and may rely on the flexibility of the protease core. Locking the enzyme into the "closed" conformation dramatically increased activity, strongly suggesting that the "closed" conformation is the active conformation. The observed resting state of the enzyme depends largely on the construct used to express the NS2B-NS3pro complex. In an "unlinked" construct, in which the NS2B and NS3 regions exist as independent, co-expressed polypeptides, the enzyme rests predominantly in a "closed", active conformation. In contrast, in a "linked" construct, in which NS2B and NS3 are attached by a nine-amino acid linker, NS2B-NS3pro adopts a more relaxed, alternative conformation. Nevertheless, even the unlinked construct samples both the "closed" and other alternative conformations. Given our findings, and the more realistic resemblance of NS2B-NS3pro to the native enzyme, these data strongly suggest that studies should focus on the "unlinked" constructs moving forward. Additionally, the results from these studies provide a more detailed understanding of the various poses of the dengue virus NS2B-NS3 protease and should help guide future drug discovery efforts aimed at this enzyme.

Silibinin, A Natural Blend In Polytherapy Formulation For Targeting Cd44v6 Expressing Colon Cancer Stem Cells

Colon cancer stem cells have been attributed to poor prognosis, therapeutic resistance and aggressive nature of the malignancy. Recent reports associated CD44v6 expression with relapse, metastasis and reduced 5-year survival of colon cancer patients, thereby making it a potential therapeutic target. Thus, in this study, comprehensive prediction and screening of CD44v6 against 1674 lead compounds was conducted. Silibinin was identified as a potential compound targeting CD44v6. Inorder to substantiate these findings, the cytotoxic effect of 5FU, Silibinin and 5FU+ Silibinin was assessed on human colon carcinoma cell line HCT116 derived CD44+ subpopulation. 5FU+ Silibinin inhibited cell proliferation of CD44+ subpopulation at lower concentration than Silibinin standalone. Further, corresponding to CD44v6 knockdown cells, 5FU+ Silibinin treatment significantly decreased CD44v6, Nanog, CTNNB1 and CDKN2A expression whereas increased E-cadherin expression in HCT116 derived CD44+ cells. Moreover, synergistic effect of these drugs suppressed sphere formation, inhibited cell migration, triggered PARP cleavage and perturbation in mitochondrial membrane potential, thereby activating intrinsic apoptotic pathways and induced autophagic cell death. Importantly, 5FU+ Silibinin could inhibit PI3K/MAPK dual activation and arrest the cell cycle at G0/G1 phase. Thus, our study suggests that inhibition of CD44v6 attenuates stemness of colon cancer stem cells and holds a prospect of potent therapeutic target.

Identification and analyses of natural compounds as potential inhibitors of TRAF6-Basigin interactions in melanoma using structure-based virtual screening and molecular dynamics simulations

The interaction of the proteins, tumor necrosis factor receptor-associated factor6 (TRAF6) and Basigin (CD147), is known to be associated with the over-expression of matrix metalloproteinases (MMPs) in melanoma cells. MMPs are known to be responsible for melanoma metastasis. Hence, the TRAF6-Basigin complex can act as a potential therapeutic target. In previous studies, amino acid residues Lys340, Lys 384, Glu417 and Glu511 of TRAF6 were identified as the most vital residues on the basis of their contributions to interaction energy, relative solvent accessibility and electrostatic interactions in the TRAF6-Basigin protein-protein interaction (PPI) scheme. In our current work, we performed structure-based virtual screenings of some natural compounds obtained from ZINC database (n = 14509) to search for molecules which can act as inhibitors against the formation of TRAF6-Basigin complex. Three potential inhibitors were identified which were observed to make intermolecular interactions with Lys384 and Glu511 of TRAF6. Molecular dynamics simulation results suggested the substantial pharmacological importance of the ligand molecules as it was observed that there was total destabilization of TRAF6-Basigin complex upon binding of the molecule ZINC02578057. From our studies, we could conclude that the ligands termed as ZINC49048033, ZINC02578057 and ZINC72320240 could have great potentials to act as inhibitors to prevent melanoma metastasis.

Conformational changes allow processing of bulky substrates by a haloalkane dehalogenase with a small and buried active site

Haloalkane dehalogenases catalyze the hydrolysis of halogen-carbon bonds in organic halogenated compounds and as such are of great utility as biocatalysts. The crystal structures of the haloalkane dehalogenase DhlA from the bacterium from Xanthobacter autotrophicus GJ10, specifically adapted for the conversion of the small 1,2-dichloroethane (DCE) molecule, display the smallest catalytic site (110 Å3) within this enzyme family. However, during a substrate-specificity screening, we noted that DhlA can catalyze the conversion of far bulkier substrates, such as the 4-(bromomethyl)-6,7-dimethoxy-coumarin (220 Å3). This large substrate cannot bind to DhlA without conformational alterations. These conformational changes have been previously inferred from kinetic analysis, but their structural basis has not been understood. Using molecular dynamic simulations, we demonstrate here the intrinsic flexibility of part of the cap domain that allows DhlA to accommodate bulky substrates. The simulations displayed two routes for transport of substrates to the active site, one of which requires the conformational change and is likely the route for bulky substrates. These results provide insights into the structure-dynamics function relationships in enzymes with deeply buried active sites. Moreover, understanding the structural basis for the molecular adaptation of DhlA to 1,2-dichloroethane introduced into the biosphere during the industrial revolution provides a valuable lesson in enzyme design by nature.

Dehalogenases: From Improved Performance to Potential Microbial Dehalogenation Applications

The variety of halogenated substances and their derivatives widely used as pesticides, herbicides and other industrial products is of great concern due to the hazardous nature of these compounds owing to their toxicity, and persistent environmental pollution. Therefore, from the viewpoint of environmental technology, the need for environmentally relevant enzymes involved in biodegradation of these pollutants has received a great boost. One result of this great deal of attention has been the identification of environmentally relevant bacteria that produce hydrolytic dehalogenases—key enzymes which are considered cost-effective and eco-friendly in the removal and detoxification of these pollutants. These group of enzymes catalyzing the cleavage of the carbon-halogen bond of organohalogen compounds have potential applications in the chemical industry and bioremediation. The dehalogenases make use of fundamentally different strategies with a common mechanism to cleave carbon-halogen bonds whereby, an active-site carboxylate group attacks the substrate C atom bound to the halogen atom to form an ester intermediate and a halide ion with subsequent hydrolysis of the intermediate. Structurally, these dehalogenases have been characterized and shown to use substitution mechanisms that proceed via a covalent aspartyl intermediate. More so, the widest dehalogenation spectrum of electron acceptors tested with bacterial strains which could dehalogenate recalcitrant organohalides has further proven the versatility of bacterial dehalogenators to be considered when determining the fate of halogenated organics at contaminated sites. In this review, the general features of most widely studied bacterial dehalogenases, their structural properties, basis of the degradation of organohalides and their derivatives and how they have been improved for various applications is discussed.

Engineering more stable proteins

The dynamic native, functional folded forms of proteins are unstable mainly because they readily unfold into flexible unstructured forms.

Temperature-dependent conformational dynamics of cytochrome c: Implications in apoptosis

Heat, electric shock, and burn injuries induce apoptosis by releasing cytochrome c (cyt-c) from mitochondria and by subsequently activating the death protease, caspases-3. During apoptosis, cyt-c undergoes changes in the secondary structure that have been suggested to increase its peroxidase activity. Information about these structural changes will provide better understanding of the apoptotic mechanism. Hence, temperature-dependent conformational dynamics of cyt-c has been investigated through molecular dynamics (MD) simulations to explain the structural changes and to correlate them with its apoptotic behavior. We observe that, at lower temperatures (223, 248, and 300K), the secondary structure of cyt-c, remains stable, while at higher temperatures (323, 373, 423, and 473K), the secondary structural regions change significantly. Further, our MD results indicate that these structural changes are mainly localized on α-helices, turns, β-sheets, and important loops that were involved in the stabilization of the heme conformation. This conformational transition between specific regions of secondary structure of cyt-c directly affects the electron tunneling properties of the proteins as observed experimentally. We quantify and compare these changes and explain that the temperature plays a vital role in assuring the structural stability of cyt-c and thus its functions. Our findings from this MD study reproduce experimental results at high temperatures and provide evidence for the alteration of the heme through the disruption of the H-bonding interactions between specific regions of cyt-c, thereby enhancing its peroxidase activity which plays a crucial role in the apoptotic process.

Comparison of Five Protein Engineering Strategies for Stabilizing an α/β-Hydrolase

A review of the previous stabilization of α/β-hydrolase fold enzymes revealed many different strategies, but no comparison of strategies on the same enzyme. For this reason, we compared five strategies to identify stabilizing mutations in a model α/β-hydrolase fold enzyme, salicylic acid binding protein 2, to reversible denaturation by urea and to irreversible denaturation by heat. The five strategies included one location agnostic approach (random mutagenesis using error-prone polymerase chain reaction), two structure-based approaches [computational design (Rosetta, FoldX) and mutation of flexible regions], and two sequence-based approaches (addition of proline at locations where a more stable homologue has proline and mutation to consensus). All strategies identified stabilizing mutations, but the best balance of success rate, degree of stabilization, and ease of implementation was mutation to consensus. A web-based automated program that predicts substitutions needed to mutate to consensus is available at http://kazlab.umn.edu .

Molecular Dynamic Simulation of Space and Earth-Grown Crystal Structures of Thermostable T1 Lipase Geobacillus zalihae Revealed a Better Structure

Less sedimentation and convection in a microgravity environment has become a well-suited condition for growing high quality protein crystals. Thermostable T1 lipase derived from bacterium Geobacilluszalihae has been crystallized using the counter diffusion method under space and earth conditions. Preliminary study using YASARA molecular modeling structure program for both structures showed differences in number of hydrogen bond, ionic interaction, and conformation. The space-grown crystal structure contains more hydrogen bonds as compared with the earth-grown crystal structure. A molecular dynamics simulation study was used to provide insight on the fluctuations and conformational changes of both T1 lipase structures. The analysis of root mean square deviation (RMSD), radius of gyration, and root mean square fluctuation (RMSF) showed that space-grown structure is more stable than the earth-grown structure. Space-structure also showed more hydrogen bonds and ion interactions compared to the earth-grown structure. Further analysis also revealed that the space-grown structure has long-lived interactions, hence it is considered as the more stable structure. This study provides the conformational dynamics of T1 lipase crystal structure grown in space and earth condition.

Destabilizing polymorphism in cervid prion protein hydrophobic core determines prion conformation and conversion efficiency

Prion diseases are infectious neurodegenerative disorders of humans and animals caused by misfolded forms of the cellular prion protein PrP^C. Prions cause disease by converting PrP^C into aggregation-prone PrP^Sc. Chronic wasting disease (CWD) is the most contagious prion disease with substantial lateral transmission, affecting free-ranging and farmed cervids. Although the PrP primary structure is highly conserved among cervids, the disease phenotype can be modulated by species-specific polymorphisms in the prion protein gene. How the resulting amino-acid substitutions impact PrP^C and PrP^Sc structure and propagation is poorly understood. We investigated the effects of the cervid 116A>G substitution, located in the most conserved PrP domain, on PrP^C structure and conversion and on 116AG-prion conformation and infectivity. Molecular dynamics simulations revealed structural de-stabilization of 116G-PrP, which enhanced its in vitro conversion efficiency when used as recombinant PrP substrate in real-time quaking-induced conversion (RT-QuIC). We demonstrate that 116AG-prions are conformationally less stable, show lower activity as a seed in RT-QuIC and exhibit reduced infectivity in vitro and in vivo. Infectivity of 116AG-prions was significantly enhanced upon secondary passage in mice, yet conformational features were retained. These findings indicate that structurally de-stabilized PrP^C is readily convertible by cervid prions of different genetic background and results in a prion conformation adaptable to cervid wild-type PrP. Conformation is an important criterion when assessing transmission barrier, and conformational variants can target a different host range. Therefore, a thorough analysis of CWD isolates and re-assessment of species-barriers is important in order to fully exclude a zoonotic potential of CWD.

Chronic wasting disease (CWD) is a prion disease which affects wild and captive cervids. Prion diseases are infectious neurodegenerative disorders, and the causative agent consists of abnormally folded prion protein termed PrP^Sc. Prions replicate without genetic information, and their three-dimensional structure is thought to encode heritable information necessary to propagate using the cellular prion protein PrP^C as a substrate for conversion. In this study, we use in vitro and in vivo techniques to analyze the effect of a polymorphism at codon 116 (A>G) of the white-tailed deer prion protein on CWD prion conformation, propagation and pathogenesis. We observed differences in conformation, infectivity and seeding activity in vitro between CWD prions isolated from white-tailed deer encoding wild-type (116AA) PrP^C or 116AG-PrP^C. In mouse bioassays conformational differences are retained, however, 116AG CWD prions resulted in significantly shortened incubation times upon passages. Molecular dynamics simulations suggest that the structure of 116G-PrP^C is more flexible, which is supported by an improved convertibility in an in vitro conversion assay. Altogether these data indicate the importance of a variation in the most conserved PrP domain, and highlight the relationship between PrP^C structural flexibility, prion conformation and conversion, and pathogenesis of prion disease in vivo.

Structural and dynamical insight into thermally induced functional inactivation of firefly luciferase

Luciferase is the key component of light production in bioluminescence process. Extensive and advantageous application of this enzyme in biotechnology is restricted due to its low thermal stability. Here we report the effect of heating up above Tm on the structure and dynamical properties of luciferase enzyme compared to temperature at 298 K. In this way we demonstrate that the number of hydrogen bonds between N- and C-domain is increased for the free enzyme at 325 K. Increased inter domain hydrogen bonds by three at 325 K suggests that inter domain contact is strengthened. The appearance of simultaneous strong salt bridge and hydrogen bond between K529 and D422 and increased existence probability between R533 and E389 could mechanistically explain stronger contact between N- and C-domain. Mutagenesis studies demonstrated the importance of K529 and D422 experimentally. Also the significant reduction in SASA for experimentally important residues K529, D422 and T343 which are involved in active site region was observed. Principle component analysis (PCA) in our study shows that the dynamical behavior of the enzyme is changed upon heating up which mainly originated from the change of motion modes and associated extent of those motions with respect to 298 K. These findings could explain why heating up of the enzyme or thermal fluctuation of protein conformation reduces luciferase activity in course of time as a possible mechanism of thermal functional inactivation. According to these results we proposed two strategies to improve thermal stability of functional luciferase.

Ancestral Haloalkane Dehalogenases Show Robustness and Unique Substrate Specificity

Ancestral sequence reconstruction (ASR) represents a powerful approach for empirical testing structure-function relationships of diverse proteins. We employed ASR to predict sequences of five ancestral haloalkane dehalogenases (HLDs) from the HLD-II subfamily. Genes encoding the inferred ancestral sequences were synthesized and expressed in Escherichia coli, and the resurrected ancestral enzymes (AncHLD1-5) were experimentally characterized. Strikingly, the ancestral HLDs exhibited significantly enhanced thermodynamic stability compared to extant enzymes (ΔT_m up to 24 °C), as well as higher specific activities with preference for short multi-substituted halogenated substrates. Moreover, multivariate statistical analysis revealed a shift in the substrate specificity profiles of AncHLD1 and AncHLD2. This is extremely difficult to achieve by rational protein engineering. The study highlights that ASR is an efficient approach for the development of novel biocatalysts and robust templates for directed evolution.

[Prediction of protein thermostability from their primary structure: the current state and development factors]

The construction of proteins and peptides with desired properties, including resistance to high temperatures, as well as optimization of their amino acid composition, is an important and complex task, which attracts much attention in various branches of the basic sciences, and also in biomedicine and biotechnology. This raises the question: what method is more relevant for the at the pilot stage of research in order to estimate the influence of the planned amino acid substitutions on the thermostability of the resultant protein construct? In this brief review we have classified existing basic practical and theoretical approaches used in studies and predicting the thermal stability of native and recombinant polypeptides. Particular attention has been paid to the predictive potential of statistical methods for studying the thermodynamic parameters of the primary protein structure and prospects of their use.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

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.

A general and efficient strategy for generating the stable enzymes

The local flexibility of an enzyme's active center plays pivotal roles in catalysis, however, little is known about how the flexibility of these flexible residues affects stability. In this study, we proposed an active center stabilization (ACS) strategy to improve the kinetic thermostability of Candida rugosa lipase1. Based on the B-factor ranking at the region ~10 Å within the catalytic Ser209, 18 residues were selected for site-saturation mutagenesis. Based on three-tier high-throughput screening and ordered recombination mutagenesis, the mutant VarB3 (F344I/F434Y/F133Y/F121Y) was shown to be the most stable, with a 40-fold longer in half-life at 60 °C and a 12.7 °C higher Tm value than that of the wild type, without a decrease in catalytic activity. Further analysis of enzymes with different structural complexities revealed that focusing mutations on the flexible residues within around 10 Å of the catalytic residue might increase the success rate for enzyme stabilization. In summary, this study identifies a panel of flexible residues within the active center that affect enzyme stability. This finding not only provides clues regarding the molecular evolution of enzyme stability but also indicates that ACS is a general and efficient strategy for exploring the functional robustness of enzymes for industrial applications.

Molecular Characterization and Enhancement of Anticancer Activity of Caffeic Acid Phenethyl Ester by γ Cyclodextrin

Caffeic Acid Phenethyl Ester (CAPE) is a key component in New Zealand propolis, known for a variety of health promoting and therapeutic potentials. We investigated the molecular mechanism of anticancer and anti-metastasis activities of CAPE. cDNA array performed on the control and CAPE-treated breast cancer cells revealed activation of DNA damage signaling involving upregulation of GADD45α and p53 tumor suppressor proteins. Molecular docking analysis revealed that CAPE is capable of disrupting mortalin-p53 complexes. We provide experimental evidence and demonstrate that CAPE induced disruption of mortalin-p53 complexes led to nuclear translocation and activation of p53 resulting in growth arrest in cancer cells. Furthermore, CAPE-treated cells exhibited downregulation of mortalin and several other key regulators of cell migration accountable for its anti-metastasis activity. Of note, we found that whereas CAPE was unstable in the culture medium (as it gets degraded into caffeic acid by secreted esterases), its complex with gamma cyclodextrin (γCD) showed high efficacy in anti-tumor and anti-metastasis assays in vitro and in vivo (when administered through either intraperitoneal or oral route). The data proposes that CAPE-γCD complex is a potent anti-cancer and anti-metastasis reagent.

Structure-guided discovery of thiazolidine-2,4-dione derivatives as a novel class of Leishmania major pteridine reductase 1 inhibitors

Leishmania major, as other protozoan parasites, plague human kind since pre-historic times but it remains a worldwide ailment for which the therapeutic arsenal remains scarce. Although L. major is pteridine- and purine-auxotroph, well-established folate biosynthesis inhibitors, such as methotrexate, have poor effect over the parasite survival. The lack of efficiency is related to an alternative biochemical pathway in which pteridine reductase 1 (PTR1) plays a major role. For this reason, this enzyme has been considered a promising target for anti-leishmanial drug development and several inhibitors that share the substrate scaffold have been reported. In order to design a novel class of PTR1 inhibitors, we employed the thiazolidinone ring as a bioisosteric replacement for pteridine/purine ring. Among seven novel thiazolidine-2,4-dione derivatives reported herein, 2d was identified as the most promising lead by thermal shift assays (ΔTm = 11 °C, p = 0,01). Kinetic assays reveal that 2d has IC50 = 44.67 ± 1.74 μM and shows a noncompetitive behavior. This information guided docking studies and molecular dynamics simulations (50 000 ps) that supports 2d putative binding profile (H-bonding to Ser-111 and Leu-66) and shall be useful to design more potent inhibitors.

Molecular dynamics-based identification of novel natural mortalin-p53 abrogators as anticancer agents

INTRODUCTION

Cancer is one of the leading causes of mortality worldwide that requires attention in terms of extensive study and research. Eradication of mortalin-p53 interaction that leads to the inhibition of transcriptional activation or blocking of p53 from functioning as a suppressor and induction of nuclear translocation of p53 can prove to be one of the useful approaches for cancer management.

RESULTS

In this study, we used structure-based approach to target the p53-binding domain of mortalin in order to prevent mortalin-p53 complex formation. We screened compounds from ZINC database against the modeled mortalin protein using Glide virtual screening. The top two compounds, DTOM (ZINC 28639308) and TTOM (ZINC 38143676) with Glide score of -12.27 and -12.16, respectively, were identified with the potential to abrogate mortalin-p53 interaction. Finally, molecular dynamics simulations were used to analyze the dynamic stability of the ligand-bound complex and it was observed that residues Tyr196, Asn198, Val264 and Thr267 were involved in intermolecular interactions in both the simulated ligand-bound complexes, and thus, these residues may have a paramount role in stabilizing the binding of the ligands with the protein.

CONCLUSION

These detailed insights can further facilitate the development of potent inhibitors against mortalin-p53 complex.

Computational investigation of curcumin, a natural polyphenol that inhibits the destabilization and the aggregation of human SOD1 mutant (Ala4Val)

Curcumin inhibits the aberrant aggregation in mutant SOD1 protein, thereby decreasing the propensity of β-sheets and the toxicity level.