Total free energy analysis of fully hydrated proteins.

The total free energy of a hydrated biomolecule and its corresponding decomposition of energy and entropy provides detailed information about regions of thermodynamic stability or instability. The free energies of four hydrated globular proteins with different net charges are calculated from a molecular dynamics simulation, with the energy coming from the system Hamiltonian and entropy using multiscale cell correlation. Water is found to be most stable around anionic residues, intermediate around cationic and polar residues, and least stable near hydrophobic residues, especially when more buried, with stability displaying moderate entropy-enthalpy compensation. Conversely, anionic residues in the proteins are energetically destabilized relative to singly solvated amino acids, while trends for other residues are less clear-cut. Almost all residues lose intraresidue entropy when in the protein, enthalpy changes are negative on average but may be positive or negative, and the resulting overall stability is moderate for some proteins and negligible for others. The free energy of water around single amino acids is found to closely match existing hydrophobicity scales. Regarding the effect of secondary structure, water is slightly more stable around loops, of intermediate stability around ß strands and turns, and least stable around helices. An interesting asymmetry observed is that cationic residues stabilize a residue when bonded to its N-terminal side but destabilize it when on the C-terminal side, with a weaker reversed trend for anionic residues.

A model for pH coupling of the SARS-CoV-2 spike protein open/closed equilibrium.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causative agent of the coronavirus disease 2019 (COVID-19) pandemic, is thought to release its RNA genome at either the cell surface or within endosomes, the balance being dependent on spike protein stability, and the complement of receptors, co-receptors and proteases. To investigate possible mediators of pH-dependence, pKa calculations have been made on a set of structures for spike protein ectodomain and fragments from SARS-CoV-2 and other coronaviruses. Dominating a heat map of the aggregated predictions, three histidine residues in S2 are consistently predicted as destabilizing in pre-fusion (all three) and post-fusion (two of the three) structures. Other predicted features include the more moderate energetics of surface salt-bridge interactions and sidechain-mainchain interactions. Two aspartic acid residues in partially buried salt-bridges (D290-R273 and R355-D398) have pKas that are calculated to be elevated and destabilizing in more open forms of the spike trimer. These aspartic acids are most stabilized in a tightly closed conformation that has been observed when linoleic acid is bound, and which also affects the interactions of D614. The D614G mutation is known to modulate the balance of closed to open trimer. It is suggested that D398 in particular contributes to a pH-dependence of the open/closed equilibrium, potentially coupled to the effects of linoleic acid binding and D614G mutation, and possibly also A570D mutation. These observations are discussed in the context of SARS-CoV-2 infection, mutagenesis studies, and other human coronaviruses.

Investigating the Orientation of an Interfacially Adsorbed Monoclonal Antibody and Its Fragments Using Neutron Reflection.

Interfacial adsorption is a molecular process occurring during the production, purification, transport, and storage of antibodies, with a direct impact on their structural stability and subsequent implications on their bioactivities. While the average conformational orientation of an adsorbed protein can be readily determined, its associated structures are more complex to characterize. Neutron reflection has been used in this work to investigate the conformational orientations of the monoclonal antibody COE-3 and its Fab and Fc fragments at the oil/water and air/water interfaces. Rigid body rotation modeling was found to be suitable for globular and relatively rigid proteins such as the Fab and Fc fragments but less so for relatively flexible proteins such as full COE-3. Fab and Fc fragments adopted a 'flat-on' orientation at the air/water interface, minimizing the thickness of the protein layer, but they adopted a substantially tilted orientation at the oil/water interface with increased layer thickness. In contrast, COE-3 was found to adsorb in tilted orientations at both interfaces, with one fragment protruding into the solution. This work demonstrates that rigid-body modeling can provide additional insights into protein layers at various interfaces relevant to bioprocess engineering.

Chemometrics in Protein Formulation: Stability Governed by Repulsion and Protein Unfolding.

