AptaMat: a matrix-based algorithm to compare single-stranded oligonucleotides secondary structures.

MOTIVATION: Comparing single-stranded nucleic acids (ssNAs) secondary structures is fundamental when investigating their function and evolution and predicting the effect of mutations on their structures. Many comparison metrics exist, although they are either too elaborate or not sensitive enough to distinguish close ssNAs structures. RESULTS: In this context, we developed AptaMat, a simple and sensitive algorithm for ssNAs secondary structures comparison based on matrices representing the ssNAs secondary structures and a metric built upon the Manhattan distance in the plane. We applied AptaMat to several examples and compared the results to those obtained by the most frequently used metrics, namely the Hamming distance and the RNAdistance, and by a recently developed image-based approach. We showed that AptaMat is able to discriminate between similar sequences, outperforming all the other here considered metrics. In addition, we showed that AptaMat was able to correctly classify 14 RFAM families within a clustering procedure. AVAILABILITY AND IMPLEMENTATION: The python code for AptaMat is available at https://github.com/GEC-git/AptaMat.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparative Study of Single-stranded Oligonucleotides Secondary Structure Prediction Tools.

BACKGROUND: Single-stranded nucleic acids (ssNAs) have important biological roles and a high biotechnological potential linked to their ability to bind to numerous molecular targets. This depends on the different spatial conformations they can assume. The first level of ssNAs spatial organisation corresponds to their base pairs pattern, i.e. their secondary structure. Many computational tools have been developed to predict the ssNAs secondary structures, making the choice of the appropriate tool difficult, and an up-to-date guide on the limits and applicability of current secondary structure prediction tools is missing. Therefore, we performed a comparative study of the performances of 9 freely available tools (mfold, RNAfold, CentroidFold, CONTRAfold, MC-Fold, LinearFold, UFold, SPOT-RNA, and MXfold2) on a dataset of 538 ssNAs with known experimental secondary structure. RESULTS: The minimum free energy-based tools, namely mfold and RNAfold, and some tools based on artificial intelligence, namely CONTRAfold and MXfold2, provided the best results, with [Formula: see text] of exact predictions, whilst MC-fold seemed to be the worst performing tool, with only [Formula: see text] of exact predictions. In addition, UFold and SPOT-RNA are the only options for pseudoknots prediction. Including in the analysis of mfold and RNAfold results 5-10 suboptimal solutions further improved the performances of these tools. Nevertheless, we could observe issues in predicting particular motifs, such as multiple-ways junctions and mini-dumbbells, or the ssNAs whose structure has been determined in complex with a protein. In addition, our benchmark shows that some effort has to be paid for ssDNA secondary structure predictions. CONCLUSIONS: In general, Mfold, RNAfold, and MXfold2 seem to currently be the best choice for the ssNAs secondary structure prediction, although they still show some limits linked to specific structural motifs. Nevertheless, actual trends suggest that artificial intelligence has a high potential to overcome these remaining issues, for example the recently developed UFold and SPOT-RNA have a high success rate in predicting pseudoknots.

Lyme borreliosis diagnosis: state of the art of improvements and innovations.

With almost 700 000 estimated cases each year in the United States and Europe, Lyme borreliosis (LB), also called Lyme disease, is the most common tick-borne illness in the world. Transmitted by ticks of the genus Ixodes and caused by bacteria Borrelia burgdorferi sensu lato, LB occurs with various symptoms, such as erythema migrans, which is characteristic, whereas others involve blurred clinical features such as fatigue, headaches, arthralgia, and myalgia. The diagnosis of Lyme borreliosis, based on a standard two-tiered serology, is the subject of many debates and controversies, since it relies on an indirect approach which suffers from a low sensitivity depending on the stage of the disease. Above all, early detection of the disease raises some issues. Inappropriate diagnosis of Lyme borreliosis leads to therapeutic wandering, inducing potential chronic infection with a strong antibody response that fails to clear the infection. Early and proper detection of Lyme disease is essential to propose an adequate treatment to patients and avoid the persistence of the pathogen. This review presents the available tests, with an emphasis on the improvements of the current diagnosis, the innovative methods and ideas which, ultimately, will allow more precise detection of LB.

Synthetic Peptide Antibodies as TNF-α Inhibitors: Molecularly Imprinted Polymer Nanogels Neutralize the Inflammatory Activity of TNF-α in THP-1 Derived Macrophages.

