Accelerated Molecular Dynamics for Peptide Folding: Benchmarking Different Combinations of Force Fields and Explicit Solvent Models.

Accelerated molecular dynamics (aMD) protocols were assessed on predicting the secondary structure of eight peptides, of which two are helical, three are ß-hairpins, and three are disordered. Protocols consisted of combinations of three force fields (ff99SB, ff14SB, ff19SB) and two explicit solvation models (TIP3P and OPC), and were evaluated in two independent aMD simulations, one starting from an extended conformation, the other starting from a misfolded conformation. The results of these analyses indicate that all three combinations performed well on helical peptides. As for ß-hairpins, ff19SB performed well with both solvation methods, with a slight preference for the TIP3P solvation model, even though performance was dependent on both peptide sequence and initial conformation. The ff19SB/OPC combination had the best performance on intrinsically disordered peptides. In general, ff14SB/TIP3P suffered the strongest helical bias.

Novel Polymyxin-Inspired Peptidomimetics Targeting the SARS-CoV-2 Spike:hACE2 Interface.

Though the bulk of the COVID-19 pandemic is behind, the search for effective and safe anti-SARS-CoV-2 drugs continues to be relevant. A highly pursued approach for antiviral drug development involves targeting the viral spike (S) protein of SARS-CoV-2 to prevent its attachment to the cellular receptor ACE2. Here, we exploited the core structure of polymyxin B, a naturally occurring antibiotic, to design and synthesize unprecedented peptidomimetics (PMs), intended to target contemporarily two defined, non-overlapping regions of the S receptor-binding domain (RBD). Monomers 1, 2, and 8, and heterodimers 7 and 10 bound to the S-RBD with micromolar affinity in cell-free surface plasmon resonance assays (KD ranging from 2.31 µM to 2.78 µM for dimers and 8.56 µM to 10.12 µM for monomers). Although the PMs were not able to fully protect cell cultures from infection with authentic live SARS-CoV-2, dimer 10 exerted a minimal but detectable inhibition of SARS-CoV-2 entry in U87.ACE2+ and A549.ACE2.TMPRSS2+ cells. These results validated a previous modeling study and provided the first proof-of-feasibility of using medium-sized heterodimeric PMs for targeting the S-RBD. Thus, heterodimers 7 and 10 may serve as a lead for the development of optimized compounds, which are structurally related to polymyxin, with improved S-RBD affinity and anti-SARS-CoV-2 potential.

In Vitro Antimalarial Activity of Inhibitors of the Human GTPase Rac1.

Malaria accounts for millions of cases and thousands of deaths every year. In the absence of an effective vaccine, drugs are still the most important tool in the fight against the disease. Plasmodium parasites developed resistance to all classes of known antimalarial drugs. Thus, the search for antimalarial drugs with novel mechanisms of action is compelling. The human GTPase Rac1 plays a role in parasite invasion of the host cell in many intracellular pathogens. Also, in Plasmodium falciparum, the involvement of Rac1 during both the invasion process and parasite intracellular development was suggested. The aim of this work is to test a panel of Rac1 inhibitors as potential antimalarial drugs. Fourteen commercially available or newly synthesized inhibitors of Rac1 were tested for antimalarial activity. Among these, EHop-016 was the most effective against P. falciparum in vitro, with nanomolar 50% inhibitory concentrations (IC50s) (138.8 ± 16.0 nM on the chloroquine-sensitive D10 strain and 321.5 ± 28.5 nM on the chloroquine-resistant W2 strain) and a selectivity index of 37.8. EHop-016 did not inhibit parasite invasion of red blood cells but affected parasite growth inside them. Among the tested Rac1 inhibitors, EHop-016 showed promising activity that raises attention to this class of molecules as potential antimalarials and deserves further investigation.

Modular synthesis of 2,4-diaminoanilines as CNS drug-like non-covalent inhibitors of asparagine endopeptidase.

Asparagine endopeptidase (AEP), also called legumain, is a pH-dependent endolysosomal cysteine protease that cleaves its substrates after asparagine residues. Recent studies showed that it possesses δ-secretase activity and that it is implicated in numerous neurological diseases such as Alzheimer's disease (AD). Following evidence of aryl-morpholines as useful asparagine endopeptidase inhibitors, a series of morpholinoanilines with diverse substituents at ortho position were synthesized in view of improving the potency and scope of this molecular scaffold, allowing to identify ethyl 2-isonipecotate-4-morpholinoaniline possessing inhibition potency in the nanomolar range. CNS MPO (CNS MultiParameter Optimization) calculations revealed that most of the compounds developed in this work show physicochemical parameters in the desirable range for CNS drug candidates.