Therapeutic proteins can be challenging to develop due to their complexity and the requirement of an acceptable formulation to ensure patient safety and efficacy. To date, there is no universal formulation development strategy that can identify optimal formulation conditions for all types of proteins in a fast and reliable manner. In this work, high-throughput characterization, employing a toolbox of five techniques, was performed on 14 structurally different proteins formulated in 6 different buffer conditions and in the presence of 4 different excipients. Multivariate data analysis and chemometrics were used to analyze the data in an unbiased way. First, observed changes in stability were primarily determined by the individual protein. Second, pH and ionic strength are the two most important factors determining the physical stability of proteins, where there exists a significant statistical interaction between protein and pH/ionic strength. Additionally, we developed prediction methods by partial least-squares regression. Colloidal stability indicators are important for prediction of real-time stability, while conformational stability indicators are important for prediction of stability under accelerated stress conditions at 40 °C. In order to predict real-time storage stability, protein-protein repulsion and the initial monomer fraction are the most important properties to monitor.

Evolution and divergence of teleost adrenergic receptors: why sometimes 'the drugs don't work' in fish.

Adrenaline and noradrenaline, released as hormones and/or neurotransmitters, exert diverse physiological functions in vertebrates, and teleost fishes are widely used as model organisms to study adrenergic regulation; however, such investigations often rely on receptor subtype-specific pharmacological agents (agonists and antagonists; see Glossary) developed and validated in mammals. Meanwhile, evolutionary (phylogenetic and comparative genomic) studies have begun to unravel the diversification of adrenergic receptors (ARs) and reveal that whole-genome duplications and pseudogenization events in fishes results in notable distinctions from mammals in their genomic repertoire of ARs, while lineage-specific gene losses within teleosts have generated significant interspecific variability. In this Review, we visit the evolutionary history of ARs (including α1-, α2- and ß-ARs) to highlight the prominent interspecific differences in teleosts, as well as between teleosts and other vertebrates. We also show that structural modelling of teleost ARs predicts differences in ligand binding affinity compared with mammalian orthologs. To emphasize the difficulty of studying the roles of different AR subtypes in fish, we collate examples from the literature of fish ARs behaving atypically compared with standard mammalian pharmacology. Thereafter, we focus on specific case studies of the liver, heart and red blood cells, where our understanding of AR expression has benefited from combining pharmacological approaches with molecular genetics. Finally, we briefly discuss the ongoing advances in 'omics' technologies that, alongside classical pharmacology, will provide abundant opportunities to further explore adrenergic signalling in teleosts.

Computational Investigation of Mechanisms for pH Modulation of Human Chloride Channels.

Many transmembrane proteins are modulated by intracellular or extracellular pH. Investigation of pH dependence generally proceeds by mutagenesis of a wide set of amino acids, guided by properties such as amino-acid conservation and structure. Prediction of pKas can streamline this process, allowing rapid and effective identification of amino acids of interest with respect to pH dependence. Commencing with the calcium-activated chloride channel bestrophin 1, the carboxylate ligand structure around calcium sites relaxes in the absence of calcium, consistent with a measured lack of pH dependence. By contrast, less relaxation in the absence of calcium in TMEM16A, and maintenance of elevated carboxylate sidechain pKas, is suggested to give rise to pH-dependent chloride channel activity. This hypothesis, modulation of calcium/proton coupling and pH-dependent activity through the extent of structural relaxation, is shown to apply to the well-characterised cytosolic proteins calmodulin (pH-independent) and calbindin D9k (pH-dependent). Further application of destabilised, ionisable charge sites, or electrostatic frustration, is made to other human chloride channels (that are not calcium-activated), ClC-2, GABAA, and GlyR. Experimentally determined sites of pH modulation are readily identified. Structure-based tools for pKa prediction are freely available, allowing users to focus on mutagenesis studies, construct hypothetical proton pathways, and derive hypotheses such as the model for control of pH-dependent calcium activation through structural flexibility. Predicting altered pH dependence for mutations in ion channel disorders can support experimentation and, ultimately, clinical intervention.

Tadalafil Rescues the p.M325T Mutant of Best1 Chloride Channel.