Tumor Necrosis Factor-α (TNF-α) is a cytokine that is normally produced by immune cells when fighting an infection. But, when too much TNF-α is produced as in autoimmune diseases, this leads to unwanted and persistent inflammation. Anti-TNF-α monoclonal antibodies have revolutionized the therapy of these disorders by blocking TNF-α and preventing its binding to TNF-α receptors, thus suppressing the inflammation. Herein, we propose an alternative in the form of molecularly imprinted polymer nanogels (MIP-NGs). MIP-NGs are synthetic antibodies obtained by nanomoulding the 3-dimensional shape and chemical functionalities of a desired target in a synthetic polymer. Using an in-house developed in silico rational approach, epitope peptides of TNF-α were generated and 'synthetic peptide antibodies' were prepared. The resultant MIP-NGs bind the template peptide and recombinant TNF-α with high affinity and selectivity, and can block the binding of TNF-α to its receptor. Consequently they were applied to neutralize pro-inflammatory TNF-α in the supernatant of human THP-1 macrophages, leading to a downregulation of the secretion of pro-inflammatory cytokines. Our results suggest that MIP-NGs, which are thermally and biochemically more stable and easier to manufacture than antibodies, and cost-effective, are very promising as next generation TNF-α inhibitors for the treatment of inflammatory diseases.

Molecularly Imprinted Polymer Nanogels for Protein Recognition: Direct Proof of Specific Binding Sites by Solution STD and WaterLOGSY NMR Spectroscopies.

Molecularly imprinted polymers (MIPs) are tailor-made synthetic antibodies possessing specific binding cavities designed for a target molecule. Currently, MIPs for protein targets are synthesized by imprinting a short surface-exposed fragment of the protein, called epitope or antigenic determinant. However, finding the epitope par excellence that will yield a peptide "synthetic antibody" cross-reacting exclusively with the protein from which it is derived, is not easy. We propose a computer-based rational approach to unambiguously identify the "best" epitope candidate. Then, using Saturation Transfer Difference (STD) and WaterLOGSY NMR spectroscopies, we prove the existence of specific binding sites created by the imprinting of this peptide epitope in the MIP nanogel. The optimized MIP nanogel could bind the epitope and cognate protein with a high affinity and selectivity. The study was performed on Hepatitis A Virus Cell Receptor-1 protein, also known as KIM-1 and TIM-1, for its ubiquitous implication in numerous pathologies.

In Silico Drug Repurposing for SARS-CoV-2 Main Proteinase and Spike Proteins.

The pandemic caused by SARS-CoV-2 is currently representing a major health and economic threat to humanity. So far, no specific treatment to this viral infection has been developed and the emergency still requires an efficient intervention. In this work, we used virtual screening to facilitate drug repurposing against SARS-CoV-2, targeting viral main proteinase and spike protein with 3000 existing drugs. We used a protocol based on a docking step followed by a short molecular dynamic simulation and rescoring by the Nwat-MMGBSA approach. Our results provide suggestions for prioritizing in vitro and/or in vivo tests of already available compounds.

Thermal Adaptation of Enzymes: Impacts of Conformational Shifts on Catalytic Activation Energy and Optimum Temperature.

Thermal adaptation of enzymes is essential for both living organism development in extreme conditions and efficient biocatalytic applications. However, the molecular mechanisms leading to a shift in catalytic activity optimum temperatures remain unclear, and there is increasing experimental evidence that thermal adaptation involves complex changes in both structural and reactive properties. Here, a combination of enhanced protein conformational sampling with an explicit chemical reaction description was applied to mesophilic and thermophilic homologues of the dihydrofolate reductase enzyme, and a quantitative description of the stability and catalytic activity shifts between homologues was obtained. The key role played by temperature-induced shifts in protein conformational distributions is revealed. In contrast with pictures focusing on protein flexibility and dynamics, it is shown that while the homologues' reaction free energies are similar, the striking discrepancy between their activation energies is caused by their different conformational changes with temperature. An analytic model is proposed that combines catalytic activity and structural stability, and which quantitatively predicts the shift in homologue optimum temperatures. It is shown that this general model provides a molecular explanation of changes in optimum temperatures for several other enzymes.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes.