Morpholino-based peptide oligomers: Synthesis and DNA binding properties.

The chemical structure of oligonucleotide analogues dictates the conformation of oligonucleotide analogue oligomers, their ability to hybridize complementary DNA and RNA, their stability to degradation and their pharmacokinetic properties. In a study aimed at investigating new analogues featuring a neutral backbone, we explored the ability of oligomers containing a morpholino-peptide backbone to bind oligonucleotides. Circular Dichroism studies revealed the ability of our oligomers to interact with DNA, molecular modelling studies revealed the interaction responsible for complex stabilization.

Central residues of FSHß (89-97) peptide are not critical for FSHR binding: Implications for peptidomimetic design.

In our previous study, we had identified a 9-mer peptide (FSHß (89-97)) derived from seat belt loop of human FSHß and demonstrated its ability to function as FSHR antagonist in vivo. Structure analysis revealed that the four central residues 91STDC94 within this peptide may not be critical for receptor binding. In the present study, 91STDC94 residues were substituted with alanine to generate ΔFSHß 89-97(91STDC94/AAAA) peptide. Analogous to the parent peptide, ΔFSHß 89-97(91STDC94/AAAA) peptide inhibited binding of iodinated FSH to rat FSHR and reduced FSH-induced cAMP production. The peptide could impede granulosa cell proliferation leading to reduction in FSH-mediated ovarian weight gain in immature female rats. In these rats, peptide administration further downregulated androgen receptor and estrogen receptor-alpha expression and upregulated estrogen receptor-beta expression. The results indicate that substitution of 91STDC94 with alanine did not significantly alter FSHR antagonist activity of FSHß (89-97) peptide implying that these residues are not critical for FSH-FSHR interaction and can be replaced with non-peptidic moieties for development of more potent peptidomimetics.

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.

Discovery of a d-pro-lys peptidomimetic inhibitor of MMP9: Addressing the gelatinase selectivity beyond S1' subsite.

Despite a high degree of structural similarity, it is known that MMP2 and MMP9 have distinct roles in the angiogenic switch and in cell migration, as they activate diverse signaling pathways. Indeed, inhibition of MMP2 and MMP9 can show beneficial or detrimental effects depending on the stage of tumor progression. Thus, the selective inhibition of gelatinases is of relevance for a successful drug lead, which has to be achieved despite the high structural similarity of the two gelatinases. Herein, the synthesis and evaluation of d-proline-derived hydroxamic acids containing amino appendages at C-4 as gelatinase inhibitors are reported. Inhibition assays enabled the identification of a > 200-fold selective MMP9 inhibitor when Lys was considered as a C-4 substituent, thus addressing gelatinase selectivity beyond the S1' subsite, which is a major driver for selectivity. Molecular docking studies revealed the basic moiety of Lys as detrimental for inhibition of MMP2 as compared to MMP9.

Rescoring Virtual Screening Results with the MM-PBSA Methods: Beware of Internal Dielectric Constants.

With the potential of improving virtual screening outcome, MM-PB/GBSA has become a disputed method that requires extensive testing and tuning to provide the optimal results. One of the tuning factors is the internal or solute dielectric constant. We have applied three test sets with receptors of different categories and libraries from different sources to investigate the underlying issue related to this constant. We discovered that increasing internal dielectric value does not improve the virtual screening enrichment qualitatively. More interestingly, nonpolar and polar calculated energies act differently in libraries with different molecular weight distributions. From this work, the performance of MM-PBSA rescoring in virtual screening is more library- than receptor-dependent.

Identification of highly potent and selective MMP2 inhibitors addressing the S1' subsite with d-proline-based compounds.

