Search

The development of models of macromolecular electrostatics capable of delivering improved fidelity to quantum mechanical calculations is an active field of research in computational chemistry. Most molecular force field development takes place in the context of models with full Cartesian coordinate degrees of freedom. Nevertheless, a number of macromolecular modeling programs use a reduced set of conformational variables limited to rotatable bonds. Efficient algorithms for minimizing the energies of macromolecular systems with torsional degrees of freedom have been developed with the assumption that all atom-atom interaction potentials are isotropic. We describe novel modifications to address the anisotropy of higher order multipole terms while retaining the efficiency of these approaches. In addition, we present a treatment for obtaining derivatives of atom-centered tensors with respect to torsional degrees of freedom. We apply these results to enable minimization of the Amoeba multipole electrostatics potential in a system with torsional degrees of freedom, and validate the correctness of the gradients by comparison to finite difference approximations. In the interest of enabling a complete model of electrostatics with implicit treatment of solvent-mediated effects, we also derive expressions for the derivative of solvent accessible surface area with respect to torsional degrees of freedom.

Background: Sequence-specific binding by transcription factors (TFs) plays a significant role in the selection and regulation of target genes. At the protein:DNA interface, amino acid side-chains construct a diverse physicochemical network of specific and non-specific interactions, and seemingly subtle changes in amino acid identity at certain positions may dramatically impact TF:DNA binding. Variation of these specificity-determining residues (SDRs) is a major mechanism of functional divergence between TFs with strong structural or sequence homology., Methods: In this study, we employed a combination of high-throughput specificity profiling by SELEX and Spec-seq, structural modeling, and evolutionary analysis to probe the binding preferences of winged helix-turn-helix TFs belonging to the OmpR sub-family in Escherichia coli ., Results: We found that E. coli OmpR paralogs recognize tandem, variably spaced repeats composed of "GT-A" or "GCT"-containing half-sites. Some divergent sequence preferences observed within the "GT-A" mode correlate with amino acid similarity; conversely, "GCT"-based motifs were observed for a subset of paralogs with low sequence homology. Direct specificity profiling of a subset of OmpR homologues (CpxR, RstA, and OmpR) as well as predicted "SDR-swap" variants revealed that individual SDRs may impact sequence preferences locally through direct contact with DNA bases or distally via the DNA backbone., Conclusions: Overall, our work provides evidence for a common structural code for sequence-specific wHTH:DNA interactions, and demonstrates that surprisingly modest residue changes can enable recognition of highly divergent sequence motifs. Further examination of SDR predictions will likely reveal additional mechanisms controlling the evolutionary divergence of this important class of transcriptional regulators., Competing Interests: COMPLIANCE WITH ETHICS GUIDELINES The authors Adam P. Joyce and James J. Havranek declare that they have no conflict of interests. This article does not contain any studies with human or animal subjects performed by any of the authors.

Cooperative binding of transcription factors is known to be important in the regulation of gene expression programs conferring cellular identities. However, current methods to measure cooperativity parameters have been laborious and therefore limited to studying only a few sequence variants at a time. We developed Coop-seq (cooperativity by sequencing) that is capable of efficiently and accurately determining the cooperativity parameters for hundreds of different DNA sequences in a single experiment. We apply Coop-seq to 12 dimer pairs from the Sox and POU families of transcription factors using 324 unique sequences with changed half-site orientation, altered spacing and discrete randomization within the binding elements. The study reveals specific dimerization profiles of different Sox factors with Oct4. By contrast, Oct4 and the three neural class III POU factors Brn2, Brn4 and Oct6 assemble with Sox2 in a surprisingly indistinguishable manner. Two novel half-site configurations can support functional Sox/Oct dimerization in addition to known composite motifs. Moreover, Coop-seq uncovers a nucleotide switch within the POU half-site when spacing is altered, which is mirrored in genomic loci bound by Sox2/Oct4 complexes., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Cre recombinase catalyzes the cleavage and religation of DNA at loxP sites. The enzyme is a homotetramer in its functional state, and the symmetry of the protein complex enforces a pseudo-palindromic symmetry upon the loxP sequence. The Cre-lox system is a powerful tool for many researchers. However, broader application of the system is limited by the fixed sequence preferences of Cre, which are determined by both the direct DNA contacts and the homotetrameric arrangement of the Cre monomers. As a first step toward achieving recombination at arbitrary asymmetric target sites, we have broken the symmetry of the Cre tetramer assembly. Using a combination of computational and rational protein design, we have engineered an alternative interface between Cre monomers that is functional yet incompatible with the wild-type interface. Wild-type and engineered interface halves can be mixed to create two distinct Cre mutants, neither of which are functional in isolation, but which can form an active heterotetramer when combined. When these distinct mutants possess different DNA specificities, control over complex assembly directly discourages recombination at unwanted half-site combinations, enhancing the specificity of asymmetric site recombination. The engineered Cre mutants exhibit this assembly pattern in a variety of contexts, including mammalian cells., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Leber congenital amaurosis 9 (LCA9) is an autosomal recessive retinal degeneration condition caused by mutations in the NAD(+) biosynthetic enzyme NMNAT1. This condition leads to early blindness but no other consistent deficits have been reported in patients with NMNAT1 mutations despite its central role in metabolism and ubiquitous expression. To study how these mutations affect NMNAT1 function and ultimately lead to the retinal degeneration phenotype, we performed detailed analysis of LCA-associated NMNAT1 mutants, including the expression, nuclear localization, enzymatic activity, secondary structure, oligomerization, and promotion of axonal and cellular integrity in response to injury. In many assays, most mutants produced results similar to wild type NMNAT1. Indeed, NAD(+) synthetic activity is unlikely to be a primary mechanism underlying retinal degeneration as most LCA-associated NMNAT1 mutants had normal enzymatic activity. In contrast, the secondary structure of many NMNAT1 mutants was relatively less stable as they lost enzymatic activity after heat shock, whereas wild type NMNAT1 retains significant activity after this stress. These results suggest that LCA-associated NMNAT1 mutants are more vulnerable to stressful conditions that lead to protein unfolding, a potential contributor to the retinal degeneration observed in this syndrome., (© 2015 by The American Society for Biochemistry and Molecular Biology, Inc.)

