Thermodynamics of Antimicrobial Lipopeptide Binding to Membranes: Origins of Affinity and Selectivity

The development of novel antibiotic drugs is one of the most pressing biomedical problems due to the increasing number of antibiotic-resistant pathogens. Antimicrobial peptides and lipopeptides are a promising category of candidates, but the molecular origins of their antimembrane activity is unclear. Here we explore a series of recently developed antimicrobial lipopeptides, using coarse-grained molecular-dynamics simulations and free energy methods to uncover the thermodynamics governing their binding to membranes. Specifically, we quantify C16-KGGK’s binding affinity to the two types of membrane by umbrella sampling. We also examined the origin of C16-KGGK’s selectivity for bacterial versus mammalian membranes by systematically varying the peptide sequence and salt concentration. Our data showed that the C16 hydrophobic tail is the main contributor to its affinity to lipid membrane, whereas the peptide portion is mainly responsible for its selectivity. Furthermore, the electrostatic interaction between the cationic peptide and anionic bacterial membrane plays a significant role in the selectivity.     

Characterization of a potent antimicrobial lipopeptide via coarse-grained molecular dynamics

The prevalence of antibiotic-resistant pathogens is a major medical concern, prompting increased interest in the development of novel antimicrobial compounds. One such set of naturally occurring compounds, known as antimicrobial peptides (AMPs), have broad-spectrum activity, but come with many limitations for clinical use. Recent work has resulted in a set of antimicrobial lipopeptides (AMLPs) with micromolar minimum inhibitory concentrations and excellent selectivity for bacterial membranes. To characterize a potent, synthetic lipopeptide, C16-KGGK, we used multi-microsecond coarse-grained simulations with the MARTINI forcefield, with a total simulation time of nearly 46μs. These simulations show rapid binding of C16-KGGK, which forms micelles in solution, to model bacterial lipid bilayers. Furthermore, upon binding to the surface of the bilayer, these lipopeptides alter the local lipid organization by recruiting negatively charged POPG lipids to the site of binding. It is likely that this drastic reorganization of the bilayer has major effects on bilayer dynamics and cellular processes that depend on specific bilayer compositions. By contrast, the simulations revealed no association between the lipopeptides and model mammalian bilayers. These simulations provide biophysical insights into lipopeptide selectivity and suggest a possible mechanism for antimicrobial action. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Dynamics of Bcl-xL in Water and Membrane: Molecular Simulations

The Bcl2 family of proteins is capable of switching the apoptotic machinery by directly controlling the release of apoptotic factors from the mitochondrial outer membrane. They have ‘pro’ and ‘anti’-apoptotic subgroups of proteins which antagonize each other’s function; however a detailed atomistic understanding of their mechanisms based on the dynamical events, particularly in the membrane, is lacking. Using molecular dynamics simulations totaling 1.6µs we outline the major differences between the conformational dynamics in water and in membrane. Using implicit models of solvent and membrane, the simulated results reveal a picture that is in agreement with the ‘hit-and run’ concept which states that BH3-only peptides displace the tail (which acts as a pseudo substrate of the protein itself) from its binding pocket; this helps the membrane association of the protein after which the BH3 peptide becomes free. From simulations, Bcl-xL appears to be auto-inhibited by its C-terminal tail that embeds into and covers the hydrophobic binding pocket. However the tail is unable to energetically compete with BH3-peptides in water. In contrast, in the membrane, neither the tail nor the BH3-peptides are stable in the binding pocket and appear to be easily dissociated off as the pocket expands in response to the hydrophobic environment. This renders the binding pocket large and open, thus receptive to interactions with other protein partners. Principal components of the motions are dramatically different in the aqueous and in the membrane environments and provide clues regarding the conformational transitions that Bcl-xL undergoes in the membrane, in agreement with the biochemical data.     

Interaction of the Antimicrobial Peptide Polymyxin B1 with Both Membranes of E. coli: A Molecular Dynamics Study

