Protein Flexibility and Enzymatic Catalysis

ATLAS: protein flexibility description from atomistic molecular dynamics simulations

Dynamical behaviour is one of the most crucial protein characteristics. Despite the advances in the field of protein structure resolution and prediction, analysis and prediction of protein dynamic properties remains a major challenge, mostly due to the low accessibility of data and its diversity and heterogeneity. To address this issue, we present ATLAS, a database of standardised all-atom molecular dynamics simulations, accompanied by their analysis in the form of interactive diagrams and trajectory visualisation. ATLAS offers a large-scale view and valuable insights on protein dynamics for a large and representative set of proteins, by combining data obtained through molecular dynamics simulations with information extracted from experimental structures. Users can easily analyse dynamic properties of functional protein regions, such as domain limits (hinge positions) and residues involved in interaction with other biological molecules. Additionally, the database enables exploration of proteins with uncommon dynamic properties conditioned by their environment such as chameleon subsequences and Dual Personality Fragments. The ATLAS database is freely available at https://www.dsimb.inserm.fr/ATLAS.

A computational analysis to evaluate deleterious SNPs of GSK3β, a multifunctional and regulatory protein, for metabolism, wound healing, and migratory processes

This study focused on GSK-3β, a critical serine/threonine kinase with diverse cellular functions. However, there is limited understanding of the impact of non-synonymous single nucleotide polymorphisms (nsSNPs) on its structure and function. Through an exhaustive in-silico investigation 12 harmful nsSNPs were predicted from a pool of 172 acquired from the NCBI dbSNP database using 12 established tools that detects deleterious SNPs. Consistently, these nsSNPs were discovered in locations with high levels of conservation. Notably, the three harmful nsSNPs F67C, A83T, and T138I were situated in the active/binding site of GSK-3β, which may affect the protein's capacity to bind to substrates and other proteins. Molecular dynamics simulations revealed that the F67C and T138I mutants had stable structures, indicating rigidness, whereas the A83T mutant was unstable. Analysis of secondary structures revealed different modifications in all mutant forms, which may affect the stability, functioning, and interactions of the protein. These mutations appear to alter the structural dynamics of GSK-3β, which may have functional ramifications, such as the formation of novel secondary structures and variations in coil-to-helix transitions. In conclusion, this study illuminates the possible structural and functional ramifications of these GSK-3 nsSNPs, revealing how protein compactness, stiffness, and interactions may affect biological activities.

Bioinformatics study of phytase from Aspergillus niger for use as feed additive in livestock feed

BACKGROUND

Phytase supplementation in rations can reduce their phytic acid composition in order to enhance their nutritional value. Aspergillus niger is a fungus that can encode phytase. This study aims to determine the characteristics of its DNA sequences and amino acid composition that encode the phytase enzyme, as well as to determine the primer designs.

METHOD

This study used gene sequence data and protein-encoding phytase from Aspergillus niger that was collected manually from NCBI and PDB. The data was analyzed using SPDBV and then be aligned using the ClustalW Multiple Alignment features. The phylogenetic tree was built by Mega11 software. Primers were designed from selected candidate sequences that were analyzed. The designed primers were then simulated for PCR using FastPCR and SnapGene software.

RESULTS

There are 18 Aspergillus niger phytases in NCBI which is 14.87% of the total Aspergillus. There are 14 Aspergillus niger phytases that have identity above 95%. Aspergillus niger 110. M94550.1 is the closest strain to the PDB template. Candidate sources of phytase genes are Aspergillus niger 110.M94550.1, 48.2.BCMY01000003.1, and 92.JQ654450.1. The primer design has 2 possibilities of self-annealing and high melting temperature on the reverse primer. PCR simulation shows that the primer design can attach completely but still has the possibility of mispriming.

CONCLUSION

This study suggests promising results for the future development of phytase enzyme production from Aspergillus niger as a feed additive using genetic engineering to enhance the quality of livestock feed in Indonesia.

MDexciteR: Enhanced Sampling Molecular Dynamics by Excited Normal Modes or Principal Components Obtained from Experiments