Bestrophin 1 (Best1) is a chloride channel that localises to the plasma membrane of retinal pigment epithelium (RPE) cells. Mutations in the BEST1 gene are associated with a group of untreatable inherited retinal dystrophies (IRDs) called bestrophinopathies, caused by protein instability and loss-of-function of the Best1 protein. 4PBA and 2-NOAA have been shown to rescue the function, expression, and localisation of Best1 mutants; however, it is of interest to find more potent analogues as the concentration of the drugs required is too high (2.5 mM) to be given therapeutically. A virtual docking model of the COPII Sec24a site, where 4PBA has been shown to bind, was generated and a library of 1416 FDA-approved compounds was screened at the site. The top binding compounds were tested in vitro in whole-cell patch-clamp experiments of HEK293T cells expressing mutant Best1. The application of 25 µM tadalafil resulted in full rescue of Cl- conductance, comparable to wild type Best1 levels, for p.M325T mutant Best1 but not for p.R141H or p.L234V mutants.

Activation of p53: How phosphorylated Ser15 triggers sequential phosphorylation of p53 at Thr18 by CK1δ.

The N-terminal transactivation domain (TAD) of p53 is a disordered region with multiple phosphorylation sites. Phosphorylation at Thr18 is crucial for the release of p53 from its negative regulator, MDM2. In stressed cells, CK1δ is responsible for phosphorylating Thr18, but requires Ser15 to be phosphorylated. To understand the mechanistic underpinnings of this sequential phosphorylation, molecular modeling and molecular dynamics simulation studies of these phosphorylation events were carried out. Our models suggest that a positively charged region on CK1δ near the adenosine triphosphate (ATP) binding pocket, which is conserved across species, sequesters the negatively charged pSer15, thereby constraining the positioning of the rest of the peptide, such that the side chain of Thr18 is positioned close to the γ-phosphate of ATP. Furthermore, our studies show that the phosphorylated p53 TAD1 (p53pSer15) peptide binds more strongly to CK1δ than does p53. p53 adopts a helical structure when bound to CK1δ, which is lost upon phosphorylation at Ser15, thus gaining higher flexibility and ability to morph into the binding site. We propose that upon phosphorylation at Ser15 the p53 TAD1 peptide binds to CK1δ through an electrostatically driven induced fit mechanism resulting in a flanking fuzzy complex.

Protein-sol pKa: prediction of electrostatic frustration, with application to coronaviruses.

MOTIVATION: Evolution couples differences in ambient pH to biological function through protonatable groups, in particular, those that switch from buried to exposed and alter protonation state in doing so. We present a tool focusing on structure-based discovery and display of these groups. RESULTS: Since prediction of buried group pKas is computationally intensive, solvent accessibility of ionizable groups is displayed, from which the user can iteratively select pKa calculation centers. Results are color-coded, with emphasis on buried groups. Utility is demonstrated with benchmarking against known pH sensing sites in influenza virus hemagglutinin and in variants of murine hepatitis virus, a coronavirus. A pair of histidine residues, which are conserved in coronavirus spike proteins, are predicted to be electrostatically frustrated at acidic pH in both pre- and post-fusion conformations. We suggest that an intermediate expanded conformation at endosomal pH could relax the frustration, allowing histidine protonation and facilitating conformational conversion of coronavirus spike protein. AVAILABILITY AND IMPLEMENTATION: This tool is available at http://www.protein-sol.manchester.ac.uk/pka/.

How Do Electrostatic Perturbations of the Protein Affect the Bifurcation Pathways of Substrate Hydroxylation versus Desaturation in the Nonheme Iron-Dependent Viomycin Biosynthesis Enzyme?

