Locuaz: an in silico platform for protein binders optimization.

MOTIVATION: Engineering high-affinity binders targeting specific antigenic determinants remains a challenging and often daunting task, requiring extensive experimental screening. Computational methods have the potential to accelerate this process, reducing costs and time, but only if they demonstrate broad applicability and efficiency in exploring mutations, evaluating affinity, and pruning unproductive mutation paths. RESULTS: In response to these challenges, we introduce a new computational platform for optimizing protein binders towards their targets. The platform is organized as a series of modules, performing mutation selection and application, molecular dynamics simulations to sample conformations around interaction poses, and mutation prioritization using suitable scoring functions. Notably, the platform supports parallel exploration of different mutation streams, enabling in silico high-throughput screening on High Performance Computing (HPC) systems. Furthermore, the platform is highly customizable, allowing users to implement their own protocols. AVAILABILITY AND IMPLEMENTATION: The source code is available at https://github.com/pgbarletta/locuaz and documentation is at https://locuaz.readthedocs.io/. The data underlying this article are available at https://github.com/pgbarletta/suppl_info_locuaz.

Design, Synthesis, and Cytotoxic Assessment of New Haloperidol Analogues as Potential Anticancer Compounds Targeting Sigma Receptors.

Sigma receptors (SRs), including SR1 and SR2 subtypes, have attracted increasing interest in recent years due to their involvement in a wide range of activities, including the modulation of opioid analgesia, neuroprotection, and potential anticancer activity. In this context, haloperidol (HAL), a commonly used antipsychotic drug, also possesses SR activity and cytotoxic effects. Herein, we describe the identification of novel SR ligands, obtained by a chemical hybridization approach. There wereendowed with pan-affinity for both SR subtypes and evaluated their potential anticancer activity against SH-SY5Y and HUH-7 cancer cell lines. Through a chemical hybridization approach, we identified novel compounds (4d, 4e, 4g, and 4j) with dual affinity for SR1 and SR2 receptors. These compounds were subjected to cytotoxicity testing using a resazurin assay. The results revealed potent cytotoxic effects against both cancer cell lines, with IC50 values comparable to HAL. Interestingly, the cytotoxic potency of the novel compounds resembled that of the SR1 antagonist HAL rather than the SR2 agonist siramesine (SRM), indicating the potential role of SR1 antagonism in their mechanism of action. The further exploration of their structure-activity relationships and their evaluation in additional cancer cell lines will elucidate their therapeutic potential and may pave the way for the development of novel anticancer agents that target SRs.

De novo KCNA6 variants with attenuated KV 1.6 channel deactivation in patients with epilepsy.

OBJECTIVE: Mutations in the genes encoding neuronal ion channels are a common cause of Mendelian neurological diseases. We sought to identify novel de novo sequence variants in cases with early infantile epileptic phenotypes and neurodevelopmental anomalies. METHODS: Following clinical diagnosis, we performed whole exome sequencing of the index cases and their parents. Identified channel variants were expressed in Xenopus oocytes and their functional properties assessed using two-electrode voltage clamp. RESULTS: We identified novel de novo variants in KCNA6 in four unrelated individuals variably affected with neurodevelopmental disorders and seizures with onset in the first year of life. Three of the four identified mutations affect the pore-lining S6 α-helix of KV 1.6. A prominent finding of functional characterization in Xenopus oocytes was that the channel variants showed only minor effects on channel activation but slowed channel closure and shifted the voltage dependence of deactivation in a hyperpolarizing direction. Channels with a mutation affecting the S6 helix display dominant effects on channel deactivation when co-expressed with wild-type KV 1.6 or KV 1.1 subunits. SIGNIFICANCE: This is the first report of de novo nonsynonymous variants in KCNA6 associated with neurological or any clinical features. Channel variants showed a consistent effect on channel deactivation, slowing the rate of channel closure following normal activation. This specific gain-of-function feature is likely to underlie the neurological phenotype in our patients. Our data highlight KCNA6 as a novel channelopathy gene associated with early infantile epileptic phenotypes and neurodevelopmental anomalies.

Cytotoxicity Profiles and Neuroprotective Properties of the Novel Ifenprodil Analogues as Sigma Ligands.