Molecular dynamics with excited normal modes (MDeNM) is an enhanced sampling method for exploring conformational changes in proteins with minimal biases. The excitation corresponds to injecting kinetic energy along normal modes describing intrinsic collective motions. Herein, we developed a new automated open-source implementation, MDexciteR (https://github.com/mcosta27/MDexciteR), enabling the integration of MDeNM with two commonly used simulation programs with GPU support. Second, we generalized the method to include the excitation of principal components calculated from experimental ensembles. Finally, we evaluated whether the use of coarse-grained normal modes calculated with elastic network representations preserved the performance and accuracy of the method. The advantages and limitations of these new approaches are discussed based on results obtained for three different protein test cases: two globular and a protein/membrane system.

Structural Insight of KSIII (β-Ketoacyl-ACP Synthase)-like Acyltransferase ChlB3 in the Biosynthesis of Chlorothricin

Chlorothricin (CHL) belongs to a spirotetronate antibiotic family produced by Streptomyces antibioticus that inhibits pyruvate carboxylase and malate dehydrogenase. For the biosynthesis of CHL, ChlB3 plays a crucial role by introducing the 6-methylsalicylic acid (6MSA) moiety to ChlB2, an acyl carrier protein (ACP). However, the structural insight and catalytic mechanism of ChlB3 was unclear. In the current study, the crystal structure of ChlB3 was solved at 3.1 Å-resolution and a catalytic mechanism was proposed on the basis of conserved residues of structurally related enzymes. ChlB3 is a dimer having the same active sites as CerJ (a structural homologous enzyme) and uses a KSIII-like fold to work as an acyltransferase. The relaxed substrate specificity of ChlB3 was defined by its catalytic efficiencies (kcat/Km) for non-ACP tethered synthetic substrates such as 6MSA-SNAC, acetyl-SNAC, and cyclohexonyl-SNAC. ChlB3 successfully detached the 6MSA moiety from 6MSA-SNAC substrate and this hydrolytic activity demonstrated that ChlB3 has the potential to catalyze non-ACP tethered substrates. Structural comparison indicated that ChlB3 belongs to FabH family and showed 0.6–2.5 Å root mean square deviation (RMSD) with structural homologous enzymes. Molecular docking and dynamics simulations were implemented to understand substrate active site and structural behavior such as the open and closed conformation of the ChlB3 protein. The resultant catalytic and substrate recognition mechanism suggested that ChlB3 has the potential to use non-native substrates and minimize the labor of expressing ACP protein. This versatile acyltransferase activity may pave the way for manufacturing CHL variants and may help to hydrolyze several thioester-based compounds.

Synthesis, Crystallographic, Quantum Chemical, Antitumor, and Molecular Docking/Dynamic Studies of 4-Hydroxycoumarin-Neurotransmitter Derivatives

In this contribution, four new compounds synthesized from 4-hydroxycoumarin and tyramine/octopamine/norepinephrine/3-methoxytyramine are characterized spectroscopically (IR and NMR), chromatographically (UHPLC-DAD), and structurally at the B3LYP/6-311++G*(d,p) level of theory. The crystal structure of the 4-hydroxycoumarin-octopamine derivative was solved and used as a starting geometry for structural optimization. Along with the previously obtained 4-hydroxycoumarin-dopamine derivative, the intramolecular interactions governing the stability of these compounds were quantified by NBO and QTAIM analyses. Condensed Fukui functions and the HOMO-LUMO gap were calculated and correlated with the number and position of OH groups in the structures. In vitro cytotoxicity experiments were performed to elucidate the possible antitumor activity of the tested substances. For this purpose, four cell lines were selected, namely human colon cancer (HCT-116), human adenocarcinoma (HeLa), human breast cancer (MDA-MB-231), and healthy human lung fibroblast (MRC-5) lines. A significant selectivity towards colorectal carcinoma cells was observed. Molecular docking and molecular dynamics studies with carbonic anhydrase, a prognostic factor in several cancers, complemented the experimental results. The calculated MD binding energies coincided well with the experimental activity, and indicated 4-hydroxycoumarin-dopamine and 4-hydroxycoumarin-3-methoxytyramine as the most active compounds. The ecotoxicology assessment proved that the obtained compounds have a low impact on the daphnia, fish, and green algae population.

Comparative MD Study of Inhibitory Activity of Opaganib and Adamantane-Isothiourea Derivatives toward COVID-19 Main Protease Mpro

In this study, the inhibitory potency of four adamantly- isothiourea derivatives (compounds 1 [4-bromobenzyl (Z)-N'-(adamantan-1-yl)-4-phenylpiperazine-1-carbothioimidate], 2 [3,5-bis(trifluoromethyl)benzyl (Z)-N'-(adamantan-1-yl)-4-phenylpiperazine-1-carbothioimidate], 3 [4-bromobenzyl (Z)-N-(adamantan-1-yl)morpholine-4-carbothioimidate] and 4 [3,5-bis(trifluoromethyl)benzyl (Z)-N-(adamantan-1-yl)morpholine-4-carbothioimidate]) was evaluated against SARS-CoV-2 targeted proteins. The investigated compounds 1-4 possess a similar structure to opaganib, which is used in studies like a potential drug for COVID-19 treatment. Since examined adamantly-isothiourea derivatives (1-4) shown broad-spectrum of antibacterial activity and significant in vitro cytotoxic effects against five human tumor cell lines and shown similarity in structure with opaganib, it was of interest to study their inhibitory potency toward some SARS-CoV-2 proteins such as SARS-CoV-2 main protease Mpro and mutation of SARS-CoV-2 Spike (S) Protein D614G. The inhibitory potency of studied compounds is examined using molecular docking and molecular dynamic simulations. The results of molecular docking simulations indicate compound 1 as the most prominent candidate of inhibition of SARS-CoV-2 main protease Mpro (▵Gbind=11.24 kcal/mol), while almost the same inhibition potency of all studied compounds is exhibited toward D614G. Regarding the results obtained by molecular dynamic simulations, compounds 1 and 4 possess similar inhibitory potency toward SARS-CoV-2 main protease Mpro as opaganib (▵Gbind ≈ 40 kcal/mol).

Reductionism in the study of enzyme adaptation

One of the principal goals of comparative biology is the elucidation of mechanisms by which organisms adapt to different environments. The study of enzyme structure, function, and stability has contributed significantly to this effort, by revealing adaptation at a molecular level. Comparative biochemistry, including enzymology, necessarily pursues a reductionist approach in describing the function and structure of biomolecules, allowing more straightforward study of molecular systems by removing much of the complexity of their biological milieu. Although this reductionism has allowed a remarkable series of discoveries linking chemical processes to metabolism and to whole-organism function in the context of the environment, it also has the potential to mislead when careful consideration is not made of the simplifying assumptions inherent to such research. In this review, a brief history of the growth of enzymology, its reliance on a reductionist philosophy, and its contributions to our understanding of biological systems is given. Examples then are provided of research techniques, based on a reductionist approach, that have advanced our knowledge about enzyme adaptation to environmental stresses, including stability assays, enzyme kinetics, and the impact of solute composition on enzyme function. In each case, the benefits of the reductionist nature of the approach is emphasized, notable advances are described, but potential drawbacks due to inherent oversimplification of the study system are also identified.

Flexibility of enzymatic transitions as a hallmark of optimized enzyme steady-state kinetics and thermodynamics

We investigate the relations between the enzyme kinetic flexibility, the rate of entropy production, and the Shannon information entropy in a steady-state enzyme reaction. All these quantities are maximized with respect to enzyme rate constants. We show that the steady-state, which is characterized by the most flexible enzymatic transitions between the enzyme conformational states, coincides with the global maxima of the Shannon information entropy and the rate of entropy production. This steady-state of an enzyme is referred to as globally optimal. This theoretical approach is then used for the analysis of the kinetic and the thermodynamic performance of the enzyme triose-phosphate isomerase. The analysis reveals that there exist well-defined maxima of the kinetic flexibility, the rate of entropy production, and the Shannon information entropy with respect to any arbitrarily chosen rate constant of the enzyme and that these maxima, calculated from the measured kinetic rate constants for the triose-phosphate isomerase are lower, however of the same order of magnitude, as the maxima of the globally optimal state of the enzyme. This suggests that the triose-phosphate isomerase could be a well, but not fully evolved enzyme, as it was previously claimed. Herein presented theoretical investigations also provide clear evidence that the flexibility of enzymatic transitions between the enzyme conformational states is a requirement for the maximal Shannon information entropy and the maximal rate of entropy production.

Insight to Functional Conformation and Noncovalent Interactions of Protein-Protein Assembly Using MALDI Mass Spectrometry

Noncovalent interactions are the keys to the structural organization of biomolecule e.g., proteins, glycans, lipids in the process of molecular recognition processes e.g., enzyme-substrate, antigen-antibody. Protein interactions lead to conformational changes, which dictate the functionality of that protein-protein complex. Besides biophysics techniques, noncovalent interaction and conformational dynamics, can be studied via mass spectrometry (MS), which represents a powerful tool, due to its low sample consumption, high sensitivity, and label-free sample. In this review, the focus will be placed on Matrix-Assisted Laser Desorption Ionization Mass Spectrometry (MALDI-MS) and its role in the analysis of protein-protein noncovalent assemblies exploring the relationship within noncovalent interaction, conformation, and biological function.

Membrane Lipids Assist Catalysis by CTP: Phosphocholine Cytidylyltransferase

While most of the articles in this issue review the workings of integral membrane enzymes, in this review, we describe the catalytic mechanism of an enzyme that contains a soluble catalytic domain but appears to catalyze its reaction on the membrane surface, anchored and assisted by a separate regulatory amphipathic helical domain and inter-domain linker. Membrane partitioning of CTP: phosphocholine cytidylyltransferase (CCT), a key regulatory enzyme of phosphatidylcholine metabolism, is regulated chiefly by changes in membrane phospholipid composition, and boosts the enzyme's catalytic efficiency >200-fold. Catalytic enhancement by membrane binding involves the displacement of an auto-inhibitory helix from the active site entrance-way and promotion of a new conformational ensemble for the inter-domain, allosteric linker that has an active role in the catalytic cycle. We describe the evidence for close contact between membrane lipid, a compact allosteric linker, and the CCT active site, and discuss potential ways that this interaction enhances catalysis.

Conformational gating in ammonia lyases

BACKGROUND

Ammonia lyases are enzymes of industrial and biomedical interest. Knowledge of structure-dynamics-function relationship in ammonia lyases is instrumental for exploiting the potential of these enzymes in industrial or biomedical applications.

METHODS

We investigated the conformational changes in the proximity of the catalytic pocket of a 3-methylaspartate ammonia lyase (MAL) as a model system. At this scope, we used microsecond all-atom molecular dynamics simulations, analyzed with dimensionality reduction techniques, as well as in terms of contact networks and correlated motions.

RESULTS

We identify two regulatory elements in the MAL structure, i.e., the β5-α2 loop and the helix-hairpin-loop subdomain. These regulatory elements undergo conformational changes switching from 'occluded' to 'open' states. The rearrangements are coupled to changes in the accessibility of the active site. The β5-α2 loop and the helix-hairpin-loop subdomain modulate the formation of tunnels from the protein surface to the catalytic site, making the active site more accessible to the substrate when they are in an open state.

CONCLUSIONS

Our work pinpoints a sequential mechanism, in which the helix-hairpin-loop subdomain of MAL needs to break a subset of intramolecular interactions first to favor the displacement of the β5-α2 loop. The coupled conformational changes of these two elements contribute to modulate the accessibility of the catalytic site.

GENERAL SIGNIFICANCE

Similar molecular mechanisms can have broad relevance in other ammonia lyases with similar regulatory loops. Our results also imply that it is important to account for protein dynamics in the design of variants of ammonia lyases for industrial and biomedical applications.

Single-Nucleotide Polymorphisms (SNP) Mining and Their Effect on the Tridimensional Protein Structure Prediction in a Set of Immunity-Related Expressed Sequence Tags (EST) in Atlantic Salmon (Salmo salar)

Single-nucleotide polymorphisms (SNPs) are single genetic code variations considered one of the most common forms of nucleotide modifications. Such SNPs can be located in genes associated to immune response and, therefore, they may have direct implications over the phenotype of susceptibility to infections affecting the productive sector. In this study, a set of immune-related genes (cc motif chemokine 19 precursor [ccl19], integrin β2 (itβ2, also named cd18), glutathione transferase omega-1 [gsto-1], heat shock 70 KDa protein [hsp70], major histocompatibility complex class I [mhc-I]) were analyzed to identify SNPs by data mining. These genes were chosen based on their previously reported expression on infectious pancreatic necrosis virus (IPNV)-infected Atlantic salmon phenotype. The available EST sequences for these genes were obtained from the Unigene database. Twenty-eight SNPs were found in the genes evaluated and identified most of them as transition base changes. The effect of the SNPs located on the 5’-untranslated region (UTR) or 3’-UTR upon transcription factor binding sites and alternative splicing regulatory motifs was assessed and ranked with a low-medium predicted FASTSNP score risk. Synonymous SNPs were found on itβ2 (c.2275G > A), gsto-1 (c.558G > A), and hsp70 (c.1950C > T) with low FASTSNP predicted score risk. The difference in the relative synonymous codon usage (RSCU) value between the variant codons and the wild-type codon (ΔRSCU) showed one negative (hsp70 c.1950C > T) and two positive ΔRSCU values (itβ2 c.2275G > A; gsto-1 c.558G > A), suggesting that these synonymous SNPs (sSNPs) may be associated to differences in the local rate of elongation. Nonsynonymous SNPs (nsSNPs) in the gsto-1 translatable gene region were ranked, using SIFT and POLYPHEN web-tools, with the second highest (c.205A > G; c484T > C) and the highest (c.499T > C; c.769A > C) predicted score risk possible. Using homology modeling to predict the effect of these nonsynonymous SNPs, the most relevant nucleotide changes for gsto-1 were observed for the nsSNPs c.205A > G, c484T > C, and c.769A > C. Molecular dynamics was assessed to analyze if these GSTO-1 variants have significant differences in their conformational dynamics, suggesting these SNPs could have allosteric effects modulating its catalysis. Altogether, these results suggest that candidate SNPs identified may play a crucial potential role in the immune response of Atlantic salmon.

Evolutionary redesign of the lysosomal enzyme arylsulfatase A increases efficacy of enzyme replacement therapy for metachromatic leukodystrophy

Protein engineering is a means to optimize protein therapeutics developed for the treatment of so far incurable diseases including cancers and genetic disorders. Here we report on an engineering approach in which we successfully increased the catalytic rate constant of an enzyme that is presently evaluated in enzyme replacement therapies (ERT) of a lysosomal storage disease (LSD). Although ERT is a treatment option for many LSDs, outcomes are lagging far behind expectations for most of them. This has been ascribed to insufficient enzyme activities accumulating in tissues difficult to target such as brain and peripheral nerves. We show for human arylsulfatase A (hARSA) that the activity of a therapeutic enzyme can be substantially increased by reversing activity-diminishing and by inserting activity-promoting amino acid substitutions that had occurred in the evolution of hominids and non-human mammals, respectively. The potential of this approach, here designated as evolutionary redesign, was highlighted by the observation that murinization of only 1 or 3 amino acid positions increased the hARSA activity 3- and 5-fold, with little impact on stability, respectively. The two kinetically optimized hARSA variants showed no immunogenic potential in ERT of a humanized ARSA knockout mouse model of metachromatic leukodystrophy (MLD) and reduced lysosomal storage of kidney, peripheral and central nervous system up to 3-fold more efficiently than wild-type hARSA. Due to their safety profile and higher therapeutic potential the engineered hARSA variants might represent major advances for future enzyme-based therapies of MLD and stimulate analogous approaches for other enzyme therapeutics.

Anti-HIV drug repurposing against SARS-CoV-2

A novel severe acute respiratory syndrome human coronavirus (SARS HCoV) was identified from respiratory illness patients (named SARS-CoV-2 by ICTV) in December 2019 and has recently emerged as a serious threat to world public health. However, no approved drugs have been found to effectively inhibit the virus. Since it has been reported that HIV protease inhibitors can be used as anti-SARS drugs by targeting SARS-CoV-1 3CLpro, we chose six approved anti-HIV drugs and investigated their binding interactions with 3CLpro to evaluate their potential to become clinical drugs for the new coronavirus pneumonia (COVID-19) caused by SARS-CoV-2 infection. The molecular docking results indicate that the 3CLpro of SARS-CoV-2 has a higher binding affinity for all the studied inhibitors than does SARS-CoV-1. Two docking complexes (indinavir and darunavir) with high docking scores were further subjected to MM-PBSA binding free energy calculations to detail the molecular interactions between these two protease inhibitors and SARS HCoV 3CLpro. Our results show that, among the inhibitors tested, darunavir has the highest binding affinity with SARS-CoV-2 and SARS-CoV-1 3CLpro, indicating that it may have the potential to be used as an anti-COVID-19 clinical drug. The mechanism behind the increased binding affinity of HIV protease inhibitors toward SARS-CoV-2 3CLpro (as compared to SARS-CoV-1) was investigated by MD simulations. Our study provides insight into the possible role of structural flexibility during interactions between SARS HCoV 3CLpro and inhibitors and sheds light on structure-based design of anti-COVID-19 drugs targeting SARS-CoV-2 3CLpro.

A novel severe acute respiratory syndrome human coronavirus (SARS HCoV) was identified from respiratory illness patients (named SARS-CoV-2 by ICTV) in December 2019 and has recently emerged as a serious threat to world public health.

Enhancing the thermostability and activity of uronate dehydrogenase from Agrobacterium tumefaciens LBA4404 by semi-rational engineering

Background

Glucaric acid, one of the aldaric acids, has been declared a “top value-added chemical from biomass”, and is especially important in the food and pharmaceutical industries. Biocatalytic production of glucaric acid from glucuronic acid is more environmentally friendly, efficient and economical than chemical synthesis. Uronate dehydrogenases (UDHs) are the key enzymes for the preparation of glucaric acid in this way, but the poor thermostability and low activity of UDH limit its industrial application. Therefore, improving the thermostability and activity of UDH, for example by semi-rational design, is a major research goal.

Results

In the present work, three UDHs were obtained from different Agrobacterium tumefaciens strains. The three UDHs have an approximate molecular weight of 32 kDa and all contain typically conserved UDH motifs. All three UDHs showed optimal activity within a pH range of 6.0–8.5 and at a temperature of 30 °C, but the UDH from A. tumefaciens (At) LBA4404 had a better catalytic efficiency than the other two UDHs (800 vs 600 and 530 s−1 mM−1). To further boost the catalytic performance of the UDH from AtLBA4404, site-directed mutagenesis based on semi-rational design was carried out. An A39P/H99Y/H234K triple mutant showed a 400-fold improvement in half-life at 59 °C, a 5 °C improvement in $$ {\text{T}}_{ 5 0}^{ 1 0} $$ T 50 10 value and a 2.5-fold improvement in specific activity at 30 °C compared to wild-type UDH.

Conclusions

In this study, we successfully obtained a triple mutant (A39P/H99Y/H234K) with simultaneously enhanced activity and thermostability, which provides a novel alternative for the industrial production of glucaric acid from glucuronic acid.

A super-Gaussian Poisson-Boltzmann model for electrostatic free energy calculation: smooth dielectric distribution for protein cavities and in both water and vacuum states

Calculations of electrostatic potential and solvation free energy of macromolecules are essential for understanding the mechanism of many biological processes. In the classical implicit solvent Poisson-Boltzmann (PB) model, the macromolecule and water are modeled as two-dielectric media with a sharp border. However, the dielectric property of interior cavities and ion-channels is difficult to model realistically in a two-dielectric setting. In fact, the detection of water molecules in a protein cavity remains to be an experimental challenge. This introduces an uncertainty, which affects the subsequent solvation free energy calculation. In order to compensate this uncertainty, a novel super-Gaussian dielectric PB model is introduced in this work, which devices an inhomogeneous dielectric distribution to represent the compactness of atoms and characterizes empty cavities via a gap dielectric value. Moreover, the minimal molecular surface level set function is adopted so that the dielectric profile remains to be smooth when the protein is transferred from water phase to vacuum. An important feature of this new model is that as the order of super-Gaussian function approaches the infinity, the dielectric distribution reduces to a piecewise constant of the two-dielectric model. Mathematically, an effective dielectric constant analysis is introduced in this work to benchmark the dielectric model and select optimal parameter values. Computationally, a pseudo-time alternative direction implicit (ADI) algorithm is utilized for solving the super-Gaussian PB equation, which is found to be unconditionally stable in a smooth dielectric setting. Solvation free energy calculation of a Kirkwood sphere and various proteins is carried out to validate the super-Gaussian model and ADI algorithm. One macromolecule with both water filled and empty cavities is employed to demonstrate how the cavity uncertainty in protein structure can be bypassed through dielectric modeling in biomolecular electrostatic analysis.

Enhancement of the thermostability of Aspergillus niger α-l-rhamnosidase based on PoPMuSiC algorithm

α-l-Rhamnosidase is a biotechnologically important enzyme in food industry and in the preparation of drugs and drug precursors. To expand the functionality of our previously cloned α-l-rhamnosidase from Aspergillus niger JMU-TS528, 14 mutants were constructed based on the changes of the folding free energy (ΔΔG), predicted by the PoPMuSiC algorithm. Among them, six single-site mutants displayed higher thermal stability than wild type (WT). The combinational mutant K573V-E631F displayed even higher thermostability than six single-site mutants. The spectra analyses displayed that the WT and K573V-E631F had almost similar secondary and tertiary structure profiles. The simulated protein structure-based interaction analysis and molecular dynamics calculation were further implemented to assess the conformational preferences of the K573V-E631F. The improved thermostability of mutant K573V-E631F may be attributed to the introduction of new cation-π and hydrophobic interactions, and the overall improvement of the enzyme conformation. PRACTICAL APPLICATIONS: The stability of enzymes, particularly with regards to thermal stability remains a critical issue in industrial biotechnology and industrial processing generally tends to higher ambient temperature to inhibit microbial growth. Most of the α-l-rhamnosidases are usually active at temperature from 30 to 60°C, which are apt to denature at temperatures over 60°C. To expand the functionality of our previously cloned α-l-rhamnosidase from Aspergillus niger JMU-TS528, we used protein engineering methods to increase the thermal stability of the α-l-rhamnosidase. Practically, conducting reactions at high temperatures enhances the solubility of substrates and products, increases the reaction rate thus reducing the reaction time, and inhibits the growth of contaminating microorganisms. Thus, the improvement on the thermostability of α-l-rhamnosidase on the one hand can increase enzyme efficacy; on the other hand, the high ambient temperature would enhance the solubility of natural substrates of α-l-rhamnosidase, such as naringin, rutin, and hesperidin, which are poorly dissolved in water at room temperature. Protein thermal resistance is an important issue beyond its obvious industrial importance. The current study also helps in the structure-function relationship study of α-l-rhamnosidase.

Cold survival strategies for bacteria, recent advancement and potential industrial applications

Microorganisms have evolved themselves to thrive under various extreme environmental conditions such as extremely high or low temperature, alkalinity, and salinity. These microorganisms adapted several metabolic processes to survive and reproduce efficiently under such extreme environments. As the major proportion of earth is covered with the cold environment and is exploited by human beings, these sites are not pristine anymore. Human interventions are a great reason for disturbing the natural biogeochemical cycles in these regions. The survival strategies of these organisms have shown great potential for helping us to restore these pristine sites and the use of isolated cold-adapted enzymes from these organisms has also revolutionized various industrial products. This review gives you the insight of psychrophilic enzyme adaptations and their industrial applications.

Flexible Peptide Linkers Enhance the Antimicrobial Activity of Surface-Immobilized Bacteriolytic Enzymes

Chemical linkers are frequently used in enzyme immobilization to improve enzyme flexibility and activity, whereas peptide linkers, although ubiquitous in protein engineering, are much less explored in enzyme immobilization. Here, we report peptide-linker-assisted noncovalent immobilization of the bacteriolytic enzyme lysostaphin (Lst) to generate anti- Staphylococcus aureus surfaces. Lst was immobilized through affinity tags onto a silica surface (glass slides) and nickel nitrilotriacetic acid (NiNTA) agarose beads via silica-binding peptides (SiBPs) or a hexahistidine tag (His-tag) fused at the C-terminus of Lst, respectively. By inserting specific peptide linkers upstream of the SiBP or His-tag, the immobilized enzymes killed >99.5% of S. aureus ATCC 6538 cells (10⁸ CFU/mL) within 3 h in buffer and could be reused multiple times without significant loss of activity. In contrast, immobilized Lst without a peptide linker was less active/stable. Molecular modeling of Lst-linker-affinity tag constructs illustrated that the presence of the peptide linkers enhanced the molecular flexibility of the proximal Lst binding domain, which interacts with the bacterial substrate, and such increased flexibility correlated with increased antimicrobial activity. We further show that Lst immobilized onto NiNTA beads retained the ability to kill ∼99% of a 10⁸ CFU/mL microbial challenge even in the presence of 1% of a commercial anionic surfactant, C12-14 alcohol EO 3:1 sodium sulfate, when the Lst construct contained a decapeptide linker containing glycine, serine, and alanine residues. This linker-assisted immobilization strategy could be extended to an unrelated lytic enzyme, the endolysin PlyPH, to target Bacillus anthracis Sterne cells either in buffer or in the presence of anionic surfactants. Our approach, therefore, provides a facile route to the use of antimicrobial enzymes on surfaces.

Solving the riddle: Unraveling the mechanisms of blocking the binding of leukotoxin by therapeutic antagonists in periodontal diseases

Aggregatibacter actinomycetemcomitans is a Gram-negative bacteria that has gained wide recognition for its causative role in the development of various immune diseases, which includes localized aggressive periodontitis. Its ability to evade host defense mechanisms is mediated by the secretion of leukotoxin (LtxA), which induces death of white blood cells (leukocytes) by specific binding to their surface-expressed leukocyte function-associated receptor (LFA-1) in its active state. Therapeutic compounds that interfere with this pathogenic process and abrogate A. actinomycetemcomitans virulence have been reported in literature. These include doxycycline, and more recently phytochemical compounds such as hamamelitanin, resveratrol, naringin, and quercetin. However, the question remains how do they work? Therefore, with the aid of computational tools, we explore the molecular mechanisms by which they possibly elicit their therapeutic functions. Molecular mechanics Poisson/Boltzmann surface area analyses revealed that these compounds bind favorably to active LFA-1 with high affinity and considerable stability, indicative of their ability to occupy the LtxA binding site (LBS) and prevent LtxA binding. The conformational transition of open LFA-1 to its closed state further describe the mechanistic activity of these compounds. In addition to notable reductions in structural mobility and flexibility, the burial of surface-exposed interactive side chains at the LBS was observed, an occurrence that could alter the complementary binding of LtxA. It is also important to mention that these occurrences were induced more prominently by the phytochemicals. We believe that these findings will enhance the scope of drug design and discovery for potent LtxA antagonists with improved activities and therapeutic efficacies in the treatment of virulent A. actinomycetemcomitans diseases.

Differential Enzyme Flexibility Probed Using Solid-State Nanopores

Enzymes and motor proteins are dynamic macromolecules that coexist in a number of conformations of similar energies. Protein function is usually accompanied by a change in structure and flexibility, often induced upon binding to ligands. However, while measuring protein flexibility changes between active and resting states is of therapeutic significance, it remains a challenge. Recently, our group has demonstrated that breadth of signal amplitudes in measured electrical signatures as an ensemble of individual protein molecules is driven through solid-state nanopores and correlates with protein conformational dynamics. Here, we extend our study to resolve subtle flexibility variation in dihydrofolate reductase mutants from unlabeled single molecules in solution. We first demonstrate using a canonical protein system, adenylate kinase, that both size and flexibility changes can be observed upon binding to a substrate that locks the protein in a closed conformation. Next, we investigate the influence of voltage bias and pore geometry on the measured electrical pulse statistics during protein transport. Finally, using the optimal experimental conditions, we systematically study a series of wild-type and mutant dihydrofolate reductase proteins, finding a good correlation between nanopore-measured protein conformational dynamics and equilibrium bulk fluorescence probe measurements. Our results unequivocally demonstrate that nanopore-based measurements reliably probe conformational diversity in native protein ensembles.

A Memetic Algorithm for 3-D Protein Structure Prediction Problem

Memetic Algorithms are population-based metaheuristics intrinsically concerned with exploiting all available knowledge about the problem under study. The incorporation of problem domain knowledge is not an optional mechanism, but a fundamental feature of the Memetic Algorithms. In this paper, we present a Memetic Algorithm to tackle the three-dimensional protein structure prediction problem. The method uses a structured population and incorporates a Simulated Annealing algorithm as a local search strategy, as well as ad-hoc crossover and mutation operators to deal with the problem. It takes advantage of structural knowledge stored in the Protein Data Bank, by using an Angle Probability List that helps to reduce the search space and to guide the search strategy. The proposed algorithm was tested on nineteen protein sequences of amino acid residues, and the results show the ability of the algorithm to find native-like protein structures. Experimental results have revealed that the proposed algorithm can find good solutions regarding root-mean-square deviation and global distance total score test in comparison with the experimental protein structures. We also show that our results are comparable in terms of folding organization with state-of-the-art prediction methods, corroborating the effectiveness of our proposal.

MD Simulations Reveal Complex Water Paths in Squalene-Hopene Cyclase: Tunnel-Obstructing Mutations Increase the Flow of Water in the Active Site

Squalene-hopene cyclase catalyzes the cyclization of squalene to hopanoids. A previous study has identified a network of tunnels in the protein, where water molecules have been indicated to move. Blocking these tunnels by site-directed mutagenesis was found to change the activation entropy of the catalytic reaction from positive to negative with a concomitant lowering of the activation enthalpy. As a consequence, some variants are faster and others are slower than the wild type (wt) in vitro under optimal reaction conditions for the wt. In this study, molecular dynamics (MD) simulations have been performed for the wt and the variants to investigate how the mutations affect the protein structure and the water flow in the enzyme, hypothetically influencing the activation parameters. Interestingly, the tunnel-obstructing variants are associated with an increased flow of water in the active site, particularly close to the catalytic residue Asp376. MD simulations with the substrate present in the active site indicate that the distance for the rate-determining proton transfer between Asp376 and the substrate is longer in the tunnel-obstructing protein variants than in the wt. On the basis of the previous experimental results and the current MD results, we propose that the tunnel-obstructing variants, at least partly, could operate by a different catalytic mechanism, where the proton transfer may have contributions from a Grotthuss-like mechanism.

Obtaining a high activity subtilisin preparation by controlled thermal stress in n-octane

The use of enzymes in organic solvents has considerably widened their repertoire of applications. Such low water containing media also offer the possibility of carrying out enzymatic reactions at higher temperatures and enhancing reaction yields. The utility of such preparations is limited by the damage caused to the protein structure during freeze-drying. This work investigates the result of exposing the proteolytic enzyme subtilisin to high temperature in low water containing n-octane on its activity in aqueous and non-aqueous media. Exposing subtilisin at 90 °C for 5 h led to 18-fold improvement in its transesterification activity even at the normal assay temperature (37 °C) when compared with the untreated enzyme. The use of n-octane as the reaction medium was important as it helped to retain the three-dimensional architecture of the enzyme and should be considered while designing strategies for obtaining high activity preparations of other enzymes. Structural analysis using differential scanning fluorimetry showed that the enzyme lost its structure after heating in aqueous medium but retained it when heated in organic solvent. The simplicity and general applicability of the strategy should make it useful for obtaining highly active preparations of other enzymes as well.

In silico screening, molecular docking, and molecular dynamics studies of SNP-derived human P5CR mutants

Pyrroline-5-carboxylate reductase (P5CR) encoded by PYCR1 gene is a housekeeping enzyme that catalyzes the reduction of P5C to proline using NAD(P)H as the cofactor. In this study, we used in silico approaches to examine the role of nonsynonymous single-nucleotide polymorphisms in the PYCR1 gene and their putative functions in the pathogenesis of Cutis Laxa. Among the 348 identified SNPs, 15 were predicted to be potentially damaging by both SIFT and PolyPhen tools; of them two SNP-derived mutations, R119G and G206W, have been previously reported to correlate with Cutis Laxa. These two mutations were therefore selected to be mapped to the wild-type (WT) P5CR structure for further structural and functional analyses. The results of comparative computational analyses using I-Mutant and Autodock reveal reductions in both stability and cofactor binding affinity of these two mutants. Comparative molecular dynamics (MD) simulations were performed to evaluate the changes in dynamic properties of P5CR upon mutations. The results reveal that the two mutations enhance the rigidity of P5CR structure, especially that of cofactor binding site, which could result in decreased kinetics of cofactor entrance and egress. Comparison between the structural properties of the WT and mutants during MD simulations shows that the enhanced rigidity of mutants results most likely from the increased number of inter-atomic interactions and the decreased number of dynamic hydrogen bonds. Our study provides novel insight into the deleterious effects of the R119G and G206W mutations on P5CR, and sheds light on the mechanisms by which these mutations mediate Cutis Laxa.

Action of a GH115 α-glucuronidase from Amphibacillus xylanus at alkaline condition promotes release of 4-O-methylglucopyranosyluronic acid from glucuronoxylan and arabinoglucuronoxylan

Glucuronic acid and/or 4-O-methyl-glucuronic acid (GlcA/MeGlcA) are substituents of the main xylans present in hardwoods, conifers, and many cereal grains. α-Glucuronidases from glycoside hydrolase family GH115 can target GlcA/MeGlcA from both internally and terminally substituted regions of xylans. The current study describes the first GH115 α-glucuronidase, AxyAgu115A, from the alkaliphilic organism Amphilbacillus xylanus. AxyAgu115A was active in a wide pH range, and demonstrated better performance in alkaline condition compared to other characterized GH115 α-glucuronidases, which generally show optimal activity in acidic conditions. Specifically, its relative activity between pH 5.0 and pH 8.5 was above 80%, and was 35% of maximum at pH 10.5; although the enzyme lost 30% and 80% relative residual activity after 24-h pre-incubation at pH 9 and pH 10, respectively. AxyAgu115A was also similarly active towards glucuronoxylan as well as comparatively complex xylans such as spruce arabinoglucurunoxylan. Accommodation of complex xylans was supported by docking analyses that predicted accessibility of AxyAgu115A to branched xylo-oligosaccharides. MeGlcA release by AxyAgu115A from each xylan sample was increased by up to 30% by performing the reaction at pH 11.0 rather than pH 4.0, revealing applied benefits of AxyAgu115A for xylan recovery and processing.

Molecular motions and free-energy landscape of serine proteinase K in relation to its cold-adaptation: a comparative molecular dynamics simulation study and the underlying mechanisms

The physicochemical bases for enzyme cold-adaptation remain elusive.

Probing the role of cation-π interaction in the thermotolerance and catalytic performance of endo-polygalacturonases

Understanding the dynamics of the key pectinase, polygalacturonase, and improving its thermotolerance and catalytic efficiency are of importance for the cost-competitive bioconversion of pectic materials. By combining structure analysis and molecular dynamics (MD) simulations, eight mutagenesis sites having the potential to form cation-π interactions were identified in the widely used fungal endo-polygalacturonase PG63. In comparison to the wild-type, three single mutants H58Y, T71Y and T304Y showed improved thermostability (the apparent Tms increased by 0.6-3.9 °C) and catalytic efficiency (by up to 32-fold). Chromatogram analysis of the hydrolysis products indicated that a larger amount of shorter sugars were released from the polygalacturonic acid by these three mutants than by the wild-type. MD analysis of the enzyme-substrate complexes illustrated that the mutants with introduced cation-π interaction have modified conformations of catalytic crevice, which provide an enviable environment for the catalytic process. Moreover, the lower plasticity of T3 loop 2 at the edge of the subsite tunnel appears to recruit the reducing ends of oligogalacturonide into the active site tunnel and initiates new hydrolysis reactions. This study demonstrates the importance of cation-π interaction in protein conformation and provides a realistic strategy to enhance the thermotolerance and catalytic performance of endo-polygalacturonases.

Effects of the linker region on the structure and function of modular GH5 cellulases

The association of glycosyl hydrolases with catalytically inactive modules is a successful evolutionary strategy that is commonly used by biomass-degrading microorganisms to digest plant cell walls. The presence of accessory domains in these enzymes is associated with properties such as higher catalytic efficiency, extension of the catalytic interface and targeting of the enzyme to the proper substrate. However, the importance of the linker region in the synergistic action of the catalytic and accessory domains remains poorly understood. Thus, this study examined how the inter-domain region affects the structure and function of modular GH5 endoglucanases, by using cellulase 5A from Bacillus subtilis (BsCel5A) as a model. BsCel5A variants featuring linkers with different stiffnesses or sizes were designed and extensively characterized, revealing that changes in flexibility or rigidity in this region differentially affect kinetic behavior. Regarding the linker length, we found that precise inter-domain spacing is required to enable efficient hydrolysis because excessively long or short linkers were equally detrimental to catalysis. Together, these findings identify molecular and structural features that may contribute to the rational design of chimeric and multimodular glycosyl hydrolases.

Variability of the Cyclin-Dependent Kinase 2 Flexibility Without Significant Change in the Initial Conformation of the Protein or Its Environment; a Computational Study

BACKGROUND

Protein flexibility, which has been referred as a dynamic behavior has various roles in proteins' functions. Furthermore, for some developed tools in bioinformatics, such as protein-protein docking software, considering the protein flexibility, causes a higher degree of accuracy. Through undertaking the present work, we have accomplished the quantification plus analysis of the variations in the human Cyclin Dependent Kinase 2 (hCDK2) protein flexibility without affecting a significant change in its initial environment or the protein per se.

OBJECTIVES

The main goal of the present research was to calculate variations in the flexibility for each residue of the hCDK2, analysis of their flexibility variations through clustering, and to investigate the functional aspects of the residues with high flexibility variations.

MATERIALS AND METHODS

Using Gromacs package (version 4.5.4), three independent molecular dynamics (MD) simulations of the hCDK2 protein (PDB ID: 1HCL) was accomplished with no significant changes in their initial environments, structures, or conformations, followed by Root Mean Square Fluctuations (RMSF) calculation of these MD trajectories. The amount of variations in these three curves of RMSF was calculated using two formulas.

RESULTS

More than 50% of the variation in the flexibility (the distance between the maximum and the minimum amount of the RMSF) was found at the region of Val-154. As well, there are other major flexibility fluctuations in other residues. These residues were mostly positioned in the vicinity of the functional residues. The subsequent works were done, as followed by clustering all hCDK2 residues into four groups considering the amount of their variability with respect to flexibility and their position in the RMSF curves.

CONCLUSIONS

This work has introduced a new class of flexibility aspect of the proteins' residues. It could also help designing and engineering proteins, with introducing a new dynamic aspect of hCDK2, and accordingly, for the other similar globular proteins. In addition, it could provide a better computational calculation of the protein flexibility, which is, especially important in the comparative studies of the proteins' flexibility.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

Saccharification of ozonated sugarcane bagasse using enzymes from Myceliophthora thermophila JCP 1-4 for sugars release and ethanol production

The saccharification of ozonated sugarcane bagasse (SCB) by enzymes from Myceliophthora thermophila JCP 1-4 was studied. Fungal enzymes provided slightly higher sugar release than commercial enzymes, working at 50°C. Sugar release increased with temperature increase. Kinetic studies showed remarkable glucose release (4.99 g/L, 3%w/w dry matter) at 60°C, 8 h of hydrolysis, using an enzyme load of 10 FPU (filter paper unit). FPase and β-glucosidase activities increased during saccharification (284% and 270%, respectively). No further significant improvement on glucose release was observed increasing the enzyme load above 7.5 FPU per g of cellulose. Higher dry matter contents increased sugars release, but not yields. The fermentation of hydrolysates by Saccharomyces cerevisiae provided glucose-to-ethanol conversions around to 63%.

Structural features determining thermal adaptation of esterases

The adaptation of microorganisms to extreme living temperatures requires the evolution of enzymes with a high catalytic efficiency under these conditions. Such extremophilic enzymes represent valuable tools to study the relationship between protein stability, dynamics and function. Nevertheless, the multiple effects of temperature on the structure and function of enzymes are still poorly understood at the molecular level. Our analysis of four homologous esterases isolated from bacteria living at temperatures ranging from 10°C to 70°C suggested an adaptation route for the modulation of protein thermal properties through the optimization of local flexibility at the protein surface. While the biochemical properties of the recombinant esterases are conserved, their thermal properties have evolved to resemble those of the respective bacterial habitats. Molecular dynamics simulations at temperatures around the optimal temperatures for enzyme catalysis revealed temperature-dependent flexibility of four surface-exposed loops. While the flexibility of some loops increased with raising the temperature and decreased with lowering the temperature, as expected for those loops contributing to the protein stability, other loops showed an increment of flexibility upon lowering and raising the temperature. Preserved flexibility in these regions seems to be important for proper enzyme function. The structural differences of these four loops, distant from the active site, are substantially larger than for the overall protein structure, indicating that amino acid exchanges within these loops occurred more frequently thereby allowing the bacteria to tune atomic interactions for different temperature requirements without interfering with the overall enzyme function.

Molecular dynamics simulations: advances and applications

Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.

Proteomics in evolutionary ecology

Evolutionary ecologists are traditionally gene-focused, as genes propagate phenotypic traits across generations and mutations and recombination in the DNA generate genetic diversity required for evolutionary processes. As a consequence, the inheritance of changed DNA provides a molecular explanation for the functional changes associated with natural selection. A direct focus on proteins on the other hand, the actual molecular agents responsible for the expression of a phenotypic trait, receives far less interest from ecologists and evolutionary biologists. This is partially due to the central dogma of molecular biology that appears to define proteins as the 'dead-end of molecular information flow' as well as technical limitations in identifying and studying proteins and their diversity in the field and in many of the more exotic genera often favored in ecological studies. Here we provide an overview of a newly forming field of research that we refer to as 'Evolutionary Proteomics'. We point out that the origins of cellular function are related to the properties of polypeptide and RNA and their interactions with the environment, rather than DNA descent, and that the critical role of horizontal gene transfer in evolution is more about coopting new proteins to impact cellular processes than it is about modifying gene function. Furthermore, post-transcriptional and post-translational processes generate a remarkable diversity of mature proteins from a single gene, and the properties of these mature proteins can also influence inheritance through genetic and perhaps epigenetic mechanisms. The influence of post-transcriptional diversification on evolutionary processes could provide a novel mechanistic underpinning for elements of rapid, directed evolutionary changes and adaptations as observed for a variety of evolutionary processes. Modern state-of the art technologies based on mass spectrometry are now available to identify and quantify peptides, proteins, protein modifications and protein interactions of interest with high accuracy and assess protein diversity and function. Therefore, proteomic technologies can be viewed as providing evolutionary biologist with exciting novel opportunities to understand very early events in functional variation of cellular molecular machinery that are acting as part of evolutionary processes.

Discovery of New Classes of Compounds that Reactivate Acetylcholinesterase Inhibited by Organophosphates

Acetylcholinesterase (AChE) that has been covalently inhibited by organophosphate compounds (OPCs), such as nerve agents and pesticides, has traditionally been reactivated by using nucleophilic oximes. There is, however, a clearly recognized need for new classes of compounds with the ability to reactivate inhibited AChE with improved in vivo efficacy. Here we describe our discovery of new functional groups--Mannich phenols and general bases--that are capable of reactivating OPC--inhibited AChE more efficiently than standard oximes and we describe the cooperative mechanism by which these functionalities are delivered to the active site. These discoveries, supported by preliminary in vivo results and crystallographic data, significantly broaden the available approaches for reactivation of AChE.

Conformational variation of an extreme drug resistant mutant of HIV protease

Molecular mechanisms leading to high level drug resistance have been analyzed for the clinical variant of HIV-1 protease bearing 20 mutations (PR20); which has several orders of magnitude worse affinity for tested drugs. Two crystal structures of ligand-free PR20 with the D25N mutation of the catalytic aspartate (PR20D25N) revealed three dimers with different flap conformations. The diverse conformations of PR20D25N included a dimer with one flap in a unique "tucked" conformation; directed into the active site. Analysis of molecular dynamics (MD) simulations of the ligand-free PR20 and wild-type enzymes showed that the mutations in PR20 alter the correlated interactions between two monomers in the dimer. The two flaps tend to fluctuate more independently in PR20 than in the wild type enzyme. Combining the results of structural analysis by X-ray crystallography and MD simulations; unusual flap conformations and weakly correlated inter-subunit motions may contribute to the high level resistance of PR20.

Substrate Specificity, Inhibitor Selectivity and Structure-Function Relationships of Aldo-Keto Reductase 1B15: A Novel Human Retinaldehyde Reductase

Human aldo-keto reductase 1B15 (AKR1B15) is a newly discovered enzyme which shares 92% amino acid sequence identity with AKR1B10. While AKR1B10 is a well characterized enzyme with high retinaldehyde reductase activity, involved in the development of several cancer types, the enzymatic activity and physiological role of AKR1B15 are still poorly known. Here, the purified recombinant enzyme has been subjected to substrate specificity characterization, kinetic analysis and inhibitor screening, combined with structural modeling. AKR1B15 is active towards a variety of carbonyl substrates, including retinoids, with lower kcat and Km values than AKR1B10. In contrast to AKR1B10, which strongly prefers all-trans-retinaldehyde, AKR1B15 exhibits superior catalytic efficiency with 9-cis-retinaldehyde, the best substrate found for this enzyme. With ketone and dicarbonyl substrates, AKR1B15 also shows higher catalytic activity than AKR1B10. Several typical AKR inhibitors do not significantly affect AKR1B15 activity. Amino acid substitutions clustered in loops A and C result in a smaller, more hydrophobic and more rigid active site in AKR1B15 compared with the AKR1B10 pocket, consistent with distinct substrate specificity and narrower inhibitor selectivity for AKR1B15.

Homologous ligands accommodated by discrete conformations of a buried cavity

Conformational change in protein-ligand complexes is widely modeled, but the protein accommodation expected on binding a congeneric series of ligands has received less attention. Given their use in medicinal chemistry, there are surprisingly few substantial series of congeneric ligand complexes in the Protein Data Bank (PDB). Here we determine the structures of eight alkyl benzenes, in single-methylene increases from benzene to n-hexylbenzene, bound to an enclosed cavity in T4 lysozyme. The volume of the apo cavity suffices to accommodate benzene but, even with toluene, larger cavity conformations become observable in the electron density, and over the series two other major conformations are observed. These involve discrete changes in main-chain conformation, expanding the site; few continuous changes in the site are observed. In most structures, two discrete protein conformations are observed simultaneously, and energetic considerations suggest that these conformations are low in energy relative to the ground state. An analysis of 121 lysozyme cavity structures in the PDB finds that these three conformations dominate the previously determined structures, largely modeled in a single conformation. An investigation of the few congeneric series in the PDB suggests that discrete changes are common adaptations to a series of growing ligands. The discrete, but relatively few, conformational states observed here, and their energetic accessibility, may have implications for anticipating protein conformational change in ligand design.