Antimicrobial peptides are small, cationic proteins that can induce lysis of bacterial cells through interaction with their membranes. Different mechanisms for cell lysis have been proposed, but these models tend to neglect the role of the chemical composition of the membrane, which differs between bacterial species and can be heterogeneous even within a single cell. Moreover, the cell envelope of Gram-negative bacteria such as E. coli contains two membranes with differing compositions. To this end, we report the first molecular dynamics simulation study of the interaction of the antimicrobial peptide, polymyxin B1 with complex models of both the inner and outer membranes of E. coli. The results of >16 microseconds of simulation predict that polymyxin B1 is likely to interact with the membranes via distinct mechanisms. The lipopeptides aggregate in the lipopolysaccharide headgroup region of the outer membrane with limited tendency for insertion within the lipid A tails. In contrast, the lipopeptides readily insert into the inner membrane core, and the concomitant increased hydration may be responsible for bilayer destabilization and antimicrobial function. Given the urgent need to develop novel, potent antibiotics, the results presented here reveal key mechanistic details that may be exploited for future rational drug development.     

Short arginine-rich lipopeptides: From self-assembly to antimicrobial activity

In this paper, we examine antimicrobial and cytotoxic activities, self-assembly and interactions with anionic and zwitterionic membranes of short arginine-rich lipopeptides: C₁₆-RRRR-NH₂, C₁₄-RRRR-NH₂, C₁₂-RRRR-NH₂, and C₁₆-PRRR-NH₂. They show a tendency to self-assembly into micelles, but it is not required for antimicrobial activity. The membrane binding of the lipopeptides can be accompanied by other factors such as: peptide aggregation, pore formation or micellization of phospholipid bilayer. The shortening of the acyl chain results in compounds with a lower haemolytic activity and a slightly improved antimicrobial activity against Gram-positive bacteria, what indicates enhanced cell specificity. Results of coarse-grained molecular dynamics simulations indicate different organization of membrane lipids upon binding of arginine-based lipopeptides and the previously studied lysine-based ones.     

Selective acylation enhances membrane charge sensitivity of the antimicrobial peptide mastoparan-x

The partitioning of the wasp venom peptide mastoparan-X (MPX) into neutral and negatively charged lipid membranes has been compared with two new synthetic analogs of MPX where the N^α-terminal of MPX was acylated with propanoic acid (PA) and octanoic acid (OA). The acylation caused a considerable change in the membrane partitioning properties of MPX and it was found that the shorter acylation with PA gave improved affinity and selectivity toward negatively charged membranes, whereas OA decreased the selectivity. Based on these findings, we hypothesize that minor differences in the embedding and positioning of the peptide in the membrane caused by either PA or OA acylation play a critical role in the fine-tuning of the effective charge of the peptide and thereby the fine-tuning of the peptide's selectivity between neutral and negatively charged lipid membranes. This finding is unique compared to previous reports where peptide acylation enhanced membrane affinity but also resulted in impaired selectivity. Our result may provide a method of enhancing selectivity of antimicrobial peptides toward bacterial membranes due to their high negative charge—a finding that should be investigated for other, more potent antimicrobial peptides in future studies.     

Contributions of the Complementarity Determining Regions to the Thermal Stability of a Single-Domain Antibody

Single domain antibodies (sdAbs) are the recombinantly-expressed variable domain from camelid (or shark) heavy chain only antibodies and provide rugged recognition elements. Many sdAbs possess excellent affinity and specificity; most refold and are able to bind antigen after thermal denaturation. The sdAb A3, specific for the toxin Staphylococcal enterotoxin B (SEB), shows both sub-nanomolar affinity for its cognate antigen (0.14 nM) and an unusually high melting point of 85°C. Understanding the source of sdAb A3’s high melting temperature could provide a route for engineering improved melting temperatures into other sdAbs. The goal of this work was to determine how much of sdAb A3’s stability is derived from its complementarity determining regions (CDRs) versus its framework. Towards answering this question we constructed a series of CDR swap mutants in which the CDRs from unrelated sdAbs were integrated into A3’s framework and where A3’s CDRs were integrated into the framework of the other sdAbs. All three CDRs from A3 were moved to the frameworks of sdAb D1 (a ricin binder that melts at 50°C) and the anti-ricin sdAb C8 (melting point of 60°C). Similarly, the CDRs from sdAb D1 and sdAb C8 were moved to the sdAb A3 framework. In addition individual CDRs of sdAb A3 and sdAb D1 were swapped. Melting temperature and binding ability were assessed for each of the CDR-exchange mutants. This work showed that CDR2 plays a critical role in sdAb A3’s binding and stability. Overall, results from the CDR swaps indicate CDR interactions play a major role in the protein stability.     

Effect of self-assembly on antimicrobial activity of double-chain short cationic lipopeptides