The viomycin biosynthesis enzyme VioC is a nonheme iron and α-ketoglutarate-dependent dioxygenase involved in the selective hydroxylation of l-arginine at the C3-position for antibiotics biosynthesis. Interestingly, experimental studies showed that using the substrate analogue, namely, l-homo-arginine, a mixture of products was obtained originating from C3-hydroxylation, C4-hydroxylation, and C3-C4-desaturation. To understand how the addition of one CH2 group to a substrate can lead to such a dramatic change in selectivity and activity, we decided to perform a computational study using quantum mechanical (QM) cluster models. We set up a large active-site cluster model of 245 atoms that includes the oxidant with its first- and second-coordination sphere influences as well as the substrate binding pocket. The model was validated against experimental work from the literature on related enzymes and previous computational studies. Thereafter, possible pathways leading to products and byproducts were investigated for a model containing l-Arg and one for l-homo-Arg as substrate. The calculated free energies of activation predict product distributions that match the experimental observation and give a low-energy C3-hydroxylation pathway for l-Arg, while for l-homo-Arg, several barriers are found to be close in energy leading to a mixture of products. We then analyzed the origins of the differences in product distributions using thermochemical, valence bond, and electrostatic models. Our studies show that the C3-H and C4-H bond strengths of l-Arg and l-homo-Arg are similar; however, external perturbations from an induced electric field of the protein affect the relative C-H bond strengths of l-Arg dramatically and make the C3-H bond the weakest and guide the reaction to a selective C3-hydroxylation channel. Therefore, the charge distribution in the protein and the induced electric dipole field of the active site of VioC guides the l-Arg substrate activation to C3-hydroxylation and disfavors the C4-hydroxylation pathway, while this does not occur for l-homo-Arg. Tight substrate positioning and electrostatic perturbations from the second-coordination sphere residues in VioC also result in a slower overall reaction for l-Arg; however, they enable a high substrate selectivity. Our studies highlight the importance of the second-coordination sphere in proteins that position the substrate and oxidant, perturb charge distributions, and enable substrate selectivity.

Predictive approaches to guide the expression of recombinant vaccine targets in Escherichia coli: a case study presentation utilising Absynth Biologics Ltd. proprietary Clostridium difficile vaccine antigens.

Bacterial expression systems remain a widely used host for recombinant protein production. However, overexpression of recombinant target proteins in bacterial systems such as Escherichia coli can result in poor solubility and the formation of insoluble aggregates. As a consequence, numerous strategies or alternative engineering approaches have been employed to increase recombinant protein production. In this case study, we present the strategies used to increase the recombinant production and solubility of 'difficult-to-express' bacterial antigens, termed Ant2 and Ant3, from Absynth Biologics Ltd.'s Clostridium difficile vaccine programme. Single recombinant antigens (Ant2 and Ant3) and fusion proteins (Ant2-3 and Ant3-2) formed insoluble aggregates (inclusion bodies) when overexpressed in bacterial cells. Further, proteolytic cleavage of Ant2-3 was observed. Optimisation of culture conditions and changes to the construct design to include N-terminal solubility tags did not improve antigen solubility. However, screening of different buffer/additives showed that the addition of 1-15 mM dithiothreitol alone decreased the formation of insoluble aggregates and improved the stability of both Ant2 and Ant3. Structural models were generated for Ant2 and Ant3, and solubility-based prediction tools were employed to determine the role of hydrophobicity and charge on protein production. The results showed that a large non-polar region (containing hydrophobic amino acids) was detected on the surface of Ant2 structures, whereas positively charged regions (containing lysine and arginine amino acids) were observed for Ant3, both of which were associated with poor protein solubility. We present a guide of strategies and predictive approaches that aim to guide the construct design, prior to expression studies, to define and engineer sequences/structures that could lead to increased expression and stability of single and potentially multi-domain (or fusion) antigens in bacterial expression systems.

How Do Metal Ions Modulate the Rate-Determining Electron-Transfer Step in Cytochrome P450 Reactions?

