Enriching productive mutational paths accelerates enzyme evolution.

Darwinian evolution has given rise to all the enzymes that enable life on Earth. Mimicking natural selection, scientists have learned to tailor these biocatalysts through recursive cycles of mutation, selection and amplification, often relying on screening large protein libraries to productively modulate the complex interplay between protein structure, dynamics and function. Here we show that by removing destabilizing mutations at the library design stage and taking advantage of recent advances in gene synthesis, we can accelerate the evolution of a computationally designed enzyme. In only five rounds of evolution, we generated a Kemp eliminase-an enzymatic model system for proton transfer from carbon-that accelerates the proton abstraction step >108-fold over the uncatalyzed reaction. Recombining the resulting variant with a previously evolved Kemp eliminase HG3.17, which exhibits similar activity but differs by 29 substitutions, allowed us to chart the topography of the designer enzyme's fitness landscape, highlighting that a given protein scaffold can accommodate several, equally viable solutions to a specific catalytic problem.

Darobactin Substrate Engineering and Computation Show Radical Stability Governs Ether versus C-C Bond Formation.

The Gram-negative selective antibiotic darobactin A has attracted interest owing to its intriguing fused bicyclic structure and unique targeting of the outer membrane protein BamA. Darobactin, a ribosomally synthesized and post-translationally modified peptide (RiPP), is produced by a radical S-adenosyl methionine (rSAM)-dependent enzyme (DarE) and contains one ether and one C-C cross-link. Herein, we analyze the substrate tolerance of DarE and describe an underlying catalytic principle of the enzyme. These efforts produced 51 enzymatically modified darobactin variants, revealing that DarE can install the ether and C-C cross-links independently and in different locations on the substrate. Notable variants with fused bicyclic structures were characterized, including darobactin W3Y, with a non-Trp residue at the twice-modified central position, and darobactin K5F, which displays a fused diether ring pattern. While lacking antibiotic activity, quantum mechanical modeling of darobactins W3Y and K5F aided in the elucidation of the requisite features for high-affinity BamA engagement. We also provide experimental evidence for ß-oxo modification, which adds support for a proposed DarE mechanism. Based on these results, ether and C-C cross-link formation was investigated computationally, and it was determined that more stable and longer-lived aromatic Cß radicals correlated with ether formation. Further, molecular docking and transition state structures based on high-level quantum mechanical calculations support the different indole connectivity observed for ether (Trp-C7) and C-C (Trp-C6) cross-links. Finally, mutational analysis and protein structural predictions identified substrate residues that govern engagement to DarE. Our work informs on darobactin scaffold engineering and further unveils the underlying principles of rSAM catalysis.

AlphaFold2 Predicts Alternative Conformation Populations in Green Fluorescent Protein Variants.

Artificial intelligence-based protein structure prediction methods such as AlphaFold2 have emerged as powerful tools for characterizing proteins sequence-structure relationship offering unprecedented opportunities for the molecular interpretation of biological and biochemical phenomena. While initially confined to providing a static representation of proteins through their global free-energy minimum, AlphaFold2 has demonstrated the ability to partially sample conformational landscapes, providing insights into protein dynamics, which is fundamental for interpreting and potentially tuning the function of natural and artificial proteins. In this study, we show that targeted column masking of AlphaFold2's multiple sequence alignment enables the characterization and estimation of the population ratio of the two main conformations of engineered green fluorescent proteins with alternative ß-strands. The possibility of quickly estimating relative populations through AlphaFold2 predictions is expected to speed-up the computational design of related systems for sensing applications.

Conformationally Restricted ß-Sheet Breaker Peptides Incorporating Cyclic α-Methylisoserine Sulfamidates.

Peptides containing variations of the ß-amyloid hydrophobic core and five-membered sulfamidates derived from ß-amino acid α-methylisoserine have been synthesized and fully characterized in the gas phase, solid state and in aqueous solution by a combination of experimental and computational techniques. The cyclic sulfamidate group effectively locks the secondary structure at the N-terminus of such hybrid peptides imposing a conformational restriction and stabilizing non-extended structures. This conformational bias, which is maintained in the gas phase, solid state and aqueous solution, is shown to be resistant to structure templating through assays of in vitro ß-amyloid aggregation, acting as ß-sheet breaker peptides with moderate activity.

