Comparing four different approaches for the determination of inter‐residue interactions provides insight for the structure prediction of helical membrane proteins

Studying inter‐residue interactions provides insight into the folding and stability of both soluble and membrane proteins and is essential for developing computational tools for protein structure prediction. As the first step, various approaches for elucidating such interactions within protein structures have been proposed and proven useful. Since different approaches may grasp different aspects of protein structural folds, it is of interest to systematically compare them. In this work, we applied four approaches for determining inter‐residue interactions to the analysis of three distinct structure datasets of helical membrane proteins and compared their correlation to the three individual quality measures of structures in these datasets. These datasets included one of 35 structures of rhodopsin receptors and bacterial rhodopsins determined at various resolutions, one derived from the HOMEP benchmark dataset previously reported, and one comprising of 139 homology models. It was found that the correlation between the average number of inter‐residue interactions obtained by applying the four approaches and the available structure quality measures varied quite significantly among them. The best correlation was achieved by the approach focusing exclusively on favorable inter‐residue interactions. These results provide interesting insight for the development of objective quality measure for the structure prediction of helical membrane proteins. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 547–556, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Structural analysis of the human cannabinoid receptor one carboxyl‐terminus identifies two amphipathic helices

Recent research has implicated the C‐terminus of G‐protein coupled receptors in key events such as receptor activation and subsequent intracellular sorting, yet obtaining structural information of the entire C‐tail has proven a formidable task. Here, a peptide corresponding to the full‐length C‐tail of the human CB1 receptor (residues 400–472) was expressed in E.coli and purified in a soluble form. Circular dichroism (CD) spectroscopy revealed that the peptide adopts an α‐helical conformation in negatively charged and zwitterionic detergents (48–51% and 36–38%, respectively), whereas it exhibited the CD signature of unordered structure at low concentration in aqueous solution. Interestingly, 27% helicity was displayed at high peptide concentration suggesting that self‐association induces helix formation in the absence of a membrane mimetic. NMR spectroscopy of the doubly labeled (¹⁵N‐ and ¹³C‐) C‐terminus in dodecylphosphocholine (DPC) identified two amphipathic α‐helical domains. The first domain, S401‐F412, corresponds to the helix 8 common to G protein‐coupled receptors while the second domain, A440‐M461, is a newly identified structural motif in the distal region of the carboxyl‐terminus of the receptor. Molecular modeling of the C‐tail in DPC indicates that both helices lie parallel to the plane of the membrane with their hydrophobic and hydrophilic faces poised for critical interactions. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 565–573, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Revisiting “reverse hydrophobic effect”: Applicable only to coil mutations at the surface

In a seminal paper, Pakula and Sauer (Nature, 1990, 344, 363–364) demonstrated that the increase in side‐chain hydrophobicity has a reverse relationship with protein stability. We have addressed this problem with several examples of mutants that span at different locations in protein structure based on secondary structure and solvent accessibility. We confirmed that the stability change upon single coil mutation at exposed region is reversely correlated with hydrophobicity with a single exception. In addition, we found the existence of such relationship in partially buried coil mutants. The stability of exposed helical mutants is governed by conformational properties. In buried and partially buried helical and strand mutants properties reflecting hydrophobicity have direct relationship with stability, whereas an opposite relationship was obtained with entropy and flexibility. The structural analysis of partially buried/exposed mutants showed that the surrounding residues are important for the stability change upon mutation. These results provide insights to understand the general behavior for the stability of proteins upon amino acid substitutions. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 591–599, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Study of the interactions between a proline‐rich protein and a flavan‐3‐ol by NMR: Residual structures in the natively unfolded protein provides anchorage points for the ligands