Cytochrome P450 (CYP450) enzymes play important roles in maintaining human health and their reaction rates are dependent on the first electron transfer from the reduction partner. Interestingly, experimental work has shown that this step is highly influenced by the addition of metal ions. To understand the effect of external perturbations on the CYP450 first reduction step, we have performed a computational study with model complexes in the presence of metal and organic ions, solvent molecules, and an electric field. The results show that these medium-range interactions affect the driving force as well as electron-transfer rates dramatically. Based on the location, distance, and direction of the ions/electric field, the catalytic reaction rates are enhanced or impaired. Calculations on a large crystal structure with bonded alkali metal ions indicated inhibition patterns of the ions. Therefore, we predict that the active forms of the natural CYP450 isozymes will not have more than one alkali metal ion bound in the second-coordination sphere. As such, this study provides an insight into the activity of CYP450 enzymes and the effects of ions and electric field perturbations on their activity.

Model for Counterion Binding and Charge Reversal on Protein Surfaces.

The structural stability and solubility of proteins in liquid therapeutic formulations is important, especially since new generations of therapeutics are designed for efficacy before consideration of stability. We introduce an electrostatic binding model to measure the net charge of proteins with bound ions in solution. The electrostatic potential on a protein surface is used to separately group together acidic and basic amino acids into patches, which are then iteratively bound with oppositely charged counterions. This model is aimed toward formulation chemists for initial screening of a range of conditions prior to lab-work. Computed results compare well with experimental zeta potential measurements from the literature covering a range of solution conditions. Importantly, the binding model reproduces the charge reversal phenomenon that is observed with polyvalent ion binding to proteins and its dependence on ion charge and concentration. Intriguingly, protein sequence can be used to give similarly good agreement with experiment as protein structure, interpreted as resulting from the close proximity of charged side chains on a protein surface. Further, application of the model to human proteins suggests that polyanion binding and overcharging, including charge reversal for cationic proteins, is a general feature. These results add to evidence that addition of polyanions to protein formulations could be a general mechanism for modulating solution stability.

Arginine to Lysine Mutations Increase the Aggregation Stability of a Single-Chain Variable Fragment through Unfolded-State Interactions.

Increased protein solubility is known to correlate with an increase in the proportion of lysine over arginine residues. Previous work has shown that the aggregation propensity of a single-chain variable fragment (scFv) does not correlate with its conformational stability or native-state protein-protein interactions. Here, we test the hypothesis that aggregation is driven by the colloidal stability of partially unfolded states, studying the behavior of scFv mutants harboring single or multiple site-specific arginine to lysine mutations in denaturing buffers. In 6 M guanidine hydrochloride (GdmCl) or 8 M urea, repulsive protein-protein interactions were measured for the wild-type and lysine-enriched (4RK) scFvs reflecting weakened short-range attractions and increased excluded volume. In contrast to the arginine-enriched mutant (7KR) scFv exhibited strong reversible association. In 3 M GdmCl, the minimum concentration at which the scFvs were unfolded, the hydrodynamic radius of 4RK remained constant but increased for the wild type and especially for 7KR. Studies of single-point arginine to lysine scFv mutants indicated that the observed aggregation propensity of arginine under denaturing conditions was nonspecific. Interestingly, one such swap generated a scFv with especially low aggregation rates under low/high ionic strengths and denaturing buffers; molecular modeling identified hydrogen bonding between the arginine side chain and main chain peptide groups, stabilizing the structure. The arginine/lysine ratio is not routinely considered in biopharmaceutical scaffold design or current amyloid prediction methods. This work therefore suggests a simple method for increasing the stability of a biopharmaceutical protein against aggregation.

Surface patches on recombinant erythropoietin predict protein solubility: engineering proteins to minimise aggregation.