Unraveling the Photophysical Characteristics, Aromaticity, and Stability of π-Extended Acene-Quinodimethyl Thioamides†.

Poly-aromatic systems that contain quinodimethyl (QDM) units are appealing for several photonic and spintronic applications owing to the unique electronic structure, aromaticity, and spin state(s) of the QDM ring. Herein, we report the synthesis and characterization of novel QDM-based chromophores 1-3, which exhibit unique photo-excited behavior and aromaticity. Extending the aromatic core with a biphenyl/phenanthryl- and a pyrrolo-fragment led to reducing the optoelectronic bandgap and modulating the photophysics QDM 1-3. Yet, QDM 2 and 3 suffer from "aromaticity imbalance" and become relatively unstable compared to the parent compound QDM 1. Further assessment of local aromaticity using computational tools revealed that the pseudo-quinoidal ring B is the main driving force allowing to easily populate the excited triplet state of these chromophores. The present study provides complementary guidelines for designing novel non-classical poly-aromatic systems.

Accurate Prediction of Enzyme Thermostabilization with Rosetta Using AlphaFold Ensembles.

Thermostability enhancement is a fundamental aspect of protein engineering as a biocatalyst's half-life is key for its industrial and biotechnological application, particularly at high temperatures and under harsh conditions. Thermostability changes upon mutation originate from modifications of the free energy of unfolding (ΔGu), making thermostabilization extremely challenging to predict with computational methods. In this contribution, we combine global conformational sampling with energy prediction using AlphaFold and Rosetta to develop a new computational protocol for the quantitative prediction of thermostability changes upon laboratory evolution of acyltransferase LovD and lipase LipA. We highlight how using an ensemble of protein conformations rather than a single three-dimensional model is mandatory for accurate thermostability predictions. By comparing our approaches with existing ones, we show that ensembles based on AlphaFold models provide more accurate and robust calculated thermostability trends than ensembles based solely on crystallographic structures as the latter introduce a strong distortion (scaffold bias) in computed thermostabilities. Eliminating this bias is critical for computer-guided enzyme design and evaluating the effect of multiple mutations on protein stability.

The NCIWEB Server: A Novel Implementation of the Noncovalent Interactions Index for Biomolecular Systems.