Protein:DNA interactions are essential to a range of processes that maintain and express the information encoded in the genome. Structural modeling is an approach that aims to understand these interactions at the physicochemical level. It has been proposed that structural modeling can lead to deeper understanding of the mechanisms of protein:DNA interactions, and that progress in this field can not only help to rationalize the observed specificities of DNA-binding proteins but also to allow researchers to engineer novel DNA site specificities. In this review we discuss recent developments in the structural description of protein:DNA interactions and specificity, as well as the challenges facing the field in the future., (© The Author 2014. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Computational protein design relies on several approximations, including the use of fixed backbones and rotamers, to reduce protein design to a computationally tractable problem. However, allowing backbone and off-rotamer flexibility leads to more accurate designs and greater conformational diversity. Exhaustive sampling of this additional conformational space is challenging, and often impossible. Here, we report a computational method that utilizes a preselected library of native interactions to direct backbone flexibility to accommodate placement of these functional contacts. Using these native interaction modules, termed motifs, improves the likelihood that the interaction can be realized, provided that suitable backbone perturbations can be identified. Furthermore, it allows a directed search of the conformational space, reducing the sampling needed to find low energy conformations. We implemented the motif-based design algorithm in Rosetta, and tested the efficacy of this method by redesigning the substrate specificity of methionine aminopeptidase. In summary, native enzymes have evolved to catalyze a wide range of chemical reactions with extraordinary specificity. Computational enzyme design seeks to generate novel chemical activities by altering the target substrates of these existing enzymes. We have implemented a novel approach to redesign the specificity of an enzyme and demonstrated its effectiveness on a model system., (© 2014 The Authors Protein Science published by Wiley Periodicals, Inc. on behalf of The Protein Society.)

The Rosetta molecular modeling software package provides experimentally tested and rapidly evolving tools for the 3D structure prediction and high-resolution design of proteins, nucleic acids, and a growing number of non-natural polymers. Despite its free availability to academic users and improving documentation, use of Rosetta has largely remained confined to developers and their immediate collaborators due to the code's difficulty of use, the requirement for large computational resources, and the unavailability of servers for most of the Rosetta applications. Here, we present a unified web framework for Rosetta applications called ROSIE (Rosetta Online Server that Includes Everyone). ROSIE provides (a) a common user interface for Rosetta protocols, (b) a stable application programming interface for developers to add additional protocols, (c) a flexible back-end to allow leveraging of computer cluster resources shared by RosettaCommons member institutions, and (d) centralized administration by the RosettaCommons to ensure continuous maintenance. This paper describes the ROSIE server infrastructure, a step-by-step 'serverification' protocol for use by Rosetta developers, and the deployment of the first nine ROSIE applications by six separate developer teams: Docking, RNA de novo, ERRASER, Antibody, Sequence Tolerance, Supercharge, Beta peptide design, NCBB design, and VIP redesign. As illustrated by the number and diversity of these applications, ROSIE offers a general and speedy paradigm for serverification of Rosetta applications that incurs negligible cost to developers and lowers barriers to Rosetta use for the broader biological community. ROSIE is available at http://rosie.rosettacommons.org.

Accurate energy functions are critical to macromolecular modeling and design. We describe new tools for identifying inaccuracies in energy functions and guiding their improvement, and illustrate the application of these tools to the improvement of the Rosetta energy function. The feature analysis tool identifies discrepancies between structures deposited in the PDB and low-energy structures generated by Rosetta; these likely arise from inaccuracies in the energy function. The optE tool optimizes the weights on the different components of the energy function by maximizing the recapitulation of a wide range of experimental observations. We use the tools to examine three proposed modifications to the Rosetta energy function: improving the unfolded state energy model (reference energies), using bicubic spline interpolation to generate knowledge-based torisonal potentials, and incorporating the recently developed Dunbrack 2010 rotamer library (Shapovalov & Dunbrack, 2011)., (Copyright © 2013 Elsevier Inc. All rights reserved.)

The ability to engineer novel protein folds, conformations, and enzymatic activities offers enormous potential for the development of new protein therapeutics and biocatalysts. However, many de novo and redesigned proteins exhibit poor hydrophobic packing in their predicted structures, leading to instability or insolubility. The general utility of rational, structure-based design would greatly benefit from an improved ability to generate well-packed conformations. Here we present an automated protocol within the RosettaDesign framework that can identify and improve poorly packed protein cores by selecting a series of stabilizing point mutations. We apply our method to previously characterized designed proteins that exhibited a decrease in stability after a full computational redesign. We further demonstrate the ability of our method to improve the thermostability of a well-behaved native protein. In each instance, biophysical characterization reveals that we were able to stabilize the original proteins against chemical and thermal denaturation. We believe our method will be a valuable tool for both improving upon designed proteins and conferring increased stability upon native proteins.

We have recently completed a full re-architecturing of the ROSETTA molecular modeling program, generalizing and expanding its existing functionality. The new architecture enables the rapid prototyping of novel protocols by providing easy-to-use interfaces to powerful tools for molecular modeling. The source code of this rearchitecturing has been released as ROSETTA3 and is freely available for academic use. At the time of its release, it contained 470,000 lines of code. Counting currently unpublished protocols at the time of this writing, the source includes 1,285,000 lines. Its rapid growth is a testament to its ease of use. This chapter describes the requirements for our new architecture, justifies the design decisions, sketches out central classes, and highlights a few of the common tasks that the new software can perform., (© 2011 Elsevier Inc. All rights reserved.)

Modeling the conformational changes that occur on binding of macromolecules is an unsolved challenge. In previous rounds of the Critical Assessment of PRediction of Interactions (CAPRI), it was demonstrated that the Rosetta approach to macromolecular modeling could capture side chain conformational changes on binding with high accuracy. In rounds 13-19 we tested the ability of various backbone remodeling strategies to capture the main-chain conformational changes observed during binding events. These approaches span a wide range of backbone motions, from limited refinement of loops to relieve clashes in homologous docking, through extensive remodeling of loop segments, to large-scale remodeling of RNA. Although the results are encouraging, major improvements in sampling and energy evaluation are clearly required for consistent high accuracy modeling. Analysis of our failures in the CAPRI challenges suggest that conformational sampling at the termini of exposed beta strands is a particularly pressing area for improvement., (© Wiley-Liss, Inc.)

A long-standing goal of computational protein design is to create proteins similar to those found in Nature. One motivation is to harness the exquisite functional capabilities of proteins for our own purposes. The extent of similarity between designed and natural proteins also reports on how faithfully our models represent the selective pressures that determine protein sequences. As the field of protein design shifts emphasis from reproducing native-like protein structure to function, it has become important that these models treat the notion of specificity in molecular interactions. Although specificity may, in some cases, be achieved by optimization of a desired protein in isolation, methods have been developed to address directly the desire for proteins that exhibit specific functions and interactions.

Site-specific homing endonucleases are capable of inducing gene conversion via homologous recombination. Reprogramming their cleavage specificities allows the targeting of specific biological sites for gene correction or conversion. We used computational protein design to alter the cleavage specificity of I-MsoI for three contiguous base pair substitutions, resulting in an endonuclease whose activity and specificity for its new site rival that of wild-type I-MsoI for the original site. Concerted design for all simultaneous substitutions was more successful than a modular approach against individual substitutions, highlighting the importance of context-dependent redesign and optimization of protein-DNA interactions. We then used computational design based on the crystal structure of the designed complex, which revealed significant unanticipated shifts in DNA conformation, to create an endonuclease that specifically cleaves a site with four contiguous base pair substitutions. Our results demonstrate that specificity switches for multiple concerted base pair substitutions can be computationally designed, and that iteration between design and structure determination provides a route to large scale reprogramming of specificity.

Enzymes use substrate-binding energy both to promote ground-state association and to stabilize the reaction transition state selectively. The monomeric homing endonuclease I-AniI cleaves with high sequence specificity in the centre of a 20-base-pair (bp) DNA target site, with the amino (N)-terminal domain of the enzyme making extensive binding interactions with the left (-) side of the target site and the similarly structured carboxy (C)-terminal domain interacting with the right (+) side. Here we show that, despite the approximate twofold symmetry of the enzyme-DNA complex, there is almost complete segregation of interactions responsible for substrate binding to the (-) side of the interface and interactions responsible for transition-state stabilization to the (+) side. Although single base-pair substitutions throughout the entire DNA target site reduce catalytic efficiency, mutations in the (-) DNA half-site almost exclusively increase the dissociation constant (K(D)) and the Michaelis constant under single-turnover conditions (K(M)*), and those in the (+) half-site primarily decrease the turnover number (k(cat)*). The reduction of activity produced by mutations on the (-) side, but not mutations on the (+) side, can be suppressed by tethering the substrate to the endonuclease displayed on the surface of yeast. This dramatic asymmetry in the use of enzyme-substrate binding energy for catalysis has direct relevance to the redesign of endonucleases to cleave genomic target sites for gene therapy and other applications. Computationally redesigned enzymes that achieve new specificities on the (-) side do so by modulating K(M)*, whereas redesigns with altered specificities on the (+) side modulate k(cat)*. Our results illustrate how classical enzymology and modern protein design can each inform the other.

Computational protein design relies on a number of approximations to efficiently search the huge sequence space available to proteins. The fixed backbone and rotamer approximations in particular are important for formulating protein design as a discrete combinatorial optimization problem. However, the resulting coarse-grained sampling of possible side-chain terminal positions is problematic for the design of protein function, which depends on precise positioning of side-chain atoms. Although backbone flexibility can greatly increase the conformation freedom of side-chain functional groups, it is not obvious which backbone movements will generate the critical constellation of atoms responsible for protein function. Here, we report an automated method for identifying protein backbone movements that can give rise to any specified set of desired side-chain atomic placements and interactions, using protein-DNA interfaces as a model system. We use a library of previously observed protein-DNA interactions (motifs) and a rotamer-based description of side-chain conformation freedom to identify placements for the protein backbone that can give rise to a favorable side-chain interaction with DNA. We describe a tree-search algorithm for identifying those combinations of interactions from the library that can be realized with minimal perturbation of the protein backbone. We compare the efficiency of this method with the alternative approach of building and screening alternate backbone conformations.

Protein tyrosine kinases are critical cell signaling enzymes. These enzymes have a highly conserved Arg residue in their catalytic loop which is present two residues or four residues downstream from an absolutely conserved Asp catalytic base. Prior studies on protein tyrosine kinases Csk and Src revealed the potential for chemical rescue of catalytically deficient mutant kinases (Arg to Ala mutations) by small diamino compounds, particularly imidazole; however, the potency and efficiency of rescue was greater for Src. This current study further examines the structural and kinetic basis of rescue for mutant Src as compared to mutant Abl tyrosine kinase. An X-ray crystal structure of R388A Src revealed the surprising finding that a histidine residue of the N-terminus of a symmetry-related kinase inserts into the active site of the adjacent Src and mimics the hydrogen-bonding pattern seen in wild-type protein tyrosine kinases. Abl R367A shows potent and efficient rescue more comparable to Src, even though its catalytic loop is more like that of Csk. Various enzyme redesigns of the active sites indicate that the degree and specificity of rescue are somewhat flexible, but the overall properties of the enzymes and rescue agents play an overarching role. The newly discovered rescue agent 2-aminoimidazole is about as efficient as imidazole in rescuing R/A Src and Abl. Rate vs pH studies with these imidazole analogues suggest that the protonated imidazolium is the preferred form for chemical rescue, consistent with structural models. The efficient rescue seen with mutant Abl points to the potential of this approach to be used effectively to analyze Abl phosphorylation pathways in cells.

We recently reported that certain mutations in the FK506-rapamycin binding (FRB) domain disrupt its stability in vitro and in vivo (Stankunas et al. Mol. Cell, 2003, 12, 1615). To determine the precise residues that cause instability, we calculated the folding free energy (Delta G) of a collection of FRB mutants by measuring their intrinsic tryptophan fluorescence during reversible chaotropic denaturation. Our results implicate the T2098L point mutation as a key determinant of instability. Further, we found that some of the mutants in this collection were destabilized by up to 6 kcal mol(-1) relative to the wild type. To investigate how these mutants behave in cells, we expressed firefly luciferase fused to FRB mutants in African green monkey kidney (COS) cell lines and mouse embryonic fibroblasts (MEFs). When unstable FRB mutants were used, we found that the protein levels and the luminescence intensities were low. However, addition of a chemical ligand for FRB, rapamycin, restored luciferase activity. Interestingly, we found a roughly linear relationship between the Delta G of the FRB mutants calculated in vitro and the relative chemical rescue in cells. Because rapamycin is capable of simultaneously binding both FRB and the chaperone, FK506-binding protein (FKBP), we next examined whether FKBP might contribute to the protection of FRB mutants. Using both in vitro experiments and a cell-based model, we found that FKBP stabilizes the mutants. These findings are consistent with recent models that suggest damage to intrinsic Delta G can be corrected by pharmacological chaperones. Further, these results provide a collection of conditionally stable fusion partners for use in controlling protein stability.

Recent efforts to design de novo or redesign the sequence and structure of proteins using computational techniques have met with significant success. Most, if not all, of these computational methodologies attempt to model atomic-level interactions, and hence high-resolution structural characterization of the designed proteins is critical for evaluating the atomic-level accuracy of the underlying design force-fields. We previously used our computational protein design protocol RosettaDesign to completely redesign the sequence of the activation domain of human procarboxypeptidase A2. With 68% of the wild-type sequence changed, the designed protein, AYEdesign, is over 10 kcal/mol more stable than the wild-type protein. Here, we describe the high-resolution crystal structure and solution NMR structure of AYEdesign, which show that the experimentally determined backbone and side-chains conformations are effectively superimposable with the computational model at atomic resolution. To isolate the origins of the remarkable stabilization, we have designed and characterized a new series of procarboxypeptidase mutants that gain significant thermodynamic stability with a minimal number of mutations; one mutant gains more than 5 kcal/mol of stability over the wild-type protein with only four amino acid changes. We explore the relationship between force-field smoothing and conformational sampling by comparing the experimentally determined free energies of the overall design and these focused subsets of mutations to those predicted using modified force-fields, and both fixed and flexible backbone sampling protocols.

Background: Growing evidence suggests that distal regulatory elements are essential for cellular function and states. The sequences within these distal elements, especially motifs for transcription factor binding, provide critical information about the underlying regulatory programs. However, cooperativities between transcription factors that recognize these motifs are nonlinear and multiplexed, rendering traditional modeling methods insufficient to capture the underlying mechanisms. Recent development of attention mechanism, which exhibit superior performance in capturing dependencies across input sequences, makes them well-suited to uncover and decipher intricate dependencies between regulatory elements. Result: We present Transcription factors cooperativity Inference Analysis with Neural Attention (TIANA), a deep learning framework that focuses on interpretability. In this study, we demonstrated that TIANA could discover biologically relevant insights into co-occurring pairs of transcription factor motifs. Compared with existing tools, TIANA showed superior interpretability and robust performance in identifying putative transcription factor cooperativities from co-occurring motifs. Conclusion: Our results suggest that TIANA can be an effective tool to decipher transcription factor cooperativities from distal sequence data. TIANA can be accessed through: https://github.com/rzzli/TIANA. [ABSTRACT FROM AUTHOR]

Recent years have witnessed the transformative impact from the integration of artificial intelligence with organic and polymer synthesis. This synergy offers innovative and intelligent solutions to a range of classic problems in synthetic chemistry. These exciting advancements include the prediction of molecular property, multi-step retrosynthetic pathway planning, elucidation of the structure-performance relationship of single-step transformation, establishment of the quantitative linkage between polymer structures and their functions, design and optimization of polymerization process, prediction of the structure and sequence of biological macromolecules, as well as automated and intelligent synthesis platforms. Chemists can now explore synthetic chemistry with unprecedented precision and efficiency, creating novel reactions, catalysts, and polymer materials under the data-driven paradigm. Despite these thrilling developments, the field of artificial intelligence (AI) synthetic chemistry is still in its infancy, facing challenges and limitations in terms of data openness, model interpretability, as well as software and hardware support. This review aims to provide an overview of the current progress, key challenges, and future development suggestions in the interdisciplinary field between AI and synthetic chemistry. It is hoped that this overview will offer readers a comprehensive understanding of this emerging field, inspiring and promoting further scientific research and development. [ABSTRACT FROM AUTHOR]

With recent methodological advances in the field of computational protein design, in particular those based on deep learning, there is an increasing need for frameworks that allow for coherent, direct integration of different models and objective functions into the generative design process. Here we demonstrate how evolutionary multiobjective optimization techniques can be adapted to provide such an approach. With the established Non-dominated Sorting Genetic Algorithm II (NSGA-II) as the optimization framework, we use AlphaFold2 and ProteinMPNN confidence metrics to define the objective space, and a mutation operator composed of ESM-1v and ProteinMPNN to rank and then redesign the least favorable positions. Using the two-state design problem of the foldswitching protein RfaH as an in-depth case study, and PapD and calmodulin as examples of higher-dimensional design problems, we show that the evolutionary multiobjective optimization approach leads to significant reduction in the bias and variance in RfaH native sequence recovery, compared to a direct application of ProteinMPNN. We suggest that this improvement is due to three factors: (i) the use of an informative mutation operator that accelerates the sequence space exploration, (ii) the parallel, iterative design process inherent to the genetic algorithm that improves upon the ProteinMPNN autoregressive sequence decoding scheme, and (iii) the explicit approximation of the Pareto front that leads to optimal design candidates representing diverse tradeoff conditions. We anticipate this approach to be readily adaptable to different models and broadly relevant for protein design tasks with complex specifications. Author summary: Proteins are the fundamental building blocks of life, and engineering them has broad applications in medicine and biotechnology. Computational methods that seek to model and predict the protein sequence-structure-function relationship have seen significant advancement from the recent development in deep learning. As more models become available, it remains an open question how to effectively combine them into a coherent computational design approach. One approach is to perform computational design with one model, and filter the design candidates with the others. In this work, we demonstrate how an optimization algorithm inspired by evolution can be adapted to provide an alternative framework that outperforms the post hoc filtering approach, especially for problems with multiple competing design specifications. Such a framework may enable researchers to more effectively integrate the strengths of different modeling approaches to tackle more complex design problems. [ABSTRACT FROM AUTHOR]

The practical application of new single molecule protein sequencing (SMPS) technologies requires accurate estimates of their associated sequencing error rates. Here, we describe the development and application of two distinct parameter estimation methods for analyzing SMPS reads produced by fluorosequencing. A Hidden Markov Model (HMM) based approach, extends whatprot, where we previously used HMMs for SMPS peptide-read matching. This extension offers a principled approach for estimating key parameters for fluorosequencing experiments, including missed amino acid cleavages, dye loss, and peptide detachment. Specifically, we adapted the Baum-Welch algorithm, a standard technique to estimate transition probabilities for an HMM using expectation maximization, but modified here to estimate a small number of parameter values directly rather than estimating every transition probability independently. We demonstrate a high degree of accuracy on simulated data, but on experimental datasets, we observed that the model needed to be augmented with an additional error type, N-terminal blocking. This, in combination with data pre-processing, results in reasonable parameterizations of experimental datasets that agree with controlled experimental perturbations. A second independent implementation using a hybrid of DIRECT and Powell's method to reduce the root mean squared error (RMSE) between simulations and the real dataset was also developed. We compare these methods on both simulated and real data, finding that our Baum-Welch based approach outperforms DIRECT and Powell's method by most, but not all, criteria. Although some discrepancies between the results exist, we also find that both approaches provide similar error rate estimates from experimental single molecule fluorosequencing datasets. Author summary: Diverse new technologies are being developed for single-molecule protein sequencing, capable of identifying and quantifying mixtures of proteins at the level of individual molecules. There are many biochemical challenges intrinsic to high-throughput studies of proteins at such high sensitivity arising from their heterogeneous chemistries, sizes, and abundances. Beyond these challenges, the technologies themselves involve complex multi-step analytical processes. Thus, in developing and optimizing these technologies, it is important to consider the accuracy of each step and to have reliable approaches for estimating these accuracies. We focus on one particular single-molecule sequencing technology known as flourosequencing. We report and validate two methods for simultaneously determining the error-rates of each of the various steps of the fluorosequencing process. These new error estimation techniques will help researchers to better interpret the effects of changes to the chemistry and sample preparation used in fluorosequencing so that these steps can be improved. Further, more accurate determination of error rates will aid in the creation of better tools for the interpretation of this data. [ABSTRACT FROM AUTHOR]

The reprogramming of DNA-binding specificity is an important challenge for computational protein design that tests current understanding of protein-DNA recognition, and has considerable practical relevance for biotechnology and medicine. Here we describe the computational redesign of the cleavage specificity of the intron-encoded homing endonuclease I-MsoI using a physically realistic atomic-level forcefield. Using an in silico screen, we identified single base-pair substitutions predicted to disrupt binding by the wild-type enzyme, and then optimized the identities and conformations of clusters of amino acids around each of these unfavourable substitutions using Monte Carlo sampling. A redesigned enzyme that was predicted to display altered target site specificity, while maintaining wild-type binding affinity, was experimentally characterized. The redesigned enzyme binds and cleaves the redesigned recognition site approximately 10,000 times more effectively than does the wild-type enzyme, with a level of target discrimination comparable to the original endonuclease. Determination of the structure of the redesigned nuclease-recognition site complex by X-ray crystallography confirms the accuracy of the computationally predicted interface. These results suggest that computational protein design methods can have an important role in the creation of novel highly specific endonucleases for gene therapy and other applications.

Protein-DNA interactions play a central role in transcriptional regulation and other biological processes. Investigating the mechanism of binding affinity and specificity in protein-DNA complexes is thus an important goal. Here we develop a simple physical energy function, which uses electrostatics, solvation, hydrogen bonds and atom-packing terms to model direct readout and sequence-specific DNA conformational energy to model indirect readout of DNA sequence by the bound protein. The predictive capability of the model is tested against another model based only on the knowledge of the consensus sequence and the number of contacts between amino acids and DNA bases. Both models are used to carry out predictions of protein-DNA binding affinities which are then compared with experimental measurements. The nearly additive nature of protein-DNA interaction energies in our model allows us to construct position-specific weight matrices by computing base pair probabilities independently for each position in the binding site. Our approach is less data intensive than knowledge-based models of protein-DNA interactions, and is not limited to any specific family of transcription factors. However, native structures of protein-DNA complexes or their close homologs are required as input to the model. Use of homology modeling can significantly increase the extent of our approach, making it a useful tool for studying regulatory pathways in many organisms and cell types.

Protein-DNA interactions are crucial for many biological processes. Attempts to model these interactions have generally taken the form of amino acid-base recognition codes or purely sequence-based profile methods, which depend on the availability of extensive sequence and structural information for specific structural families, neglect side-chain conformational variability, and lack generality beyond the structural family used to train the model. Here, we take advantage of recent advances in rotamer-based protein design and the large number of structurally characterized protein-DNA complexes to develop and parameterize a simple physical model for protein-DNA interactions. The model shows considerable promise for redesigning amino acids at protein-DNA interfaces, as design calculations recover the amino acid residue identities and conformations at these interfaces with accuracies comparable to sequence recovery in globular proteins. The model shows promise also for predicting DNA-binding specificity for fixed protein sequences: native DNA sequences are selected correctly from pools of competing DNA substrates; however, incorporation of backbone movement will likely be required to improve performance in homology modeling applications. Interestingly, optimization of zinc finger protein amino acid sequences for high-affinity binding to specific DNA sequences results in proteins with little or no predicted specificity, suggesting that naturally occurring DNA-binding proteins are optimized for specificity rather than affinity. When combined with algorithms that optimize specificity directly, the simple computational model developed here should be useful for the engineering of proteins with novel DNA-binding specificities.

Specific protein-protein interactions are crucial in signaling networks and for the assembly of multi-protein complexes, and represent a challenging goal for protein design. Optimizing interaction specificity requires both positive design, the stabilization of a desired interaction, and negative design, the destabilization of undesired interactions. Currently, no automated protein-design algorithms use explicit negative design to guide a sequence search. We describe a multi-state framework for engineering specificity that selects sequences maximizing the transfer free energy of a protein from a target conformation to a set of undesired competitor conformations. To test the multi-state framework, we engineered coiled-coil interfaces that direct the formation of either homodimers or heterodimers. The algorithm identified three specificity motifs that have not been observed in naturally occurring coiled coils. In all cases, experimental results confirm the predicted specificities.

Solvent plays a significant role in determining the electrostatic potential energy of proteins, most notably through its favorable interactions with charged residues and its screening of electrostatic interactions. These energetic contributions are frequently ignored in computational protein design and protein modeling methodologies because they are difficult to evaluate rapidly and accurately. To address this deficiency, we report a revised form of the original Tanford-Kirkwood continuum electrostatic model [Tanford, C. & Kirkwood, J. G. (1957) J. Am. Chem. Soc. 79, 5333-5339], which accounts for the effects of solvent polarization on charged atoms in proteins. The Tanford-Kirkwood model was modified to increase its speed and to improve its sensitivity to the details of protein structure. For the 37 electrostatic self-energies of the polar side-chains in bovine pancreatic trypsin inhibitor, and their 666 interaction energies, the modified Tanford-Kirkwood potential of mean force differs from a computationally intensive numerical potential (DelPhi) by root-mean-square errors of 0.6 kcal/mol and 0.08 kcal/mol, respectively. The Tanford-Kirkwood approach makes possible a realistic treatment of electrostatics in computationally demanding protein modeling calculations. For example, pH titration calculations for ovomucoid third domain that model polar side-chain relaxation (including >2 x 10(23) rotamer conformations of the protein) provide pKa values of unprecedented accuracy.

Efficient and accurate recognition of protein–DNA interactions is vital for understanding the molecular mechanisms of related biological processes and further guiding drug discovery. Although the current experimental protocols are the most precise way to determine protein–DNA binding sites, they tend to be labor-intensive and time-consuming. There is an immediate need to design efficient computational approaches for predicting DNA-binding sites. Here, we proposed ULDNA, a new deep-learning model, to deduce DNA-binding sites from protein sequences. This model leverages an LSTM-attention architecture, embedded with three unsupervised language models that are pre-trained on large-scale sequences from multiple database sources. To prove its effectiveness, ULDNA was tested on 229 protein chains with experimental annotation of DNA-binding sites. Results from computational experiments revealed that ULDNA significantly improves the accuracy of DNA-binding site prediction in comparison with 17 state-of-the-art methods. In-depth data analyses showed that the major strength of ULDNA stems from employing three transformer language models. Specifically, these language models capture complementary feature embeddings with evolution diversity, in which the complex DNA-binding patterns are buried. Meanwhile, the specially crafted LSTM-attention network effectively decodes evolution diversity-based embeddings as DNA-binding results at the residue level. Our findings demonstrated a new pipeline for predicting DNA-binding sites on a large scale with high accuracy from protein sequence alone. [ABSTRACT FROM AUTHOR]

Objective: Glutamate transporters play a key role in central nervous system physiology by maintaining excitatory neurotransmitter homeostasis. Biological assemblies of the transporters, consisting of cyclic homotrimers, emerge as a crucial aspect of glutamate transporter modulation. Hence targeting heteromerization promises an effective approach for modulator design. On the other hand, the dynamic nature of transcription allows for the generation of transporter isoforms in structurally distinct manners., Methods: The potential isoforms were identified through the analysis of computationally generated gene-centric isoform maps. The conserved features of isoform sequences were revealed by computational chemistry methods and subsequent structural analysis of AlphaFold2 predictions. Truncated isoforms were further subjected to a wide range of docking analyses, 50ns molecular dynamics simulations, and evolutionary coupling analyses., Results: Energetic landscapes of isoform-canonical transporter complexes suggested an inhibitory potential of truncated isoforms on glutamate transporter bio-assembly. Moreover, isoforms that mimic the trimerization domain (in particular, TM2 helices) exhibited stronger interactions with canonical transporters, underscoring the role of transmembrane helices in isoform interactions. Additionally, self-assembly dynamics observed in truncated isoforms mimicking canonical TM5 helices indicate a potential protective role against unwanted interactions with canonical transporters., Conclusion: Our computational studies on glutamate transporters offer insights into the roles of alternative splicing on protein interactions and identifies potential drug targets for physiological or pathological processes., (© 2024. The Author(s).)

Amino acids (AAs) with a noncanonical backbone would be a valuable tool for protein engineering, enabling new structural motifs and building blocks. To incorporate them into an expanded genetic code, the first, key step is to obtain an appropriate aminoacyl‐tRNA synthetase. Currently, directed evolution is not available to optimize AAs with noncanonical backbones, since an appropriate selective pressure has not been discovered. Computational protein design (CPD) is an alternative. We used a new CPD method to redesign MetRS and increase its activity towards β‐Met, which has an extra backbone methylene. The new method considered a few active site positions for design and used a Monte Carlo exploration of the corresponding sequence space. During the exploration, a bias energy was adaptively learned, such that the free energy landscape of the apo enzyme was flattened. Enzyme variants could then be sampled, in the presence of the ligand and the bias energy, according to their β‐Met binding affinities. Eighteen predicted variants were chosen for experimental testing; 10 exhibited detectable activity for β‐Met adenylation. Top predicted hits were characterized experimentally in detail. Dissociation constants, catalytic rates, and Michaelis constants for both α‐Met and β‐Met were measured. The best mutant retained a preference for α‐Met over β‐Met; however, the preference was reduced, compared to the wildtype, by a factor of 29. For this mutant, high resolution crystal structures were obtained in complex with both α‐Met and β‐Met, indicating that the predicted, active conformation of β‐Met in the active site was retained. [ABSTRACT FROM AUTHOR]

Objectives: Although Geobacillus are significant thermophilic bacteria source, there are no reports of thermostable esterase gene in Geobacillus jurassicus or rational design strategies to increase the thermal stability of esterases., Results: Gene gju768 showed a highest similarity of 15.20% to esterases from Geobacillus sp. with detail enzymatic properties. Using a combination of Gibbs Unfolding Free Energy (∆∆G) calculator and the distance from the mutation site to the catalytic site (Ds Cα-Cα ) to screen suitable mutation sites with elimination of negative surface charge, the mutants (D24N, E221Q, and E253Q) displayed stable mutants with higher thermal stability than the wild-type (WT). Mutant E253Q exhibited the best thermal stability, with a half-life (T 1/2 ) at 65 °C of 32.4 min, which was 1.8-fold of the WT (17.9 min)., Conclusion: Cloning of gene gju768 and rational design based on surface charge engineering contributed to the identification of thermostable esterase from Geobacillus sp. and the exploration of evolutionary strategies for thermal stability., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

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

The critical first step in Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)–associated (CRISPR–Cas) protein-mediated gene editing is recognizing a preferred protospacer adjacent motif (PAM) on target DNAs by the protein's PAM-interacting amino acids (PIAAs). Thus, accurate computational modeling of PAM recognition is useful in assisting CRISPR–Cas engineering to relax or tighten PAM requirements for subsequent applications. Here, we describe a universal computational protein design framework (UniDesign) for designing protein–nucleic acid interactions. As a proof of concept, we applied UniDesign to decode the PAM–PIAA interactions for eight Cas9 and two Cas12a proteins. We show that, given native PIAAs, the UniDesign-predicted PAMs are largely identical to the natural PAMs of all Cas proteins. In turn, given natural PAMs, the computationally redesigned PIAA residues largely recapitulated the native PIAAs (74% and 86% in terms of identity and similarity, respectively). These results demonstrate that UniDesign faithfully captures the mutual preference between natural PAMs and native PIAAs, suggesting it is a useful tool for engineering CRISPR–Cas and other nucleic acid-interacting proteins. UniDesign is open-sourced at https://github.com/tommyhuangthu/UniDesign. [ABSTRACT FROM AUTHOR]

