\documentclass[10pt,letterpaper]{article}
%DIF LATEXDIFF DIFFERENCE FILE
%DIF DEL protrusion_ms_old.tex   Fri May  4 15:25:02 2018
%DIF ADD protrusion_ms.tex       Sun May  6 20:54:11 2018
\usepackage[top=0.85in,left=2.75in,footskip=0.75in]{geometry}
\usepackage{amsmath,amssymb}
\usepackage{changepage}
\usepackage[utf8x]{inputenc}
\usepackage{textcomp,marvosym}
\usepackage{cite}
\usepackage{nameref,hyperref}
\usepackage[right]{lineno}
\usepackage{microtype}
\DisableLigatures[f]{encoding = *, family = * }
% color can be used to apply background shading to table cells only
\usepackage[table]{xcolor}
\usepackage{array}
% create "+" rule type for thick vertical lines
\newcolumntype{+}{!{\vrule width 2pt}}

% create \thickcline for thick horizontal lines of variable length
\newlength\savedwidth
\newcommand\thickcline[1]{%
  \noalign{\global\savedwidth\arrayrulewidth\global\arrayrulewidth 2pt}%
  \cline{#1}%
  \noalign{\vskip\arrayrulewidth}%
  \noalign{\global\arrayrulewidth\savedwidth}%
}

% \thickhline command for thick horizontal lines that span the table
\newcommand\thickhline{\noalign{\global\savedwidth\arrayrulewidth\global\arrayrulewidth 2pt}%
\hline
\noalign{\global\arrayrulewidth\savedwidth}}


% Remove comment for double spacing
%\usepackage{setspace} 
%\doublespacing

% Text layout
\raggedright
\setlength{\parindent}{0.5cm}
\textwidth 5.25in 
\textheight 8.75in

% Bold the 'Figure #' in the caption and separate it from the title/caption with a period
% Captions will be left justified
\usepackage[aboveskip=1pt,labelfont=bf,labelsep=period,justification=raggedright,singlelinecheck=off]{caption}
\renewcommand{\figurename}{Fig}

% Use the PLoS provided BiBTeX style
\bibliographystyle{plos2015}

% Remove brackets from numbering in List of References
\makeatletter
\renewcommand{\@biblabel}[1]{\quad#1.}
\makeatother

% Leave date blank
\date{}

% Header and Footer with logo
\usepackage{lastpage,fancyhdr,graphicx}
\usepackage{epstopdf}
\pagestyle{myheadings}
\pagestyle{fancy}
\fancyhf{}
\setlength{\headheight}{27.023pt}
\lhead{\includegraphics[width=2.0in]{PLOS-submission.eps}}
\rfoot{\thepage/\pageref{LastPage}}
\renewcommand{\footrule}{\hrule height 2pt \vspace{2mm}}
\fancyheadoffset[L]{2.25in}
\fancyfootoffset[L]{2.25in}
\lfoot{\sf PLOS}

\usepackage{gensymb}
\usepackage{textgreek}
%% Include all macros below
\DeclareMathOperator{\R}{R}
\newcommand\OR[1]{\R{\left(#1\right)}}
\newcommand{\vv}[1]{\mathbf{#1}}

\DeclareUnicodeCharacter{2212}{-}

\def\inc/{\note{cite}}
\def\inr/{\note{ref}}
\def\check/{\note{check}}
\def\ea/{\textit{et al.}}
\def\ca/{C\textsubscript{\textalpha}}
\def\cb/{C\textsubscript{\textbeta}}
\def\matmet/{\emph{Materials and methods}}
%% END MACROS SECTION


%remove before prod
\usepackage{color}
\newcommand{\note}[1]{\textcolor{red}{\{#1\}}}
\usepackage{graphicx}


%
% Figures
% arrange differently
%
\def\periphcolor/{blue}
\def\tmcolor/{red}

\def\colwidth{6.4cm}
\def\halfwidth{4.35cm}
\def\quarterwidth{2.7cm}
%\def\bigfigside{11.4cm}


\newcommand{\FigModelIllustration}{
	\begin{figure}[ht]
		\centering
		\includegraphics[width=\colwidth]{Fig1.pdf}
%DIF 115c115
%DIF < 	\caption{The definitions of \emph{protrusions} and \emph{co-insertable protruding hydrophobes}. Panel A shows a cartoon representation of the C2 domain of human phospholipase A\textsubscript{2} (PDB ID: 1RLW), and panel B shows the convex hull for the same protein. All \ca/- and \cb/-atoms are shown as spheres. Hydrophobes are coloured orange. The convex hull for the \ca/- and \cb/-atomic coordinates is shown in blue. All spheres visible on the convex-hull representation are vertex residues. \emph{Protrusions} are defined as vertex residues with low local protein density, and shown as large spheres. \emph{Co-insertable protruding hydrophobes} are protruding hydrophobes that are adjacent vertices of the convex hull, they are shown connected by orange lines. Small black spheres are vertex residues that have high local density, and do therefore not meet the criteria for protrusions.}
%DIF -------
	\caption{The definitions of \emph{protrusions} and \emph{co-insertable protruding hydrophobes}. Panel A shows a cartoon representation of the C2 domain of human phospholipase A\textsubscript{2} (PDB ID: 1RLW), and panel B shows the convex hull for the same protein. Panel~C and~D shows the same for an example with an aliphatic helix in the ENTH-domain (PDB ID: 1H0A). All \ca/- and \cb/-atoms are shown as spheres. Hydrophobes are coloured orange. The convex hull for the \ca/- and \cb/-atomic coordinates is shown in blue. All spheres visible on the convex-hull representation are vertex residues. \emph{Protrusions} are defined as vertex residues with low local protein density, and shown as large spheres. \emph{Co-insertable protruding hydrophobes} are protruding hydrophobes that are adjacent vertices of the convex hull, they are shown connected by orange lines. Small black spheres are at vertex residues that have high local density, and do therefore not meet the criteria for protrusions.} %DIF > 
%DIF -------
	\label{FigModelIllustration}
	\end{figure}
}

\newcommand{\FigTmPeriphAaComparison}{
	\begin{figure}[ht]
		\centering
			\includegraphics[width=\colwidth]{Fig10.pdf}
	\caption{Large aliphatic and aromatic side chains are particularly over-represented on protrusion on peripheral proteins. Panel A shows the weighted fractions (Eq~\ref{eq:weighted_frec}) of hydrophobic amino acids on protrusions from peripheral proteins (\periphcolor/) and from proteins in the reference set (\tmcolor/). In panel B, the contrast between the two sets is quantified by the odds ratio (Eq~\ref{eq:odds_ratio}), so that positive values reflect higher frequencies in the set of peripheral proteins than in the reference set. More precisely the vertical axis denote $\ln\OR{\mathrm{peripheral}, \mathrm{reference}, \hat{F}_{aa, \mathrm{protrusion}}}$, with $aa$ representing each of the standard amino acids. Error bars are 95\% confidence intervals.}
	\label{FigTmPeriphAaComparison}
	\end{figure}
}

\newcommand{\FigProtrusionDensityHistograms}{
	\begin{figure}[ht]
		\centering
		
		\includegraphics[width=\colwidth]{Fig2.pdf}

%DIF 135c135
%DIF < 	\caption{Hydrophobes are more common on protruding positions in peripheral proteins, than in the reference set. The plots show frequencies of hydrophobes on surface amino acids, both on protrusions (A and C) and among all solvent exposed amino acids (B and D). Compare peripheral proteins (\periphcolor/) and the reference set (\tmcolor/). The horizontal axes show the mean fraction (Eq~\ref{eq:mean_frac_family}) of protrusions or solvent exposed amino-acids that are hydrophobic. The vertical axis shows the fraction of protein families.}
%DIF -------
	\caption{Hydrophobes are more common on protruding positions in peripheral proteins, than in the reference set. The plots show frequencies of hydrophobes on surface amino acids, both on protrusions (A, and C) and among all solvent exposed amino acids (B and D). Compare peripheral proteins (\periphcolor/) and the non-binding surfaces (\tmcolor/). The horizontal axes show the mean fraction (Eq~\ref{eq:mean_frac_family}) of protrusions or solvent exposed amino-acids that are hydrophobic. The vertical axis shows the fraction of protein families for each set. A-D shows the comparison between the data sets \emph{Peripheral Proteins} and \emph{Non-binding surfaces}, and E-H the comparison between the \emph{Peripheral PISA} and \emph{Reference proteins}.} %DIF > 
%DIF -------

	\label{FigProtrusionDensityHistograms}
	
	\end{figure}
}

\newcommand{\FigContainingProtrusions}{
	\begin{figure}[ht]
		\centering
		\includegraphics[width=\colwidth]{Fig5.pdf}
%DIF 146c146
%DIF < 		\caption{\emph{Co-insertable protruding hydrophobes} are common in peripheral proteins and rare in the reference set. The plot shows the occurrence of \emph{co-insertable protruding hydrophobes} on protein surfaces. Panels~A and~B show the weighted fraction (Eq~\ref{eq:weighted_existance_frec}) of proteins that have protruding hydrophobes, in the peripheral proteins (\periphcolor/) and the reference set (\tmcolor/). We have differentiated here between protrusions that have at least one co-insertable protruding hydrophobe (labeled ``Co-ins.''), and those that have not (labeled ``isolated''). The analysis is done separately for two groups of proteins according to the total number of protrusions on the protein surface ($[0,25\rangle$ in panel A, $[25,50\rangle$ in panel B). Panel~C shows the frequency distribution of the total number of protruding residues (``\# protrusions'') for all proteins. The selections analysed in panel~A and~B are found between the dashed lines in panel~C. Error bars in panel~A and~B are 95\% confidence intervals.}
%DIF -------
		\caption{\emph{Co-insertable protruding hydrophobes} are common in peripheral proteins and rare in the reference set. The plot shows the occurrence of \emph{co-insertable protruding hydrophobes} on protein surfaces. Panels~A,~B,~D, and~E show the weighted fraction (Eq~\ref{eq:weighted_existance_frec}) of proteins that have protruding hydrophobes, in the peripheral proteins (\periphcolor/) and the reference sets (\tmcolor/). We have differentiated here between protrusions that have at least one co-insertable protruding hydrophobe (labeled ``Co-ins.''), and those that have not (labeled ``isolated''). The analysis is done separately for two groups of proteins according to the total number of protrusions on the protein surface ($[0,25\rangle$ in panels~A and~D, $[25,50\rangle$ in panels~B and~E). Panels~C and~F shows the frequency distribution of the total number of protruding residues (``\# protrusions'') for all proteins. The selections analysed in panels~A,~B,~D, and~E are found between the dashed lines in panels~C andF. Error bars in panels~A,~B,~D, and~E are 95\% confidence intervals. Panels~A-C shows the comparison between \emph{Peripheral Proteins} and \emph{Non-binding surfaces}, panels~D-F, the comparison between \emph{Peripheral PISA} and \emph{Reference proteins}.} %DIF > 
%DIF -------
\label{FigContainingProtrusions}
	\end{figure}
}

\newcommand{\FigNeighboursVsOddsRatio}{
\begin{figure}[ht]
		\centering
		\includegraphics[width=\colwidth]{Fig3.pdf}
%DIF 155c155
%DIF < 		\caption{On peripheral proteins, protrusions in low density regions are more often hydrophobes, compared to the reference set. The plot shows the logarithm of the odds-ratio (Eq~\ref{eq:odds_ratio}) comparing the frequency of hydrophobes on \emph{vertex} residues in peripheral proteins and the reference set. Positive values reflect higher frequencies in the peripheral proteins. The horizontal axis shows the protein density $d$ around the protrusion, measured as the number of \ca/ and \cb/ atoms within $1 nm$. Vertex residues are all on the convex hull, but only the vertex residues with $d<22$ are protrusions. The leftmost bar with $d<7$ corresponds mostly to chain terminals. More precisely, the vertical axis shows $\OR{A, B,\hat{F}_{\mathrm{hydrophobe}|\mathrm{vertex}\cap l<d\leq u }}$ where $A$ denotes the peripheral proteins, $B$ the reference set, $l$ and $u$ denote the lower and upper limits of the ranges given on the vertical axis, and $d$ is the local protein density defined in~\matmet/. Error bars are 95\% confidence intervals.}
%DIF -------
		\caption{On \emph{Peripheral proteins}, protrusions in low density regions are more often hydrophobes, compared to the \emph{Non-binding surfaces}. The plot shows the logarithm of the odds-ratio (Eq~\ref{eq:odds_ratio}) comparing the frequency of hydrophobes on \emph{vertex} residues in peripheral proteins and the reference set. Positive values reflect higher frequencies in the peripheral proteins. The horizontal axis shows the protein density $d$ around the protrusion, measured as the number of \ca/ and \cb/ atoms within $1 nm$. Vertex residues are all on the convex hull, but only the vertex residues with $d<22$ are protrusions. The leftmost bar with $d<7$ corresponds mostly to chain terminals. More precisely, the vertical axis shows $\OR{A, B,\hat{F}_{\mathrm{hydrophobe}|\mathrm{vertex}\cap l<d\leq u }}$ where $A$ denotes the \emph{peripheral proteins}, $B$ the {Non-binding surfaces}, $l$ and $u$ denote the lower and upper limits of the ranges given on the vertical axis, and $d$ is the local protein density defined in~\matmet/. Error bars are 95\% confidence intervals.} %DIF > 
%DIF -------
		\label{FigNeighboursVsOddsRatio}
\end{figure}
}

\newcommand{\FigSSE}{
	\begin{figure}[ht]
		\centering
		\includegraphics[width=\colwidth]{Fig9.pdf}
	\caption{In peripheral proteins, hydrophobic protrusions are more frequent on turns, bends and $\alpha$-helices, compared to the reference set. Panel A shows the weighted number (Eq~\ref{eq:weighted_count}) of \emph{protruding hydrophobes} associated with the different types of secondary structure elements. We have differentiated between protrusions that have at least one co-insertable protruding hydrophobe (right, labeled ``Co-ins.''), and those that have not (left, labeled ``Isolated''). Panel B compares the weighted frequencies (Eq~\ref{eq:weighted_frec}) of hydrophobes on protruding secondary structures between the peripheral membrane proteins and the reference set, using the odds ratio (Eq~\ref{eq:odds_ratio}). Positive values reflect higher frequencies in the peripheral proteins. More precisely, panel A show the values $N_{\mathrm{hydrophobe}|\mathrm{protrusion}\cap sse }$, and panel B the comparisons $\OR{A,B, \hat{F}_{\mathrm{hydrophobe}|\mathrm{protrusion}\cap sse }}$ where $A$ denote the peripheral proteins, $B$ the reference set, and $sse$ specifies the secondary structures given in the color legend. Error bars in panel B are 95\% confidence intervals.}
	\label{FigSSE}
	\end{figure}
}

\newcommand{\FigHfConditionals}{
	\begin{figure}
		\centering
		\includegraphics[width=\colwidth]{Fig4.pdf}
%DIF 173c173
%DIF < 	\caption{The \emph{protruding hydrophobes} tend to be \emph{co-insertable} in the peripheral proteins. The tendency for protrusions to be co-insertable is quantified by the weighted frequency of co-insertion (Eq~\ref{eq:coins_freq}), and is compared between each data set and a null model using the odds ratio (Eq~\ref{eq:odds_ratio}). Positive values reflect higher frequencies of co-insertion than in the null model. More precisely, we show the comparisons $\OR{set,null,\hat{F}^{\mathrm{pair}}_{\mathrm{one},\mathrm{both}}}$, where $set$ represents the set of peripheral proteins (\periphcolor/) and the reference set (\tmcolor/), and $null$ represent their respective null models where hydrophobes have been relocated randomly among protrusions as described in~\matmet/. Error bars are 95\% confidence intervals.}
%DIF -------
	\caption{The \emph{protruding hydrophobes} tend to be \emph{co-insertable} in the peripheral proteins. The tendency for protrusions to be co-insertable is quantified by the weighted frequency of co-insertion (Eq~\ref{eq:coins_freq}), and is compared between each data set and a null model using the odds ratio (Eq~\ref{eq:odds_ratio}). Positive values reflect higher frequencies of co-insertion than in the null model. More precisely, we show the comparisons $\OR{set,null,\hat{F}^{\mathrm{pair}}_{\mathrm{one},\mathrm{both}}}$, where $set$ represents the set of peripheral proteins (\periphcolor/) and the reference set (\tmcolor/), and $null$ represent their respective null models where hydrophobes have been relocated randomly among protrusions as described in~\matmet/. Error bars are 95\% confidence intervals. Panel A shows the comparison between the data sets \emph{Peripheral Proteins} and \emph{Non-binding surfaces}, and B the comparison between the \emph{Peripheral PISA} and \emph{Reference proteins}.} %DIF > 
%DIF -------
	\label{FigHfConditionals}
	\end{figure}
}

\newcommand{\FigLargestFacetLocation}{
	\begin{figure}[ht]
		\centering
		\includegraphics[width=\colwidth]{Fig8.pdf}

	\caption{Comparing predictions based on protruding hydrophobes with the predicted IBS in the Orientation of Proteins in Membranes (OPM) database. The plots show the distributions of the median \emph{insertion coordinate} from OPM for \emph{the Likely Inserted Hydrophobe} in each family (measured at the \ca/-atom). Values greater than or equal to zero correspond to atoms positioned in the hydrophobic core or at the boundary. Hence insertion coordinate values close to zero indicate agreement with OPM. Panel~A (C) show data for the Likely Inserted Hydrophobes and panel~B (D) for a null model of randomly selected \emph{protruding} residues. Panel~C and~D show cumulative histograms (accumulated with decreasing insertion coordinates).}
	\label{FigLargestFacetLocation}
	\end{figure}
}

\newcommand{\FigCompExperiment}{
	\begin{figure}[ht]
		\centering
		\includegraphics[width=\colwidth]{Fig7.pdf}
%DIF 192c192
%DIF < 		\caption{Protruding hydrophobes predict experimentally verified binding sites. The figure shows comparisons of predicted binding residues (\emph{the Likely Inserted Hydrophobe}) with experimentally verified binding sites for a manually curated dataset of 24 proteins (listed in Table S2). The vertical axis corresponds to values of the angle (Eq~\ref{eq:orientation}) comparing the two vectors connecting the center of the protein with either the predicted or known binding sites. Smaller angles imply better agreement between prediction and experiment. Asterisks ($*$) mark proteins where the Likely Inserted Hydrophobe is an amino acid experimentally identified to be interacting with the membrane. The grey boxplots show the distribution of angles when the known binding site residues are compared to all protruding amino acids on the protein. 1iaz is analysed in its soluble monomeric state, while it forms a transmembrane pore upon oligomerisation. The structure of the C-type lysozyme (PDBID 1f6s) has no identified protruding hydrophobes and is marked with a cross at 180\degree. Interestingly, while our analysis is performed on its crystallised form, it is known to bind membranes in a molten globule state.}
%DIF -------
		\caption{Protruding hydrophobes predict experimentally verified binding sites. The figure shows comparisons of predicted binding residues (\emph{the Likely Inserted Hydrophobe}) with experimentally verified binding sites for a manually curated dataset of 24 proteins (listed in Table S2). The vertical axis corresponds to values of the angle (Eq~\ref{eq:orientation}) comparing the two vectors connecting the center of the protein with either the predicted or known binding sites. Smaller angles imply better agreement between prediction and experiment. Asterisks ($*$) mark proteins where the Likely Inserted Hydrophobe is an amino acid experimentally identified to be interacting with the membrane. The grey boxplots show the distribution of angles when the known binding site residues are compared to all protruding amino acids on the protein. 1iaz is analysed in its soluble monomeric state, while it forms a transmembrane pore upon oligomerisation. The structure of the Bovine $\alpha$-lactalbumin (PDBID: 1F6S) has no identified protruding hydrophobes and is marked with a cross at 180\degree.} %DIF > 
%DIF -------
		\label{FigCompExperiment}
	\end{figure}
}

\newcommand{\FigDomainExamples}{
\begin{figure}[t]
		\centering
		\includegraphics[width=\colwidth]{Fig6.pdf}
%DIF 201c201
%DIF < 		\caption{Protruding hydrophobes are found on the membrane binding sites of well known membrane binding domains. The figure shows the convex hull (in blue) of the \ca/ and \cb/-atoms of selected peripheral membrane binding domains. The \cb/-atoms of \emph{the Likely Inserted Hydrophobe} are shown as orange spheres and \cb/-atoms of experimentally identified membrane-binding residues as gray spheres. The Likely Inserted Hydrophobe is an amino acid that has been experimentally verified to be a membrane binding residue for A, B, D and F. For C and E the Likely Inserted Hydrophobe is located in the same area as the residues identified by experiments. \textbf{A}: C2  domain of human phospholipase A2 (1RLW~\cite{Perisic:1998tv});  \textbf{B}:  PX domain of P40PHOX (1H6H~\cite{Bravo:2001ua}); \textbf{C}: snake phospholipase A2 (1POA~\cite{Scott:1990wf}); \textbf{D} : C1 domain of protein kinase C delta (1PTR~\cite{Zhang:1995vz}); \textbf{E}: Epsin ENTH domain (1H0A~\cite{Ford:2002if});\textbf{F}: FYVE domain of yeast vacuolar protein sorting-associated protein 27 (1VFY~\cite{Misra:1999uo}).}
%DIF -------
		\caption{Protruding hydrophobes are found on the membrane binding sites of well known membrane binding domains. The figure shows the convex hull (in blue) of the \ca/ and \cb/-atoms of selected peripheral membrane binding domains. The \cb/-atoms of \emph{the Likely Inserted Hydrophobe} are shown as orange spheres and \cb/-atoms of experimentally identified membrane-binding residues as gray spheres. The Likely Inserted Hydrophobe is an amino acid that has been experimentally verified to be a membrane binding residue for A, B, D and F. For C and E the Likely Inserted Hydrophobe is located in the same area as the residues identified by experiments. \textbf{A}: C2  domain of human phospholipase A2 (PDBID: 1RLW~\cite{Perisic:1998tv});  \textbf{B}:  PX domain of P40PHOX (PDBID: 1H6H~\cite{Bravo:2001ua}); \textbf{C}: snake phospholipase A2 (PDBID: 1POA~\cite{Scott:1990wf}); \textbf{D} : C1 domain of protein kinase C delta (PDBID: 1PTR~\cite{Zhang:1995vz}); \textbf{E}: Epsin ENTH domain (PDBID: 1H0A~\cite{Ford:2002if}); \textbf{F}: FYVE domain of yeast vacuolar protein sorting-associated protein 27 (PDBID: 1VFY~\cite{Misra:1999uo}).} %DIF > 
%DIF -------
		\label{FigDomainExamples}
%DIF 203a203-211
	\end{figure} %DIF > 
} %DIF > 
 %DIF > 
\newcommand{\FigOligomerComparison}{ %DIF > 
	\begin{figure}[t] %DIF > 
		\centering %DIF > 
		\includegraphics[width=\colwidth]{Fig11.pdf} %DIF > 
	\caption{Differences in number of polypeptide chains between the protein models present in the \emph{Peripheral Proteins} (quaternary structure model from OPM) and the models in \emph{Peripheral PISA} (quaternary structure model predicted by PISA). The difference is calculated for each of the PDB IDs occurring in both datasets. When more chains are present in the PISA models, The difference (horisontal axis) is negative.} %DIF > 
	\label{FigOligomerComparison} %DIF > 
%DIF -------
	\end{figure}
}
%DIF PREAMBLE EXTENSION ADDED BY LATEXDIFF
%DIF UNDERLINE PREAMBLE %DIF PREAMBLE
\RequirePackage[normalem]{ulem} %DIF PREAMBLE
\RequirePackage{color}\definecolor{RED}{rgb}{1,0,0}\definecolor{BLUE}{rgb}{0,0,1} %DIF PREAMBLE
\providecommand{\DIFaddtex}[1]{{\protect\color{blue}\uwave{#1}}} %DIF PREAMBLE
\providecommand{\DIFdeltex}[1]{{\protect\color{red}\sout{#1}}}                      %DIF PREAMBLE
%DIF SAFE PREAMBLE %DIF PREAMBLE
\providecommand{\DIFaddbegin}{} %DIF PREAMBLE
\providecommand{\DIFaddend}{} %DIF PREAMBLE
\providecommand{\DIFdelbegin}{} %DIF PREAMBLE
\providecommand{\DIFdelend}{} %DIF PREAMBLE
%DIF FLOATSAFE PREAMBLE %DIF PREAMBLE
\providecommand{\DIFaddFL}[1]{\DIFadd{#1}} %DIF PREAMBLE
\providecommand{\DIFdelFL}[1]{\DIFdel{#1}} %DIF PREAMBLE
\providecommand{\DIFaddbeginFL}{} %DIF PREAMBLE
\providecommand{\DIFaddendFL}{} %DIF PREAMBLE
\providecommand{\DIFdelbeginFL}{} %DIF PREAMBLE
\providecommand{\DIFdelendFL}{} %DIF PREAMBLE
%DIF END PREAMBLE EXTENSION ADDED BY LATEXDIFF
%DIF PREAMBLE EXTENSION ADDED BY LATEXDIFF
%DIF HYPERREF PREAMBLE %DIF PREAMBLE
\providecommand{\DIFadd}[1]{\texorpdfstring{\DIFaddtex{#1}}{#1}} %DIF PREAMBLE
\providecommand{\DIFdel}[1]{\texorpdfstring{\DIFdeltex{#1}}{}} %DIF PREAMBLE
%DIF END PREAMBLE EXTENSION ADDED BY LATEXDIFF

\begin{document}
\vspace*{0.2in}

% Title must be 250 characters or less.
\begin{flushleft}
{\Large
\textbf\newline{A model for hydrophobic protrusions on peripheral membrane proteins} % Please use "sentence case" for title and headings (capitalize only the first word in a title (or heading), the first word in a subtitle (or subheading), and any proper nouns).
}
\newline

Edvin Fuglebakk\textsuperscript{1,2\textcurrency*},
Nathalie Reuter\textsuperscript{1,2**}
\\
\bigskip
\textbf{1} Computational Biology Unit, University of Bergen, Pb7803, 5020 Bergen, Norway
\\
\textbf{2} Department of Molecular Biology, University of Bergen, Pb7803, 5020 Bergen, Norway
\\
\bigskip


% Current address notes
\textcurrency Current Address: Institute of Marine Research, Norway % change symbol to "\textcurrency a" if more than one current address note
% \textcurrency b Insert second current address 
% \textcurrency c Insert third current address


% Use the asterisk to denote corresponding authorship and provide email address in note below.
* Edvin.Fuglebakk@hi.no
\DIFaddbegin 

\DIFaddend ** Nathalie.Reuter@cbu.uib.no

\end{flushleft}
\newpage

%\pdfsuppresswarningpagegroup1

\section*{Abstract}
With remarkable spatial and temporal specificities, peripheral membrane proteins bind to biological membranes. They do this without compromising solubility of the protein, and their binding sites are not easily distinguished. Prototypical peripheral membrane binding sites display a combination of patches of basic and hydrophobic amino acids that are also frequently present on other protein surfaces. The purpose of this contribution is to identify simple but essential components for membrane binding, through structural criteria that distinguish exposed hydrophobes at membrane binding sites from those that are frequently found on any protein surface. We formulate the concepts of \textit{protruding hydrophobes} and \textit{co-insertability} and have analysed more than 300 families of proteins that are classified as peripheral membrane binders. We find that this structural motif strongly discriminates the surfaces of membrane-binding and non-binding proteins. Our model constitutes a novel formulation of a structural pattern for membrane recognition and emphasizes the importance of subtle structural properties of hydrophobic membrane binding sites.

\section*{Author summary}
Peripheral membrane proteins bind cellular membranes transiently, and are otherwise soluble proteins. As the interaction between proteins and membranes happen at cellular interfaces, they are naturally involved in important interfacial processes, such as recognintion, signalling and trafficking. Commonly their binding sites are also soluble, and their binding mechanisms are subtle and hard to recognize. This complicates our conceptual and quantitative models of how they function, and makes binding site detection difficult. It is therefore of great interest to discover traits that are common between these binding sites, and that distinguishes them from other protein surfaces. In this work we identify simple and general structural features that facilitates membrane recognition by soluble proteins. We find that these motifs are highly over-represented on peripheral membrane proteins.

\linenumbers

\section*{Introduction}
Biological membranes are ancient and crucial components in the organisation of life. Not only do they define the boundaries of cells and organelles, but they are central to a myriad \DIFdelbegin \DIFdel{of }\DIFdelend protein-protein and protein-lipid interactions instrumental in numerous pathways~\cite{Kutateladze:2010jz, Vogler:2008ch,Cullen:2008kv,Inaba:2016cb, Itoh:2005cp}. Besides the embedded transmembrane proteins and receptors, a number of soluble proteins interact transiently with the surface of cellular and organellar membranes achieving remarkable spatial and temporal specificities. These proteins are referred to as peripheral proteins and their membrane-binding site as interfacial binding site or IBS. \DIFdelbegin \DIFdel{Peripheral proteins include well-known }\DIFdelend \DIFaddbegin \DIFadd{Some peripheral proteins contain }\DIFaddend lipid-binding domains \DIFdelbegin \DIFdel{such as C2 and FYVE }\DIFdelend \DIFaddbegin \DIFadd{of well-known classes, such as C2-domains and FYVE-domains }\DIFaddend that confer larger proteins the ability to bind membranes~\cite{Lemmon:2008kt, Cho:2005bq}. Many other proteins such as lipid-processing enzymes, endogenous or secreted by pathogens are also included in this definition.

Unlike protein-protein or protein-ligand interactions, interfacial binding sites of peripheral proteins are poorly characterized in terms of amino acid composition and structural patterns. Embedded and transmembrane proteins contain well defined regions of hydrophobic surface, clearly identifying their membrane interacting segments. This is \DIFdelbegin \DIFdel{seldomly }\DIFdelend \DIFaddbegin \DIFadd{seldom }\DIFaddend the case for peripheral membrane proteins. Currently the prototypical peripheral membrane binding site is described as displaying a combination of basic and hydrophobic amino acids\DIFaddbegin \DIFadd{~\mbox{%DIFAUXCMD
\cite{Cho:2005bq, Johnson:1999wc} }%DIFAUXCMD
}\DIFaddend . Attempts to characterize the energetics of membrane binding has mostly focused on electrostatic complementarity of peripheral proteins with the charged surfaces of membrane~\cite{MulgrewNesbitt:2006dy}, rather than on the desolvation of hydrophobes which is more difficult to isolate in theoretical treatments. Nevertheless the predictive power of implicit membrane models in the prediction of membrane binding sites is a strong indication of the importance of the hydrophobic effect \cite{Lazaridis:2003wu} in peripheral membrane binding. For example, Lomize \DIFdelbegin \DIFdel{et al. }\DIFdelend \DIFaddbegin \ea\DIFadd{/ }\DIFaddend could correctly identify the experimentally known IBS of 53 peripheral peptides and proteins using a model that includes only hydrophobic, desolvation and ionization energy terms \cite{Lomize:2007gm}. Yet in order to assert the generality of a protein-membrane binding mechanism, it is not enough to demonstrate its validity for a selected set of true positives, but it is also important to evaluate it on a control dataset. 

As both small hydrophobic patches and charged residues are frequently present on protein surfaces it is challenging to distinguish membrane binding sites from the rest of the peripheral membrane proteins surface solely relying on amino acid composition. There are indications that structural considerations may allow signatures of membrane interacting hydrophobes to be defined. Terms like  \textit{hydrophobic spikes}~\cite{Gilbert:2002ja, Gamsjaeger:2005jd} and \textit{protruding loops}~\cite{Lomize:2007gm} have been used to describe membrane binding sites, prompting the idea of hydrophobes protruding from the protein globule. A close look at amphipathic helices, also motivates the concept of protruding hydrophobes. Amphipathic helices are characteristic of membrane-binding peptides and proteins. When such membrane binding helices exist, they are often found lining a protein, forming a cylindrical protrusion from the globule (e.g. ENTH domain of Epsin, PDBID: 1H0A~\cite{Ford:2002if}\DIFaddbegin \DIFadd{, shown in Fig~\ref{FigModelIllustration}~C and~D}\DIFaddend ). Yet, no generalization of protruding membrane binding sites has been proposed for peripheral membrane proteins.

The purpose of this contribution is to identify structural characteristics that distinguish exposed hydrophobes at membrane binding sites from those that are frequently found on any protein surface. We propose a simple definition that formalizes the concept of protruding hydrophobes, and which can be easily computed from the protein structure. This definition allows us to systematically investigate to what extent protruding hydrophobes are found on both binding and non-membrane-binding surfaces, and to identify structural criteria for recognizing exposed hydrophobes that are likely to be important for membrane binding.

A major obstacle in developing general association models for peripheral membrane proteins is the scarcity of experimentally verified binding sites, and detailed descriptions of binding orientations. So far, computational studies on the role of hydrophobes on membrane binding sites have been based on relatively small sets of proteins with known binding sites~\cite{Lomize:2007gm, BalaliMood:2009cv, Lazaridis:2003wu}. To get around this problem and to leverage the large number of proteins for which membrane binding has been identified without a detailed characterisation of the IBS, we perform a comparative statistical analysis of \DIFdelbegin \DIFdel{overall }\DIFdelend protein surfaces. Given \DIFdelbegin \DIFdel{a classification }\DIFdelend \DIFaddbegin \DIFadd{classifications }\DIFaddend of proteins that \DIFdelbegin \DIFdel{separates membrane bindersfrom non-binders}\DIFdelend \DIFaddbegin \DIFadd{identifies membrane binders}\DIFaddend , we compare peripheral membrane proteins with \DIFdelbegin \DIFdel{non-membrane-binding reference surfaces }\DIFdelend \DIFaddbegin \DIFadd{non-binding surfaces and reference proteins}\DIFaddend . With this we can extend our analysis to hundreds of protein families, rather than the few dozens for which binding sites have been partially identified by experiments.

With our simple definition of structural protrusions, we perform a statistical analysis of protruding hydrophobes in a large protein structure dataset and our results support their general role in membrane association. We find that protruding hydrophobes can be used to strongly discriminate protein surfaces involved in membrane binding from those that are not. Hydrophobes are much more frequent on protruding sites of peripheral membrane proteins than in the reference dataset, and they have a strong tendency to cluster on positions that can simultaneously interact with the membrane. 

\section*{Results and discussion}
\FigModelIllustration

Our formalisation of the concept of protruding amino acids is illustrated in Fig~\ref{FigModelIllustration} and described in details in~\matmet/. In short, it relies on firstly identifying the convex hull (in blue in Fig~\ref{FigModelIllustration}) of a coarse-grained protein model consisting of only \DIFaddbegin \DIFadd{its }\DIFaddend \ca/- and \cb/-atoms. We then identify amino acids located at vertices of the convex hull which intuitively are good candidates to be inserted into a membrane without inserting other residues, and without deforming the protein backbone. The model thus implicitly assumes that (1) proteins interact with the membrane without appreciable conformational change, or prior to such change and (2) that the membrane is locally flat, which is a valid approximation in most cases\DIFaddbegin \DIFadd{~\mbox{%DIFAUXCMD
\cite{Yang:2015fd}}%DIFAUXCMD
}\DIFaddend . In order to single out the amino acids that are most exposed to solvent, we \DIFdelbegin \DIFdel{single out }\DIFdelend \DIFaddbegin \DIFadd{identify }\DIFaddend amino acids (vertices) in regions of low protein density, characterized as having a low number of neighboring atoms. Solvent accessibility is a necessary condition for the hydrophobic effect to contribute to binding. In addition, regions of low local protein density are also likely to cause less disruption of lipid packing upon membrane insertion. \DIFaddbegin \DIFadd{The model was formulated based on inspection of the eight proteins included as supporting information (Table~S1).
}\DIFaddend 

In what follows, we present results of the application of this model to characterise %hydropathic 
hydrophobic properties of protrusions in peripheral membrane proteins. We do this by comparing \DIFdelbegin \DIFdel{a dataset of }\DIFdelend peripheral membrane proteins \DIFaddbegin \DIFadd{to a reference set of non-binding protein surface segments, }\DIFaddend and a reference set of \DIFdelbegin \DIFdel{protein surfaces not interacting }\DIFdelend \DIFaddbegin \DIFadd{typical protein structures. The reference set of non-binding surface segments (}\emph{\DIFadd{Non-binding surfaces}}\DIFadd{) is constructed from the solvent exposed regions of trans-membrane proteins and is intended to represent structures that do not interact }\DIFaddend with the membrane\DIFdelbegin \DIFdel{, as described in }\DIFdelend \DIFaddbegin \DIFadd{. The reference set of typical proteins (}\emph{\DIFadd{Reference proteins}}\DIFadd{) is constructed from a protein classification excluding proteins that are classified as membrane-interacting, and is intended to represent typical protein structures, with some unknown frequency of peripheral membrane binders that are not classified as such. For consistent treatment of non-redundancy and quaternary structure modeling, these sets are compared with two variants of a set of peripheral membrane proteins (}\emph{\DIFadd{Peripheral proteins}} \DIFadd{and }\emph{\DIFadd{Peripheral PISA}}\DIFadd{). These data sets are described in detail in }\DIFaddend \matmet/. \DIFaddbegin \DIFadd{Apart from some difference in classification, the main difference between the two set of peripheral membrane proteins is the modeling of quaternary structure, which our analysis might be sensitive to. The difference is illustrated in Fig~\ref{FigOligomerComparison}.
}\DIFaddend 

\subsection*{Protruding hydrophobes in a dataset of peripheral membrane proteins}

\FigProtrusionDensityHistograms
First we calculated the frequency of hydrophobes on protrusions in peripheral \DIFdelbegin \DIFdel{proteins }\DIFdelend \DIFaddbegin \DIFadd{protein }\DIFaddend families and compared it to the reference \DIFdelbegin \DIFdel{dataset}\DIFdelend \DIFaddbegin \DIFadd{datasets}\DIFaddend . In Fig~\ref{FigProtrusionDensityHistograms}, we observe a stark contrast between the set of peripheral proteins and the \DIFdelbegin \DIFdel{reference set (Figures}\DIFdelend \DIFaddbegin \DIFadd{non-binding surfaces (Compare figures}\DIFaddend ~\ref{FigProtrusionDensityHistograms}~A \DIFdelbegin \DIFdel{and}\DIFdelend \DIFaddbegin \DIFadd{to}\DIFaddend ~\ref{FigProtrusionDensityHistograms}~C). Hydrophobes occur with high frequency and in almost all families on protrusions of peripheral proteins. In the reference set on the other hand, hydrophobes on protrusions are much less tolerated, reflected by a histogram mode of zero. \DIFdelbegin \DIFdel{This }\DIFdelend \DIFaddbegin \DIFadd{While less pronounced, the distinction is also clear for the comparison with reference proteins (Compare figures~\ref{FigProtrusionDensityHistograms}~E to~\ref{FigProtrusionDensityHistograms}~G). Qualitatively, the frequency of hydrophobes on protrusions are similar in the two reference sets, but the sets of peripheral proteins differ somewhat, suggesting some sensitivity to quaternary structure modeling. For both comparisons, however, this }\DIFaddend trend is specific for protruding positions, and does not reflect a general difference in composition of surface exposed amino-acids between the data sets as shown by plots in Figures~\ref{FigProtrusionDensityHistograms}~B\DIFdelbegin \DIFdel{and}\DIFdelend \DIFaddbegin \DIFadd{,}\DIFaddend ~\ref{FigProtrusionDensityHistograms}~D\DIFdelbegin \DIFdel{. Indeed }\DIFdelend \DIFaddbegin \DIFadd{,~\ref{FigProtrusionDensityHistograms}~F and~\ref{FigProtrusionDensityHistograms}~H. Indeed, }\DIFaddend if we consider the frequency of hydrophobes on all solvent exposed residues, the distributions look quite similar with both sets having histogram modes close to $0.2$. This value is in agreement with the fraction of the surface of globular proteins typically reported to be hydrophobic (for instance 0.19 in Ref.~\cite{Miller:1987un}). The \DIFdelbegin \DIFdel{surfaces of the references set }\DIFdelend \DIFaddbegin \emph{\DIFadd{Non-binding surfaces}} \DIFaddend are in some cases very small, due to the way we ensure that these surfaces are not interacting with the membrane (see \matmet/). While these small surfaces are relevant samples for calculating average frequencies, the fraction of hydrophobes on such surfaces can take more extreme values (close to zero or $1$). For this reason the tails of the histograms for \DIFdelbegin \DIFdel{the }\DIFdelend \DIFaddbegin \DIFadd{this }\DIFaddend reference set are somewhat fatter than those for the peripheral membrane proteins.

\FigNeighboursVsOddsRatio
%As we define protrusions to be vertex residues that have at most $n=22$  within $c=1~\mathrm{nm}$.} 
Given the nature of our model the differences presented in Fig~\ref{FigProtrusionDensityHistograms} are naturally ascribed to two factors; the accessibility of amino acids compared to other regions of the protein (they are vertices of the convex hull) and their low local protein density $d$ defined as the number of \DIFdelbegin \DIFdel{neighbouring }\DIFdelend \DIFaddbegin \DIFadd{neighboring }\DIFaddend \ca/- or \cb/-atoms (Cf. definition in \matmet/). We here explore the dependence of this difference on $d$. In Fig~\ref{FigNeighboursVsOddsRatio} we show the difference between frequencies of hydrophobes in peripherals and \DIFdelbegin \DIFdel{reference data sets }\DIFdelend \DIFaddbegin \DIFadd{the non-binding surfaces }\DIFaddend for different ranges of the local protein density $d$. The leftmost bar ($0 \le d \le 6$) corresponds to chain terminals. The other bars corresponding to ranges covered by our definition of protruding residues ($ 7 \le d < 22$) show that hydrophobic residues are more frequently found at vertex residues with low local protein density in the peripheral proteins. This also serves as an \emph{a posteriori} justification for constricting our definition of protrusions to amino-acids with $d<22$.

\FigHfConditionals
Assuming that the over-representation of hydrophobes on protrusions in peripheral membrane proteins stems from actual membrane binding sites, \DIFdelbegin \DIFdel{one can }\DIFdelend \DIFaddbegin \DIFadd{we }\DIFaddend expect those proteins to have more than one hydrophobic protrusion\DIFaddbegin \DIFadd{.
}\DIFaddend %frequently occur on positions that can be simultaneously inserted in the membrane.
We estimated the tendency of \DIFdelbegin \DIFdel{each hydrophobic protrusion }\DIFdelend \DIFaddbegin \DIFadd{hydrophobic protrusions }\DIFaddend to be \emph{co-insertable} by calculating the weighted frequency of co-insertion (Eq~\ref{eq:coins_freq}) (Cf \matmet/) for \DIFdelbegin \DIFdel{both }\DIFdelend \DIFaddbegin \DIFadd{all }\DIFaddend datasets (Fig~\ref{FigHfConditionals}). We note that peripheral membrane proteins do indeed tend to have hydrophobes on co-insertable protrusions to a significantly larger extent than what would be expected from randomly scattering hydrophobes among protruding positions. This tendency is much lower for the \DIFdelbegin \DIFdel{reference set}\DIFdelend \DIFaddbegin \emph{\DIFadd{Non-binding surfaces}}\DIFaddend , even when considering the extremities of the error bars, which are wide precisely because there are very few protruding hydrophobes in this set. \DIFaddbegin \DIFadd{In the }\emph{\DIFadd{Reference proteins}}\DIFadd{, the analysis indicate that co-insertability is more common than in the null model, but far less so than in the Peripheral proteins.
}\DIFaddend 

\FigContainingProtrusions
We further explore the degree of co-insertability of the hydrophobic protrusions present in our \DIFdelbegin \DIFdel{dataset of peripheral proteins and in the reference dataset}\DIFdelend \DIFaddbegin \DIFadd{datasets}\DIFaddend . We seek to evaluate to what extent co-insertable hydrophobic protrusions can be used to discriminate likely peripheral membrane binders from other proteins. Fig~\ref{FigContainingProtrusions} shows the fraction of proteins in each dataset that have at least one pair of co-insertable hydrophobic protrusions (labelled \emph{Co-ins.}) and the fraction of proteins that have at least one \emph{isolated} hydrophobic protrusion (i.e a protrusion that does not satisfy the criteria that define \emph{co-insertability}). While we do see some discrimination between the data sets in the case of isolated protruding hydrophobes, the co-insertable ones prove to be very strong indicators of which proteins surfaces are membrane binding. As the coincidental occurrence of such properties increase with the size of the protein surface, we have grouped the proteins by total number of surface protrusions (regardless of hydropathic properties). We do however see no appreciable difference between the proteins of size $0-25$ and those of size $25-50$. We consider the fraction in the reference \DIFdelbegin \DIFdel{set }\DIFdelend \DIFaddbegin \DIFadd{sets }\DIFaddend to be a reasonable estimate of a false positive rate for predicting membrane binding function based on the presence of \DIFaddbegin \DIFadd{co-instertable }\DIFaddend protruding hydrophobes. \DIFdelbegin \DIFdel{We find this }\DIFdelend \DIFaddbegin \DIFadd{The reference proteins (Fig~\ref{FigContainingProtrusions}~D-F), indicate a }\DIFaddend false positive rate \DIFdelbegin \DIFdel{to be }\DIFdelend \DIFaddbegin \DIFadd{in the range of $20\%-30\%$. The lack of membrane interaction is not asserted for this set, and we do expect it to contain some proteins with undetected or unclassified membrane binding. The false positive rate is }\DIFaddend around $12\%$ for \DIFdelbegin \DIFdel{co-insertable protrusions, in both the size ranges we have analysed (see }\DIFdelend \DIFaddbegin \DIFadd{the non-binding surfaces (}\DIFaddend Fig~\ref{FigContainingProtrusions}\DIFdelbegin \DIFdel{)}\DIFdelend \DIFaddbegin \DIFadd{~A-C), but with a smaller sample size this estimate comes with somewhat wider error bars}\DIFaddend . For the peripheral membrane proteins, \DIFdelbegin \DIFdel{we estimate that }\DIFdelend \DIFaddbegin \DIFadd{around }\DIFaddend $64\%$ and $75\%$ of the peripheral membrane proteins in the respective \DIFdelbegin \DIFdel{groups }\DIFdelend \DIFaddbegin \DIFadd{size-groups }\DIFaddend have co-insertable protruding hydrophobes. \DIFdelbegin \DIFdel{Assuming that such motifs occur by chance at a rate no higher than it does in the reference set , and that the over-representation is due to the membrane binding function defining the data sets, we conservatively estimate that }\DIFdelend \DIFaddbegin \DIFadd{As for previous presented analysis, we see that the predictive power is weaker for the set }\emph{\DIFadd{Peripheral PISA}}\DIFadd{. We have interpreted this as a dependence on quaternary-structure modeling, which is corroborated by Fig~\ref{FigOligomerComparison}. Considering the manually curated oligomeric states two be more reliable, we expect the peripheral proteins presented in Fig~\ref{FigContainingProtrusions}~A-C to be a better approximation to real proteins. To get an idea of how commonly }\DIFaddend co-insertable \DIFaddbegin \DIFadd{protruding hydrophobes occur as membrane-interacting motifs, we will assume the rate of occurrence in the set }\emph{\DIFadd{Peripheral}}\DIFadd{, and conservatively assume a frequency of occurrence on non-membrane interacting sites around $20\%$, which is consistent with both the extremes of the 95$\%$-confidence intervals in the non-binding surfaces (Fig~\ref{FigContainingProtrusions}~A-C), and the estimate from the reference proteins (Fig~\ref{FigContainingProtrusions}~D-F). Even when considering that as much as $20\%$ of co-insertable }\DIFaddend protruding hydrophobes \DIFdelbegin \DIFdel{are membrane-interacting motifs for more than }\DIFdelend \DIFaddbegin \DIFadd{might not be membrane interacting, we would expect a rough estimate of around }\DIFaddend half of the \DIFdelbegin \DIFdel{peripheral membrane proteins we have analysed}\DIFdelend \DIFaddbegin \DIFadd{analysed membrane binders to have this motif at membrane-interacting sites}\DIFaddend .

\subsection*{Protruding hydrophobes vs. experimentally verified membrane-binding sites}
\FigDomainExamples
\FigCompExperiment
The analysis presented in Figures~\ref{FigNeighboursVsOddsRatio} and~\ref{FigContainingProtrusions} suggests that the concepts of protruding hydrophobes and co-insertability can be used to identify membrane binding residues. Based on these results we seek to define a predictor of membrane binding sites. We define \emph{the Likely Inserted Hydrophobe} as the protruding hydrophobe with the highest number of co-insertable protruding hydrophobes and lowest local protein density, as defined in~\matmet/. Fig~\ref{FigDomainExamples} illustrates that this simple definition is able to identify binding sites on modular membrane-binding domains: C1, C2, PX, ENTH, PLA2 and FYVE. For most of these cases, the Likely Inserted Hydrophobe has in fact been experimentally \DIFdelbegin \DIFdel{identified }\DIFdelend \DIFaddbegin \DIFadd{indicated }\DIFaddend to contribute to membrane binding. For the other examples, it is clearly positioned close to the experimentally identified binding site. A more quantitative comparison between predicted and verified membrane interacting residues is complicated by the \DIFdelbegin \DIFdel{absence }\DIFdelend \DIFaddbegin \DIFadd{lack }\DIFaddend of negative assertions from either methods. Experiments aiming at identifying membrane-binding sites will usually only target some of the amino acids suspected to belong to the membrane binding residues, and usually not conclude on other amino acids. \DIFaddbegin \DIFadd{To the extent non-binding amino-acids are investigated, interpretation of results in this context is also less straightforward, as non-interaction of an amino-acid does not strictly preclude it from being located close to a binding site. }\DIFaddend Similarly, the Likely Inserted Hydrophobe is by definition only one residue, and provide no negative prediction of which amino acids do not bind the membrane. We can however make a rough, but well defined, comparison by computing the angle between the vectors connecting the protein center with respectively; the mean position of the membrane interacting residues identified in experiments ($\vv{t}_{I_{e}}$), and the Likely Inserted Hydrophobe ($\vv{t}_{I_{p}}$, See Eq~\ref{eq:orientation}). While this comparison does not provide a quantitative evaluation of whether experimentally determined IBS and predicted residues match exactly, it allows us to separate proteins where the predicted and verified residues are ``on the same side'' of the protein ($\angle\vv{t}_{I_{e}}\vv{t}_{I_{p}}<~90$\degree) from those where they are not. We show on Fig~\ref{FigCompExperiment} such a comparison for proteins whose binding sites are experimentally determined. This is a coarse approximation to the protein orientation, which is sensitive to both protein shape, the selection of residues included in the partial biding sites, and any difference in backbone conformation between bound and unbound protein. Even so, we do expect that wrong binding site predictions should provide angles in the entire range from $0$\degree to $180$\degree with roughly uniform probability. But, we observe that almost all angles are sharper than $90$\degree, indicating a reasonable agreement with experimental data. We also observe a similar range of angles for cases where the membrane interaction of the Likely Inserted Hydrophobe has been experimentally verified (marked with asterisks ($*$) in Fig~\ref{FigCompExperiment}) and the cases where it has not. We would like to emphasise at this point that the Likely Inserted Hydrophobes that are not yet found to be membrane interacting might very well never have been tested. We also calculated all angles between the set of experimentally identified residues and protruding amino acids of all kinds. These results are displayed as box-plots in Fig~\ref{FigCompExperiment}. While they vary a bit between families, we note that all medians are close to $90$\degree, confirming that the statistical expectation for protrusions in general is to have roughly equally many observations larger than and smaller than $90$\degree. \DIFaddbegin \DIFadd{Interestingly, the Bovine $\alpha$-lactalbumin, for which we find no protruding hydrophobes, is analysed in its crystallised form while it is known to bind membranes in a molten globule state~\mbox{%DIFAUXCMD
\cite{Agasoster:2003bja}}%DIFAUXCMD
.
}\DIFaddend 

We provide as Supporting Information the complete list of amino acids experimentally identified as being part of membrane binding sites (\nameref{S1_File}, Table~S2). It overlaps with the list provided by Lomize \ea/ \cite{Lomize:2007gm}, but sometimes differ in exactly which amino acids are included, as we include \DIFaddbegin \DIFadd{indicated }\DIFaddend membrane interacting residues even when they are not inserted in the hydrophobic core of the membrane.

\subsection*{Protruding hydrophobes on predicted membrane binding sites}
\FigLargestFacetLocation
The continuum-model presented by Lomize~\ea/~\cite{Lomize:2011hl} forms the basis for a systematic effort to predict binding orientations for peripheral membrane proteins. The OPM database~\cite{Lomize:2011iz} provides prediction of spatial arrangements of membrane proteins with respect to the lipid bilayer for a selection of peripheral membrane proteins. We here investigate to what extent protruding hydrophobes are captured by the model proposed by Lomize~\ea/ \DIFdelbegin \DIFdel{. }\DIFdelend We identify The Likely Inserted Hydrophobe for each of the proteins in our dataset, and \DIFdelbegin \DIFdel{extracts }\DIFdelend \DIFaddbegin \DIFadd{extract }\DIFaddend the OPM predicted insertion coordinate of its \ca/-atom. The \emph{insertion coordinate} of an atom measures its depth of insertion into the hydrocarbon region of the membrane model, and is thus positive for atoms located in the hydrocarbon core and negative for atoms located on either side of the membrane including the interfacial region (Cf. \matmet/). Fig~\ref{FigLargestFacetLocation} shows histograms of the median insertion coordinate of the Likely Inserted Hydrophobes identified in each family.  A clear majority of those residues are located close to the interface of the membrane model in the OPM-predictions (Fig~\ref{FigLargestFacetLocation}~A) and $75\%$ of the families in the set of peripheral membrane proteins have the median insertion coordinate for the Likely Inserted Hydrophobe within a margin of $0.5$ nm from the membrane. This fraction is similar to the estimated fraction of proteins that have co-insertable protruding hydrophobes (Figures~\ref{FigContainingProtrusions}~A and~B). We allow this margin of $0.5$ nm to compensate for the assumptions of rigid protein, flat membrane, and the distance between \ca/-atoms and side-chain atoms. Fractions for other margins can be read from the cumulative histogram shown in Fig~\ref{FigLargestFacetLocation}~C. By representing position with the insertion coordinate, we effectively project residue coordinates onto the membrane normal. We therefore do not expect surface amino acids to be uniformly distributed along the insertion coordinate axis and present control statistics for randomly chosen protruding amino acids of all hydropathic properties (Fig~\ref{FigLargestFacetLocation}~B and~\ref{FigLargestFacetLocation}~D). It appears clearly that the high number of Likely Inserted Hydrophobes close to the membrane model is not an effect of it simply being more protein there.

\subsection*{Structure and amino acid composition at hydrophobic protrusions}

\FigSSE
The analysis presented in Fig~\ref{FigNeighboursVsOddsRatio} indicates that the ability to discriminate the data sets based on the frequency of hydrophobes on protrusions gets lower as the local protein density gets higher. Local protein density of a protrusion is dependent on secondary structure elements with loops, turns and bends being those that intuitively favor low local protein density. These secondary structures typically mark a clear change in direction of the backbone trace, where the neighbouring residues 'make way' for the protruding hydrophobe. Fig~\ref{FigSSE}~A shows which secondary structure elements the protruding hydrophobes are associated \DIFdelbegin \DIFdel{to }\DIFdelend \DIFaddbegin \DIFadd{with }\DIFaddend in the set of peripheral proteins. We note that loops, turns and bends are indeed abundant, but also helices and not beta-strands. Fig~\ref{FigSSE}~B shows a comparison with the reference data set. We see that protruding hydrophobes on turns and bends are not only common in the peripheral membrane proteins as we saw in Fig~\ref{FigSSE}~A, but they are also significantly more frequent than in the reference set. Interestingly, this is not the case for loops. \DIFdelbegin \DIFdel{A reason for this might be that turns and bends provide a rigid scaffold for exposing the hydrophobes, which would otherwise rearrange to desolvate when exposed to solvent, and thereby likely reduce the free-energy gain of membrane insertion. As the definition of }\emph{\DIFdel{loop}} %DIFAUXCMD
\DIFdel{here is simply }\DIFdelend \DIFaddbegin \DIFadd{Turns and bends are by definition structural elements with restricted flexibility~\mbox{%DIFAUXCMD
\cite{Petsko:2004wf} }%DIFAUXCMD
compared to }\emph{\DIFadd{loops}}\DIFadd{, which are here defined as the }\DIFaddend absence of any of the other secondary structure definitions \DIFdelbegin \DIFdel{, we would expect this }\DIFdelend \DIFaddbegin \DIFadd{(equivalent to }\emph{\DIFadd{coil}}\DIFadd{). We expect the latter }\DIFaddend category to contain less regular, more flexible structures. We \DIFaddbegin \DIFadd{speculate that turns and bends provide rigid scaffolds for exposing hydrophobic side chains, which might otherwise rearrange to desolvate when exposed to solvent. We }\DIFaddend also expect this property of rigid scaffolding from amphipathic helices, which is an established motif for membrane association. Fig~\ref{FigSSE} illustrates however that protrusions are not dominantly helices, confirming that the concept of protruding hydrophobes provides a useful generalisation for the shapes of membrane-binding sites.

\FigTmPeriphAaComparison
For purposes of isolating the structural component of hydrophobic membrane association, we have until now used a dichotomous definition of hydrophobicity based on \DIFdelbegin \DIFdel{the Wimley-White scale for interfacial insertion~\mbox{%DIFAUXCMD
\cite{Wimley:1996vpa} }%DIFAUXCMD
}\DIFdelend \DIFaddbegin \DIFadd{signs of free energy of transfer determined by Wimley and White~\mbox{%DIFAUXCMD
\cite{Wimley:1996vpa} }%DIFAUXCMD
(leucine, isoleucine, phenylalanine, tyrosine, tryptophan, cysteine and methionine have been considered to be hydrophobic)}\DIFaddend . Yet, we do expect different amino acids to have varying contributions to the free energy of binding. We have therefore also assessed the relative importance of different amino acids for discriminating between our \DIFdelbegin \DIFdel{two data }\DIFdelend sets. Fig~\ref{FigTmPeriphAaComparison}~B shows the comparison of the frequencies of different hydrophobic amino acids on protrusions in the \DIFdelbegin \DIFdel{two data sets }\DIFdelend \DIFaddbegin \emph{\DIFadd{Peripheral proteins}} \DIFadd{and the }\emph{\DIFadd{Non-binding surfaces}}\DIFadd{. Analysis of the other two sets can be found as supporting information (}\nameref{S1_File}\DIFadd{)}\DIFaddend . As expected, we find non-polar residues with large aliphatic or aromatic side chains to be much more frequent at the protrusions of \DIFdelbegin \DIFdel{peripheral proteins }\DIFdelend \DIFaddbegin \emph{\DIFadd{Peripheral proteins}} \DIFaddend than in the \DIFdelbegin \DIFdel{reference data set }\DIFdelend \DIFaddbegin \DIFadd{set of }\emph{\DIFadd{non-binding surfaces}}\DIFaddend . While the error bars in Fig~\ref{FigTmPeriphAaComparison}~B are not corrected for multiple testing, the signal for the hydrophobes as a group is quite clear. They all occur as over-represented in the set of \DIFdelbegin \DIFdel{peripheral proteins }\DIFdelend \DIFaddbegin \emph{\DIFadd{Peripheral proteins}} \DIFaddend and the odds-ratio is much larger for phenylalanine, leucine and tryptophan than for any of the amino-acids that are over-represented in the \DIFdelbegin \DIFdel{reference set. }\DIFdelend \DIFaddbegin \emph{\DIFadd{Non-binding surfaces}}\DIFadd{. Analysis of the other two sets can be found as supporting information (}\nameref{S1_File}\DIFadd{). }\DIFaddend Recall that $\ln{R}$ (Eq~\ref{eq:odds_ratio}) is symmetric around 0, so the magnitude of the bar representing phenylalanine on one end, can be directly compared to that of the bar representing threonine in the negative direction. Tyrosine on the other hand discriminates the sets poorly compared to its high hydrophobicity score in the Wimley-White scale. We consider this a possible consequence of the orientational restrictions on the binding sites of peripheral membrane proteins. The typical orientations consistent with shallow binding, has the residue anchored above the membrane. This probably allows less freedom for the hydroxyl group of tyrosine to orient towards regions of higher water density, than it has in the peptides used for the Wimley-White experiments or in transmembrane proteins. We also note with interest that proline is among the residues that are somewhat over-represented in the set of peripheral proteins. In general, prolines are conformationally important protein components, that restricts the backbone with respect to its immediate neighbours along the peptide chain, and are therefore likely to promote local rigidity. They also serve to induce sharp changes in the backbone trace, which would facilitate solvent exposure of neighbouring side-chains, as discussed above. Specifically, they are in general frequently found on turns~\cite{Wilmot:1988vb}.

\DIFaddbegin \subsection*{\DIFadd{Comments on the protein model}}
\DIFadd{The convex-hull representation presents a useful abstraction of proteins for investigating surface properties of approximately rigid protein conformers interacting shallowly with an approximately flat membrane. The model enables statistical analysis of protein structures, which is prohibited by high-resolution models where model parameters and quality controls typically have to be made subjectively for individual protein-membrane systems. We have employed this abstraction specifically to quantify and understand aspects of hydrophobes in peripheral membrane binding. In order to isolate components contributing to membrane binding, we have purposefully avoided complicating the interpretation with other known important factors, such as electrostatics, conformational flexibility, and even relative hydrophobicity. For the purposes of understanding the competition and complementarity between different generic models in explaining membrane-binding, it will be necessary to take these other binding factors into account in ways that allows decomposition of their contribution. In the framework of non-energetic structural analysis presented here, it is natural to do that in terms of comparing presence, absence and location of predicted binding sites between protein models. Particularly, models of electrostatic binding are well developed and readily applicable to surface representations of rigid protein conformers. For comparison of partial binding-site indications we have developed a measure of agreement that can be applied to models of low resolution, and to compare models of different resolution. We expect this to be useful in analysing the complementarity of hydrophobic and electrostatic binding modes. While complex energetic models, or machine learning approaches, can be expected to yield high performance in predicting membrane-binding properties of proteins, the transparent kind of model presented here provides a clear interpretation of why a specific protein is predicted to bind. Apart from the obvious academic interest, this connection to expert knowledge is invaluable for interpreting automated classifications where the models can not be reliably parameterised against negative data, that is; definitely non-binding proteins. Co-localized binding-site indicators based on different generic binding models, such as hydrophobic and electrostatic models, can provide a much improved performance in such prediction while maintaining interpretability. Such decomposition of binding modes are also useful for inference or interpretation about the specificity towards interacting membranes, as electrostatic binding modes are expected in interactions with highly charged membranes, while predominantly hydrophobic interactions are more likely to be specific to membrane or lipid morphologies.
}

\DIFaddend \section*{Conclusion}

Protein-membrane interactions are typically studied \emph{in vitro} or \emph{in silico} and inference to their biological context have to carry over from greatly simplified membrane models. To make sense of such experiments and simulations, it is essential to formulate general models that explain protein association in terms of factors that are present in both model systems and the relevant \emph{in vivo} counterpart. In pursuit of such general models for membrane recognition, we have formulated the concepts of protruding hydrophobes and co-insertability. We have analysed more than 300 families of proteins that are classified as peripheral membrane binders and identified this model to be a good fit to \DIFdelbegin \DIFdel{more than }\DIFdelend \DIFaddbegin \DIFadd{at least }\DIFaddend half of them, after \DIFaddbegin \DIFadd{cautiously }\DIFaddend correcting for the \DIFdelbegin \DIFdel{small false positive rate }\DIFdelend \DIFaddbegin \DIFadd{false positive rates }\DIFaddend estimated from the reference \DIFdelbegin \DIFdel{set }\DIFdelend \DIFaddbegin \DIFadd{sets }\DIFaddend (Fig~\ref{FigContainingProtrusions}). The generality of the model is corroborated by three important points. Hydrophobes are clearly over-represented on the protrusions of peripheral membrane proteins (compare Fig~\ref{FigProtrusionDensityHistograms}~A and~\ref{FigProtrusionDensityHistograms}~C, and see Fig~\ref{FigNeighboursVsOddsRatio}), they tend to locate on co-insertable protrusions (see Fig~\ref{FigHfConditionals} and Fig~\ref{FigContainingProtrusions}), and protruding hydrophobes are generally positioned consistent with experimentally identified binding sites (Fig~\ref{FigDomainExamples} and Fig~\ref{FigCompExperiment}). Amphipathic helices are already well known membrane binding motifs which our definition of protrusion is well suited to capture, whenever these are stably folded and exposed. We do however find that the majority of identified protruding hydrophobes are not helices (Fig~\ref{FigSSE}~A) and that hydrophobes are also highly over-represented on protruding turns and bends (Fig~\ref{FigSSE}~B). We therefore propose the concept of protruding hydrophobes as a useful generalisation upon binding motifs that are identified in terms of secondary structure.

\DIFdelbegin \DIFdel{Both the choice of reference set, and the choice of quaternary structure modelling comes with some assumptions. We have elaborated on these in “materials and methods”. We have also performed some checks on how sensitive our analyses are to violations of these assumptions, and found that our conclusions are robust. We present details of these analyses as Supporting Information (}%DIFDELCMD < \nameref{S1_File}%%%
\DIFdel{).
}%DIFDELCMD < 

%DIFDELCMD < %%%
\DIFdelend Investigation of the interfacial binding sites of numerous peripheral membrane proteins has revealed the presence of hydrophobic amino acids, and of basic amino acids such as arginines and lysines. This reflect the two almost universal traits of biological membranes; their hydrophobic core and anionic surface. Yet the focus on the electrostatic component of the free energy of transfer from water to membrane - often referred to as being long-range - has overshadowed the importance of hydrophobic contribution which is sometimes referred to as being short-range. The focus on electrostatic interaction is at least in part to be attributed to the difficulties in evaluating the hydrophobic contribution as opposed to for example, the computational tractability of continuum electrostatic models. In principle the contribution of hydrophobes to membrane binding can only be determined with a rigorous treatment of the hydrophobic effect, which requires very accurate treatment of large systems involving both protein, membrane and solvent. The mere presence of hydrophobes on the protein surface is to a large extent tolerated by non-membrane-binding proteins as well, and for both hydrophobes and basic amino acids, it is challenging to determine when their presence on protein surfaces are coincidental, and when they are important for membrane binding. Moreover, amino acids on membrane binding sites are not typically strongly conserved~\cite{Park:2016kw}, so modeling their generic binding modes is important both for relating binding sites between homologs, and for understanding how additional factors determine differences in membrane specificities. Fortunately, as evident from the results presented in this contribution, the role of hydrophobes can often be understood in much simpler terms than what is required for an exact estimate of the energetics of the hydrophobic effect, and their importance for membrane-binding can be inferred from comparative statistical analyses. The subtle considerations of protein structure encoded in our definition of protrusions, strongly distinguishes the small hydrophobic patches on peripheral membrane proteins from those on other protein surfaces. This provides good reason to assume their importance for binding.

\section*{Materials and methods}
\subsection*{Data sets}
We \DIFdelbegin \DIFdel{obtained data setsfrom the collection of proteins in }\DIFdelend \DIFaddbegin \DIFadd{have compiled four data sets, two versions of a set of peripheral proteins, and two different reference sets:
} \begin{itemize} 
	\item \emph{\DIFadd{Peripheral Proteins}}\DIFadd{: A set of peripheral membrane binders obtained from }\DIFaddend the OPM-database~\cite{Lomize:2011iz} \DIFaddbegin \DIFadd{using the OPM quaternary structure models.
	}\item \emph{\DIFadd{Peripheral PISA}}\DIFadd{: A subset of }\emph{\DIFadd{Peripheral Proteins}} \DIFadd{where no protein overlap in terms of their SCOPe-family classification ~\mbox{%DIFAUXCMD
\cite{Fox:2014vx} }%DIFAUXCMD
and with quaternary structure predicted by PISA~\mbox{%DIFAUXCMD
\cite{Krissinel:2007wq}}%DIFAUXCMD
.
	}\item \emph{\DIFadd{Non-binding Surfaces}}\DIFadd{: A set of protein surfaces obtained from the solvent exposed regions of transmembrane proteins}\DIFaddend .
	\DIFdelbegin \DIFdel{Our set of peripheral proteinsare all }\DIFdelend \DIFaddbegin \item \emph{\DIFadd{Reference Proteins}}\DIFadd{: A non-redundant set of proteins from 5 SCOPe-classes, obeying the following conditions: (1) none of these proteins have a domain represented in OPM, and (2) none of the proteins in the dataset have domains belonging to the same SCOPe-family (the same restriction as for }\emph{\DIFadd{Peripheral PISA}}\DIFadd{).
} \end{itemize} 
\DIFadd{In our analysis }\emph{\DIFadd{Peripheral Proteins}} \DIFadd{is always compared to }\emph{\DIFadd{Non-binding surfaces}}\DIFadd{, and }\emph{\DIFadd{Peripheral PISA}} \DIFadd{to }\emph{\DIFadd{Reference Proteins}}\DIFadd{.
}

\emph{\DIFadd{Peripheral Proteins}} \DIFadd{are all the }\DIFaddend proteins in OPM classified as \emph{type}: \emph{Monotopic/peripheral}. While the OPM has strict criteria for inclusion, membrane binding is not asserted by experiment in all cases, and the set might contain false positives. \DIFaddbegin \DIFadd{This data set is provided as supporting information (}\nameref{S1_Set}\DIFadd{).
}\DIFaddend 

The \DIFdelbegin \DIFdel{reference set }\DIFdelend \DIFaddbegin \DIFadd{set }\emph{\DIFadd{Non-binding surfaces}} \DIFaddend consist of fragments of transmembrane complexes. We obtained these protein fragments from all proteins classified as \emph{type}: \emph{Transmembrane} in OPM. The fragments analysed are composed of all amino acids whose \ca/-coordinates are at least $1.5$~nm from the hydrocarbon region of the membrane model (The parameter $Z_{HDC}$ in the OPM model~\cite{Lomize:2011gj}). We rely here on membrane models positioned by the OPM, which we deem very reliable for transmembrane proteins. While the entire protein complex was considered when calculating structural properties, only the fragments meeting this distance criteria were considered in the statistical analyses. When these proteins interact with secondary membranes or interact with membranes of extremely high curvature, it is not captured by the OPM model, and the assumption that these surfaces are not interacting with membrane may be violated. We have assumed that such issues are exceptional. \DIFaddbegin \DIFadd{This data set is provided as supporting information (}\nameref{S2_Set}\DIFadd{).
}\DIFaddend 

We do consider the assumptions mentioned above to be conservative. Inclusion of non-membrane-binding proteins in our set of peripheral membrane proteins would likely weaken any general signal from membrane binding proteins, and inclusion of secondary membrane interactions sites in the reference set would probably inflate the number of hydrophobes on protrusions in that set. 

All protein structures \DIFaddbegin \DIFadd{in these two sets }\DIFaddend are obtained by X-ray crystallography and NMR spectroscopy and we have assumed that at least the backbone coordinates are representative of the solvated state of the proteins. As the source of structural information for this database is the Protein Data Bank (PDB)\cite{Berman:2000fma} the relevant oligomeric state is not always determined, The curators of the OPM-database have decided on oligomer models, upon which we have relied \DIFaddbegin \DIFadd{for the sets }\emph{\DIFadd{Peripheral Proteins}} \DIFadd{and }\emph{\DIFadd{Non-binding surfaces}}\DIFaddend . These are taken from PDBe~\cite{Velankar:2009he}, generated by PISA~\cite{Krissinel:2007wq} or obtained from literature as described by Lomize \ea/~\cite{Lomize:2011iz}. 
\DIFdelbegin \DIFdel{As weak protein-protein interaction interfaces may also contain exposed hydrophobic patches, we }\DIFdelend \DIFaddbegin 

\DIFadd{Even if the solvent exposed regions of the proteins in the set }\emph{\DIFadd{Non-binding surfaces}} \DIFadd{are extracted after relevant properties for potential membrane interaction was calculated, we cannot exclude totally that the surface constructed reflect artifacts of the extration of fragments from complete protein models. In addition we }\DIFaddend expect our analysis to be sensitive to \DIFdelbegin \DIFdel{how protein quaternary structure is modelled}\DIFdelend \DIFaddbegin \DIFadd{quaternary structure modeling, as oligomeric protein-protein interfaces may also contain exposed hydrophobic patches~\mbox{%DIFAUXCMD
\cite{Jones:2000vm, TonddastNavaei:2015bm}}%DIFAUXCMD
}\DIFaddend . As a quality control, we therefore also performed our analysis relying solely on computationally predicted quaternary structures, \DIFdelbegin \DIFdel{which we provide in the Supporting Information (}%DIFDELCMD < \nameref{S1_File}%%%
\DIFdel{)}\DIFdelend \DIFaddbegin \DIFadd{and complete protein structures. This is implemented in the comparison of }\emph{\DIFadd{Peripheral PISA}} \DIFadd{and }\emph{\DIFadd{Reference Proteins}}\DIFaddend . 
\DIFdelbegin \DIFdel{This control reproduced qualitatively all observations that we have interpreted. In the Supporting Information we also report analysis on the sensitivity of the results to how the reference setis obtained, using a reference }\DIFdelend \DIFaddbegin 

\DIFadd{The set }\emph{\DIFadd{Reference Proteins}} \DIFadd{is constructed as a subset of all PDB IDs determined by X-ray crystallography, with some domain classified in SCOPe~\mbox{%DIFAUXCMD
\cite{Fox:2014vx} }%DIFAUXCMD
in the classes }\emph{\DIFadd{All alpha proteins}} \DIFadd{(sunid: 46456), }\emph{\DIFadd{All beta proteins}} \DIFadd{(sunid: 48724), }\emph{\DIFadd{Alpha and beta proteins (a+b)}} \DIFadd{(sunid: 51349), }\emph{\DIFadd{Alpha and beta proteins (a/b)}} \DIFadd{(sunid: 53931) or }\emph{\DIFadd{Multi domain proteins}} \DIFadd{(sunid: 56572). The exclusion of structures not determined by X-ray crystallography was done to get consistent quaternary structure predictions. All PDB IDs that has one or more domains classified in the same SCOPe-family as any domain in the OPM-database~\mbox{%DIFAUXCMD
\cite{Lomize:2011iz} }%DIFAUXCMD
was excluded from the set. This excludes not only the peripheral membrane binders, but also any transmembrane protein found in the reference set used for our primary analysis. In order to avoid redundancy, we iteratively removed proteins with domains that share SCOPe-family classification with any other domain in the }\DIFaddend set\DIFdelbegin \DIFdel{based on the SCOPe-classification~\mbox{%DIFAUXCMD
\cite{Fox:2014vx}}%DIFAUXCMD
. }\DIFdelend \DIFaddbegin \DIFadd{, until there were no such shared classifications left. This process ensures that there is at most one representative for each SCOPe family in the set. We generated quaternary structure models using PISA~\mbox{%DIFAUXCMD
\cite{Krissinel:2007wq} }%DIFAUXCMD
for all members of this set. While this data set consist of more complete protein surfaces, it is only intended to be a reference for typical protein surfaces, we do expect it to be a mix of both membrane interacting and non-interacting proteins. This data set is provided as supporting information (}\nameref{S4_Set}\DIFadd{).
}\DIFaddend 

\DIFaddbegin \DIFadd{The set }\emph{\DIFadd{Peripheral PISA}} \DIFadd{was derived from }\emph{\DIFadd{Peripheral Proteins}}\DIFadd{, for comparability with }\emph{\DIFadd{Reference Proteins}}\DIFadd{. All structures not determined by X-ray crystallography was excluded, and proteins with domains that share SCOPe-family classification with any other domain in the set was iteratively removed until no such shared classifications were left. Quaternary structure models were predicted using PISA. This data set is provided as supporting information (}\nameref{S3_Set}\DIFadd{).
}

\DIFaddend A few structures meeting the above criteria, were not included in the analysis for technical reasons, such as issues with formatting of PDB files. After exclusion of these cases, the final set of \DIFdelbegin \DIFdel{peripheral proteins }\DIFdelend \DIFaddbegin \emph{\DIFadd{Peripheral proteins}} \DIFaddend contains 1012 protein structures classified into 326 families. The final \DIFdelbegin \DIFdel{reference set }\DIFdelend \DIFaddbegin \DIFadd{set of }\emph{\DIFadd{Non-binding surfaces}} \DIFaddend contains 495 protein structures classified into 158 families. \DIFaddbegin \DIFadd{After filtering out the few structure files that did not exactly comply to the expected PDB format, the final set of }\emph{\DIFadd{Peripheral PISA}} \DIFadd{binders contained $170$ proteins, and the set }\emph{\DIFadd{Reference Proteins}} \DIFadd{contained $2250$ proteins.
}\DIFaddend 

\DIFaddbegin \FigOligomerComparison
\DIFadd{The two set of peripheral proteins are both derived from OPM, but }\emph{\DIFadd{Peripheral PISA}} \DIFadd{is organized in a different classification from }\emph{\DIFadd{Peripheral Proteins}}\DIFadd{, and retains fewer structures. In addition their quaternary structures, which are not completely determined by X-ray crystallography, are modeled differently. In Figure~\ref{FigOligomerComparison}, we illustrate this difference in quartnery structure by showing the difference in the number of polypetide chains present in these protein-models in the two sets.
}

\DIFaddend Based on experiments reported in available literature~\cite{Hedin:2002eu, Grauffel:2013uf, Malmberg:2003uk, Stahelin:2003cv, Stahelin:uj, Wang:2001jw, Stahelin:2003is, Stahelin:2001cr, Frazier:2003wf, Gerber:2002fp, Rufener:2005ie, CorbalanGarcia:2003ba, Kim:2000em, Gilbert:2002ja, Feng:2002fc, Feng:2003ik, Grauffel:2013uf, Walther:2004ut, Kohout:2003bs, Goh:2016hv, BCampos:1998fa, Isas:2004fw, Kutateladze:2001ir, Stahelin:uj, Wang:2001jw, Anderluh:2005kk, Shenkarev:2006fg, Lin:1998uz, Canaan:2002dp, Lathrop:2001wu, Chen:2000vf, SekinoSuzuki:1996uy, Phillips:2005gn, Thennarasu:2005jw, Tatulian:2005ja, Oldham:2005fb, Mathias:2009bj, Agasoster:2003bja, Jian:2015do}, we made a data set of partially identified membrane binding sites on proteins with resolved structures. This set contains membrane interacting residues of 34 protein structures, classified into 22 families. A detailed description is provided in the Supporting Information (\nameref{S1_File}, Table~S2).

\subsection*{Definitions}

\subsubsection*{Structural characteristics of protein surfaces}

We characterise the surface of proteins with different criteria designed to capture solvent\emph {exposed} residues, \emph{protruding} residues and \emph{co-insertable} protruding residues. The two latter are illustrated in Fig~\ref{FigModelIllustration}. 

\emph{Exposed} amino acids are defined as all amino acids that have a solvent accessible side-chain area greater than $0.2~\textrm{nm}^{2}$, as calculated with a probe with a radius of $0.14$ nm, following the procedure described in Eisenhaber \ea/~\cite{Eisenhaber:1995bd} using van der Waals radi reported by Bondi~\cite{Bondi:1964un}.

We identify a \emph{protrusion} or a \emph{protruding} residue via the calculation of the convex hull of the \ca/- and \cb/-coordinates of the protein. The convex hull of a set of points \emph{S} is the smallest possible convex set containing \emph{S}. We define \emph{vertex} residues as residues whose \cb/-atom is a vertex of this convex hull. A \emph{protrusion} or a \emph{protruding} residue, is defined as a \emph{vertex} residue that also has low local protein density. For the purposes of this work, we will define the local protein density $d$ of a residue, as the number of \ca/- or \cb/-atoms within a distance $c$ of its \cb/-atom. We will designate a local protein density as low, if $d<n$, with $n=22$ and $c=1~\mathrm{nm}$. These parameters were manually chosen based on a set of six different families of peripheral membrane proteins (C2-domain, PX-domain, Discodin domain, ENTH domain, Lipoxygenases and a Bacterial Phospholipase C). A list of these proteins are provided as Supporting Information (\nameref{S1_File}, Table~S1).

We define two protrusions to be \emph{co-insertable} or a \emph{co-insertable pair}, if the straight line connecting them is an edge of the convex hull polygon.

\subsubsection*{Hydrophobic residues}
An amino acid is defined to be \emph{hydrophobic}, or a \emph{hydrophobe}, if it contributes favourably to membrane interface partitioning of peptides, as determined in the Wimley-White scale for interfacial insertion~\cite{Wimley:1996vpa}. These amino acids are: leucine, isoleucine, phenylalanine, tyrosine, tryptophan, cysteine and methionine.

\subsubsection*{Secondary structure}
We use DSSP definitions~\cite{Kabsch:1983bp} for protein secondary structure. DSSP codes H, G or I are reported as \emph{helix}, DSSP codes B or E as $\beta$, DSSP code T as \emph{bend} and DSSP code S as \emph{turn}. All other residues are considered to be in \emph{loops}.

\subsubsection*{Likely Inserted Hydrophobe}
The \emph{Likely Inserted Hydrophobe} is defined as the protruding hydrophobe with the largest number of co-insertable protruding hydrophobes in a protein. Ties are resolved by choosing the likely inserted hydrophobe with the smallest local protein density $d$. Further ties are resolved by random selection, so that each protein has exactly one Likely Inserted Hydrophobe, unless it has no protruding hydrophobes at all.

\subsubsection*{Insertion coordinate}
For comparisons with OPM predictions, we define the \emph{insertion coordinate} of atoms. This coordinate measures how deeply into the OPM membrane model an atom is inserted, and is therefore negative on the solvated side of the membrane. The membrane perimeter, where the insertion coordinate is $0$, is the end of the hydrocarbon region. We identify this boundary as it is done in the model used to predict the OPM orientations, namely the planes where the volume fraction of total hydrocarbon is equal to $0.5$. See Eq~2 in~\cite{Lomize:2011gj}.

\subsection*{Measures}
\subsubsection*{Averages of residues}
We compare protein surfaces with respect to structural and hydropathic properties, reflected in different selection criteria and averaged over families or the entire data sets.

The mean fraction of residues having property $s$ with respect to a reference property $r$ in a family is:
\begin{equation}
	\label{eq:mean_frac_family}
	\hat{f}_{s|r} = \frac{1}{\left|C\right|}\sum_{G \in C} \frac{\left|G_{s} \cap G_{r}\right|}{\left|G_{r}\right|}
\end{equation}
where $C$ is the set of proteins in a family, $G$ is a protein, and, $G_{s}$ is the set of residues on a protein meeting criteria $s$. Vertical bars denote size of sets. We will specify $s$ and $r$ according to the definitions above, using intersect notation to combine criteria when necessary. $\hat{f}_{\mathrm{hydrophobe}|\mathrm{protrusion} \cap \mathrm{helix}}$, for instance, should be interpreted as the mean fraction of hydrophobes out of all protruding amino acids that are in helices.

We estimate weighted data set counts of amino acids with property $s$ as:
\begin{equation}
	\label{eq:weighted_count}
	\hat{N}_{s} = \sum_{C \in D} \left(\frac{1}{\left|C\right|}\sum_{G \in C} \left| G_{s} \right| \right)
\end{equation}
where $D$ is a data set, such as the set of peripheral proteins or the reference set. Similarly we quantify the weighted count of proteins that have at least one amino acid with property $s$ as:
\begin{equation}
	\label{eq:weighted_existence}
	\hat{M}_{s} = \sum_{C \in D} \left(\frac{1}{\left|C\right|}\sum_{G \in C} \mathrm{H}\left(\left| G_{s}\right| \right)\right)
\end{equation}
where $\mathrm{H}$ is the Heaviside step function. Given a property $s$ and reference property $r$, we estimate the weighted fraction in a data set, $\hat{F}_{s|r}$:
\begin{equation}
	\label{eq:weighted_frec}
	\hat{F}_{s|r} = \frac{\hat{N}_{s \cap r}}{\hat{N}_{r}}
\end{equation}
or the weighted fraction of proteins that have at least one residue with the given property $s$:
\begin{equation}
	\label{eq:weighted_existance_frec}
	\hat{E}_{s} = \frac{\hat{M}_{s}}{\left|D\right|}
\end{equation}
With $\left|D\right|$ being the number of families in the data set. When such fractions (Eqs.~\ref{eq:weighted_frec} or~\ref{eq:weighted_existance_frec}) are reported, we estimate $95\%$-confidence intervals using a normal approximation to the binomial distribution, with $|D|$ the total number of trials (Eq~\ref{eq:weighted_existance_frec}), or $\hat{N}_{r}$ serving as a real-number analog to the total number of trials (Eq~\ref{eq:weighted_frec}).

\subsubsection*{Averages of co-insertable pairs}
To analyse co-insertable residues, we estimate weighted data set counts of co-insertable pairs of residues with property $s$, as:
\begin{equation}
		\hat{N}^{\mathrm{pair}}_{s} = \sum_{C \in D} \left(\frac{1}{\left|C\right|}\sum_{G \in C} \left| G^{\mathrm{pair}}_{s}\right| \right)
\end{equation}
where $\left|G^{\mathrm{pair}}_{s}\right|$ are the number of co-insertable amino acids pairs with property $s$. For quantification of the weighted count of proteins that have at least one co-insertable pair with property $s$, we calculate:
\begin{equation}
	\hat{M}^{\mathrm{pair}}_{s} = \sum_{C \in D} \left(\frac{1}{\left|C\right|}\sum_{G \in C} H\left(\left| G^{\mathrm{pair}}_{s}\right|\right) \right)
\end{equation}
Considering the set of co-insertable amino acid pairs in a protein, $G^{\mathrm{pair}}$, we will denote the set of pairs where at least one of the amino acids is a protruding hydrophobe as $G^{\mathrm{pair}}_{\mathrm{one}}$, and the set where both are protruding hydrophobes as $G^{\mathrm{pair}}_{\mathrm{both}}$. We will report the weighted fraction of proteins that have co-insertable protruding hydrophobes as:
\begin{equation}
	\label{eq:coins_existance_freq}
	\hat{E}^{\mathrm{pair}}_{\mathrm{both}} = \frac{\hat{M}^{\mathrm{pair}}_{\mathrm{both}}}{\left|D\right|}
\end{equation}
and the weighted frequency of co-insertion of protruding hydrophobes as:
\begin{equation}
	\label{eq:coins_freq}
	\hat{F}^{\mathrm{pair}}_{\mathrm{both}|\mathrm{one}} = \frac{\hat{N}^{\mathrm{pair}}_{\mathrm{both}}}{\hat{N}^{\mathrm{pair}}_{\mathrm{one}}}
\end{equation}
Note that $\hat{F}^{\mathrm{pair}}_{\mathrm{both}|\mathrm{one}}$ estimates the conditional probability that both amino acids of a co-insertable pair are protruding hydrophobes, given that one of them is. The tendency for protruding hydrophobes to be located at co-insertable positions can then be quantified by comparing with a null model for each set. We obtain these null models by randomly reassigning the hydrophobic amino acids to other protruding locations in the same protein.

\subsubsection*{Comparison between data sets}
The frequency of properties in different data sets, are compared via weighted fractions. For two data sets, $A$ and $B$, we compare a certain weighted fraction $\hat{F}$ using the odds ratio, $\OR{A, B, \hat{F}}$:
\begin{equation}
	\label{eq:odds_ratio}
	\OR{A, B, \hat{F}}=\frac{\hat{F}^{A}\left(1-\hat{F}^{B}\right)}{\hat{F}^{B}\left(1-\hat{F}^{A}\right)}
\end{equation}
where $\hat{F}^{A}$ denotes the fraction $\hat{F}_{s|r}$ obtained for data set $A$. We will report $\ln{\R}$, which is symmetric around 0, so that $\ln{\OR{A,B, \hat{F}}}=-\ln{\OR{B,A,\hat{F}}}$. Wald $95\%$-confidence intervals for $\ln{\R}$ are calculated with $\hat{N}_{s \cap r}$ and $\left(\hat{N}_{r}-\hat{N}_{s \cap r}\right)$ serving as real number analogs for the count of successes and failures in the data sets compared. When $\hat{F}^{\mathrm{pair}}_{\mathrm{both}|\mathrm{one}}$ is compared, the corresponding counts of successes and failures are $\hat{N}^{\mathrm{pair}}_{\mathrm{both}}$ and $\hat{N}^{\mathrm{pair}}_{\mathrm{both}} - \hat{N}^{\mathrm{pair}}_{\mathrm{one}}$, respectively.

\subsubsection*{Comparison of experimentally verified and predicted binding sites}
We define two vectors which we then compare to evaluate the distance between experimentally verified and predicted membrane binding residues. The \ca/-coordinate of experimentally verified membrane binding residues functions as a proxy for the membrane, and the vector defined by the latter residues and the center of mass (COM) of the protein is used as a reference to which we compare the vector defined by the protein COM and the Likely Inserted Hydrophobe. Given a set of identified or predicted membrane interacting resides, $I$, we compute the vector, $\vv{t}_{I}$:
\begin{equation}
	\label{eq:orientation}
	\vv{t}_{I} = \frac{1}{\left|I\right|}\sum_{a \in I}\vv{V}_{a} - \frac{1}{\left|G_{*}\right|}\sum_{a \in G_{*}}\vv{V}_{a}
\end{equation}
where $\vv{V}_{a}$ denotes the \ca/-coordinates of residue $a$, and $G_{*}$ is the set of all residues in the protein. We will denote vectors obtained for experimentally identified membrane binding residues as $\vv{t}_{I_{e}}$, and those obtained for a Likely Inserted Hydrophobe as $\vv{t}_{I_{p}}$. We then measure the angle $\angle\vv{t}_{I_{e}}\vv{t}_{I_{p}}$ between the two vectors for each protein in the dataset of known binding sites.

\subsection*{Implementation}
The solvent accessible area was calculated with MMTK~\cite{Hinsen:2000wia} (version 2.9.0), and the convex hull was calculated with Qhull~\cite{Barber:1996iv} via scipy~\cite{SciPyOpenSourceS:g1eMLzGE} (version 0.13.3). Proportion test confidence intervals were calculated with R~\cite{Team:2008wf} (Version 2.12.0), odds ratios and corresponding confidence intervals were calculated with the R-package epitools~\cite{epitoolsEpidemiolo:dQp6gsqn} (version 0.5-6). Secondary structure annotations were computed with the CMBI DSSP implementation~\cite{Touw:2015ff} (version 2.0.4). \DIFaddbegin \DIFadd{For construction of the set }\emph{\DIFadd{Peripheral PISA}} \DIFadd{and }\emph{\DIFadd{Reference proteins}} \DIFadd{SCOPe version 2.06 was used. PISA predictions were obtained through the ``Protein interfaces, surfaces and assemblies'' service PISA at the European Bioinformatics Institute. (}\url{http://www.ebi.ac.uk/pdbe/prot_int/pistart.html}\DIFadd{). Where PISA predicted that the asymmetric unit represents the most stable quaternary structure in solution, we obtained structures from the protein data bank (}\url{http://www.rcsb.org/}\DIFadd{)\mbox{%DIFAUXCMD
\cite{Berman:2000fma}}%DIFAUXCMD
. }\DIFaddend Otherwise the analyses were implemented by us, using Python and R. Plots were produced with R, and other visualisations using VMD (Visual Molecular Dynamics)~\cite{Humphrey:1996bk}. \DIFaddbegin \DIFadd{Data sets of peripheral membrane proteins were generated on a snapshot of the OPM-database extracted the 23. Dec. 2013.
}\DIFaddend 

\section*{Supporting information}

\paragraph*{S1 Text.}
\label{S1_File}
{\bf Supplementary tables and analysis.} Analysis to assess the roboustness of \DIFaddbegin \DIFadd{some }\DIFaddend results to quaternary structure modelling\DIFdelbegin \DIFdel{and the definition of the reference set}\DIFdelend \DIFaddbegin \DIFadd{, and specification of proteins and binding sites compared with experiment}\DIFaddend .

\paragraph*{S1 Data Set.}
\label{S1_Set}
{\bf The set of \DIFdelbegin \DIFdel{peripheral membrane protein structures}\DIFdelend \DIFaddbegin \DIFadd{Peripheral proteins}\DIFaddend .} Comma-separated file with PDB IDs for the set and the OPM classification of these at the time of analysis. The PDB IDs also serve to identify entries in the OPM database.

\paragraph*{S2 Data Set.}
\label{S2_Set}
{\bf The set of \DIFdelbegin \DIFdel{structures from which reference surfaceswas extracted}\DIFdelend \DIFaddbegin \DIFadd{Non-binding surfaces}\DIFaddend .} Comma-separated file with PDB IDs for the set and the OPM classification of these at the time of analysis. The PDB IDs also serve to identify entries in the OPM database\DIFaddbegin \DIFadd{.
}

\paragraph*{\DIFadd{S3 Data Set.}}
\label{S3_Set}
{\bf \DIFadd{The set: Peripheral PISA.}} \DIFadd{PDB IDs for the set and the OPM classification of these at the time of analysis. The PDB IDs also serve to identify entries in the OPM database.
}

\paragraph*{\DIFadd{S4 Data Set.}}
\label{S4_Set}
{\bf \DIFadd{The set: Reference proteins.}} \DIFadd{PDB IDs for the set.
}

\paragraph*{\DIFadd{S5 Data Set.}}
\label{S5_Set}
{\bf \DIFadd{Surface properties for all sets.}} \DIFadd{Calculated properties for exposed residues in all sets. For the }\emph{\DIFadd{Non-binding surfaces}} \DIFadd{we have only included residues in the analysed fragments. PDB IDs, chain ids and residue ids refer to OPM quaternary structure models for }\emph{\DIFadd{Peripheral proteins}} \DIFadd{and }\emph{\DIFadd{Non-binding surfaces}}\DIFadd{, and to the PISA generated quaternary structure models for }\emph{\DIFadd{Peripheral PISA}} \DIFadd{and }\emph{\DIFadd{Reference proteins}}\DIFadd{. The local density parameter $d$ is identified as }\emph{\DIFadd{local density}}\DIFadd{, the number of hydrophobes co-insertable to a residue is identified as }\emph{\DIFadd{co insertables}}\DIFadd{, the solvent accessible surface area of a side chain ($(nm)^{2}$) is identified as }\emph{\DIFadd{sidechain sasa}}\DIFadd{, the }\emph{\DIFadd{Likely inserted hydrophobe}} \DIFadd{is identified as }\emph{\DIFadd{LIH}}\DIFadd{, and other column names are self-explanatory}\DIFaddend .

%\paragraph*{S1 Table.}
%\label{S1_Table}
%{\bf Protein structures used for parameterisation.} Table of protein structures inspected for formulation of the model, and for choice of parameters.

%\paragraph*{S2 Table.}
%\label{S2_Table}
%{\bf Binding sites.} Table of binding residues indicated by experiment, used for comparison with model predictions in Fig~\ref{FigDomainExamples} and~\ref{FigCompExperiment}.


\section*{Acknowledgments}
We thank Anne Gershenson at the University of Massachusetts Amherst, and Angèle Abboud at the University of Bergen for their valuable comments. This work was supported by grants from the Norwegian Research Council (FRIMEDBIO 214167 and FRIMEDBIO 251247).

\nolinenumbers

\bibliography{library}

\end{document}