A key aspect of life's evolution on Earth is the adaptation of proteins to be stable and work in a very wide range of temperature conditions. A detailed understanding of the associated molecular mechanisms would also help to design enzymes optimized for biotechnological processes. Despite important advances, a comprehensive picture of how thermophilic enzymes succeed in functioning under extreme temperatures remains incomplete. Here, we examine the temperature dependence of stability and of flexibility in the mesophilic monomeric Escherichia coli (Ec) and thermophilic dimeric Thermotoga maritima (Tm) homologs of the paradigm dihydrofolate reductase (DHFR) enzyme. We use all-atom molecular dynamics simulations and a replica-exchange scheme that allows to enhance the conformational sampling while providing at the same time a detailed understanding of the enzymes' behavior at increasing temperatures. We show that this approach reproduces the stability shift between the two homologs, and provides a molecular description of the denaturation mechanism by identifying the sequence of secondary structure elements melting as temperature increases, which is not straightforwardly obtained in the experiments. By repeating our approach on the hypothetical TmDHFR monomer, we further determine the respective effects of sequence and oligomerization in the exceptional stability of TmDFHR. We show that the intuitive expectation that protein flexibility and thermal stability are correlated is not verified. Finally, our simulations reveal that significant conformational fluctuations already take place much below the melting temperature. While the difference between the TmDHFR and EcDHFR catalytic activities is often interpreted via a simplified two-state picture involving the open and closed conformations of the key M20 loop, our simulations suggest that the two homologs' markedly different activity temperature dependences are caused by changes in the ligand-cofactor distance distributions in response to these conformational changes.

Improved Computation of Protein-Protein Relative Binding Energies with the Nwat-MMGBSA Method.

A MMGBSA variant (here referred to as Nwat-MMGBSA), based on the inclusion of a certain number of explicit water molecules (Nwat) during the calculations, has been tested on a set of 20 protein-protein complexes, using the correlation between predicted and experimental binding energy as the evaluation metric. Besides the Nwat parameter, the effect of the force field, the molecular dynamics simulation length, and the implicit solvent model used in the MMGBSA analysis have been also evaluated. We found that considering 30 interfacial water molecules improved the correlation between predicted and experimental binding energies by up to 30%, compared to the standard approach. Moreover, the correlation resulted in being rather sensitive to the force field and, to a minor extent, to the implicit solvent model and to the length of the MD simulation.

Biocatalytic Dynamic Kinetic Resolution for the Synthesis of Atropisomeric Biaryl N-Oxide Lewis Base Catalysts.

Atropisomeric biaryl pyridine and isoquinoline N-oxides were synthesized enantioselectively by dynamic kinetic resolution (DKR) of rapidly racemizing precursors exhibiting free bond rotation. The DKR was achieved by ketoreductase (KRED) catalyzed reduction of an aldehyde to form a configurationally stable atropisomeric alcohol, with the substantial increase in rotational barrier arising from the loss of a bonding interaction between the N-oxide and the aldehyde. Use of different KREDs allowed either the M or P enantiomer to be synthesized in excellent enantiopurity. The enantioenriched biaryl N-oxide compounds catalyze the asymmetric allylation of benzaldehyde derivatives with allyltrichlorosilane.

The functionality of a therapeutic antibody candidate restored by a single mutation from proline to threonine in the variable region.

mAbs play an essential role in the therapeutic arsenal. Our laboratory has patented the Rendomab-B49 mAb targeting the endothelin B receptor (ETB). This G protein-coupled receptor plays a driving role in the progression of numerous cancers. We chimerized our mAb (xiRB49) to evaluate its preclinical therapeutic efficacy in different ETB+ tumor models with an antibody drug conjugate approach. As previously reported, the chimerization process of an antibody can alter its functionality. In this article, we present the chimerization of RB49. xiRB49 purified by Protein A remained perfectly soluble and did not aggregate, but it lost all its ability to recognize ETB. A detailed analysis of its variable region using IMGT tools allowed us to identify an unusual proline at position 125. In silico mAb modeling and in vitro experiments were performed for a better understanding of xiRB49 structure-function relationships. Our results show that the proline in position 125 on the heavy chain alters the xiRB49 CDR3 light chain conformation and its mutation to threonine allows complete functional recovery.

Molecularly imprinted polymer nanogels targeting the HAV motif in cadherins inhibit cell-cell adhesion and migration.