Neurodegeneration is a slow and progressive loss of neuronal cells or their function in specific regions of the brain or in the peripheral system. Among several causes responsible for the most common neurodegenerative diseases (NDDs), cholinergic/dopaminergic pathways, but also some endogenous receptors, are often involved. In this context, sigma 1 receptor (S1R) modulators can be used as neuroprotective and antiamnesic agents. Herein, we describe the identification of novel S1R ligands endowed with antioxidant properties, potentially useful as neuroprotective agents. We also computationally assessed how the most promising compounds might interact with the S1R protein's binding sites. The in silico predicted ADME properties suggested that they could be able to cross the brain-blood-barrier (BBB), and to reach the targets. Finally, the observation that at least two novel ifenprodil analogues (5d and 5i) induce an increase of the mRNA levels of the antioxidant NRF2 and SOD1 genes in SH-SY5Y cells suggests that they might be effective agents for protecting neurons against oxidative damage.

Mutations in the Neuronal Vesicular SNARE VAMP2 Affect Synaptic Membrane Fusion and Impair Human Neurodevelopment.

VAMP2 encodes the vesicular SNARE protein VAMP2 (also called synaptobrevin-2). Together with its partners syntaxin-1A and synaptosomal-associated protein 25 (SNAP25), VAMP2 mediates fusion of synaptic vesicles to release neurotransmitters. VAMP2 is essential for vesicular exocytosis and activity-dependent neurotransmitter release. Here, we report five heterozygous de novo mutations in VAMP2 in unrelated individuals presenting with a neurodevelopmental disorder characterized by axial hypotonia (which had been present since birth), intellectual disability, and autistic features. In total, we identified two single-amino-acid deletions and three non-synonymous variants affecting conserved residues within the C terminus of the VAMP2 SNARE motif. Affected individuals carrying de novo non-synonymous variants involving the C-terminal region presented a more severe phenotype with additional neurological features, including central visual impairment, hyperkinetic movement disorder, and epilepsy or electroencephalography abnormalities. Reconstituted fusion involving a lipid-mixing assay indicated impairment in vesicle fusion as one of the possible associated disease mechanisms. The genetic synaptopathy caused by VAMP2 de novo mutations highlights the key roles of this gene in human brain development and function.

A homozygous MED11 C-terminal variant causes a lethal neurodegenerative disease.

PURPOSE: The mediator (MED) multisubunit-complex modulates the activity of the transcriptional machinery, and genetic defects in different MED subunits (17, 20, 27) have been implicated in neurologic diseases. In this study, we identified a recurrent homozygous variant in MED11 (c.325C>T; p.Arg109Ter) in 7 affected individuals from 5 unrelated families. METHODS: To investigate the genetic cause of the disease, exome or genome sequencing were performed in 5 unrelated families identified via different research networks and Matchmaker Exchange. Deep clinical and brain imaging evaluations were performed by clinical pediatric neurologists and neuroradiologists. The functional effect of the candidate variant on both MED11 RNA and protein was assessed using reverse transcriptase polymerase chain reaction and western blotting using fibroblast cell lines derived from 1 affected individual and controls and through computational approaches. Knockouts in zebrafish were generated using clustered regularly interspaced short palindromic repeats/Cas9. RESULTS: The disease was characterized by microcephaly, profound neurodevelopmental impairment, exaggerated startle response, myoclonic seizures, progressive widespread neurodegeneration, and premature death. Functional studies on patient-derived fibroblasts did not show a loss of protein function but rather disruption of the C-terminal of MED11, likely impairing binding to other MED subunits. A zebrafish knockout model recapitulates key clinical phenotypes. CONCLUSION: Loss of the C-terminal of MED subunit 11 may affect its binding efficiency to other MED subunits, thus implicating the MED-complex stability in brain development and neurodegeneration.

Design, synthesis and biological evaluation of novel aminopropylcarboxamide derivatives as sigma ligands.

In our continuing effort to develop novel sigma receptor (SR) ligands, we present the design, synthesis and binding studies of a small library of aminopropylcarboxamide derivatives, obtained from a deconstruction of the piperidine ring of previously synthesized piperidine-based compounds. The best results were achieved with benzofuran (5c, 5g) and quinoline (5a, 5e) derivatives. These compounds revealed the highest affinity for both receptor subtypes. In particular, the 3,4-dimethoxyphenyl derivatives 5e and 5g showed the highest selectivity profile for S2R, especially the quinoline derivative 5e exhibited a 35-fold higher affinity for S2R subtype. The cytotoxic activity of aforementioned compounds was evaluated against SKBR3 and MCF7 cell lines, widely used for breast cancer studies. Whereas the potency of 5g was similar that of Siramesine and Haloperidol in both cell lines, compounds 5a, 5c and 5e exhibited a potency at least comparable to that of Haloperidol in SKBR3 cells. A molecular modelling evaluation towards the S2R binding site, confirmed the strong interaction of compound 5e thus justifying its highest S2R affinity.