Short cationic antimicrobial lipopeptides with surfactant-like structure are promising antibiotic candidates that preferentially target microbial membranes. Therefore, we focused our study on double-chain lipopeptides, (C_10-16)₂Dab-KKK-NH₂ and (C_10-16)₂Dap-KKK-NH₂, where Dab and Dap are 2,4-diaminobutyric and 2,3-diaminopropionic acids, respectively. We tried to answer a question how the self-assembly behaviour affects biological activities of the tested compounds. The subject compounds were synthesized by solid-phase method and screened for their antimicrobial and haemolytic activities. Cytotoxicity tests on human keratinocytes were carried out for the most promising lipopeptides. Self-assembly properties were evaluated by both experimental and theoretical methods. Interactions with membrane models were examined using the ITC and FTIR techniques. All the lipopeptides studied showed the tendency to self-assembly in solution, and this behaviour was affected by the length of the hydrocarbon chains. Acyl chain elongation supported the formation of the bilayer structure and deprived the lipopeptides of antimicrobial activity. A multi-step mechanism of interaction with a negatively charged membrane was observed for the short-chain lipopeptides, indicating other processes accompanying the binding process. Short-chain lipopeptides were able to penetrate into the liposome’s interior and/or cause the rupture of the liposome, this being compatible with their high antimicrobial activity.     

Crystal structures of N-myristoylated lipopeptide-bound HLA class I complexes indicate reorganization of B-pocket architecture upon ligand binding

Rhesus monkeys have evolved MHC-encoded class I allomorphs such as Mamu-B∗098 that are capable of binding N-myristoylated short lipopeptides rather than conventional long peptides; however, it remains unknown whether such antigen-binding molecules exist in other species, including humans. We herein demonstrate that human leukocyte antigen (HLA)-A∗24:02 and HLA-C∗14:02 proteins, which are known to bind conventional long peptides, also have the potential to bind N-myristoylated short lipopeptides. These HLA class I molecules shared a serine at position 9 (Ser9) with Mamu-B∗098, in contrast to most MHC class I molecules that harbor a larger amino acid residue, such as tyrosine, at this position. High resolution X-ray crystallographic analyses of lipopeptide-bound HLA-A∗24:02 and HLA-C∗14:02 complexes indicated that Ser9 was at the bottom of the B pocket with its small hydroxymethyl side chain directed away from the B-pocket cavity, thereby contributing to the formation of a deep hydrophobic cavity suitable for accommodating the long-chain fatty acid moiety of lipopeptide ligands. Upon peptide binding, however, we found the hydrogen-bond network involving Ser9 was reorganized, and the remodeled B pocket was able to capture the second amino acid residue (P2) of peptide ligands. Apart from the B pocket, virtually no marked alterations were observed for the A and F pockets upon peptide and lipopeptide binding. Thus, we concluded that the structural flexibility of the large B pocket of HLA-A∗2402 and HLA-C∗1402 primarily accounted for their previously unrecognized capacity to bind such chemically distinct ligands as conventional peptides and N-myristoylated lipopeptides.     

Biophysical Investigation of the Membrane-Disrupting Mechanism of the Antimicrobial and Amyloid-Like Peptide Dermaseptin S9

Dermaseptin S9 (Drs S9) is an atypical cationic antimicrobial peptide with a long hydrophobic core and with a propensity to form amyloid-like fibrils. Here we investigated its membrane interaction using a variety of biophysical techniques. Rather surprisingly, we found that Drs S9 induces efficient permeabilisation in zwitterionic phosphatidylcholine (PC) vesicles, but not in anionic phosphatidylglycerol (PG) vesicles. We also found that the peptide inserts more efficiently in PC than in PG monolayers. Therefore, electrostatic interactions between the cationic Drs S9 and anionic membranes cannot explain the selectivity of the peptide towards bacterial membranes. CD spectroscopy, electron microscopy and ThT fluorescence experiments showed that the peptide adopts slightly more β-sheet and has a higher tendency to form amyloid-like fibrils in the presence of PC membranes as compared to PG membranes. Thus, induction of leakage may be related to peptide aggregation. The use of a pre-incorporation protocol to reduce peptide/peptide interactions characteristic of aggregates in solution resulted in more α-helix formation and a more pronounced effect on the cooperativity of the gel-fluid lipid phase transition in all lipid systems tested. Calorimetric data together with ²H- and ³¹P-NMR experiments indicated that the peptide has a significant impact on the dynamic organization of lipid bilayers, albeit slightly less for zwitterionic than for anionic membranes. Taken together, our data suggest that in particular in membranes of zwitterionic lipids the peptide binds in an aggregated state resulting in membrane leakage. We propose that also the antimicrobial activity of Drs S9 may be a result of binding of the peptide in an aggregated state, but that specific binding and aggregation to bacterial membranes is regulated not by anionic lipids but by as yet unknown factors.     

