Modeling beta-sheet peptide-protein interactions: Rosetta FlexPepDock in CAPRI rounds 38-45

Peptide-protein docking is challenging due to the considerable conformational freedom of the peptide. CAPRI rounds 38-45 included two peptide-protein interactions, both characterized by a peptide forming an additional beta strand of a beta sheet in the receptor. Using the Rosetta FlexPepDock peptide docking protocol we generated top-performing, high-accuracy models for targets 134 and 135, involving an interaction between a peptide derived from L-MAG with DLC8. In addition, we were able to generate the only medium-accuracy models for a particularly challenging target, T121. In contrast to the classical peptide-mediated interaction, in which receptor side chains contact both peptide backbone and side chains, beta-sheet complementation involves a major contribution to binding by hydrogen bonds between main chain atoms. To establish how binding affinity and specificity are established in this special class of peptide-protein interactions, we extracted PeptiDBeta, a benchmark of solved structures of different protein domains that are bound by peptides via beta-sheet complementation, and tested our protocol for global peptide-docking PIPER-FlexPepDock on this dataset. We find that the beta-strand part of the peptide is sufficient to generate approximate and even high resolution models of many interactions, but inclusion of adjacent motif residues often provides additional information necessary to achieve high resolution model quality.     

Interplay between hydrophobic cluster and loop propensity in beta-hairpin formation

Autonomously folding beta-hairpins have recently emerged as powerful tools for elucidating the origins of antiparallel beta-sheet folding preferences. Analysis of such model systems has suggested four potential sources of beta-sheet stability: (1) the conformational propensity of the loop segment that connects adjacent strands; (2) favorable contacts between side-chains on adjacent strands; (3) interstrand hydrogen bonds; and (4) the intrinsic beta-sheet propensities of the strand residues. We describe the design and analysis of a series of isomeric 20 residue peptides in which factors (1)-(4) are identical. Differences in beta-hairpin formation within this series demonstrate that these four factors, individually, are not sufficient to explain beta-sheet stability. In agreement with the prediction of a simple statistical mechanical model for beta-hairpin formation, our results show that the separation between the loop segment and an interstrand cluster of hydrophobic side-chains strongly influences beta-hairpin size and stability, with a smaller separation leading to greater stability.     

Side chain interactions determine the amyloid organization: a single layer beta-sheet molecular structure of the calcitonin peptide segment 15-19

In this paper we present a detailed atomic model for a protofilament, the most basic organization level, of the amyloid fibre formed by the peptide DFNKF. This pentapeptide is a segment derived from the human calcitonin, a natural amyloidogenic protein. Our model, which represents the outcome of extensive explicit solvent molecular dynamics (MD) simulations of different strand/sheet organizations, is a single beta-sheet filament largely without a hydrophobic core. Nevertheless, this structure is capable of reproducing the main features of the characteristic amyloid fibril organization and provides clues to the molecular basis of its experimental aggregation behaviour. Our results show that the side chains' chemical diversity induces the formation of a complex network of interactions that finally determine the microscopic arrangement of the strands at the protofilament level. This network of interactions, consisting of both side chain-side chain and backbone-side chain interactions, confers on the final single beta-sheet arrangement an unexpected stability, both by enhancing the association of related chemical groups and, at the same time, by shielding the hydrophobic segments from the polar solvent. The chemical physical characterization of this protofilament provides hints to the possible thermodynamical basis of the supra molecular organization that allows the formation of the filaments by lateral association of the preformed protofibrils. Its regular, highly polarized structure shows how other protofilaments can assemble. In terms of structural biology, our results clearly indicate that an amyloid organization implies a degree of complexity far beyond a simple nonspecific association of peptide strands via amide hydrogen bonds.     

Sequence dependence of beta-hairpin structure: comparison of a salt bridge and an aromatic interaction