Cadherins are cell-surface proteins that mediate cell-cell adhesion. By regulating their grip formation and strength, cadherins play a pivotal role during normal tissue morphogenesis and homeostasis of multicellular organisms. However, their dysfunction is associated with cell migration and proliferation, cancer progression and metastasis. The conserved amino acid sequence His-Ala-Val (HAV) in the extracellular domain of cadherins is implicated in cadherin-mediated adhesion and migration. Antagonists of cadherin adhesion such as monoclonal antibodies and small molecule inhibitors based on HAV peptides, are of high therapeutic value in cancer treatment. However, antibodies are not stable outside their natural environment and are expensive to produce, while peptides have certain limitations as a drug as they are prone to proteolysis. Herein, we propose as alternative, a synthetic antibody based on molecularly imprinted polymer nanogels (MIP-NGs) to target the HAV domain. The MIP-NGs are biocompatible, have high affinity for N-cadherin and inhibit cell adhesion and migration of human cervical adenocarcinoma (HeLa) cells, as demonstrated by cell aggregation and Matrigel invasion assays, respectively. The emergence of MIPs as therapeutics for fighting cancer is still in its infancy and this novel demonstration reinforces the fact that they have a rightful place in cancer treatment.

Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations.

BACKGROUND: The inclusion of direct effects mediated by water during the ligandreceptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions. OBJECTIVE: In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered. RESULTS: Here, we analyse the software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules. CONCLUSIONS: Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon.

Fluoro-Aryl Substituted α,ß2,3-Peptides in the Development of Foldameric Antiparallel ß-Sheets: A Conformational Study.

α,ß2,3-Disteroisomeric foldamers of general formula Boc(S-Ala-ß-2R,3R-Fpg)nOMe or Boc(S-Ala-ß-2S,3S-Fpg)nOMe were prepared from both enantiomers of syn H-2-(2-F-Phe)-h-PheGly-OH (named ß-Fpg) and S-alanine. Our peptides show two appealing features for biomedical applications: the presence of fluorine, attractive for non-covalent interactions, and aryl groups, crucial for π-stacking. A conformational study was performed, using IR, NMR and computational studies of diastereoisomeric tetra- and hexapeptides containing the ß2,3-amino acid in the R,R- and S,S-stereochemistry, respectively. We found that the stability of peptide conformation is dependent on the stereochemistry of the ß-amino acid. Combining S-Ala with ß-2R,3R-Fpg, a stable extended ß-strand conformation was obtained. Furthermore, ß-2R,3R-Fpg containing hexapeptide self-assembles to form antiparallel ß-sheet structure stabilized by intermolecular H-bonds and π,π-interactions. These features make peptides containing the ß2,3-fluoro amino acid very appealing for the development of bioactive proteolytically stable foldameric ß-sheets as modulators of protein-protein interaction (PPI).

Genetic background and immunological status influence B cell repertoire diversity in mice.

The relationship between the immune repertoire and the physiopathological status of individuals is essential to apprehend the genesis and the evolution of numerous pathologies. Nevertheless, the methodological approaches to understand these complex interactions are challenging. We performed a study evaluating the diversity harbored by different immune repertoires as a function of their physiopathological status. In this study, we base our analysis on a murine scFv library previously described and representing four different immune repertoires: i) healthy and naïve, ii) healthy and immunized, iii) autoimmune prone and naïve, and iv) autoimmune prone and immunized. This library, 2.6 × 109 in size, is submitted to high throughput sequencing (Next Generation Sequencing, NGS) in order to analyze the gene subgroups encoding for immunoglobulins. A comparative study of the distribution of immunoglobulin gene subgroups present in the four libraries has revealed shifts in the B cell repertoire originating from differences in genetic background and immunological status of mice.

An Efficient Implementation of the Nwat-MMGBSA Method to Rescore Docking Results in Medium-Throughput Virtual Screenings.

Nwat-MMGBSA is a variant of MM-PB/GBSA based on the inclusion of a number of explicit water molecules that are the closest to the ligand in each frame of a molecular dynamics trajectory. This method demonstrated improved correlations between calculated and experimental binding energies in both protein-protein interactions and ligand-receptor complexes, in comparison to the standard MM-GBSA. A protocol optimization, aimed to maximize efficacy and efficiency, is discussed here considering penicillopepsin, HIV1-protease, and BCL-XL as test cases. Calculations were performed in triplicates on both classic HPC environments and on standard workstations equipped by a GPU card, evidencing no statistical differences in the results. No relevant differences in correlation to experiments were also observed when performing Nwat-MMGBSA calculations on 4 or 1 ns long trajectories. A fully automatic workflow for structure-based virtual screening, performing from library set-up to docking and Nwat-MMGBSA rescoring, has then been developed. The protocol has been tested against no rescoring or standard MM-GBSA rescoring within a retrospective virtual screening of inhibitors of AmpC ß-lactamase and of the Rac1-Tiam1 protein-protein interaction. In both cases, Nwat-MMGBSA rescoring provided a statistically significant increase in the ROC AUCs of between 20 and 30%, compared to docking scoring or to standard MM-GBSA rescoring.