Membrane fluidity,composition, and charge affect the activity and selectivity of the AMP ascaphin-8

Ascaphins are cationic antimicrobial peptides that have been shown to have potential in the treatment of infectious diseases caused by multidrug-resistant pathogens (MDR). However, to date, their principal molecular target and mechanism of action are unknown. Results from peptide prediction software and molecular dynamics simulations confirmed that ascaphin-8 is an alpha-helical peptide. For the first time, the peptide was described as membranotrophic using biophysical approaches including calcein liposome leakage, Laurdan general polarization, and dynamic light scattering. Ascaphin-8’s activity and selectivity were modulated by rearranging the spatial distribution of lysine (Var-K5), aspartic acid (Var-D4) residues, or substitution of phenylalanine with tyrosine (Var-Y). The parental peptide and its variants presented high affinity toward the bacterial membrane model (≤2 μM), but lost activity in sterol-enriched membranes (mammal and fungal models, with cholesterol and ergosterol, respectively). The peptide-induced pore size was estimated to be >20 nm in the bacterial model, with no difference among peptides. The same pattern was observed in membrane fluidity (general polarization) assays, where all peptides reduced membrane fluidity of the bacterial model but not in the models containing sterols. The peptides also showed high activity toward MDR bacteria. Moreover, peptide sensitivity of the artificial membrane models compared with pathogenic bacterial isolates were in good agreement.     

Thermodynamics of Micelle Formation and Membrane Fusion Modulate Antimicrobial Lipopeptide Activity

Antimicrobial lipopeptides (AMLPs) are antimicrobial drug candidates that preferentially target microbial membranes. One class of AMLPs, composed of cationic tetrapeptides attached to an acyl chain, have minimal inhibitory concentrations in the micromolar range against a range of bacteria and fungi. Previously, we used coarse-grained molecular dynamics simulations and free energy methods to study the thermodynamics of their interaction with membranes in their monomeric state. Here, we extended the study to the biologically relevant micellar state, using, to our knowledge, a novel reaction coordinate based on hydrophobic contacts. Using umbrella sampling along this reaction coordinate, we identified the critical transition states when micelles insert into membranes. The results indicate that the binding of these AMLP micelles to membranes is thermodynamically favorable, but in contrast to the monomeric case, there are significant free energy barriers. The height of these free energy barriers depends on the membrane composition, suggesting that the AMLPs’ ability to selectively target bacterial membranes may be as much kinetic as thermodynamic. This mechanism highlights the importance of considering oligomeric state in solution as criterion when optimizing peptides or lipopeptides as antibiotic leads.     

Molecular Basis of S100A1 Activation at Saturating and Subsaturating Calcium Concentrations

The S100A1 protein mediates a wide variety of physiological processes through its binding of calcium (Ca²⁺) and endogenous target proteins. S100A1 presents two Ca²⁺-binding domains: a high-affinity “canonical” EF (cEF) hand and a low-affinity “pseudo” EF (pEF) hand. Accumulating evidence suggests that both Ca²⁺-binding sites must be saturated to stabilize an open state conducive to peptide recognition, yet the pEF hand’s low affinity limits Ca²⁺ binding at normal physiological concentrations. To understand the molecular basis of Ca²⁺ binding and open-state stabilization, we performed 100 ns molecular dynamics simulations of S100A1 in the apo/holo (Ca²⁺-free/bound) states and a half-saturated state, for which only the cEF sites are Ca²⁺-bound. Our simulations indicate that the pattern of oxygen coordination about Ca²⁺ in the cEF relative to the pEF site contributes to the former’s higher affinity, whereas Ca²⁺ binding strongly reshapes the protein’s conformational dynamics by disrupting β-sheet coupling between EF hands. Moreover, modeling of the half-saturated configuration suggests that the open state is unstable and reverts toward a closed state in the absence of the pEF Ca²⁺ ion. These findings indicate that Ca²⁺ binding at the cEF site alone is insufficient to stabilize opening; thus, posttranslational modification of the protein may be required for target peptide binding at subsaturating intracellular Ca²⁺ levels.     