Astringency is one of the major organoleptic properties of food and beverages that are made from plants, such as tea, chocolate, beer, or red wine. This sensation is thought to be due to interactions between tannins and salivary proline‐rich proteins, which are natively unfolded proteins. A human salivary proline‐rich protein, namely IB‐5, was produced by the recombinant method. Its interactions with a model tannin, epigallocatechin gallate (EGCG), the major flavan‐3‐ol in green tea, were studied here. Circular dichroism experiments showed that IB‐5 presents residual structures (PPII helices) when the ionic strength is close to that in saliva. In the presence of these residual structures, IB‐5 undergoes an increase in structural content upon binding to EGCG. NMR data corroborated the presence of preformed structural elements within the protein prior to binding and a partial assignment was proposed, showing partial structuration. TOCSY experiments showed that amino acids that are involved in PPII helices are more likely to interact with EGCG than those in random coil regions, as if they were anchorage points for the ligand. The signal from IB‐5 in the DOSY NMR spectrum revealed an increase in polydispersity upon addition of EGCG while the mean hydrodynamic radius remained unchanged. This strongly suggests the formation of IB‐5/EGCG aggregates. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 745–756, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Stability of loops in the structure of lactose permease

Structural analysis of peptide fragments has provided useful information on the secondary structure of integral membrane proteins built from a helical bundle (up to seven transmembrane segments). Comparison of those results to recent X-ray crystallographic results showed agreement between the structures of the fragments and the structures of the intact proteins. Lactose permease of Escherichia coli (lac Y) offers an opportunity to test that hypothesis on a substantially larger integral membrane protein. Lac Y contains a bundle of 12 transmembrane segments connected by 11 loops. Eleven segments, each corresponding to one of the loops in this protein, were studied. Five of these segments form defined structures in solution as determined by multidimensional nuclear magnetic resonance. Four peptides form turns, and one peptide reveals the end of one of the transmembrane helices. These results suggest that some loops in helical bundles are stabilized by short-range interactions, particularly in smaller bundles, and such intrinsically stable loops may contribute to protein stability and influence the pathway of folding. Greater conformational flexibility may be found in large integral membrane proteins.     

Sequence and conformational preferences at termini of α‐helices in membrane proteins: Role of the helix environment

α‐helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These α‐helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C‐termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze α‐helices in a high‐resolution dataset of integral α‐helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C‐termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near‐helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins. Proteins 2014; 82:3420–3436. © 2014 Wiley Periodicals, Inc.     

3(10)-Helix adjoining alpha-helix and beta-strand: sequence and structural features and their conservation

Does the amino acid use at the terminal positions of an α‐helix become altered depending on the context—more specifically, when there is an adjoining 3₁₀‐helix, and can a single helical cylinder encompass the resultant composite helix? An analysis of 138 and 107 cases of 3₁₀–α and α–3₁₀ composite helices, respectively, found in known protein structures indicate that the secondary structural element occurring first imposes its characteristics on the sequence of the structural element coming next. Thus, when preceded by a 3₁₀‐helix, the preference of proline to occur at the N1 position of an α‐helix is shifted to the N2 position, a typical characteristic of the C‐terminal capping of the 3₁₀‐helix. When an α‐ or a 3₁₀‐helix leads into a helix of the other type, there is a bend at the junction, especially for the 3₁₀–α composite, with the two junction residues facing inward and buried within the structure. Thus a single helical cylinder may not properly represent a composite helix, the bend providing a means for the tertiary structure to assume a globular shape, very much akin to what a proline‐induced kink does to an α‐helix. The tertiary structural context in which β–3₁₀ and 3₁₀–β composites occurs can be different, causing the angle between the secondary structural elements in the two cases to be different. Composites of 3₁₀‐helices and β‐strands are much more conserved among members in families of homologous structures than those between two types of helices; in many of the former instances, the 3₁₀‐helix constitutes the loops in β‐hairpin or β–β‐corner motifs. The overall fold of the chain may be more conserved than the actual identify of the secondary structure elements in a composite. © 2005 Wiley Periodicals, Inc. Biopolymers 78: 147–162, 2005 This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Structural studies of SNARE complex and its interaction with complexin by molecular dynamics simulation

