Examination of pH dependency and orientation differences of membrane spanning alpha helices carrying a single or pair of buried histidine residues

The Physical Basis for pH Sensitivity in Biomolecular Structure and Function, With Application to the Spike Protein of SARS-CoV-2

Since pH sensitivity has a fundamental role in biology, much effort has been committed to establishing physical models to rationalize and predict pH dependence from molecular structures. Two of the key challenges are to accurately calculate ionizable group solvation and hydration and then to apply this modeling to all conformations relevant to the process in question. Explicit solvent methods coupled to molecular dynamics simulation are increasingly complementing lower resolution implicit solvent techniques, but equally, the scale of biological data acquisition leaves a role for high-throughput modeling. Additionally, determination of ranges of structures for a system allows sampling of key stages in solvation. In a review of the area, it is emphasized that pH sensors in biology beyond the most obvious candidate (histidine side chain, with an unshifted pK a near neutral pH) should be considered; that modeling can benefit from other concepts in bioinformatics, in particular modulation of interactions and function in families of homologs; and that it can also be beneficial to incorporate as many experimental structures as possible, to mitigate against small variations in conformation and to analyze larger, functional, conformational changes. These aspects are then demonstrated with new work on the spike protein of SARS-CoV-2, looking at the pH dependence of variants, including prediction of a change in the balance of locked, closed, and open forms at neutral pH for the Omicron variant spike protein.

Illuminating Disorder Induced by Glu in a Stable Arg-Anchored Transmembrane Helix

Membrane proteins are vital for biological function and are complex to study. Even in model peptide-lipid systems, the combined influence or interaction of pairs of chemical groups still is not well understood. Disordered proteins, whether in solution or near lipid membranes, are an emerging paradigm for the initiation and control of biological function. The disorder can involve molecular orientation as well as molecular folding. This paper reports an astonishing induction of disorder when one Glu residue is introduced into a highly stable 23-residue transmembrane helix. The parent helix is anchored by a single Arg residue, tilted at a well-defined angle with respect to the DOPC bilayer normal and undergoes rapid cone precession. When Glu is introduced two residues away from Arg, with 200° (or 160°) radial separation, the helix properties change radically to exhibit a multiplicity of three or more disordered states. The helix characteristics have been monitored by deuterium (²H) NMR spectroscopy as functions of the pH and lipid bilayer composition. The disordered multistate behavior of the (Glu, Arg)-containing helix varies with the lipid bilayer thickness and pH. The results highlight a fundamental induction of protein multistate properties by a single Glu residue in a lipid membrane environment.

Thermodynamics-Based Molecular Modeling of α-Helices in Membranes and Micelles

The Folding of Membrane-Associated Peptides (FMAP) method was developed for modeling α-helix formation by linear peptides in micelles and lipid bilayers. FMAP 2.0 identifies locations of α-helices in the amino acid sequence, generates their three-dimensional models in planar bilayers or spherical micelles, and estimates their thermodynamic stabilities and tilt angles, depending on temperature and pH. The method was tested for 723 peptides (926 data points) experimentally studied in different environments and for 170 single-pass transmembrane (TM) proteins with available crystal structures. FMAP 2.0 detected more than 95% of experimentally observed α-helices with an average error in helix end determination of around 2, 3, 4, and 5 residues per helix for peptides in water, micelles, bilayers, and TM proteins, respectively. Helical and nonhelical residue states were predicted with an accuracy from 0.86 to 0.96, and the Matthews correlation coefficient was from 0.64 to 0.88 depending on the environment. Experimental micelle- and membrane-binding energies and tilt angles of peptides were reproduced with a root-mean-square deviation of around 2 kcal/mol and 7°, respectively. The TM and non-TM states of hydrophobic and pH-triggered α-helical peptides in various lipid bilayers were reproduced in more than 95% of cases. The FMAP 2.0 web server (https://membranome.org/fmap) is publicly available to explore the structural polymorphism of antimicrobial, cell-penetrating, fusion, and other membrane-binding peptides, which is important for understanding the mechanisms of their biological activities.