A novel anti-aldolase C antibody specifically interacts with residues 85–102 of the protein

Aldolase C is a brain-specific glycolytic isozyme whose complete repertoire of functions are obscure. This lack of knowledge can be addressed using molecular tools that discriminate the protein from the homologous, ubiquitous paralog aldolase A. The anti-aldolase C antibodies currently available are polyclonal and not highly specific. We obtained the novel monoclonal antibody 9F against human aldolase C, characterized its isoform specificity and tested its performance. First, we investigated the specificity of 9F for aldolase C. Then, using bioinformatic tools coupled to molecular cloning and chemical synthesis approaches, we produced truncated human aldolase C fragments, and assessed 9F binding to these fragments by western blot and ELISA assays. This strategy revealed that residues 85–102 harbor the epitope-containing region recognized by 9F. The efficiency of 9F was demonstrated also for immunoprecipitation assays. Finally, surface plasmon resonance revealed that the protein has a high affinity toward the epitope-containing peptide. Taken together, our findings show that epitope recognition is sequence-driven and is independent of the three-dimensional structure. In conclusion, given its specific molecular interaction, 9F is a novel and powerful tool to investigate aldolase C’s functions in the brain.     

Non-extensitivity and criticality of atomic hydropathicity around a voltage-gated sodium channel’s pore: a modeling study

Voltage-gated sodium channels (NavChs) are pore-forming membrane proteins that regulate the transport of sodium ions through the cell membrane. Understanding the structure and function of NavChs is of major biophysical, as well as clinical, importance given their key role in cellular pathophysiology. In this work, we provide a computational framework for modeling system-size-dependent, i.e., cumulative, atomic properties around a NavCh’s pore. We illustrate our methodologies on the bacterial NavAb channel captured in a closed-pore state where we demonstrate that the atomic environment around its pore exhibits a bi-phasic spatial organization dictated by the structural separation of the pore domains (PDs) from the voltage-sensing domains (VSDs). Accordingly, a mathematical model describing packing of atoms around NavAb’s pore is constructed that allows—under certain conservation conditions—for a power-law approximation of the cumulative hydropathic dipole field effect acting along NavAb’s pore. This verified the non-extensitivity hypothesis for the closed-pore NavAb channel and revealed a long-range hydropathic interactions law regulating atom-packing around the NavAb’s selectivity filter. Our model predicts a PDs-VSDs coupling energy of ∼282.1 kcal/mol corresponding to a global maximum of the atom-packing energy profile. Crucially, we demonstrate for the first time how critical phenomena can emerge in a single-channel structure as a consequence of the non-extensive character of its atomic porous environment.Supplementary InformationThe online version contains supplementary material available at 10.1007/s10867-021-09565-w.     

Integrating Solid-State NMR and Computational Modeling to Investigate the Structure and Dynamics of Membrane-Associated Ghrelin

The peptide hormone ghrelin activates the growth hormone secretagogue receptor 1a, also known as the ghrelin receptor. This 28-residue peptide is acylated at Ser3 and is the only peptide hormone in the human body that is lipid-modified by an octanoyl group. Little is known about the structure and dynamics of membrane-associated ghrelin. We carried out solid-state NMR studies of ghrelin in lipid vesicles, followed by computational modeling of the peptide using Rosetta. Isotropic chemical shift data of isotopically labeled ghrelin provide information about the peptide’s secondary structure. Spin diffusion experiments indicate that ghrelin binds to membranes via its lipidated Ser3. Further, Phe4, as well as electrostatics involving the peptide’s positively charged residues and lipid polar headgroups, contribute to the binding energy. Other than the lipid anchor, ghrelin is highly flexible and mobile at the membrane surface. This observation is supported by our predicted model ensemble, which is in good agreement with experimentally determined chemical shifts. In the final ensemble of models, residues 8–17 form an α-helix, while residues 21–23 and 26–27 often adopt a polyproline II helical conformation. These helices appear to assist the peptide in forming an amphipathic conformation so that it can bind to the membrane.     

Designing effective anticancer-radiopeptides. A Molecular Dynamics study of their interaction with model tumor and healthy cell membranes