A tendril perversion in a helical oligomer: trapping and characterizing a mobile screw-sense reversal.

Helical oligomers of achiral monomers adopt domains of uniform screw sense, which are occasionally interrupted by screw-sense reversals. These rare, elusive, and fast-moving features have eluded detailed characterization. We now describe the structure and habits of a screw-sense reversal trapped within a fragment of a helical oligoamide foldamer of the achiral quaternary amino acid 2-aminoisobutyric acid (Aib). The reversal was enforced by compelling the amide oligomer to adopt a right-handed screw sense at one end and a left-handed screw sense at the other. The trapped reversal was characterized by X-ray crystallography, and its dynamic properties were monitored by NMR and circular dichroism, and modelled computationally. Raman spectroscopy indicated that a predominantly helical architecture was maintained despite the reversal. NMR and computational results indicated a stepwise shift from one screw sense to another on moving along the helical chain, indicating that in solution the reversal is not localised at a specific location, but is free to migrate across a number of residues. Analogous unconstrained screw-sense reversals that are free to move within a helical structure are likely to provide the mechanism by which comparable helical polymers and foldamers undergo screw-sense inversion.

An Updated Test of AMBER Force Fields and Implicit Solvent Models in Predicting the Secondary Structure of Helical, ß-Hairpin, and Intrinsically Disordered Peptides.

Replica exchange molecular dynamics simulations were performed to test the ability of six AMBER force fields and three implicit solvent models of predicting the native conformation of two helical peptides, three ß-hairpins, and three intrinsically disordered peptides. Although a combination of the force field and implicit solvation models able to accurately predict the native structure of all the considered peptides was not identified, we found that the GB-Neck2 model seems to well compensate for some of the conformational biases showed by ff96 and ff99SB/ildn/ildn-φ. Indeed, the force fields of the ff99SB series coupled with GB-Neck2 reasonably discriminated helices from disordered peptides, while a good prediction of ß-hairpin conformations was only achieved by performing two independent simulations: one with the ff96/GB-Neck2 combination and the other with GB-Neck2 coupled with any of the ff99SB/ildn/ildn-φ force fields.

Origin of Helical Screw Sense Selectivity Induced by Chiral Constrained Cα-Tetrasubstituted α-Amino Acids in Aib-based Peptides.

The mechanisms behind the propensity of chiral constrained Cα-tetrasubstituted amino acids (cCTAAs) to induce one particular helical screw sense, when included in the Ac-Aib2-cCTAA-Aib2-NHMe peptide model, were studied through replica exchange molecular dynamics, potential of mean force, and quantum theory of atoms in molecules calculations. We observed that cCTAAs exert their effect on helical screw sense selectivity through the positioning of the side chain to generate steric hindrance in either the (-x, +y, +z) or (+x, +y, -z) sectors of a right-handed 3D Cartesian space, where the z axis corresponds to the axis of the helix and the Cα lies on the +y semiaxis (0, +y, 0). The different strengthening of the noncovalent interactions, also comprising C-H···O interactions, exerted by the cCTAA in the right-handed or left-handed helix was also found important to define the preference of a cCTAA for a particular helix screw sense.

Mechanism of stabilization of helix secondary structure by constrained Cα-tetrasubstituted α-amino acids.

The theoretical basis behind the ability of constrained Cα-tetrasubstituted amino acids (CTAAs) to induce stable helical conformations has been studied through Replica Exchange Molecular Dynamics Potential of Mean Force Quantum Theory of Atoms In Molecules calculations on Ac-l-Ala-CTAA-l-Ala-Aib-l-Ala-NHMe peptide models. We found that the origin of helix stabilization by CTAAs can be ascribed to at least two complementary mechanisms limiting the backbone conformational freedom: steric hindrance predominantly in the (+x,+y,-z) sector of a right-handed 3D Cartesian space, where the z axis coincides with the helical axis and the Cα of the CTAA lies on the +y axis (0,+y,0), and the establishment of additional and relatively strong C-H···O interactions involving the CTAA.