A comparison of the contributions and position dependence of cross-strand electrostatic and aromatic side-chain interactions to beta-sheet stability has been performed by using nuclear magnetic resonance in a well-folded beta-hairpin peptide of the general sequence XRTVXVdPGOXITQX. Phe-Phe and Glu-Lys pairs were varied at the internal and terminal non-hydrogen-bonded position, and the resulting stability was measured by the effects on alpha-hydrogen and aromatic hydrogen chemical shifts. It was determined that the introduction of a Phe-Phe pair resulted in a more folded peptide, regardless of position, and a more tightly folded core. Substitution of the Glu-Lys pair at the internal position results in a less folded peptide and increased fraying at the terminal residues. Upfield shifting of the aromatic hydrogens provided evidence for an edge-face aromatic interaction, regardless of position of the Phe-Phe pair. In peptides with two Phe-Phe pairs, substitution with Glu-Lys at either position resulted in a weakening of the aromatic interaction and a subsequent decrease in peptide stability. Thermal denaturation of the peptides containing Phe-Phe indicates that the aromatic interaction is enthalpically favored, whereas the folding of hairpins with cross-strand Glu-Lys pairs was less enthalpically favorable but entropically more favorable.     

Insights into stabilizing weak interactions in designed peptide beta-hairpins

beta-Hairpin peptides (two anti-parallel strands linked by a reverse beta-turn) have emerged as the simplest systems for probing weak interactions in beta-sheet folding. We describe a model 16-residue hairpin system (peptide beta1: KKYTVSINGKKITVSI) designed around the anti-parallel beta-sheet DNA binding motif of the Met repressor dimer in which two beta-strand sequences are linked through an Asn-Gly type I' beta-turn. The peptide is significantly folded in aqueous solution and has a well-defined conformation as evident from an abundance of NOE data. We review a number of analogues of beta1 designed to estimate the energetic contribution of electrostatic (ion pairing) interactions to hairpin stability, to examine effects of cooperativity and preorganization in determining the energetics of weak interactions, and examine the effects on stability and conformation of incorporation of a three-histidine motif on one face of the hairpin capable of zinc complexation.     

De novo design of a monomeric three-stranded antiparallel beta-sheet

Here we describe the NMR conformational study of a 20-residue linear peptide designed to fold into a monomeric three-stranded antiparallel beta-sheet in aqueous solution. Experimental and statistical data on amino acid beta-turn and beta-sheet propensities, cross-strand side-chain interactions, solubility criteria, and our previous experience with beta-hairpins were considered for a rational selection of the peptide sequence. Sedimentation equilibrium measurements and NMR dilution experiments provide evidence that the peptide is monomeric. Analysis of 1H and 13C-NMR parameters of the peptide, in particular NOEs and chemical shifts, and comparison with data obtained for two 12-residue peptides encompassing the N- and C-segments of the designed sequence indicates that the 20-residue peptide folds into the expected conformation. Assuming a two-state model, the exchange kinetics between the beta-sheet and the unfolded peptide molecules is in a suitable range to estimate the folding rate on the basis of the NMR linewidths of several resonances. The time constant for the coil-beta-sheet transition is of the order of several microseconds in the designed peptide. Future designs based on this peptide system are expected to contribute greatly to our knowledge of the many factors involved in beta-sheet formation and stability.     

Position effect of cross-strand side-chain interactions on beta-hairpin formation

Previous conformational analysis of 10-residue linear peptides enabled us to identify some cross-strand side-chain interactions that stabilize beta-hairpin conformations. The stabilizing influence of these interactions appeared to be greatly reduced when the interaction was located at the N- and C-termini of these 10-residue peptides. To investigate the effect of the position relative to the turn of favorable interactions on beta-hairpin formation, we have designed two 15-residue beta-hairpin forming peptides with the same residue composition and differing only in the location of two residues within the strand region. The conformational properties of these two peptides in aqueous solution were studied by 1H and 13C NMR. Differences in the conformational behavior of the two designed 15-residue peptides suggest that the influence of stabilizing factors for beta-hairpin formation, in particular, cross-strand side-chain interactions, depends on their proximity to the turn. Residues adjacent to the turn are most efficient in that concern. This result agrees with the proposal that the turn region acts as the driving force in beta-hairpin folding.     

Design and characterization of a model αβ peptide