One of the greatest merit of the use of radiopeptides in oncology is their selectivity which, however, brings about the drawback that each radiopeptide is specific for a given tumor type. To overcome this problem the direction currently taken in drug design is that of radiolabelling peptide hormones (or their analogues), relying on their intrinsic ability to bind to specific receptors in precise areas of the human body, at the cost, however, of a poor selectivity against healthy cells. We present here an extensive Molecular Dynamics study of a promising alternative inspired by the mechanism through which antimicrobial peptides interact with the negatively charged bacterial membranes. Appropriately modifying the human antimicrobial peptide, LL-37, we designed a functionalized radionuclide carrier capable of binding more strongly to the negatively charged (model) tumor membranes than to the neutral healthy ones. The mechanism behind this behaviour relies on the fact that at the slight acidic pH surrounding tumor tissues the histidines belonging to the peptide get protonated thus making it positively charged. We have investigated by an extended numerical study the way in which this artificial peptide interacts with models of tumor and healthy cell membranes, proving by Potential Mean Force calculations that the affinity of the peptide to model tumor membranes is significantly larger than to healthy ones. These features (high affinity and generic tumor selectivity) recommend antimicrobial derived customized carriers as promising theranostic constructs in cancer diagnostic and therapy.     

Synthesis of thiol-reactive lipopeptide adjuvants. Incorporation into liposomes and study of their mitogenic effect on mouse splenocytes

Synthetic analogues of triacylated and diacylated lipopeptides derived from the N-terminal domain of respectively bacterial and mycoplasmal lipoproteins are highly potent immunoadjuvants when administered either in combination with protein antigens or covalently linked to small peptide epitopes. Because of their amphipathic properties, lipopeptides, such as S-[2,3-bis(palmitoyloxy)-(2RS)-propyl]-N-palmitoyl-(R)-cysteinyl-alanyl-glycine (Pam(3)CAG), can be conveniently incorporated into liposomes and serve as anchors for antigens that are linked to them. To design vaccination constructs based on synthetic peptides and liposomes as vectors. we have accordingly synthesized a series of lipopeptides that differ by the number (Pam(3)C vs Pam(2)C) and nature of the acyl chains (palmitoyl vs oleoyl) and by the presence at their C-terminus of thiol-reactive functions, such as maleimide or bromoacetyl. When incorporated into liposomes, these latter functionalized lipopeptides allow, in aqueous media, a well controlled chemoselective conjugation of HS-peptides to the surface of the vesicles. Using a BALB/c mice splenocyte proliferation assay ([(3)H]thymidine incorporation), we have measured the lymphocyte activation potency of the different lipopeptides. We found that, compared to their free (emulsified) forms, the liposomal lipopeptides were endowed with enhanced mitogenic activities; i.e., up to 2 orders of magnitude for Pam(3)CAG which was more potent than Pam(2)CAG. The impact of functionalization on the cellular activity of Pam(3)CAG was dependent on the thiol-reactive group introduced: whereas the bromoacetyl derivative retained its full activity, the presence of a maleimide group virtually abolished the lymphocyte activation of the lipopeptide. Finally, the substitution of saturated palmitoyl chains by unsaturated oleoyl chains was inhibitory. Thus, thiol-reactive Ol(3)CAG derivatives were the least active mitogens in our assay. Taken together, our findings are of importance for the further optimization of antigen-specific liposomal-based synthetic vaccines; the bromoacetyl derivative of Pam(3)CAG should be a promising lipopeptide derivative serving as an anchor for peptide epitopes while retaining its lymphocyte activation activity.     

Differential binding of Escherichia coli McrA protein to DNA sequences that contain the dinucleotide m5CpG

The Escherichia coli McrA protein, a putative C⁵-methylcytosine/C⁵-hydroxyl methylcytosine-specific nuclease, binds DNA with symmetrically methylated HpaII sequences (Cm5CGG), but its precise recognition sequence remains undefined. To determine McrA’s binding specificity, we cloned and expressed recombinant McrA with a C-terminal StrepII tag (rMcrA-S) to facilitate protein purification and affinity capture of human DNA fragments with m5C residues. Sequence analysis of a subset of these fragments and electrophoretic mobility shift assays with model methylated and unmethylated oligonucleotides suggest that N(Y > R) m5CGR is the canonical binding site for rMcrA-S. In addition to binding HpaII-methylated double-stranded DNA, rMcrA-S binds DNA containing a single, hemimethylated HpaII site; however, it does not bind if A, C, T or U is placed across from the m5C residue, but does if I is opposite the m5C. These results provide the first systematic analysis of McrA’s in vitro binding specificity.