BACKGROUND: Protein solubility characteristics are important determinants of success for recombinant proteins in relation to expression, purification, storage and administration. Escherichia coli offers a cost-efficient expression system. An important limitation, whether for biophysical studies or industrial-scale production, is the formation of insoluble protein aggregates in the cytoplasm. Several strategies have been implemented to improve soluble expression, ranging from modification of culture conditions to inclusion of solubility-enhancing tags. RESULTS: Surface patch analysis has been applied to predict amino acid changes that can alter the solubility of expressed recombinant human erythropoietin (rHuEPO) in E. coli, a factor that has importance for both yield and subsequent downstream processing of recombinant proteins. A set of rHuEPO proteins (rHuEPO E13K, F48D, R150D, and F48D/R150D) was designed (from the framework of wild-type protein, rHuEPO WT, via amino acid mutations) that varied in terms of positively-charged patches. A variant predicted to promote aggregation (rHuEPO E13K) decreased solubility significantly compared to rHuEPO WT. In contrast, variants predicted to diminish aggregation (rHuEPO F48D, R150D, and F48D/R150D) increased solubility up to 60% in relation to rHuEPO WT. CONCLUSIONS: These findings are discussed in the wider context of biophysical calculations applied to the family of EPO orthologues, yielding a diverse range of calculated values. It is suggested that combining such calculations with naturally-occurring sequence variation, and 3D model generation, could lead to a valuable tool for protein solubility design.

Determination of Protein-Protein Interactions in a Mixture of Two Monoclonal Antibodies.

The coformulation of monoclonal antibody (mAb) mixtures provides an attractive route to achieving therapeutic efficacy where the targeting of multiple epitopes is necessary. Controlling and predicting the behavior of such mixtures requires elucidating the molecular basis for the self- and cross-protein-protein interactions and how they depend on solution variables. While self-interactions are now beginning to be well understood, systematic studies of cross-interactions between mAbs in solution do not exist. Here, we have used static light scattering to measure the set of self- and cross-osmotic second virial coefficients in a solution containing a mixture of two mAbs, mAbA and mAbB, as a function of ionic strength and pH. mAbB exhibits strong association at a low ionic strength, which is attributed to an electrostatic attraction that is enhanced by the presence of a strong short-ranged attraction of nonelectrostatic origin. Under all solution conditions, the measured cross-interactions are intermediate self-interactions and follow similar patterns of behavior. There is a strong electrostatic attraction at higher pH values, reflecting the behavior of mAbB. Protein-protein interactions become more attractive with an increasing pH due to reducing the overall protein net charges, an effect that is attenuated with an increasing ionic strength due to the screening of electrostatic interactions. Under moderate ionic strength conditions, the reduced cross-virial coefficient, which reflects only the energetic contribution to protein-protein interactions, is given by a geometric average of the corresponding self-coefficients. We show the relationship can be rationalized using a patchy sphere model, where the interaction energy between sites i and j is given by the arithmetic mean of the i-i and j-j interactions. The geometric mean does not necessarily apply to all mAb mixtures and is expected to break down at a lower ionic strength due to the nonadditivity of electrostatic interactions.

Energetic Fingerprinting of Ligand Binding to Paralogous Proteins: The Case of the Apoptotic Pathway.

Networks of biological molecules are key to cellular function, governing processes ranging from signal cascade propagation to metabolic pathway regulation. Genetic duplication processes give rise to sets of regulatory proteins that have evolved from a common ancestor, leading to interactomes whose dysregulation is often associated with disease. A better understanding of the determinants of specificity at interfaces shared by functionally related proteins is crucial to the rational design of novel pharmacotherapeutic agents. To this end, a comprehensive data set of drug and drug-like binders was assembled for the Bcl-xL and Bcl-2 antiapoptotic proteins-archetypal examples of regulatory systems governed by evolutionarily conserved protein-protein interactions. These were first used to derive a two-dimensional quantitative structure-activity relationship (2D QSAR) model, predicting ligand specificity for these homologous proteins. The strengths and weaknesses of high-throughput 2D QSAR were then compared and contrasted to those of theoretically rigorous thermodynamic integration calculations performed on 14 complexes of Bcl-xL-specific, Bcl-2-specific, and potent dual binders bound to the Bcl-xL and Bcl-2 proteins. We demonstrate that free energy calculations provide an added layer of essential information, which traditional QSAR cannot capture. Moreover, we show that protein energetic responses to different ligands, expressed as per-residue energy values, can be used to fingerprint the protein-ligand interaction, extending the framework of four-dimensional molecular dynamics/quantitative structure-activity relationships (4D-MD/QSAR) toward the facilitation of future drug design strategies.

Protein-Sol: a web tool for predicting protein solubility from sequence.