CD and nmr spectroscopy were used to compare the conformational properties of two related peptides. One of the peptides, Model AB, was designed to adopt a helix-turn-extended strand (αβ) tertiary structure in water that might be stabilized by hydrophobic interactions between two leucine residues in the amino-terminal segment and two methionine residues in the carboxyl terminal segment. The other peptide, AB Helix, has the same amino acid sequence as Model AB except that it lacks the-Pro-Met-Thr-Met-Thr-Gly segment at the carboxyl-terminus. Although the carboxyl-terminal segment of Model AB was found to be unstructured, its presence increases the number of residues in a helical conformation, shifts the pKas of three ionizable side chains by 1 pH unit or more compared to an unstructured peptide, stabilizes the peptide as a monomer in high concentrations of ammonium sulfate, increases the conformational stability of residues at the terminal ends of the helix, and results in many slowly exchanging amide protons throughout the entire backbone of the peptide. These results suggest that interactions between adjacent segments in a small peptide can have significant structure organizing effects. Similar kinds of interactions may be important in determining the structure of early intermediates in protein folding and may be useful in the de novo design of independently folding peptides. © 1995 John Wiley & Sons, Inc.     

Design and NMR conformational study of a beta-sheet peptide based on Betanova and WW domains

A good approach to test our current knowledge on formation of protein beta-sheets is de novo protein design. To obtain a three-stranded beta-sheet mini-protein, we have built a series of chimeric peptides by taking as a template a previously designed beta-sheet peptide, Betanova-LLM, and incorporating N- and/or C-terminal extensions taken from WW domains, the smallest natural beta-sheet domain that is stable in absence of disulfide bridges. Some Betanova-LLM strand residues were also substituted by those of a prototype WW domain. The designed peptides were cloned and expressed in Escherichia coli. The ability of the purified peptides to adopt beta-sheet structures was examined by circular dichroism (CD). Then, the peptide showing the highest beta-sheet population according to the CD spectra, named 3SBWW-2, was further investigated by 1H and 13C NMR. Based on NOE and chemical shift data, peptide 3SBWW-2 adopts a well defined three-stranded antiparallel beta-sheet structure with a disordered C-terminal tail. To discern between the contributions to beta-sheet stability of strand residues and the C-terminal extension, the structural behavior of a control peptide with the same strand residues as 3SBWW-2 but lacking the C-terminal extension, named Betanova-LYYL, was also investigated. beta-Sheet stability in these two peptides, in the parent Betanova-LLM and in WW-P, a prototype WW domain, decreased in the order WW-P > 3SBWW-2 > Betanova-LYYL > Betanova-LLM. Conclusions about the contributions to beta-sheet stability were drawn by comparing structural properties of these four peptides.     

Correlations among morphology, beta-sheet stability, and molecular structure in prion peptide aggregates

The misfolding of proteins into beta-sheets and the subsequent aggregation of these sheets into fibrous networks underlies many diseases. In this paper, the role of peptide structure in determining the ordering of beta-sheet aggregates and the morphology of fibrils and protofibrils is dissected. Using a series of peptides based on residues 109-122 of the Syrian hamster prion protein (H1) with a range of substitutions at position 117, the link between side chain interactions and beta-sheet thermal stability has been investigated. The thermal stability of beta-sheets is associated with the peptides' ability to adopt the same alignment as wild-type H1, with residue 117 in register across all beta-strands [Silva, R. A. G. D., Barber-Armstrong, W., and Decatur, S. M. (2003) J. Am. Chem. Soc. 125, 13674-13675]. These aligned strands are capable of forming long, rigid, and twisted fibrils (as visualized by atomic force microscopy) which are thermostable. Peptides which do not adopt this strand alignment aggregate to form thin, flexible, and smooth protofibrils. The ability to form ordered aggregates, and thus to form twisted fibrils, is modulated by the structure of the side chain of residue 117.     

The turn sequence directs beta-strand alignment in designed beta-hairpins

A previous NMR investigation of model decapeptides with identical beta-strand sequences and different turn sequences demonstrated that, in these peptide systems, the turn residues played a more predominant role in defining the type of beta-hairpin adopted than cross-strand side-chain interactions. This result needed to be tested in longer beta-hairpin forming peptides, containing more potentially stabilizing cross-strand hydrogen bonds and side-chain interactions that might counterbalance the influence of the turn sequence. In that direction, we report here on the design and 1H NMR conformational study of three beta-hairpin forming pentadecapeptides. The design consists of adding two and three residues at the N- and C-termini, respectively, of the previously studied decapeptides. One of the designed pentadecapeptides includes a potentially stabilizing R-E salt bridge to investigate the influence of this interaction on beta-hairpin stability. We suggest that this peptide self-associates by forming intermolecular salt bridges. The other two pentadecapeptides behave as monomers. A conformational analysis of their 1H NMR spectra reveals that they adopt different types of beta-hairpin structure despite having identical strand sequences. Hence, the beta-turn sequence drives beta-hairpin formation in the investigated pentadecapeptides that adopt beta-hairpins that are longer than the average protein beta-hairpins. These results reinforce our previous suggestion concerning the key role played by the turn sequence in directing the kind of beta-hairpin formed by designed peptides.     