MMP2 and MMP9, also called gelatinases, play a primary role in the angiogenic switch, as a fundamental step of tumor progression, and show high degree of structural similarity. Clinically successful gelatinase inhibitors need to be highly selective as opposite effects have been found for the two enzymes, and the S1' subsite is the major driver to attain selective and potent inhibitors. The synthesis of d-proline-derived hydroxamic acids containing diverse appendages at the amino group, varying in length and decoration allowed to give insight on the MMP2/MMP9 selectivity around the S1' subsite, resulting in the identification of sub-nanomolar compounds with high selectivity up to 730. Molecular docking studies revealed the existence of an additional hydrophobic channel at the bottom of S1' subsite for MMP2 enzyme useful to drive selectivity towards such gelatinase.

Tetrahydro-4 H-(pyrrolo[3,4- d]isoxazol-3-yl)methanamine: A Bicyclic Diamino Scaffold Stabilizing Parallel Turn Conformations.

A tetrahydro-4 H-(pyrrolo[3,4- d]isoxazol-3-yl)methanamine scaffold was designed as a diamino derivative to stabilize parallel turn conformations. Its synthesis took advantage of a [1,3]-dipolar cycloaddition reaction between the nitrile oxide derived from the inexpensive enantiopure l -phenylalanine and N-benzyl-3-pyrroline. Two diastereoisomers were formed, whose distribution depends on the selected base. 3a R,6a S-Isomer is favored in organic bases, which formation is driven by π-interactions. However, the above interactions were significantly prevented using an inorganic base due to the chaotropic effect of the cation, decreasing the amount of the above isomer. Finally, we demonstrated that this isomer is able to stabilize parallel turn conformations when inserted in short peptide sequences.

Stereodivergent synthesis of 5-aminopipecolic acids and application in the preparation of a cyclic RGD peptidomimetic as a nanomolar αVß3 integrin ligand.

A stereodivergent strategy was devised to obtain enantiopure cis and trans 5-aminopipecolic acids (5-APAs) in suitably protected forms to be employed in peptide synthesis as conformationally constrained α- and δ-amino acids. The cis isomer was used as a δ-amino acid to construct a cyclic RGD-containing peptidomimetic, the ability of which to compete with biotinylated vitronectin for binding with the isolated αVß3 integrin was measured (IC50 = 4.2 ± 0.9 nM). A complete 1H NMR and computational conformational analysis was performed to elucidate the reasons for the high affinity of this cyclic peptidomimetic in comparison with cilengitide.

Tandem Tetrahydroisoquinoline-4-carboxylic Acid/ß-Alanine as a New Construct Able To Induce a Flexible Turn.

Tetrahydroisoquinoline-4-carboxylic acid, a constrained ß2 -amino acid named ß-TIC, was synthesised for the first time in enantiopure form. The biocatalytic route applied herein represents one of the few successful examples of enzymatic resolution of ß2 -amino acids. Model tetrapeptides, namely, Fmoc-l-Ala-ß-TIC-ß-Ala-l-Val-OBn (Fmoc=fluorenylmethyloxycarbonyl, Bn=benzyl), containing both isomers of ß-TIC, were prepared. Both computational and NMR spectroscopy studies were performed. A reverse-turn conformation was observed in the case of (R)-ß-TIC enantiomer that was obtained in 99 % enantiomeric excess by enzymatic resolution. The ß-TIC/ß-Ala construct represents the first example of a flexible turn mimetic containing a cyclic and an acyclic ß-amino acid. Furthermore, the presence of an aromatic ring of ß-TIC could facilitate non-covalent interactions to increase the potential of this scaffold for the preparation of protein-protein interaction modulators.

Peptide modulators of Rac1/Tiam1 protein-protein interaction: An alternative approach for cardiovascular diseases.

Rac1 GTPase interaction with guanine nucleotide exchange factor Tiam1 is involved in several cancer types and cardiovascular diseases. Although small molecules interfering with their protein-protein interaction (PPI) were identified and studied, the ability of small peptides and peptide mimics acting as Rac1/Tiam1 PPI inhibitors has not been yet explored. Using computational alanine scanning (CAS), the "hot" interfacial residues have been determined allowing the design of a small library of putative PPI inhibitors. In particular, the insertion of an unnatural alpha, alpha disubstituted amino acid, that is norbornane amino acid, and the side chain stapling have been evaluated regarding both conformational stability and biological activity. REMD calculations and CD studies have indicated that one single norbornane amino acid at the N-terminus is not sufficient to stabilize the helix structure, while the side-chain stapling is a more efficient strategy. Furthermore, both engineered peptides have been found able to reduce Rac1-GTP levels in cultured human smooth muscle cells, while wild type sequence is not active.