Leber congenital amaurosis (LCA) is the most severe form of inherited retinal dystrophies and the most frequent cause of congenital blindness in children. To date, 25 genes have been implicated in the pathogenesis of this rare disorder. Performing an accurate molecular diagnosis is crucial as gene therapy is becoming available. This study aimed to report the molecular basis of Leber congenital amaurosis, especially novel and rare variants in 27 Polish families with a clinical diagnosis of LCA fully confirmed by molecular analyses. Whole exome sequencing or targeted next-generation sequencing (NGS) of inherited retinal dystrophies-associated (IRD) genes was applied to identify potentially pathogenic variants. Bidirectional Sanger sequencing and quantitative PCR (qPCR) were carried out for validation and segregation analysis of the variants identified within the families. We identified 28 potentially pathogenic variants, including 11 novel, in 8 LCA genes: CEP290, CRB1, GUCY2D, NMNAT1, RPGRIP1, CRX, LRAT1, and LCA5. This study expands the mutational spectrum of the LCA genes. Moreover, these results, together with the conclusions from our previous studies, allow us to point to the most frequently mutated genes and variants in the Polish cohort of LCA patients. [ABSTRACT FROM AUTHOR]

Binding of transcription factor (TF) proteins to regulatory DNA sites is key to accurate control of gene expression in response to environmental stimuli. Theoretical modeling of transcription regulation is often focused on a limited set of genes of interest, while binding of the TF to other genomic sites is seldom considered. The total number of TF binding sites (TFBSs) affects the availability of TF protein molecules and sequestration of a TF by TFBSs can promote bistability. For many signaling systems where a graded response is desirable for continuous control over the input range, biochemical parameters of the regulatory proteins need be tuned to avoid bistability. Here we analyze the mono-/bistable parameter range for positively autoregulated two-component systems (TCSs) in the presence of different numbers of competing TFBSs. TCS signaling, one of the major bacterial signaling strategies, couples signal perception with output responses via protein phosphorylation. For bistability, competition for TF proteins by TFBSs lowers the requirement for high fold change of the autoregulated transcription but demands high phosphorylation activities of TCS proteins. We show that bistability can be avoided with a low phosphorylation capacity of TCSs, a high TF affinity for the autoregulated promoter or a low fold change in signaling protein levels upon induction. These may represent general design rules for TCSs to ensure uniform graded responses. Examining the mono-/bistability parameter range allows qualitative prediction of steady-state responses, which are experimentally validated in the E. coli CusRS system. Author summary: Cell survival in an ever-changing environment depends on appropriate responses to stimuli of different strengths. Bistability, i.e., two different stable states in responses to otherwise identical environments, can be beneficial in some systems but may need to be avoided by many signaling systems. Promoting or preventing bistable responses requires specific architectures of gene regulatory networks as well as proper abundance and activities for the network building blocks. Here we use a mathematical model to study how the requirement for bistable or monostable responses places constraints on biochemical properties of autoregulated bacterial signaling systems. A particular focus is on how the relative abundance of transcription factor proteins to the number of their DNA target sites constrains the system design because competition for the limited TF protein molecules by DNA sites can promote bistability. We find that a strong binding affinity to the autoregulated TF promoter and low phosphorylation activities are preferred for monostability. We used an E. coli signaling system as an example and experimentally validated predictions of the model. Our results can help rationalize the regulatory features observed in naturally occurring systems as well as inform engineering of novel biological circuits for diverse signaling tasks. [ABSTRACT FROM AUTHOR]