Synthesis, Biological Evaluation and Computational Studies of New Hydrazide Derivatives Containing 1,3,4-Oxadiazole as Antitubercular Agents.

To extend our screening for novel antimycobacterial molecules, we have designed, synthesized, and biologically evaluated a library of 14 new hydrazide derivatives containing 1,3,4-oxadiazole core. A variety of mycobacterial strains, including some drug-resistant strains, were tested for antimycobacterial activity. Among the compounds tested, five showed high antimycobacterial activity (MIC values of 8 µg/mL) against M. tuberculosis H37Ra attenuated strain, and two derivatives were effective (MIC of 4 µg/mL) against pyrazinamide-resistant strains. Furthermore, the novel compounds were tested against the fungal C. albicans strain, showing no antimycotic activity, and thus demonstrating a good selectivity profile. Notably, they also exhibited low cytotoxicity against human SH-SY5Y cells. The molecular modeling carried out suggested a plausible mechanism of action towards the active site of the InhA enzyme, which confirmed our hypothesis. In conclusion, the active compounds were predicted in silico for ADME properties, and all proved to be potentially orally absorbed in humans.

Insights on peptide topology in the computational design of protein ligands: the example of lysozyme binding peptides.

Herein, we compared the ability of linear and cyclic peptides generated in silico to target different protein sites: internal pockets and solvent-exposed sites. We selected human lysozyme (HuL) as a model target protein combined with the computational evolution of linear and cyclic peptides. The sequence evolution of these peptides was based on the PARCE algorithm. The generated peptides were screened based on their aqueous solubility and HuL binding affinity. The latter was evaluated by means of scoring functions and atomistic molecular dynamics (MD) trajectories in water, which allowed prediction of the structural features of the protein-peptide complexes. The computational results demonstrated that cyclic peptides constitute the optimal choice for solvent exposed sites, while both linear and cyclic peptides are capable of targeting the HuL pocket effectively. The most promising binders found in silico were investigated experimentally by surface plasmon resonance (SPR), nuclear magnetic resonance (NMR), and electrospray ionization mass spectrometry (ESI-MS) techniques. All tested peptides displayed dissociation constants in the micromolar range, as assessed by SPR; however, both NMR and ESI-MS suggested multiple binding modes, at least for the pocket binding peptides. A detailed NMR analysis confirmed that both linear and cyclic pocket peptides correctly target the binding site they were designed for.

Computational Evolution of Beta-2-Microglobulin Binding Peptides for Nanopatterned Surface Sensors.

The bottom-up design of smart nanodevices largely depends on the accuracy by which each of the inherent nanometric components can be functionally designed with predictive methods. Here, we present a rationally designed, self-assembled nanochip capable of capturing a target protein by means of pre-selected binding sites. The sensing elements comprise computationally evolved peptides, designed to target an arbitrarily selected binding site on the surface of beta-2-Microglobulin (ß2m), a globular protein that lacks well-defined pockets. The nanopatterned surface was generated by an atomic force microscopy (AFM)-based, tip force-driven nanolithography technique termed nanografting to construct laterally confined self-assembled nanopatches of single stranded (ss)DNA. These were subsequently associated with an ssDNA-peptide conjugate by means of DNA-directed immobilization, therefore allowing control of the peptide's spatial orientation. We characterized the sensitivity of such peptide-containing systems against ß2m in solution by means of AFM-based differential topographic imaging and surface plasmon resonance (SPR) spectroscopy. Our results show that the confined peptides are capable of specifically capturing ß2m from the surface-liquid interface with micromolar affinity, hence providing a viable proof-of-concept for our approach to peptide design.

Exploring the use of leucine zippers for the generation of a new class of inclusion bodies for pharma and biotechnological applications.