Since neurotransmitter releasing into the synaptic space delivers electrical signals from presynaptic neural cell to the postsynaptic cell, neurotransmitter secretion must be much orchestrated. Crowded intracellular vesicles involving neurotransmitters present a question of the how secretory vesicles fuse onto the plasma membrane in a fast synchronized fashion. Complexin is one of the most experimentally studied proteins that regulate assembly of fusogenic four‐helix SNARE complex to synchronized neurotransmitter secretion. We used MD simulation to investigate the interaction of complexin with the neural SNARE complex in detail. Our results show that the SNARE complex interacts with the complexin central helix by forming salt bridges and hydrogen bonds. Complexin also can interact with the Q‐SNARE complex instead of synaptobrevin to decrease the Q‐SNARE flexibility. The complexin alpha‐accessory helix and the C‐terminal region of synaptobrevin can interact with the same region of syntaxin. Although the alpha‐accessory helix aids the tight binding of the central helix to the SNARE complex, its proximity with synaptobrevin causes the destabilization of syntaxin and Sn1 helices. This study suggests that the alpha‐accessory helix of complexin can be an inhibiting factor for membrane fusion by competing with synaptobrevin for binding to the Q‐SNARE complex. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 560–570, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Sequence context and modified hydrophobic moment plots help identify 'horizontal' surface helices in transmembrane protein structure prediction

Transmembrane proteins make up at least one-fifth of the genome of most organisms and are critical components of key pathways for cell survival and interactions with the environment. The function of helices found at the membrane surface in transmembrane proteins has not been greatly explored, but it is likely that they play an ancillary role to membrane spanning helices and are analogous to the surface active helices of peripheral membrane proteins, being involved in: lipid association, membrane perturbation, transmembrane signal transduction and regulation, and transmembrane helical bundle formation. Due to the difficulties in obtaining high-resolution structural data for this class of proteins, structure-from-sequence predictive methods continue to be developed as a means to obtain structural models for these largely intractable systems. A simple but effective variant of the hydrophobic moment analysis of amino acid sequences is described here as part of a protocol for distinguishing helical sequences that are parallel to or 'horizontal' at the membrane bilayer/aqueous phase interface from helices that are membrane-embedded or located in extra-membranous domains. This protocol when tested on transmembrane spanning protein amino acid sequences not used in its development, was found to be 84-91% accurate when the results were compared to the partition locations in the corresponding structures determined by X-ray crystallography, and 72% accurate in determining which helices lie horizontal or near horizontal at the lipid interface.     

Structural clues in the sequences of the aquaporins

The large number of sequences available for the aquaporin family represents a valuable source of information to incorporate into three-dimensional structure determination. Phylogenetic analysis was used to define type sequences to avoid extreme over-representation of some subfamilies, and as a measure of the quality of multiple sequence alignment. Inspection of the sequence alignment suggested eight conserved segments that define the core architecture of six transmembrane helices and two functional loops, B and E, projecting into the plane of the membrane. The sum of the core segments and the minimum lengths of the interlinking loops constitute the 208 residues necessary to satisfy the aquaporin architecture. Analysis of hydrophobic and conservation periodicity and of correlated mutations across the alignment indicated the likely assignment and orientation of the helices in the bilayer. This assignment is examined with respect to the structure of the erythrocyte aquaporin 1 determined by electron crystallography. The aquaporin 1 tetramer is described as three rings of helices, each ring with a different exposure to the lipid environment. The sequence analysis clearly suggests that two helices are exposed along their whole lengths, two helices are exposed only at their N termini, and two helices are not exposed to lipid. It is further proposed that, besides loops B and E, the highly conserved motifs on helices 1 and 4, ExxxTxxF/L, could line the water channel.     

Interplay between protein homeostasis networks in protein aggregation and proteotoxicity

The misfolding and aggregation of disease proteins is characteristic of numerous neurodegenerative diseases. Particular neuronal populations are more vulnerable to proteotoxicity while others are more apt to tolerate the misfolding and aggregation of disease proteins. Thus, the cellular environment must play a significant role in determining whether disease proteins are converted into toxic or benign forms. The endomembrane network of eukaryotes divides the cell into different subcellular compartments that possess distinct sets of molecular chaperones and protein interaction networks. Chaperones act as agonists and antagonists of disease protein aggregation to prevent the accumulation of toxic intermediates in the aggregation pathway. Interacting partners can also modulate the conformation and localization of disease proteins and thereby influence proteotoxicity. Thus, interplay between these protein homeostasis network components can modulate the self‐association of disease proteins and determine whether they elicit a toxic or benign outcome. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 229–236, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Exploring the interactions of gliadins with model membranes: Effect of confined geometry and interfaces