Cross-strand pairing and amyloid assembly

Amino acid cross-strand pairing interactions along a beta-sheet surface have been implicated in protein beta-structural assembly and stability, yet the relative contributions have been difficult to evaluate directly. Here we develop the central core sequence of the Abeta peptide associated with Alzheimer's disease, Abeta(16-22), as an experimental system for evaluating these interactions. The peptide allows for internal comparisons between electrostatic and steric interactions within the beta-sheet and an evaluation of these cross-strand pair contributions to beta-sheet registry. A morphological transition from fibers to hollow nanotubes arises from changes in beta-sheet surface complementarity and provides a convenient indicator of the beta-strand strand registry. The intrinsic beta-sequence and pair correlations are critical to regulate secondary assembly. These studies provide evidence for a critical desolvation step that is not present in most models of the nucleation-dependent pathway for amyloid assembly.     

Phage-display as a tool for quantifying protein stability determinants.

To address questions of protein stability, researchers have increasingly turned to combinatorial approaches that permit the rapid analysis of libraries of protein variants. Phage-display has proved to be a powerful tool for analyzing protein stability due to the large library size and the robustness of the phage particle to a variety of denaturing conditions. With the B1 domain of protein G (GB1) and a camelid heavy chain antibody as model systems, we are using phage-display libraries to experimentally address questions that have generally been addressed in silico, either through computational studies or statistical analysis of known protein structures. One effort has focused on identifying novel solutions to repacking the hydrophobic core of GB1, while maintaining stability comparable to the wild type protein. In a second study, a small set of substitutions in complimentarity-determining region 3 was found to stabilize the framework of the camelid antibody. Another major focus has been to obtain quantitative data on beta-sheet stability determinants. We have successfully adapted a phage-display method for quantitating affinities of protein variants (shotgun alanine scanning) to analysis of GB1 stability. Using this method, we have analyzed the energetic contributions of cross-strand side chain-side chain interactions. Finally, we discuss parameters to consider in using phage-display to discriminate subtle stability differences among fully folded variants. Overall, this method provides a fast approach for quantitatively addressing biophysical questions.     

Peptide hairpins with strand segments containing alpha- and beta-amino acid residues: cross-strand aromatic interactions of facing Phe residues

The incporation of beta-amino acid residues into the strand segments of designed beta-hairpin leads to the formation of polar sheets, since in the case of beta-peptide strands, all adjacent carbonyl groups point in one direction and the amide groups orient in the opposite direction. The conformational analysis of two designed peptide hairpins composed of alpha/beta-hybrid segments are described: Boc-Leu-betaPhe-Val-(D)-Pro-Gly-Leu-betaPhe-Val-OMe (1) and Boc-betaLeu-Phe-betaVal-D-Pro-Gly-betaLeu-Phe-betaVal-OMe (2). A 500-MHz 1H-NMR (nuclear magnetic resonance) analysis in methanol supports a significant population of hairpin conformations in both peptides. Diagnostic nuclear Overhauser effects (NOEs) are observed in both cases. X-ray diffraction studies on single crystals of peptide 1 reveal a beta-hairpin conformation in both the molecules, which constitute the crystallographic asymmetric unit. Three cross-strand hydrogen bonds and a nucleating type II' beta-turn at the D-Pro-Gly segment are observed in the two independent molecules. In peptide 1, the betaPhe residues at positions 2 and 7 occur at the nonhydrogen-bonding position, with the benzyl side chains pointing on opposite faces of the beta-sheet. The observed aromatic centroid-to-centroid distances are 8.92 A (molecule A) and 8.94 A (molecule B). In peptide 2, the aromatic rings must occupy facing positions in antiparallel strands, in the NMR-derived structure.Peptide 1 yields a normal "hairpin-like" CD spectrum in methanol with a minimum at 224 nm. The CD spectrum of peptide 2 reveals a negative band at 234 nm and a positive band at 221 nm, suggestive of an exciton split doublet. Modeling of the facing Phe side chains at the hydrogen-bonding position of a canonical beta-hairpin suggests that interring separation is approximately 4.78 A for the gauche+ gauche- (g+ g-) rotamer. A previously reported peptide beta-hairpin composed of only alpha-amino acids, Boc-Leu-Phe-Val-D-Pro-Gly-Leu-Phe-Val-OMe also exhibited an anomalous far-UV (ultraviolet) CD (circular dichroism) spectrum, which was interpreted in terms of interactions between facing aromatic chromophores, Phe 2 and Phe 7 (C. Zhao, P. L. Polavarapu, C. Das, and P. Balaram, Journal of the American Chemical Society, 2000, Vol 122, pp. 8228-8231).     