Background: The overprescription and misuse of classical antimicrobial compounds to treat gastrointestinal or systemic salmonellosis have been accelerating the surge of antibiotic-recalcitrant bacterial populations, posing a major public health challenge. Therefore, alternative therapeutic approaches to treat Salmonella infections are urgently required.Objectives: To identify and characterize actinobacterial secreted compounds with inhibitory properties against the Salmonella enterica PhoP/PhoQ signal transduction system, crucial for virulence regulation.Methods: The methodology was based on a combination of the measurement of the activity of PhoP/PhoQ-dependent and -independent reporter genes and bioguided assays to screen for bioactive inhibitory metabolites present in culture supernatants obtained from a collection of actinobacterial isolates. Analogues of azomycin were used to analyse the functional groups required for the detected bioactivity and Salmonella mutants and complemented strains helped to dissect the azomycin mechanism of action. The tetrazolium dye colorimetric assay was used to investigate azomycin potential cytotoxicity on cultured macrophages. Salmonella intramacrophage replication capacity upon azomycin treatment was assessed using the gentamicin protection assay.Results: Sublethal concentrations of azomycin, a nitroheterocyclic compound naturally produced by Streptomyces eurocidicus, repressed the Salmonella PhoP/PhoQ system activity by targeting PhoP and inhibiting its transcriptional activity in a PhoQ- and aspartate phosphorylation-independent manner. Sublethal, non-cytotoxic concentrations of azomycin prevented Salmonella intramacrophage replication.Conclusions: Azomycin selectively inhibits the activity of the Salmonella virulence regulator PhoP, a new activity described for this nitroheterocyclic compound that can be repurposed to develop novel anti-Salmonella therapeutic approaches. [ABSTRACT FROM AUTHOR]