MOTIVATION: Protein solubility is an important property in industrial and therapeutic applications. Prediction is a challenge, despite a growing understanding of the relevant physicochemical properties. RESULTS: Protein-Sol is a web server for predicting protein solubility. Using available data for Escherichia coli protein solubility in a cell-free expression system, 35 sequence-based properties are calculated. Feature weights are determined from separation of low and high solubility subsets. The model returns a predicted solubility and an indication of the features which deviate most from average values. Two other properties are profiled in windowed calculation along the sequence: fold propensity, and net segment charge. The utility of these additional features is demonstrated with the example of thioredoxin. AVAILABILITY AND IMPLEMENTATION: The Protein-Sol webserver is available at http://protein-sol.manchester.ac.uk. CONTACT: jim.warwicker@manchester.ac.uk.

Charged residues next to transmembrane regions revisited: "Positive-inside rule" is complemented by the "negative inside depletion/outside enrichment rule".

BACKGROUND: Transmembrane helices (TMHs) frequently occur amongst protein architectures as means for proteins to attach to or embed into biological membranes. Physical constraints such as the membrane's hydrophobicity and electrostatic potential apply uniform requirements to TMHs and their flanking regions; consequently, they are mirrored in their sequence patterns (in addition to TMHs being a span of generally hydrophobic residues) on top of variations enforced by the specific protein's biological functions. RESULTS: With statistics derived from a large body of protein sequences, we demonstrate that, in addition to the positive charge preference at the cytoplasmic inside (positive-inside rule), negatively charged residues preferentially occur or are even enriched at the non-cytoplasmic flank or, at least, they are suppressed at the cytoplasmic flank (negative-not-inside/negative-outside (NNI/NO) rule). As negative residues are generally rare within or near TMHs, the statistical significance is sensitive with regard to details of TMH alignment and residue frequency normalisation and also to dataset size; therefore, this trend was obscured in previous work. We observe variations amongst taxa as well as for organelles along the secretory pathway. The effect is most pronounced for TMHs from single-pass transmembrane (bitopic) proteins compared to those with multiple TMHs (polytopic proteins) and especially for the class of simple TMHs that evolved for the sole role as membrane anchors. CONCLUSIONS: The charged-residue flank bias is only one of the TMH sequence features with a role in the anchorage mechanisms, others apparently being the leucine intra-helix propensity skew towards the cytoplasmic side, tryptophan flanking as well as the cysteine and tyrosine inside preference. These observations will stimulate new prediction methods for TMHs and protein topology from a sequence as well as new engineering designs for artificial membrane proteins.

Proofreading of substrate structure by the Twin-Arginine Translocase is highly dependent on substrate conformational flexibility but surprisingly tolerant of surface charge and hydrophobicity changes.

The Tat system transports folded proteins across the bacterial plasma membrane, and in Escherichia coli preferentially transports correctly-folded proteins. Little is known of the mechanism by which Tat proofreads a substrate's conformational state, and in this study we have addressed this question using a heterologous single-chain variable fragment (scFv) with a defined structure. We introduced mutations to surface residues while leaving the folded structure intact, and also tested the importance of conformational flexibility. We show that while the scFv is stably folded and active in the reduced form, formation of the 2 intra-domain disulphide bonds enhances Tat-dependent export 10-fold, indicating Tat senses the conformational flexibility and preferentially exports the more rigid structure. We further show that a 26-residue unstructured tail at the C-terminus blocks export, suggesting that even this short sequence can be sensed by the proofreading system. In contrast, the Tat system can tolerate significant changes in charge or hydrophobicity on the scFv surface; substitution of uncharged residues by up to 3 Lys-Glu pairs has little effect, as has the introduction of up to 5 Lys or Glu residues in a confined domain, or the introduction of a patch of 4 to 6 Leu residues in a hydrophilic region. We propose that the proofreading system has evolved to sense conformational flexibility and detect even very transiently-exposed internal regions, or the presence of unfolded peptide sections. In contrast, it tolerates major changes in surface charge or hydrophobicity.