The roles of side chain and backbone in protein structure probed with glycine- and sarcosine-rich synthetic leucine zipper peptides

The protein folding problem has long been a formidable challenge. Here we present a synthetic natural motif approach that exploits small preexisting structural models for the dissection of forces important in protein folding. An example for this approach is shown in the modification of a 31-residue leucine zipper peptide with the helix-breaking amino acid glycine and the hydrogen bond-breaking imino acid sarcosine. Circular dichroism and NMR experiments have shown that the glycine-modified leucine zipper peptide adopts a stable helical conformation similar to the native conformation while the sarcosine-modified leucine zipper peptide adopts a random coil conformation. These results provide valuable insight into the current controversy over the relative importance of long-range side chain-side chain interactions versus local backbone interactions in protein structure and suggest that the natural motif strategy may represent a useful model to study protein folding.     

Ester-linked hen egg white lysozyme shows a compact fold in a molecular dynamics simulation - possible causes and sensitivity of experimentally observable quantities to structural changes maintaining this compact fold

Prediction and understanding of the folding and stability of the 3D structure of proteins is still a challenge. The different atomic interactions, such as non polar contacts and hydrogen bonding, are known but their exact relative weights and roles when contributing to protein folding and stability are not identified. Initiated by a previous molecular dynamics simulation of fully ester-linked hen egg white lysozyme (HEWL), which showed a more compact fold of the ester-linked molecule compared to the native one, three variants of this protein are analyzed in the present study. These are 129-residue native HEWL, partly ester-linked HEWL, in which only 34 peptide linkages that are not involved in the helical or β-strand parts of native HEWL were replaced by ester linkages, and fully (126 residues) ester-linked HEWL. Native and partly ester-linked HEWL showed comparable behaviour, whereas fully ester-linked HEWL could not maintain the native secondary structure of HEWL in the simulation and adopted a more compact fold. The conformational changes were analyzed by comparing simulation averaged values of quantities that can be measured by NMR, such as (1)H-(15)N backbone order parameters, residual dipolar couplings, proton-proton NOE distances and (3)J-couplings with the corresponding values derived from experimental NMR data for native HEWL. The information content of the latter appeared to be insufficient to detect the local conformational rearrangements upon esterification of the loop regions of the protein. For fully ester-linked HEWL, a significantly reduced agreement was observed. Upon esterification, the backbone-side chain and side chain-side chain hydrogen-bonding pattern of HEWL changes to maintain its compactness and thus the structural stability of the ester-linked lysozymes.     

Isomerization of the Xaa-Pro peptide bond induced by ionic interactions of arginine

Inclusion of Arg or Pro residues in proteins and peptides has been proved to play an essential role in biochemical functions through ionic interactions, conformational transitions, and formation of turns as well. In this study we present the conformational properties of the Ac-Arg-Ala-Pro (1), Ac-Arg-Ala-Pro-NH₂ (2), Ac-Arg-Pro-Asp-NH₂ (3), and Ac-Arg-Pro-Asp (4) tripeptides, using ¹H-nmr spectroscopy and molecular dynamics. These peptides were modeled with the aim of studying the role of the Arg-guanidinium to carboxylate ionic interactions on the Xaa-Pro peptide bond isomerization. It was found with 1 and 4 that arginine preferentially interacts with the C-terminal carboxylate group, even though the β-carboxylate is also accessible. This tendency of the Arg moiety was found to induce the cis disposition of the Ala-Pro peptide bond in 1. It was also confirmed that the Arg…Asp side chain-side chain ionic interaction in 3 plays a key role in backbone folding and structural stabilization through a type I β-turn. The nmr pattern for 3 showed a remarkable similarity with that for various Arg-Tyr-Asp containing peptides, a sequence that is crucial for the adhesion properties of the Leishmania gp63 glycoprotein. © 1996 John Wiley & Sons, Inc.     