The aim of this work was to determine the concentrations of selected mutagenic elements (As, Cd, Co, Cr, Cu, Fe, Mn, Ni, Pb, and Se) in raw ewe milk from undisturbed environment in Orava region, northern Slovakia. There are possible sources of some of the analyzed metals which may be distributed from the metallurgical plants located in the Ostrava region, Czech Republic, and Katowice, Poland. In total, forty milk samples were analyzed in June and August using an inductively coupled plasma optical emission spectrometry. The differences in elements concentrations between the seasonal periods were not significant except of iron (p < 0.0001). The concentrations of most of the metals in ewe milk were low and under the permissible or recommended limits. However, arsenic and selenium concentrations were elevated and could pose a risk of the mutagenic effect, particularly in children. The frequency of element occurrence in June was as follows: Se > Fe > As > Cu > Mn > Ni > Co > Pb > Cr > Cd, and in August: Se > Fe = As > Cu > Mn > Pb > Co > Ni > Cr > Cd. The correlation analysis revealed very strong positive correlation between Cu:Pb (p < 0.05), very strong negative correlation between Fe:Se (p < 0.05). The strong correlations were also found between other elements. The present study showed that milk produced in the relatively undisturbed environment might contain various mutagenic elements. The relationships between the elements might result in the additive or synergistic effects of elements and increase the risk of their mutagenic effects even in low concentrations. Therefore, attention must be paid to the monitoring of metals in the areas where food sources destined especially for child nutrition are produced. [ABSTRACT FROM AUTHOR]