Mechanisms leading to the assembly of wheat storage proteins into proteins bodies within the endoplasmic reticulum (ER) of endosperm cells are unresolved today. In this work, physical chemistry parameters which could be involved in these processes were explored. To model the confined environment of proteins within the ER, the dynamic behavior of γ‐gliadins inserted inside lyotropic lamellar phases was studied using FRAP experiments. The evolution of the diffusion coefficient as a function of the lamellar periodicity enabled to propose the hypothesis of an interaction between γ‐gliadins and membranes. This interaction was further studied with the help of phospholipid Langmuir monolayers. γ‐ and ω‐gliadins were injected under DMPC and DMPG monolayers and the two‐dimensional (2D) systems were studied by Brewster angle microscopy (BAM), polarization modulation infrared reflection‐absorption spectroscopy (PM‐IRRAS), and surface tension measurements. Results showed that both gliadins adsorbed under phospholipid monolayers, considered as biological membrane models, and formed micrometer‐sized domains at equilibrium. However, their thicknesses, probed by reflectance measurements, were different: ω‐gliadins aggregates displayed a constant thickness, consistent with a monolayer, while the thickness of γ‐gliadins aggregates increased with the quantity of protein injected. These different behaviors could find some explanations in the difference of aminoacid sequence distribution: an alternate repeated ‐ unrepeated domain within γ‐gliadin sequence, while one unique repeated domain was present within ω‐gliadin sequence. All these findings enabled to propose a model of gliadins self‐assembly via a membrane interface and to highlight the predominant role of wheat prolamin repeated domain in the membrane interaction. In the biological context, these results would mean that the repeated domain could be considered as an anchor for the interaction with the ER membrane and a nucleus point for the formation and growth of protein bodies within endosperm cells. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 610–622, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Peptidic modulators of protein‐protein interactions: Progress and challenges in computational design

With the decline in productivity of drug‐development efforts, novel approaches to rational drug design are being introduced and developed. Naturally occurring and synthetic peptides are emerging as novel promising compounds that can specifically and efficiently modulate signaling pathways in vitro and in vivo. We describe sequence‐based approaches that use peptides to mimic proteins in order to inhibit the interaction of the mimicked protein with its partners. We then discuss a structure‐based approach, in which protein‐peptide complex structures are used to rationally design and optimize peptidic inhibitors. We survey flexible peptide docking techniques and discuss current challenges and future directions in the rational design of peptidic inhibitors. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 505–513, 2009. This article was originally published online as an accepted preprint. The “Published Online”date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Ion channel activity of transmembrane segment 6 of Escherichia coli proton‐dependent manganese transporter

Synthetic peptides corresponding to the sixth transmembrane segment (TMS6) of secondary‐active transporter MntH (Proton‐dependent Manganese Transporter) from Escherichia coli and its two mutations in the functionally important conserved histidine residue were used as a model for structure–function study of MntH. The secondary structure of the peptides was estimated in different environments using circular dichroism spectroscopy. These peptides interacted with and adopted helical conformations in lipid membranes. Electrophysiological experiments demonstrated that TMS6 was able to form multi‐state ion channels in model biological membranes. Electrophysiological properties of these weakly cation‐selective ion channels were strongly dependent on the surrounding pH. Manganese ion, as a physiological substrate of MntH, enhanced the conductivity of TMS6 channels, influenced the transition between closed and open states, and affected the peptide conformations. Moreover, functional properties of peptides carrying two different mutations of His²¹¹ were analogous to in vivo functional characteristics of Nramp/MntH proteins mutated at homologous residues. Hence, a single functionally important TMS can retain some of the functional properties of the full‐length protein. These findings could contribute to understanding the structure–function relationship at the molecular level. However it remains unclear to what extent the peptide‐specific channel activity represents a functional aspect of the full‐length membrane carrier protein. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 718–726, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Importance of the spatial display of charged residues in heparin–peptide interactions