BACKGROUND: Inclusion bodies (IBs) are biologically active protein aggregates forming natural nanoparticles with a high stability and a slow-release behavior. Because of their nature, IBs have been explored to be used as biocatalysts, in tissue engineering, and also for human and animal therapies. To improve the production and biological efficiency of this nanomaterial, a wide range of aggregation tags have been evaluated. However, so far, the presence in the IBs of bacterial impurities such as lipids and other proteins coexisting with the recombinant product has been poorly studied. These impurities could strongly limit the potential of IB applications, being necessary to control the composition of these bacterial nanoparticles. Thus, we have explored the use of leucine zippers as alternative tags to promote not only aggregation but also the generation of a new type of IB-like protein nanoparticles with improved physicochemical properties. RESULTS: Three different protein constructs, named GFP, J-GFP-F and J/F-GFP were engineered. J-GFP-F corresponded to a GFP flanked by two leucine zippers (Jun and Fos); J/F-GFP was formed coexpressing a GFP fused to Jun leucine zipper (J-GFP) and a GFP fused to a Fos leucine zipper (F-GFP); and, finally, GFP was used as a control without any tag. All of them were expressed in Escherichia coli and formed IBs, where the aggregation tendency was especially high for J/F-GFP. Moreover, those IBs formed by J-GFP-F and J/F-GFP constructs were smaller, rougher, and more amorphous than GFP ones, increasing surface/mass ratio and, therefore, surface for protein release. Although the lipid and carbohydrate content were not reduced with the addition of leucine zippers, interesting differences were observed in the protein specific activity and conformation with the addition of Jun and Fos. Moreover, J-GFP-F and J/F-GFP nanoparticles were purer than GFP IBs in terms of protein content. CONCLUSIONS: This study proved that the use of leucine zippers strategy allows the formation of IBs with an increased aggregation ratio and protein purity, as we observed with the J/F-GFP approach, and the formation of IBs with a higher specific activity, in the case of J-GFP-F IBs. Thus, overall, the use of leucine zippers seems to be a good system for the production of IBs with more promising characteristics useful for pharma or biotech applications.

Biallelic mutations in neurofascin cause neurodevelopmental impairment and peripheral demyelination.

Axon pathfinding and synapse formation are essential processes for nervous system development and function. The assembly of myelinated fibres and nodes of Ranvier is mediated by a number of cell adhesion molecules of the immunoglobulin superfamily including neurofascin, encoded by the NFASC gene, and its alternative isoforms Nfasc186 and Nfasc140 (located in the axonal membrane at the node of Ranvier) and Nfasc155 (a glial component of the paranodal axoglial junction). We identified 10 individuals from six unrelated families, exhibiting a neurodevelopmental disorder characterized with a spectrum of central (intellectual disability, developmental delay, motor impairment, speech difficulties) and peripheral (early onset demyelinating neuropathy) neurological involvement, who were found by exome or genome sequencing to carry one frameshift and four different homozygous non-synonymous variants in NFASC. Expression studies using immunostaining-based techniques identified absent expression of the Nfasc155 isoform as a consequence of the frameshift variant and a significant reduction of expression was also observed in association with two non-synonymous variants affecting the fibronectin type III domain. Cell aggregation studies revealed a severely impaired Nfasc155-CNTN1/CASPR1 complex interaction as a result of the identified variants. Immunofluorescence staining of myelinated fibres from two affected individuals showed a severe loss of myelinated fibres and abnormalities in the paranodal junction morphology. Our results establish that recessive variants affecting the Nfasc155 isoform can affect the formation of paranodal axoglial junctions at the nodes of Ranvier. The genetic disease caused by biallelic NFASC variants includes neurodevelopmental impairment and a spectrum of central and peripheral demyelination as part of its core clinical phenotype. Our findings support possible overlapping molecular mechanisms of paranodal damage at peripheral nerves in both the immune-mediated and the genetic disease, but the observation of prominent central neurological involvement in NFASC biallelic variant carriers highlights the importance of this gene in human brain development and function.

Design, synthesis and antimycobacterial activity of benzoxazinone derivatives and open-ring analogues: Preliminary data and computational analysis.