It is well-known that the activity and function of proteins is strictly correlated with their secondary, tertiary, and quaternary structures. Their biological role is regulated by their conformational flexibility and global fold, which, in turn, is largely governed by complex noncovalent interaction networks. Because of the large size of proteins, the analysis of their noncovalent interaction networks is challenging, but can provide insights into the energetics of conformational changes or protein-protein and protein-ligand interactions. The noncovalent interaction (NCI) index, based on the reduced density gradient, is a well-established tool for the detection of weak contacts in biological systems. In this work, we present a web-based application to expand the use of this index to proteins, which only requires a molecular structure as input and provides a mapping of the number, type, and strength of noncovalent interactions. Structure preparation is automated and allows direct importing from the PDB database, making this server (https://nciweb.dsi.upmc.fr) accessible to scientists with limited experience in bioinformatics. A quick overview of this tool and concise instructions are presented, together with an illustrative application.

Iridium-Catalyzed ortho-Selective Borylation of Aromatic Amides Enabled by 5-Trifluoromethylated Bipyridine Ligands.

Iridium-catalyzed borylations of aromatic C-H bonds are highly attractive transformations because of the diversification possibilities offered by the resulting boronates. These transformations are best carried out using bidentate bipyridine or phenanthroline ligands, and tend to be governed by steric factors, therefore resulting in the competitive functionalization of meta and/or para positions. We have now discovered that a subtle change in the bipyridine ligand, namely, the introduction of a CF3 substituent at position 5, enables a complete change of regioselectivity in the borylation of aromatic amides, allowing the synthesis of a wide variety of ortho-borylated derivatives. Importantly, thorough computational studies suggest that the exquisite regio- and chemoselectivity stems from unusual outer-sphere interactions between the amide group of the substrate and the CF3 -substituted aryl ring of the bipyridine ligand.

Regioselective Synthesis and Molecular Docking Studies of 1,5-Disubstituted 1,2,3-Triazole Derivatives of Pyrimidine Nucleobases.

1,2,3-triazoles are versatile building blocks with growing interest in medicinal chemistry. For this reason, organic chemistry focuses on the development of new synthetic pathways to obtain 1,2,3-triazole derivatives, especially with pyridine moieties. In this work, a novel series of 1,5-disubstituted-1,2,3-triazoles functionalized with pyrimidine nucleobases were prepared via 1,3-dipolar cycloaddition reaction in a regioselective manner for the first time. The N1-propargyl nucleobases, used as an alkyne intermediate, were obtained in high yields (87-92%) with a new two-step procedure that selectively led to the monoalkylated compounds. Then, FeCl3 was employed as an efficient Lewis acid catalyst for 1,3-dipolar cycloaddition between different aryl and benzyl azides and the N1-propargyl nucleobases previously synthesized. This new protocol allows the synthesis of a series of new 1,2,3-triazole derivatives with good to excellent yields (82-92%). The ADME (Absorption, Distribution, Metabolism, and Excretion) analysis showed good pharmacokinetic properties and no violations of Lipinsky's rules, suggesting an appropriate drug likeness for these new compounds. Molecular docking simulations, conducted on different targets, revealed that two of these new hybrids could be potential ligands for viral and bacterial protein receptors such as human norovirus capsid protein, SARS-CoV-2 NSP13 helicase, and metallo-ß-lactamase.

Cosolute modulation of protein oligomerization reactions in the homeostatic timescale.

Protein oligomerization processes are widespread and of crucial importance to understand degenerative diseases and healthy regulatory pathways. One particular case is the homo-oligomerization of folded domains involving domain swapping, often found as a part of the protein homeostasis in the crowded cytosol, composed of a complex mixture of cosolutes. Here, we have investigated the effect of a plethora of cosolutes of very diverse nature on the kinetics of a protein dimerization by domain swapping. In the absence of cosolutes, our system exhibits slow interconversion rates, with the reaction reaching the equilibrium within the average protein homeostasis timescale (24-48 h). In the presence of crowders, though, the oligomerization reaction in the same time frame will, depending on the protein's initial oligomeric state, either reach a pure equilibrium state or get kinetically trapped into an apparent equilibrium. Specifically, when the reaction is initiated from a large excess of dimer, it becomes unsensitive to the effect of cosolutes and reaches the same equilibrium populations as in the absence of cosolute. Conversely, when the reaction starts from a large excess of monomer, the reaction during the homeostatic timescale occurs under kinetic control, and it is exquisitely sensitive to the presence and nature of the cosolute. In this scenario (the most habitual case in intracellular oligomerization processes), the effect of cosolutes on the intermediate conformation and diffusion-mediated encounters will dictate how the cellular milieu affects the domain-swapping reaction.

Impact of Cu(II) and Al(III) on the conformational landscape of amyloidß1-42.

Metal ions have been found to play an important role in the formation of extracellular ß-amyloid plaques, a major hallmark of Alzheimer's disease. In the present study, the conformational landscape of Aß42 with Al(iii) and Cu(ii) has been explored using Gaussian accelerated molecular dynamics. Both metals reduce the flexibility of the peptide and entail a higher structural organization, although to different degrees. As a general trend, Cu(ii) binding leads to an increased α-helix content and to the formation of two α-helices that tend to organize in a U-shape. By contrast, most Al(iii) complexes induce a decrease in helical content, leading to more extended structures that favor the appearance of transitory ß-strands.

Atomistic insights into the structure of heptapeptide nanofibers.

Artificial amyloid-like nanofibers formed from short peptides are emerging as new supramolecular structures for catalysis and advanced materials. In this work, we analyze, by means of computational approaches, the preferred atomistic fibrillar architectures that result from the self-assembly of polar NY7, NF7, SY7, SF7, and GY7 peptides into steric zippers formed by two ß-sheets (describing an individual steric zipper) and by four ß-sheets. For all heptapeptides, except GY7, parallel ß-sheet organizations with polar residues packed at the steric zipper appear to be the preferred assemblies for the two ß-sheets system due to the formation of a strong network of hydrogen bonds. For GY7, however, an antiparallel organization with glycine at the steric zipper is the most stable one. The preferred architecture is mostly conserved when enlarging our model from two to four ß-sheets. The present work shows that the relative stability of different architectures results from a delicate balance between peptide composition, side chain hydrophobicity, and non-covalent interactions at the interface and provides the basis for a rational design of new improved artificial prion-inspired materials.

Small molecule inhibits α-synuclein aggregation, disrupts amyloid fibrils, and prevents degeneration of dopaminergic neurons.

Parkinson's disease (PD) is characterized by a progressive loss of dopaminergic neurons, a process that current therapeutic approaches cannot prevent. In PD, the typical pathological hallmark is the accumulation of intracellular protein inclusions, known as Lewy bodies and Lewy neurites, which are mainly composed of α-synuclein. Here, we exploited a high-throughput screening methodology to identify a small molecule (SynuClean-D) able to inhibit α-synuclein aggregation. SynuClean-D significantly reduces the in vitro aggregation of wild-type α-synuclein and the familiar A30P and H50Q variants in a substoichiometric molar ratio. This compound prevents fibril propagation in protein-misfolding cyclic amplification assays and decreases the number of α-synuclein inclusions in human neuroglioma cells. Computational analysis suggests that SynuClean-D can bind to cavities in mature α-synuclein fibrils and, indeed, it displays a strong fibril disaggregation activity. The treatment with SynuClean-D of two PD Caenorhabditis elegans models, expressing α-synuclein either in muscle or in dopaminergic neurons, significantly reduces the toxicity exerted by α-synuclein. SynuClean-D-treated worms show decreased α-synuclein aggregation in muscle and a concomitant motility recovery. More importantly, this compound is able to rescue dopaminergic neurons from α-synuclein-induced degeneration. Overall, SynuClean-D appears to be a promising molecule for therapeutic intervention in Parkinson's disease.

NCIPLOT4 Guide for Biomolecules: An Analysis Tool for Noncovalent Interactions.

This Application Note is a guide for the use of the NCIPLOT4 code in the analysis of noncovalent interactions in biomacromolecular systems. Through a series of examples, the reader is walked through the process of calculating and interpreting noncovalent interaction density integrals corresponding to attractive, van der Waals, and repulsive terms in protein-ligand and protein-protein interaction problems. These integrals are robust and powerful tools to quickly obtain a semiquantitative picture of noncovalent interactions in complex systems with a temporal resolution (along a trajectory) and allow a flexible factorization of the noncovalent interactions network that permits the simultaneous evaluation of a wide range of contacts. The NCIPLOT4 code is publicly available at https://github.com/juliacontrerasgarcia/nciplot .

Structural Characterization of N-Linked Glycans in the Receptor Binding Domain of the SARS-CoV-2 Spike Protein and their Interactions with Human Lectins.

The glycan structures of the receptor binding domain of the SARS-CoV2 spike glycoprotein expressed in human HEK293F cells have been studied by using NMR. The different possible interacting epitopes have been deeply analysed and characterized, providing evidence of the presence of glycan structures not found in previous MS-based analyses. The interaction of the RBD 13 C-labelled glycans with different human lectins, which are expressed in different organs and tissues that may be affected during the infection process, has also been evaluated by NMR. In particular, 15 N-labelled galectins (galectins-3, -7 and -8 N-terminal), Siglecs (Siglec-8, Siglec-10), and C-type lectins (DC-SIGN, MGL) have been employed. Complementary experiments from the glycoprotein perspective or from the lectin's point of view have permitted to disentangle the specific interacting epitopes in each case. Based on these findings, 3D models of the interacting complexes have been proposed.

Fluorous l-Carbidopa Precursors: Highly Enantioselective Synthesis and Computational Prediction of Bioactivity.

New fluorous enantiopure (S)-α-aminated ß-keto esters were prepared through a highly enantioselective electrophilic α-amination step in the presence of europium triflate and (R,R)-phenyl-pybox. These compounds are precursors of fluorinated analogues of l-carbidopa, which is known to inhibit DOPA decarboxylase (DDC), a key protein in Parkinson's disease. Fluorination provides better stability for biological applications, which could possibly lead to DDC inhibitors better than l-carbidopa itself. Induced fit docking computational simulations performed on the new structures interacting with DDC highlight that for an efficient binding at the DDC site, at least one hydroxyl substituent must be present at the aromatic ring of the l-carbidopa analogues and show that the presence of fluorine can further fix the position of the ligand in the active site.

Visible-Light Photocatalytic Intramolecular Cyclopropane Ring Expansion.

Described herein is a new visible-light photocatalytic strategy for the synthesis of enantioenriched dihydrofurans and cyclopentenes by an intramolecular nitro cyclopropane ring expansion reaction. Mechanistic studies and DFT calculations are used to elucidate the key factors in this new ring expansion reaction, and the need for the nitro group on the cyclopropane.

Computational study on donor-acceptor optical markers for Alzheimer's disease: a game of charge transfer and electron delocalization.

According to the amyloid cascade hypothesis, amyloid-ß (Aß) deposition is a central event in the Alzheimer's disease and thus, detection of Aß deposits is crucial to monitor the progression of the pathology. Despite its low tissue penetration, fluorescence imaging may become an alternative technique for identifying these deposits because it is less toxic and less costly than positron emission tomography. Suitable dyes, however, should emit in the near infrared (NIR) region, cross the blood-brain barrier and target Aß aggregates. In this work, we use TD-DFT, AIMD simulations and protein energy landscape exploration (PELE) to analyze the photophysical properties of a family of donor-acceptor markers and their binding to amyloid fibrils. These markers are formed by N,N-dimethylaniline donor and propanedinitrile acceptor groups separated by a spacer consisting of one, two or three conjugated double bonds. The smallest compound has a low emission wavelength, can deactivate through a non-radiative process involving a conical intersection and binds weakly to Aß fibrils. In contrast, the largest dye is a suitable compound as it shows an emission wavelength in the NIR region, does not seem to relax through conical intersection processes and binds to Aß fibrils strongly entering hydrophobic voids. Analysis of electronic excitations shows that the transition has an important charge transfer character that increases with the length of the spacer, the π bridge being an active participant in the transition. Therefore, adding double bonds to the dye skeleton has two beneficial effects: (i) it increases the emission wavelength as it enlarges the π system and (ii) it increases the charge transfer character of the transition, which increases the red-shifting of the emission wavelength in polar solvents.

Disaggregation-induced fluorescence enhancement of NIAD-4 for the optical imaging of amyloid-ß fibrils.

The main hallmark of Alzheimer's disease is the deposition of amyloid-ß (Aß) aggregates in the brain. An early diagnosis of the disease requires a fast and accurate detection of such aggregates in vivo. NIAD-4 is one of the most promising in vivo markers developed due to its high emission at λ > 600 nm and its ability to rapidly cross the blood-brain barrier (BBB) and target Aß deposits. Furthermore, it shows a dramatic fluorescence enhancement upon binding to amyloid fibrils, which is essential for attaining good imaging contrast. Aiming at establishing novel design concepts for the preparation of optimized optical probes, the current work rationalizes the excellent performance of NIAD-4 by using a pool of computational (TD-DFT and CASPT2 calculations, ab initio molecular dynamics and protein energy landscape exploration) and spectroscopic techniques. Unlike other markers operating as molecular rotors or polarity-sensitive dyes, we uncover herein that the high fluorescence imaging contrast observed upon NIAD-4 binding to amyloid fibrils results from reversible aggregation. NIAD-4 forms non-emissive assemblies in aqueous solution already at very low concentrations, which convert into the highly fluorescent monomeric species by diffusion into the hydrophobic voids of Aß deposits. This result paves the way to exploit aggregation-induced processes as a new strategy towards advanced fluorescence markers for amyloid detection.

A Protein Misfolding Shaking Amplification-based method for the spontaneous generation of hundreds of bona fide prions.

Prion diseases are a group of rapidly progressing neurodegenerative disorders caused by the misfolding of the endogenous prion protein (PrPC) into a pathogenic form (PrPSc). This process, despite being the central event underlying these disorders, remains largely unknown at a molecular level, precluding the prediction of new potential outbreaks or interspecies transmission incidents. In this work, we present a method to generate bona fide recombinant prions de novo, allowing a comprehensive analysis of protein misfolding across a wide range of prion proteins from mammalian species. We study more than 380 different prion proteins from mammals and classify them according to their spontaneous misfolding propensity and their conformational variability. This study aims to address fundamental questions in the prion research field such as defining infectivity determinants, interspecies transmission barriers or the structural influence of specific amino acids and provide invaluable information for future diagnosis and therapy applications.