Many studies have examined consensus sequences required for protein‐glycosaminoglycan interactions. Through the synthesis of helical heparin binding peptides, this study probes the relationship between spatial arrangement of positive charge and heparin binding affinity. Peptides with a linear distribution of positive charge along one face of the α‐helix had the highest affinity for heparin. Moving the basic residues away from a single face resulted in drastic changes in heparin binding affinity of up to three orders of magnitude. These findings demonstrate that amino acid sequences, different from the known heparin binding consensus sequences, will form high affinity protein‐heparin binding interactions when the charged residues are aligned linearly. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 290–298, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Membrane-integration Characteristics of Two ABC Transporters, CFTR and P-glycoprotein

To what extent do corresponding transmembrane helices in related integral membrane proteins have different membrane-insertion characteristics? Here, we compare, side-by-side, the membrane insertion characteristics of the 12 transmembrane helices in the adenosine triphosphate-binding cassette (ABC) transporters, P-glycoprotein (P-gp) and the cystic fibrosis transmembrane conductance regulator (CFTR). Our results show that 10 of the 12 CFTR transmembrane segments can insert independently into the ER membrane. In contrast, only three of the P-gp transmembrane segments are independently stable in the membrane, while the majority depend on the presence of neighboring loops and/or transmembrane segments for efficient insertion. Membrane-insertion characteristics can thus vary widely between related proteins.     

Membrane Insertion of Marginally Hydrophobic Transmembrane Helices Depends on Sequence Context

In mammalian cells, most integral membrane proteins are initially inserted into the endoplasmic reticulum membrane by the so-called Sec61 translocon. However, recent predictions suggest that many transmembrane helices (TMHs) in multispanning membrane proteins are not sufficiently hydrophobic to be recognized as such by the translocon. In this study, we have screened 16 marginally hydrophobic TMHs from membrane proteins of known three-dimensional structure. Indeed, most of these TMHs do not insert efficiently into the endoplasmic reticulum membrane by themselves. To test if loops or TMHs immediately upstream or downstream of a marginally hydrophobic helix might influence the insertion efficiency, insertion of marginally hydrophobic helices was also studied in the presence of their neighboring loops and helices. The results show that flanking loops and nearest-neighbor TMHs are sufficient to ensure the insertion of many marginally hydrophobic helices. However, for at least two of the marginally hydrophobic helices, the local interactions are not enough, indicating that post-insertional rearrangements are involved in the folding of these proteins.     

Stable intermediates determine proteins' primary unfolding sites in the presence of surfactants

Despite detailed knowledge of the overall structural changes and stoichiometries of surfactant binding, little is known about which protein regions constitute the preferred sites of attack for initial unfolding. Here we have exposed three proteins to limited proteolysis at anionic (SDS) and cationic (DTAC) surfactant concentrations corresponding to specific conformational transitions, using the surfactant‐robust broad‐specificity proteases Savinase and Alcalase. Cleavage sites are identified by SDS‐PAGE and N‐terminal sequencing. We observe well‐defined cleavage fragments, which suggest that flexibility is limited to certain regions of the protein. Cleavage sites for α‐lactalbumin and myoglobin correspond to regions identified in other studies as partially unfolded at low pH or in the presence of organic solvents. For Tnfn3, which does not form partially folded structures under other conditions, cleavage sites can be rationalized from the structure of the protein's folding transition state and the position of loops in the native state. Nevertheless, they are more sensitive to choice of surfactant and protease, probably reflecting a heterogeneous and fluctuating ensemble of partially unfolded structures. Thus, for proteins accumulating stable intermediates on the folding pathway, surfactants encourage the formation of these states, while the situation is more complex for proteins that do not form these intermediates. © 2008 Wiley Periodicals, Inc. Biopolymers 91: 221–231, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com