The stabilization of natural proteins is a long‐standing desired goal in protein engineering. Optimizing the hydrophobicity of the protein core often results in extensive stability enhancements. However, the presence of totally or partially buried catalytic charged residues, essential for protein function, has limited the applicability of this strategy. Here, focusing on the thioredoxin, we aimed to augment protein stability by removing buried charged residues in the active site without loss of catalytic activity. To this end, we performed a charged‐to‐hydrophobic substitution of a buried and functional group, resulting in a significant stability increase yet abolishing catalytic activity. Then, to simulate the catalytic role of the buried ionizable group, we designed a combinatorial library of variants targeting a set of seven surface residues adjacent to the active site. Notably, more than 50% of the library variants restored, to some extent, the catalytic activity. The combination of experimental study of 2% of the library with the prediction of the whole mutational space by partial least squares regression revealed that a single point mutation at the protein surface is sufficient to fully restore the catalytic activity without thermostability cost. As a result, we engineered one of the highest thermal stabilities reported for a protein with a natural occurring fold (137°C). Further, our hyperstable variant preserves the catalytic activity both in vitro and in vivo. [ABSTRACT FROM AUTHOR]

Background: Essential proteins are indispensable to the development and survival of cells. The identification of essential proteins not only is helpful for the understanding of the minimal requirements for cell survival, but also has practical significance in disease diagnosis, drug design and medical treatment. With the rapidly amassing of protein–protein interaction (PPI) data, computationally identifying essential proteins from protein–protein interaction networks (PINs) becomes more and more popular. Up to now, a number of various approaches for essential protein identification based on PINs have been developed. Results: In this paper, we propose a new and effective approach called iMEPP to identify essential proteins from PINs by fusing multiple types of biological data and applying the influence maximization mechanism to the PINs. Concretely, we first integrate PPI data, gene expression data and Gene Ontology to construct weighted PINs, to alleviate the impact of high false-positives in the raw PPI data. Then, we define the influence scores of nodes in PINs with both orthological data and PIN topological information. Finally, we develop an influence discount algorithm to identify essential proteins based on the influence maximization mechanism. Conclusions: We applied our method to identifying essential proteins from saccharomyces cerevisiae PIN. Experiments show that our iMEPP method outperforms the existing methods, which validates its effectiveness and advantage. [ABSTRACT FROM AUTHOR]