The role of hydrophobic interactions in amyloidogenesis: example of prion-related polypeptides

Conversion of the non-infectious, cellular form of the prion protein (PrP(C)) to the infectious form (PrP(Sc)) is thought to be driven by an alpha-helical to beta-sheet conformational transition. To reveal the sequence determinants which encourage the transition to beta-fold, we study the synthetic peptides associated with hydrophobic conserved fragments of the N-terminal region of the prion protein. The structure of peptides in solution was probed under various thermodynamic conditions employing circular dichroism and steady state fluorescence spectroscopy as well as dye binding assays. The fluorescence methods utilized pyrene residues covalently attached to the end of the model peptides. In aqueous solutions, the structure assessments indicate the formation of metastable peptide aggregates; the molecular conformations within the peptide micelles are largely coiled. This stage in molecular assembly exists without significant beta-strand formation, i.e., before the appearance of any ordered secondary structure detectable by circular dichroism. At moderate concentrations of trifluoroethanol and/or acetonitrile, the conformational ensemble shifts towards beta-strand formation, and the population of the amorphous aggregates decreases significantly. Overall, the present data indicate that hydrophobic interactions between side chains of the peptide variants prevent, in fact, the formation of the rigid beta-sheet structures. Encouragement of beta-folds requires the destabilization of local interactions in the peptide chain, which in vivo might be possible within cell membranes as well as within partly folded molecular forms.     

Motion of spin label side chains in cellular retinol-binding protein: correlation with structure and nearest-neighbor interactions in an antiparallel beta-sheet

A goal in the development of site-directed spin labeling in proteins is to correlate the motion of a nitroxide side chain with local structure, interactions, and dynamics. Significant progress toward this goal has been made using alpha-helical proteins of known structure, and the present study is the first step in a similar exploration of a beta-sheet protein, cellular retinol-binding protein (CRBP). Nitroxide side chains were introduced along both interior and edge strands. At sites in interior strands, the side-chain motion is strongly influenced by interactions with side chains of neighboring strands, giving rise to a rich variety of dynamic modes (weakly ordered, strongly ordered, immobilized) and complex electron paramagnetic resonance spectra that are modulated by strand twist. The interactions giving rise to the dynamic modes are explored using mutagenesis, and the results demonstrate the particular importance of the non-hydrogen-bonded neighbor residue in giving rise to highly ordered states. Along edge strands of the beta-sheet, the motion of the side chain is simple and weakly ordered, resembling that at solvent-exposed surfaces of an alpha-helix. A simple working model is proposed that can account for the wide variety of dynamic modes encountered. Collectively, the results suggest that the nitroxide side chain is an effective probe of side-chain interactions, and that site-directed spin labeling should be a powerful means of monitoring conformational changes that involve changes in beta-sheet topology.     

Modulating protein folding rates in vivo and in vitro by side-chain interactions between the parallel beta strands of green fluorescent protein

We have identified pairs of residues across the two parallel beta strands of green fluorescent protein that facilitate native strand register of the surface-exposed beta barrel. After constructing a suitable host environment around two guest residues, minimizing interactions of the guest residues with surrounding side-chains yet maintaining the wild-type protein structure and the chromophore environment, we introduced a library of cross-strand pairings by cassette mutagenesis. Colonies of Escherichia coli transformed with the library differ in intracellular fluorescence. Most of the fluorescent pairs have predominantly charged and polar guest site residues. The magnitude and the rate of fluorescence acquisition in vivo from transformed E. coli cells varies among the mutants despite comparable levels of protein expression. Spectroscopic measurements of purified mutants show that the native protein structure is maintained. Kinetic studies using purified protein with fully matured chromophores demonstrate that the mutants span a 10-fold range in folding rates with undetectable differences in unfolding rates. Thus, green fluorescent protein provides an ideal system for monitoring determinants of in vivo protein folding. Cross-strand pairings affect both protein stability and folding kinetics by favoring the formation of native strand register preferentially to non-native strand alignments.