This study examines in depth benzoxazine nucleus for antimycobacterial property. We synthesized some benzoxazin-2-one and benzoxazin-3-one derivatives, which were tested for activity against a panel of Mycobacterium tuberculosis (Mtb) strains, including H37Ra, H37Rv and some resistant strains. Several compounds displayed a high antimycobacterial activity and the three isoniazid analogue derivatives 8a-c exhibited a MIC range of 0.125-0.250 µg/mL (0.37-0.75 µM) against strain H37Ra, therefore lower than the isoniazid reference drug. Two benzoxazin-2-one derivatives, 1c and 5j, together with isoniazid-analogue compound 8a, also revealed low MIC values against resistant strains and proved highly selective for mycobacterial cells, compared to mammalian Vero cells. To predict whether molecule 8a is able to interact with the active site of InhA, we docked it into the crystal structure; indeed, during the molecular dynamic simulation the compound never left the protein pocket. The more active compounds were predicted for ADME properties and all proved to be potentially orally active in humans.

Binding affinity prediction of nanobody-protein complexes by scoring of molecular dynamics trajectories.

Nanobodies offer a viable alternative to antibodies for engineering high affinity binders. Their small size has an additional advantage: it allows exploiting computational protocols for optimizing their biophysical features, such as the binding affinity. The efficient prediction of this quantity is still considered a daunting task especially for modelled complexes. We show how molecular dynamics can successfully assist in the binding affinity prediction of modelled nanobody-protein complexes. The approximate initial configurations obtained by in silico design must undergo large rearrangements before achieving a stable conformation, in which the binding affinity can be meaningfully estimated. The scoring functions developed for the affinity evaluation of crystal structures will provide accurate estimates for modelled binding complexes if the scores are averaged over long finite temperature molecular dynamics simulations.

Computational design of cyclic peptides for the customized oriented immobilization of globular proteins.

The oriented immobilization of proteins, key for the development of novel responsive biomaterials, relies on the availability of effective probes. These are generally provided by standard approaches based on in vivo maturation and in vitro selection of antibodies and/or aptamers. These techniques can suffer technical problems when a non-immunogenic epitope needs to be targeted. Here we propose a strategy to circumvent this issue by in silico design. In our method molecular binders, in the form of cyclic peptides, are computationally evolved by stochastically exploring their sequence and structure space to identify high-affinity peptides for a chosen epitope of a target globular protein: here a solvent-exposed site of ß2-microglobulin (ß2m). Designed sequences were screened by explicit solvent molecular dynamics simulations (MD) followed by experimental validation. Five candidates gave dose-response surface plasmon resonance signals with dissociation constants in the micromolar range. One of them was further analyzed by means of isothermal titration calorimetry, nuclear magnetic resonance, and 250 ns of MD. Atomic-force microscopy imaging showed that this peptide is able to immobilize ß2m on a gold surface. In short, we have shown by a variety of experimental techniques that it is possible to capture a protein through an epitope of choice by computational design.

Phase behaviour of self-assembled monolayers controlled by tuning physisorbed and chemisorbed states: A lattice-model view.

The self-assembly of molecules on surfaces into 2D structures is important for the bottom-up fabrication of functional nanomaterials, and the self-assembled structure depends on the interplay between molecule-molecule interactions and molecule-surface interactions. Halogenated benzene derivatives on platinum have been shown to have two distinct adsorption states: a physisorbed state and a chemisorbed state, and the interplay between the two can be expected to have a profound effect on the self-assembly and phase behaviour of these systems. We developed a lattice model that explicitly includes both adsorption states, with representative interactions parameterised using density functional theory calculations. This model was used in Monte Carlo simulations to investigate pattern formation of hexahalogenated benzene molecules on the platinum surface. Molecules that prefer the physisorbed state were found to self-assemble with ease, depending on the interactions between physisorbed molecules. In contrast, molecules that preferentially chemisorb tend to get arrested in disordered phases. However, changing the interactions between chemisorbed and physisorbed molecules affects the phase behaviour. We propose functionalising molecules in order to tune their adsorption states, as an innovative way to control monolayer structure, leading to a promising avenue for directed assembly of novel 2D structures.

Computational Mutagenesis of Antibody Fragments: Disentangling Side Chains from ΔΔG Predictions.