Synthesis and conformational analysis of peptides embodying 2,3-methanopipecolic acids.

The conformational analysis of linear and cyclic peptides incorporating 2,3-methanopipecolic acids (or Cyclopropane Pipecolic Acids, CPAs) as conformationally constrained α-amino acids is reported. Compared to peptides containing proline or pipecolic acid, a striking increase of the cis isomer (42-92%) around the CPA amide bond is observed, both in water and organic solvents, when these unnatural amino acids are embodied in linear amino acid sequences. The rotational barrier around the same bond in water was calculated, giving results comparable to that for the prolyl cis/trans isomerization. In organic solvents, CPAs at the i + 2 position of a peptide induce the formation of a type VIa ß-turn secondary structure. When incorporated into a cyclic peptide, the cis geometry around the 2,3-methanopipecolic amide bond still prevails and, in the example studied herein (a cyclic RGD-containing ligand of αVß3 integrin mimicking Cilengitide), conservation of the backbone geometry and side chain spatial orientation of the native peptide is also found. Given the importance of the proline cis/trans isomerism in many biological processes, CPAs could be useful as proline mimetics for probing protein-ligand interactions and generating novel bioactive compounds.

Synthesis and Biological Evaluation of New Natural Phenolic (2E,4E,6E)-Octa-2,4,6-trienoic Esters.

In the present study the esterification of the OH groups of resveratrol, caffeic acid, ferulic acid, and ß-sitosterol with an antioxidant polyconjugated fatty acid, (2E,4E,6E)-octa-2,4,6-trienoic acid, was achieved. As the selective esterification of OH groups of natural compounds can affect their biological activity, a selective esterification of resveratrol and caffeic acid was performed by an enzymatic approach. The new resulting compounds were characterized spectroscopically (FT-IR, NMR mono, and bidimensional techniques); when necessary the experimental data were integrated by quantum chemical calculations. The antioxidant, anti-inflammatory and proliferative activity was evaluated. The good results encourage the use of these molecules as antioxidant and/or anti-inflammatory agents in dermocosmetic application.

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.

Molecular insights into dimerization inhibition of c-Maf transcription factor.

The Maf protein family belongs to the activator protein 1 (AP-1) superfamily of transcription factors that bind specific DNA target sequences through a basic region and exploit a leucine zipper (LZ) motif for protein-protein interactions leading to homo- or hetero-dimerization. Mafs unique DNA-binding domain contains a highly conserved extended homology region (EHR) that allows to recognize longer DNA sequences than other basic leucine zipper (bZIP) transcription factors. Inspired by the fact that overexpression of Mafs is observed in about 50% of cases of multiple myeloma, a hematological malignant disorder, we undertook a peptide inhibitor approach. The LZ domain of c-Maf, one of large Mafs, was produced by solid phase peptide synthesis. We characterized its secondary structure and dimerization properties, and found that dimerization and folding events are strictly coupled. Moreover, potential peptidic c-Maf dimerization inhibitors were computationally designed and synthesized. These compounds were demonstrated by circular dichroism (CD) spectroscopy and MALDI-TOF mass spectrometry to bind to c-Maf LZ monomers, to drive folding of their partially disordered structure and to efficiently compete with dimerization, suggesting a way for interfering with the function of c-Maf and, more generally, of intrinsically disordered proteins, till now considered undruggable targets.

1H-Azepine-2-oxo-5-amino-5-carboxylic Acid: A 310 Helix Inducer and an Effective Tool for Functionalized Gold Nanoparticles.

A new α,α-disubstituted constrained glutamine analogue has been designed to decorate gold nanoparticles and to induce a 310-helix when inserted in peptides. Using an efficient "one-pot" asymmetric Schmidt reaction between 4-disubstituted-cyclohexanone and hydroxyalkylazides, 1H-azepine-2-oxo-5-amino-5-carboxylic acid was prepared. The main (R) isomer was inserted at the N-terminus in a very short peptide sequence (i.e., PhCO-(R)-Oxo-Azn-L-Ala-Aib-L-AlaNHMe) and a stable 310-helix conformation was obtained, as verified by both NMR experiments and molecular dynamics (MD) simulations. Finally, the presence of the hydroxyl chain at the nitrogen atom of the ring allowed for the preparation of covered chiral gold nanoparticles.