The molecular interaction properties and aggregation capabilities disclosed by PrP-E200K, a pathogenic mutant of the human prion protein, were investigated in detail using multilayered computational approaches. In a previous work, we reported that the electrostatic complementarity between region1 (negative) and region3 (positive) has been assumed to lead to a head-to tail interaction between 120 and 231 PrP-E200K units and to initiation of the aggregation process. In this work, we extended the PrP-E200K structure by including the unstructured 90–120 segment which was found to assume different conformations. Plausible models of 90–231 PrP-E200K dimers were calculated and analyzed in depth to identify the nature of the involved protein–protein interactions. The unstructured 90–120 segment was found to extend the positively charged region3 involved in the association of PrP-E200K units which resulted to be driven by hydrophobic interactions. The combination of molecular dynamics, protein–protein docking, grid-based mapping, and fragment molecular orbital approaches allowed us to provide a plausible mechanism of the early state of 90–231 PrP-E200K aggregation, considered a preliminary step of amyloid conversion. [ABSTRACT FROM AUTHOR]

The enzyme nicotinamide mononucleotide adenylyltransferase (NMNAT) catalyzes a reaction central to all known NAD biosynthetic routes. In mammals, three isoforms with distinct molecular and catalytic properties, different subcellular and tissue distribution have been characterized. Each isoform is essential for cell survival, with a critical role in modulating NAD levels in a compartment‐specific manner. Each isoform supplies NAD to specific NAD‐dependent enzymes, thus regulating their activity with impact on several biological processes, including DNA repair, proteostasis, cell differentiation, and neuronal maintenance. The nuclear NMNAT1 and the cytoplasmic NMNAT2 are also emerging as relevant targets in specific types of cancers and NMNAT2 has a key role in the activation of antineoplastic compounds. This review recapitulates the biochemical properties of the three isoforms and focuses on recent advances on their protective function, involvement in human diseases and role as druggable targets. [ABSTRACT FROM AUTHOR]

Leber congenital amaurosis (LCA), although rare, is one of the most severe forms of early-onset inherited retinal dystrophy (IRD). Here, we review the molecular genetics and phenotypic characteristics of patients with NMNAT1-associated IRD. The longitudinal clinical and molecular findings of a Japanese girl diagnosed with LCA associated with pathogenic variants in NMNAT1 c.648delG, (p.Trp216Ter*) and c.709C>T (p.Arg237Cys) have been described to highlight the salient clinical features of NMNAT1-associated IRD. [ABSTRACT FROM AUTHOR]