The development of highly potent antibodies and antibody fragments as binding agents holds significant implications in fields such as biosensing and biotherapeutics. Their binding strength is intricately linked to the arrangement and composition of residues at the binding interface. Computational techniques offer a robust means to predict the three-dimensional structure of these complexes and to assess the affinity changes resulting from mutations. Given the interdependence of structure and affinity prediction, our objective here is to disentangle their roles. We aim to evaluate independently six side-chain reconstruction methods and ten binding affinity estimation techniques. This evaluation was pivotal in predicting affinity alterations due to single mutations, a key step in computational affinity maturation protocols. Our analysis focuses on a data set comprising 27 distinct antibody/hen egg white lysozyme complexes, each with crystal structures and experimentally determined binding affinities. Using six different side-chain reconstruction methods, we transformed each structure into its corresponding mutant via in silico single-point mutations. Subsequently, these structures undergo minimization and molecular dynamics simulation. We therefore estimate ΔΔG values based on the original crystal structure, its energy-minimized form, and the ensuing molecular dynamics trajectories. Our research underscores the critical importance of selecting reliable side-chain reconstruction methods and conducting thorough molecular dynamics simulations to accurately predict the impact of mutations. In summary, our study demonstrates that the integration of conformational sampling and scoring is a potent approach to precisely characterizing mutation processes in single-point mutagenesis protocols and crucial for computational antibody design.

PEGylated SOCS3 Mimetics Encapsulated into PLGA-NPs as Selective Inhibitors of JAK/STAT Pathway in TNBC Cells.

Introduction: SOCS3 (suppressor of cytokine signaling 3) protein is a crucial regulator of cytokine-induced inflammation, and its administration has been shown to have therapeutic effects. Recently, we designed a chimeric proteomimetic of SOCS3, mimicking the interfacing regions of a ternary complex composed of SOCS3, JAK2 (Janus kinase 2) and gp130 (glycoprotein 130) proteins. The derived chimeric peptide, KIRCONG chim, demonstrated limited mimetic function owing to its poor water solubility. Methods: We report investigations concerning a PEGylated variant of KIRCONG mimetic, named KIRCONG chim, bearing a PEG (Polyethylene glycol) moiety as a linker of noncontiguous SOCS3 regions. Its ability to bind to the catalytic domain of JAK2 was evaluated through MST (MicroScale Thermophoresis), as well as its stability in biological serum assays. The structural features of the cyclic compounds were investigated by CD (circular dichroism), nuclear magnetic resonance (NMR), and molecular dynamic (MD) studies. To evaluate the cellular effects, we employed a PLGA-nanoparticle as a delivery system after characterization using DLS and SEM techniques. Results: KIRCONG chim PEG-revealed selective penetration into triple-negative breast cancer (TNBC) MDA-MB-231 cells with respect to the human breast epithelial cell line (MCF10A), acting as a potent inhibitor of STAT3 phosphorylation. Discussion: Overall, the data indicated that miniaturization of the SOCS3 protein is a promising therapeutic approach for aberrant dysregulation of JAK/STAT during cancer progression.

Allelic heterogeneity and abnormal vesicle recycling in PLAA-related neurodevelopmental disorders.

The human PLAA gene encodes Phospholipase-A2-Activating-Protein (PLAA) involved in trafficking of membrane proteins. Through its PUL domain (PLAP, Ufd3p, and Lub1p), PLAA interacts with p97/VCP modulating synaptic vesicles recycling. Although few families carrying biallelic PLAA variants were reported with progressive neurodegeneration, consequences of monoallelic PLAA variants have not been elucidated. Using exome or genome sequencing we identified PLAA de-novo missense variants, affecting conserved residues within the PUL domain, in children affected with neurodevelopmental disorders (NDDs), including psychomotor regression, intellectual disability (ID) and autism spectrum disorders (ASDs). Computational and in-vitro studies of the identified variants revealed abnormal chain arrangements at C-terminal and reduced PLAA-p97/VCP interaction, respectively. These findings expand both allelic and phenotypic heterogeneity associated to PLAA-related neurological disorders, highlighting perturbed vesicle recycling as a potential disease mechanism in NDDs due to genetic defects of PLAA.

Replica-exchange optimization of antibody fragments.

In the framework of the rational design of macromolecules capable of binding to a specific target for biosensing applications, we here further develop an evolutionary protocol designed to optimize the binding affinity of protein binders. In particular we focus on the optimization of the binding portion of small antibody fragments known as nanobodies (or VHH) and choose the hen egg white lysozyme (HEWL) as our target. By implementing a replica exchange scheme for this optimization, we show that an initial hit is not needed and similar solutions can be found by either optimizing an already known anti-HEWL VHH or a randomly selected binder (here a VHH selective towards another macromolecule). While we believe that exhaustive searches of the mutation space are most appropriate when only few key residues have to be optimized, in case a lead binder is not available the proposed evolutionary algorithm should be instead the method of choice.