StickWRLD was made possible in part through funding provided to Dr. Ray by the Research Institute at Nationwide Children's Hospital, and by NSF grant DBI-1262457.

1. Software Download &#38; Installation
<ol>
	<li>Use a computer has an Intel i5 or better processor with at least 4 Gb of RAM, and is running Mac OS X or GNU/Linux (e.g., Ubuntu) OS. In addition, Python 2.7.610 and the wxPython 2.811, SciPy12, and PyOpenGL13 python libraries are required - download and install each from their respective repositories.</li>
	<li>Download <a href="http://www.stickwrld.org" target="_blank">StickWRLD</a> as a zip archive containing all of the relevant Python scripts. Download the &#8220;fasta2stick.sh&#8221; script for converting standard FASTA DNA/protein sequence alignments to StickWRLD format.</li>
	<li>Extract the archive and put the resulting StickWRLD folder on your Desktop. Place the &#8220;fasta2stick.sh&#8221; script on the desktop as well.</li>
</ol>
2. Prepare the Alignment
<ol>
	<li>Create an alignment of the protein sequences using any standard alignment software (e.g., ClustalX14). Save the alignment on the desktop in FASTA format.</li>
	<li>Open the terminal application on the Mac or GNU/Linux computer and navigate to the desktop (the location of the &#8220;fasta2stick.sh&#8221; shell script) by typing cd ~/Desktop and pressing return. Execute the &#8220;fasta2stick.sh&#8221; script by typing ./fasta2stick.sh in the terminal. If the script does not execute, ensure that it is executable &#8211; in the terminal type chmod +x fasta2stick.sh to make the script executable.</li>
	<li>Follow the onscreen instructions provided by the script to specify the input file name (the file created in 1.2 above) and the desired output name. Save the output file (which is now in the correct format for StickWRLD) on the desktop.</li>
</ol>
3. Launching StickWRLD
<ol>
	<li>Navigate into the StickWRLD executables folder using the terminal application of the Mac or GNU/Linux computer. For example, if the StickWRLD folder is on the desktop, type cd ~/Desktop/StickWRLD/exec in the terminal.</li>
	<li>Launch StickWRLD by typing python-32 stickwrld_demo.py in the terminal.</li>
	<li>Verify that the StickWRLD Data Loader panel is visible on the screen (Figure 3).</li>
</ol>
4. Loading the Data
<ol>
	<li>Load the converted protein sequence alignment by pressing the &#8220;Load Protein&#8230;&#8221; button.</li>
	<li>Select the file created in step 3 above and press &#8220;Open&#8221;. StickWRLD will open several new windows, including &#8220;StickWRLD Control&#8221; (Figure 4) and &#8220;StickWRLD - OpenGL&#8221; (Figure 5).</li>
	<li>Select the &#8220;StickWRLD &#8211; OpenGL&#8221; window. Choose &#8220;Reset View&#8221; from the &#8220;OpenGL&#8221; menu to display the default StickWRLD visualization in a &#8220;top-down&#8221; view through the cylinder representing the data in the resizable OpenGL windows..</li>
</ol>
5. View Options
<ol>
	<li>Select the boxes for &#8220;Column Labels&#8221; and &#8220;Ball Labels&#8221; in the &#8220;StickWRLD Control&#8221; pane (Figure 4) to display values for columns and balls.</li>
	<li>Deselect the box for &#8220;Column Edges&#8221; in the &#8220;StickWRLD Control&#8221; pane to hide the column edge lines.</li>
	<li>Set the &#8220;Column Thickness&#8221; to 0.1 in the &#8220;StickWRLD Control&#8221; pane to draw a thin line through the columns, making it easier to navigate the 3D view. Press return to accept the change.</li>
	<li>Reset the view in the &#8220;StickWRLD &#8211; OpenGL&#8221; window as in step 5.3 above, then press the &#8220;full screen&#8221; button to maximize the view.</li>
</ol>
6. Navigation
<ol>
	<li>Rotate the 3D StickWRLD display by holding down the left mouse button while moving the mouse in any direction.</li>
	<li>Zoom the 3D StickWRLD display by holding down the right mouse button while moving the mouse up or down.</li>
</ol>
7. Finding Interpositional Dependencies (IPDs)
<ol>
	<li>Browse the view by panning and zooming as described in step 6. Coevolving residues exceeding the threshold requirements of both p and residual are connected via edge lines as seen in Figure 6. If there are too many or too few edges connecting residues, change the Residual threshold (on the &#8220;StickWRLD Control&#8221; pane) to show fewer, or more, edges.</li>
	<li>Increase the residual threshold on the StickWRLD Control Pane until no IPD edge lines are shown and slowly ramp down until relationships appear. Continue increasing the residual until you have a sufficient number of relationships to examine.</li>
	<li>Identify relationships that involve either residues of known interest (e.g., within a motif or binding/functional site) or residues that are distal to one another within the alignment (suggesting that they are proximal in the folded protein)</li>
</ol>
8. Selecting and Saving Findings
<ol>
	<li>Using command + left click on any edges of interest. The StickWRLD Control pane will indicate the columns and connect specific residues, e.g., &#8220;(124|G) (136|H)&#8221; (Figure 7). Solid lines represent positive associations; dashed lines represent negative associations.</li>
	<li>Press the &#8220;Output Edges&#8221; button on the &#8220;StickWRLD Control&#8221; panel to save a plain-text formatted file (edge_residual.csv) of all of the visible edges, including the joined residues and their actual residual values, in the /StickWRLD/exec/ directory.</li>
</ol>

<ol>
	<li>Ray, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Addressing+the+unmet+need+for+visualizing+conditional+random+fields+in+biological+data.">Addressing the unmet need for visualizing conditional random fields in biological data.</a> BMC. 15, 202 (2014).</li><li>Sullivan, B. J., Durani, V., Magliery, T. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Triosephosphate+isomerase+by+consensus+design:+dramatic+differences+in+physical+properties+and+activity+of+related+variants.">Triosephosphate isomerase by consensus design: dramatic differences in physical properties and activity of related variants.</a> Journal of molecular biology. 413, 195-208 (2011).</li><li>Smith, M. A., Bedbrook, C. N., Wu, T., Arnold, F. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Hypocrea+jecorina+cellobiohydrolase+I+stabilizing+mutations+identified+using+noncontiguous+recombination.">Hypocrea jecorina cellobiohydrolase I stabilizing mutations identified using noncontiguous recombination.</a> ACS synthetic biology. 2, 690-696 (2013).</li><li>Ray, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Understanding+the+sequence+requirements+of+protein+families:+insights+from+the+BioVis+2013+contests.">Understanding the sequence requirements of protein families: insights from the BioVis 2013 contests.</a> BMC proceedings. 8, S1 (2014).</li><li>Eddy, S. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=What+is+a+hidden+Markov+model?.">What is a hidden Markov model?.</a> Nature biotechnology. 22, 1315-1316 (2004).</li><li>Ray, W. C., Ozer, H. G., Armbruster, D. W., Daniels, C. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Beyond+identity+-+when+classical+homology+searching+fails,+why,+and+what+you+can+do+about+it.">Beyond identity - when classical homology searching fails, why, and what you can do about it.</a> Proceedings of the 4th Ohio Collaborative Conference on Bioinformatics. , 51-56 (2009).</li><li>Ovchinnikov, S., Kamisetty, H., Baker, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Robust+and+accurate+prediction+of+residue-residue+interactions+across+protein+interfaces+using+evolutionary+information.">Robust and accurate prediction of residue-residue interactions across protein interfaces using evolutionary information.</a> eLife. 3, e02030 (2014).</li><li>Trudeau, D. L., Lee, T. M., Arnold, F. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Engineered+thermostable+fungal+cellulases+exhibit+efficient+synergistic+cellulose+hydrolysis+at+elevated+temperatures.">Engineered thermostable fungal cellulases exhibit efficient synergistic cellulose hydrolysis at elevated temperatures.</a> Biotechnology and bioengineering. 111, 2390-2397 (2014).</li><li>Ray, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=MAVL+and+StickWRLD:+visually+exploring+relationships+in+nucleic+acid+sequence+alignments.">MAVL and StickWRLD: visually exploring relationships in nucleic acid sequence alignments.</a> Nucleic acids research. 32, W59-W63 (2004).</li><li>. Python Language Reference v.2.7.6 Available from: <a target="_blank" href="https://www.python.org/download/releases/2.7.6/">https://www.python.org/download/releases/2.7.6/</a> (2014)</li><li>. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> PyOpenGL The Python OpenGL Binding. , (2014).</li><li>Larkin, M. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Clustal+W+and+Clustal+X+version+2.0.">Clustal W and Clustal X version 2.0.</a> Bioinformatics. 23, 2947-2948 (2007).</li><li>Ozer, H. G., Ray, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=MAVL/StickWRLD:+analyzing+structural+constraints+using+interpositional+dependencies+in+biomolecular+sequence+alignments.">MAVL/StickWRLD: analyzing structural constraints using interpositional dependencies in biomolecular sequence alignments.</a> Nucleic acids research. 34, W133-W136 (2006).</li><li>Ray, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=MAVL/StickWRLD+for+protein:+visualizing+protein+sequence+families+to+detect+non-consensus+features.">MAVL/StickWRLD for protein: visualizing protein sequence families to detect non-consensus features.</a> Nucleic acids research. 33, W315-W319 (2005).</li><li>Ray, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+Visual+Analytics+approach+to+identifying+protein+structural+constraints.">A Visual Analytics approach to identifying protein structural constraints.</a> IEEE. , 249-250 (2010).</li><li>Wheeler, T. J., Clements, J., Finn, R. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Skylign:+a+tool+for+creating+informative,+interactive+logos+representing+sequence+alignments+and+profile+hidden+Markov+models.">Skylign: a tool for creating informative, interactive logos representing sequence alignments and profile hidden Markov models.</a> BMC bioinformatics. 15, 7 (2014).</li><li>Perrier, V., Burlacu-Miron, S., Bourgeois, S., Surewicz, W. K., Gilles, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genetically+engineered+zinc-chelating+adenylate+kinase+from+Escherichia+coli+with+enhanced+thermal+stability.">Genetically engineered zinc-chelating adenylate kinase from Escherichia coli with enhanced thermal stability.</a> The Journal of biological chemistry. 273, 19097-19101 (1998).</li></ol>

The authors declare that they have no competing financial interests.

StickWRLD has been successfully used to identify such IPDs in the Adenylate Kinase Lid domain16, as well as associated DNA bases in rho-dependent terminators9, and a novel splice-site specificity in Archaeal tRNA intron endonuclease6 target sites. These IPDs were not detectable via a direct examination of the alignments.
StickWRLD displays each position of an alignment as a column of 20 &#8220;spheres&#8221;, where each sphere represents one of the 20 amino acid residues and the size of the sphere indicates the frequency of occurrence of that particular residue within that column (Figure 4). Columns are arranged in a cylinder, with edge lines connecting residues in different columns (indicating an IPD). These edge lines are only drawn if the corresponding residues are covarying at a frequency surpassing both the p-value (significance) and residual (expected - observed) thresholds.
Detection of co-occurring interdependent residues, or IPDs, in distal regions of a DNA or protein sequence alignment is difficult using standard sequence alignment tools6. While such tools generate a consensus, or motif, sequence, this consensus is in many cases a simple majority-rule averaging and does not convey covariation relationships that may form one or more sub-motifs &#8211; groups of residues that tend to co-evolve. Even HMM models, which are capable of detecting neighboring dependencies, cannot accurately model sequence motifs with distal IPDs5. The result is that the calculated consensus may in fact be a &#8220;synthetic&#8221; sequence not found in nature &#8211; and engineered proteins based on such computational consensus may not, in fact, be optimal. In fact, the Pfam HMM for ADK would suggest that a chimeric protein containing half of the tetracysteine motif, and half of the H,S,D,T motif, is functionally just as acceptable as any actually existing ADK. This is not the case, as such chimeras (and many other blendings of these motifs) Are catalytically dead4,19.
When looking for correlations, it is critical that the residual threshold be adjusted to allow for the discovery of relevant correlations by setting the threshold above the level at which any edges are seen and then gradually ramping the threshold back down. This ensures that only the most significant edges are considered initially.
An alternate approach is to start with the residual threshold set extremely low. This results in the display of all significant edges. From here the residual threshold can slowly be increased, allowing edges to drop out until patterns emerge. While this approach is less useful when looking for the inclusion of specific nodes (e.g., application of domain knowledge), it allows for the discovery of unexpected relationships using StickWRLD as a visual analytical tool to discover emerging patterns in the data visualization.
StickWRLD is limited primarily by the available memory of the system on which it is run as well as the resolution of the display device. While there is no theoretical limit to the number of data points StickWRLD can examine, and sequences up to 20,000 positions have been tested, in practice StickWRLD performs best with sequences up to around 1,000 positions.
The primary advantage of StickWRLD lies in its ability to identify groups of residues that covary with one another. This is a significant advantage over the traditional approach of the statistical consensus sequence, which is a simple statistical averaging and does not take coevolution into account. While in some cases covarying residues may simply be an artifact of phylogeny, even these residues have withstood the &#8220;test of selection&#8221;, and as such are unlikely to detract from the functionality of any protein engineered to include them.
While using StickWRLD to identify IPDs in a canonical DNA or protein sequence consensus/motif prior to engineering synthetic variants will reduce the potential for error and support rapid optimization of function, it should be noted that StickWRLD can be used as a generalized correlation identification tool and is not limited exclusively to protein data. StickWRLD can be used to visually discover the co-occurance of any variables in any properly encoded data set.

Protein alignments have long been used to evaluate the similarity of residues in a protein family. Frequently the most interesting features of a protein (e.g., catalytic or other binding sites) are the result of protein folding bringing distal regions of the linear sequence into contact, and as a result these apparently unrelated regions in the alignment tend to evolve and change in a coordinated fashion. In other cases, the function of a protein can be dependent on its electrostatic signature, and mutations that affect the electronic dipole are compensated for by changes to distant charged residues. Allosteric effects can also induce long-range sequential and spatial dependencies between residue identities. Regardless of their origin, these functionally-required covariations of residues - inter-positional dependencies (IPDs) - may not be obvious with visual examination of the alignment (Figure 1). Identification of IPDs - as well as of which specific residues within those positions tend to covary as a unit - can reveal important clues about the processes of protein folding and the formation of functional sites. This information can then be used to optimize synthetic (engineered) proteins in terms of thermostability and activity. It has long been known that not all point mutations toward consensus provide improved stability or activity. More recently, proteins designed to take advantage of known IPDs in their sequence have been shown to result in greater activity than the same protein designed strictly from consensus1,2 (manuscript in preparation), similar to the idea of stabilizing point mutations3.
Unfortunately, traditional consensus-building models (e.g., majority rule) only capture IPDs by accident. Consensus and Position Specific Scoring Matrix methods are ignorant of IPDs and only &#39;correctly&#39; include them in models, when the dependent residues are also the most popular residues for those positions in the family. Markov Chain models can capture IPDs when they are sequentially proximal, but their typical implementation ignores everything except immediate sequential neighbors, and even at their best, Hidden Markov Model calculations (see Figure 2) become intractable when dependencies are separated in the sequence by more than a dozen or so positions4. Since these IPDs essentially form &#8220;sub-motifs&#8221; which cannot be predicted by simple &#8220;majority rule&#8221; or even HMM-based consensus models5,6 the result can be a biologically invalid &#8220;consensus&#8221; which is not only never seen in nature but is less viable than any extant protein. Systems based on Markov Random Fields, such as GREMLIN7, attempt to overcome these problems. Additionally while sophisticated biological/biochemical techniques such as noncontiguous recombination3,8 can be used to identify essential protein elements by region, they require considerable time and bench work for single-base-pair precision to be achieved.
StickWRLD9 is a Python based program that creates an interactive 3D representation of a protein alignment that makes IPDs clear and easy to understand. Each position in the alignment is represented as a column in the display, where each column is comprised of a stack of spheres, one for each of the 20 amino acids that could be present in that position within the alignment. The sphere size is dependent on the frequency of occurrence of the amino acid, such that the user can immediately glean the consensus residue or the relative distribution of amino acids within that position by simply looking at the size of the spheres. The columns representing each position are wrapped around a cylinder. This gives every sphere representing a possible amino acid at each position in the alignment, a clear &#39;line of sight&#39; to every other amino acid possibility at every other position. Prior to visualization, StickWRLD calculates the the correlation strength between all possible combinations of residues to identify the IPDs9. To represent IPDs, lines are drawn between residues which are coevolving at a higher, or lower than would be expected if the residues present in the positions were independent(IPDs).
Not only does this visualization show which sequence positions interact evolutionarily, but as the IPD edge lines are drawn between the amino acid spheres in each column, the user can quickly determine which specific amino acids tend to be coevolving at each position. The user has the ability to rotate and explore the visualized IPD structure, as well as dynamically change the statistical thresholds controlling the display of correlations, making StickWRLD a powerful discovery tool for IPDs.
Applications such as GREMLIN7 similarly display complex relational information between residues &#8211; but these relationships are computed via more traditional Markov models, which are not designed to determine any conditional relationships. As such, these are capable of being displayed as 2D projections. By contrast, StickWRLD can compute and display multi-node conditional dependencies, which may be obfuscated if rendered as a 2D graph (a phenomenon known as edge occlusion).
StickWRLD&#8217;s 3D view also has several other advantages. By allowing users to manipulate the visual &#8211; panning, rotating, and zooming &#8211; features that may be obfuscated or unintuitive in a 2D representation can be more easily seen in the 3D cylinder of StickWRLD. StickWRLD is essentially a visual analytics tool, harnessing the power of the human brain&#8217;s pattern recognition ability to see patterns and trends, and the ability to explore the data from various perspectives lends itself to this.

<table border="0" cellpadding="0" cellspacing="0">
	<tbody>
		<tr>
			<td>Mac or Ubuntu OS computer</td>
			<td>Various</td>
			<td>NA</td>
			<td>Any Mac or linux (e.g. Ubuntu) computer capable of running python &#38; associated shell scripts</td>
		</tr>
		<tr>
			<td>Python programming language</td>
			<td>python.org</td>
			<td>NA</td>
			<td>Python version 2.7.6 or greater recommended</td>
		</tr>
		<tr>
			<td>wxPython library</td>
			<td>wxpython.org</td>
			<td>NA</td>
			<td>Latest version recommended</td>
		</tr>
		<tr>
			<td>SciPy library</td>
			<td>scipy.org</td>
			<td>NA</td>
			<td>Latest version recommended</td>
		</tr>
		<tr>
			<td>PyOpenGL library</td>
			<td>pyopengl.sourceforge.net</td>
			<td>NA</td>
			<td>Latest version recommended</td>
		</tr>
		<tr>
			<td>StickWRLD Python scripts</td>
			<td>NCH BCCM</td>
			<td>NA</td>
			<td>Available from http://www.stickwrld.org</td>
		</tr>
		<tr>
			<td>fasta2stick.sh file converter</td>
			<td>NCH BCCM</td>
			<td>NA</td>
			<td>Available from http://www.stickwrld.org</td>
		</tr>
		<tr>
			<td>Protein and/or DNA sequence data</td>
			<td>NA</td>
			<td>NA</td>
			<td>Samples available at http://www.stickwrld.org</td>
		</tr>
	</tbody>
</table>

optimization of synthetic proteins: identification of interpositional dependencies indicating structurally and/or functionally linked residues

Protein alignments are commonly used to evaluate the similarity of protein residues, and the derived consensus sequence used for identifying functional units (e.g., domains). Traditional consensus-building models fail to account for interpositional dependencies &#8211; functionally required covariation of residues that tend to appear simultaneously throughout evolution and across the phylogentic tree. These relationships can reveal important clues about the processes of protein folding, thermostability, and the formation of functional sites, which in turn can be used to inform the engineering of synthetic proteins. Unfortunately, these relationships essentially form sub-motifs which cannot be predicted by simple &#8220;majority rule&#8221; or even HMM-based consensus models, and the result can be a biologically invalid &#8220;consensus&#8221; which is not only never seen in nature but is less viable than any extant protein. We have developed a visual analytics tool, StickWRLD, which creates an interactive 3D representation of a protein alignment and clearly displays covarying residues. The user has the ability to pan and zoom, as well as dynamically change the statistical threshold underlying the identification of covariants. StickWRLD has previously been successfully used to identify functionally-required covarying residues in proteins such as Adenylate Kinase and in DNA sequences such as endonuclease target sites.

StickWRLD has been used previously to detect interpositional dependencies (IPDs) between residues in both DNA3 and protein15-17 alignments. These co-evolving residues, while often distal from one another in the sequence alignment, are often proximal to one another in the folded protein. StickWRLD allows rapid discovery of residue-specific co-occurrence at such sites, e.g., an alanine at position &#8220;x&#8221; is strongly correlated to a threonine at position &#8220;y&#8221;. Such correlations can be indicative of provable structural relationships, and typically are sites that, by necessity, co-evolve. StickWRLD is able to detect these relationships even when more &#8220;traditional&#8221; approaches using HMMs to describe motifs fail. For example, analysis of the PFAM alignment of the ADK lid domain using StickWRLD reveals a strong positive correlation between cysteines (C) at positions 4 and 8 and a coordinated pair of C at positions 35 and 38. At the same time, StickWRLD showed a similar strong positive relationship between histidine (H) and serine (S) at 4 and 8, with a strong negative relationships between these and the C quartet at 4, 8, 35, and 38, and a strong positive relationship with aspartic acid (D) and threonine (T) at positions 35 and 38 respectively. Additional IPDs exist between the H,S,D,T motif and a T and G at position **** 10 and 29 in b subtilis **** highlighting the conditional nature of these IPDs - the tetracysteine motif does not &#39;care&#39; about the identities at these two positions, while the hydrophilic H,S,D,T triad requires specific residues in these positions almost absolutely. These two completely different position-dependent residue motifs can fulfill the same role the ADK lid. As can be seen in Figure 6, a large cluster of IPDs, including a 3-node association between G (glycine) at position 132, Y (tyrosine) at position 135, and a P (proline) at position 141, is visible in the foreground (Figure 6A). In Figure 6B, the view has been skewed to position the user slightly above the cylinder, revealing an IPD between an H (histidine) at position 136 and an M (methionine) at position 29, 107 residues distant. A PFAM HMM-derived motif of the same domain (Figure 2), meanwhile, not only does not detect these as specifically co-occuring motif variants, but also defines the overall groupings in a biologically unsupported scheme16.
<img alt="Figure 1" src="/files/ftp_upload/52878/52878fig1highres.jpg" width="500" /> 
Figure 1. &#8220;Subway Map&#8221; representation of the B. subtilis Adenosine Kinase (ADK) Lid domain structure. Arrows indicate IPDs identified in the PFAM alignment of ADK Lid domain by StickWRLD. StickWRLD is able to correctly identify IPDs within a cluster of residues which are in close proximity in the folded protein. Of particular interest are the T and G pair at positions 9 and 29, which only form an IPD when the tetrad of residues at 4, 7, 24, and 27 is not C,C,C,C). Residue numbers displayed represents B. subtilis position and not PFAM alignment positions. <a href="https://jove.unicartagenaproxy.elogim.com/files/ftp_upload/52878/52878fig1highres.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" src="/files/ftp_upload/52878/52878fig2highres.jpg" width="700" /> 
Figure 2. Skylign18 Hidden Markov Model (HMM) Sequence Logo for the ADK lid domain. While HMMs are powerful tools for determining probabilities at each position as well as the contribution of each site to the overall model, the positional independence of HMMs makes them unsuitable for detecting IPDs. This model does not suggest any of the dependencies seen in the StickWRLD representations (Figure 6). <a href="https://jove.unicartagenaproxy.elogim.com/files/ftp_upload/52878/52878fig2highres.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" src="/files/ftp_upload/52878/52878fig3.jpg" /> 
Figure 3. The StickWRLD Data Loader. Users can choose from existing demo data or load their own data in the form of DNA or Protein sequence alignments.
<img alt="Figure 4" src="/files/ftp_upload/52878/52878fig4highres.jpg" width="600" /> 
Figure 4. The StickWRLD Control window. The Control pane allows the user to change various view properties as well as regulate the thresholds controlling the display of edge lines indicating relationships between residues (IPDs). Circled in red are the defaults that typically need to be adjusted for best viewing of any dataset. The Residual value sets the threshold of (observed-expected) for which connector/association lines are drawn. The controls for Column and Ball labels control whether or not the column position and residue values (e.g., &#8220;A&#8221; for arginine) are displayed. The Column Edge Line control toggles on and off the display of edge lines connecting columns &#8211; for dense data sets this is better turned off. The Column Thickness controls whether or not the column itself is displayed &#8211; setting this to a very small value (e.g., 0.1) will draw a line through the spheres in the column, making it easy to distinguish the columns from one another. <a href="https://jove.unicartagenaproxy.elogim.com/files/ftp_upload/52878/52878fig4highres.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 5" src="/files/ftp_upload/52878/52878fig5highres.jpg" width="600" /> 
Figure 5. Initial view of the StickWRLD OpenGL window with the Adenylate Kinase lid domain protein data set loaded. The initial perspective looks &#8220;down&#8221; through the cylinder comprised of the sequence alignment positions. The user can rotate the cylinder using left-mouse-click-drag, and zoom in/out using right-mouse-click-drag. The initial view is quite dense because the default display shows even small rates of co-evolution. For many proteins, at this setting, distinct modules can be detected, but even in densely co-evolving proteins the display can be rapidly and interactively simplified to find the most important IPDs using the StickWRLD interface. <a href="https://jove.unicartagenaproxy.elogim.com/files/ftp_upload/52878/52878fig5highres.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 6" src="/files/ftp_upload/52878/52878fig6highres.jpg" width="700" /> 
Figure 6. Closeup view of a StickWRLD visualization of the Adenylate Kinase lid domain protein. Here we have changed the default Residual to 0.2. This increases the threshold for display of inter-residue edges, showing fewer edges. The edges that remain indicate strongly associated IPDs. In addition the view has been rotated and zoomed to allow for easier viewing of the edges. (A) A large cluster of IPDs is visible in the foreground, including a 3-node association between G (glycine) at position 132, Y (tyrosine) at position 135, and a P (proline) at position 141. (B) The view has been skewed to position the user slightly above the cylinder, revealing an IPD between an H (histidine) at position 136 and an M (methionine) at position 29, 107 residues distant. <a href="https://jove.unicartagenaproxy.elogim.com/files/ftp_upload/52878/52878fig6highres.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 7" src="/files/ftp_upload/52878/52878fig7.jpg" /> 
Figure 7. StickWRLD Control window lower-right information view. CTRL+Left clicking on an object (e.g., sphere or edge) in the OpenGL window displays the information for the object in the lower right of the StickWLRD Control window. Here we see the information for an IPD edge between a methionine at position 29 and a histidine at position 136.

Watch this Scientific Journal Video about Optimization of Synthetic Proteins: Identification of Interpositional Dependencies Indicating Structurally and/or Functionally Linked Residues at JoVE.com

Optimization of Synthetic Proteins: Identification of Interpositional Dependencies Indicating Structurally and/or Functionally Linked Residues

Synthetic protein sequences based on consensus motifs typically ignore co-evolving residues, that imply interpositional dependencies (IPDs). IPDs can be essential to activity, and designs that disregard them may result in suboptimal results. This protocol uses StickWRLD to identify IPDs and help inform rational protein design, resulting in more efficient results.

Protein alignments are commonly used to evaluate the similarity of protein residues, and the derived consensus sequence used for identifying functional ...

optimization-synthetic-proteins-identification-interpositional

Universidad de Cartagena UDEC

Research

JoVE Journal

Chemistry

7.2K Views.  The Research Institute at Nationwide Children's Hospital. Synthetic protein sequences based on consensus motifs typically ignore co-evolving residues, that imply interpositional dependencies (IPDs). IPDs can be essential to activity, and designs that disregard them may result in suboptimal results. This protocol uses StickWRLD to identify IPDs and help inform rational protein design, resulting in more efficient results.

Video: Optimization of Synthetic Proteins: Identification of Interpositional Dependencies Indicating Structurally and/or Functionally Linked Residues

Qualitative Identification of Carboxylic Acids, Boronic Acids, and Amines Using Cruciform Fluorophores

Molecular cruciforms are X-shaped systems in which two conjugation axes intersect at a central core. If one axis of these molecules is substituted with electron-donors, and the other with electron-acceptors, cruciforms&#39; HOMO will localize along the electron-rich and LUMO along the electron-poor axis. This spatial isolation of cruciforms&#39; frontier molecular orbitals (FMOs) is essential to their use as sensors, since analyte binding to the cruciform invariably changes its HOMO-LUMO gap and the associated optical properties. Using this principle, Bunz and Miljani&#263; groups developed 1,4-distyryl-2,5-bis(arylethynyl)benzene and benzobisoxazole cruciforms, respectively, which act as fluorescent sensors for metal ions, carboxylic acids, boronic acids, phenols, amines, and anions. The emission colors observed when these cruciform are mixed with analytes are highly sensitive to the details of analyte&#39;s structure and - because of cruciforms&#39; charge-separated excited states - to the solvent in which emission is observed. Structurally closely related species can be qualitatively distinguished within several analyte classes: (a) carboxylic acids; (b) boronic acids, and (c) metals. Using a hybrid sensing system composed from benzobisoxazole cruciforms and boronic acid additives, we were also able to discern among structurally similar: (d) small organic and inorganic anions, (e) amines, and (f) phenols. The method used for this qualitative distinction is exceedingly simple. Dilute solutions (typically 10-6 M) of cruciforms in several off-the-shelf solvents are placed in UV/Vis vials. Then, analytes of interest are added, either directly as solids or in concentrated solution. Fluorescence changes occur virtually instantaneously and can be recorded through standard digital photography using a semi-professional digital camera in a dark room. With minimal graphic manipulation, representative cut-outs of emission color photographs can be arranged into panels which permit quick naked-eye distinction among analytes. For quantification purposes, Red/Green/Blue values can be extracted from these photographs and the obtained numeric data can be statistically processed.

Cross-conjugated cruciform fluorophores based on 1,4-distyryl-2,5-bis(arylethynyl)benzene and benzobisoxazole nuclei can be used to qualitatively identify diverse Lewis acidic and Lewis basic analytes. This method relies on the differences in emission colors of the cruciforms that are observed upon analyte addition. Structurally closely related species can be distinguished from each other.

Molecular cruciforms are X-shaped systems in which two conjugation axes intersect at a central core. If one axis of these molecules is substituted with ...

Towards Biomimicking Wood: Fabricated Free-standing Films of Nanocellulose, Lignin, and a Synthetic Polycation

Woody materials are comprised of plant cell walls that contain a layered secondary cell wall composed of structural polymers of polysaccharides and lignin. Layer-by-layer (LbL) assembly process which relies on the assembly of oppositely charged molecules from aqueous solutions was used to build a freestanding composite film of isolated wood polymers of lignin and oxidized nanofibril cellulose (NFC). To facilitate the assembly of these negatively charged polymers, a positively charged polyelectrolyte, poly(diallyldimethylammomium chloride) (PDDA), was used as a linking layer to create this simplified model cell wall. The layered adsorption process was studied quantitatively using quartz crystal microbalance with dissipation monitoring (QCM-D) and ellipsometry. The results showed that layer mass/thickness per adsorbed layer increased as a function of total number of layers. The surface coverage of the adsorbed layers was studied with atomic force microscopy (AFM). Complete coverage of the surface with lignin in all the deposition cycles was found for the system, however, surface coverage by NFC increased with the number of layers. The adsorption process was carried out for 250 cycles (500 bilayers) on a cellulose acetate (CA) substrate. Transparent free-standing LBL assembled nanocomposite films were obtained when the CA substrate was later dissolved in acetone. Scanning electron microscopy (SEM) of the fractured cross-sections showed a lamellar structure, and the thickness per adsorption cycle (PDDA-Lignin-PDDA-NC) was estimated to be 17 nm for two different lignin types used in the study. The data indicates a film with highly controlled architecture where nanocellulose and lignin are spatially deposited on the nanoscale (a polymer-polymer nanocomposites), similar to what is observed in the native cell wall.

The objective of this research was to form synthetic plant cell wall tissue using layer-by-layer assembly of nanocellulose fibrils and isolated lignin assembled from dilute aqueous suspensions.&#160; Surface measurement techniques of quartz crystal microbalance and atomic force microscopy were used to monitor the formation of the polymer-polymer nanocomposite material.

Woody materials are comprised of plant cell walls that contain a layered secondary cell wall composed of structural polymers of polysaccharides and ...

Stabilizing Hepatocellular Phenotype Using Optimized Synthetic Surfaces

Currently, one of the major limitations in cell biology is maintaining differentiated cell phenotype. Biological matrices are commonly used for culturing and maintaining primary and pluripotent stem cell derived hepatocytes. While biological matrices are useful, they permit short term culture of hepatocytes, limiting their widespread application. We have attempted to overcome the limitations using a synthetic polymer coating. Polymers represent one of the broadest classes of biomaterials and possess a wide range of mechanical, physical and chemical properties, which can be fine-tuned for purpose. Importantly, such materials can be scaled to quality assured standards and display batch-to-batch consistency. This is essential if cells are to be expanded for high through-put screening in the pharmaceutical testing industry or for cellular based therapy. Polyurethanes (PUs) are one group of materials that have shown promise in cell culture. Our recent progress in optimizing a polyurethane coated surface, for long-term culture of human hepatocytes displaying stable phenotype, is presented and discussed.

This article will focus on developing polymer coated surfaces for long-term, stable culture of stem cell derived human hepatocytes.

Currently, one of the major limitations in cell biology is maintaining differentiated cell phenotype. Biological matrices are commonly used for culturing ...

Diffuse Reflectance Infrared Spectroscopic Identification of Dispersant/Particle Bonding Mechanisms in Functional Inks

In additive manufacturing, or 3D printing, material is deposited drop by drop, to create micron to macroscale layers. A typical inkjet ink is a colloidal dispersion containing approximately ten components including solvent, the nano to micron scale particles which will comprise the printed layer, polymeric dispersants to stabilize the particles, and polymers to tune layer strength, surface tension and viscosity. To rationally and efficiently formulate such an ink, it is crucial to know how the components interact. Specifically, which polymers bond to the particle surfaces and how are they attached? Answering this question requires an experimental procedure that discriminates between polymer adsorbed on the particles and free polymer. Further, the method must provide details about how the functional groups of the polymer interact with the particle. In this protocol, we show how to employ centrifugation to separate particles with adsorbed polymer from the rest of the ink, prepare the separated samples for spectroscopic measurement, and use Diffuse Reflectance Fourier Transform Infrared Spectroscopy (DRIFTS) for accurate determination of dispersant/particle bonding mechanisms. A significant advantage of this methodology is that it provides high level mechanistic detail using only simple, commonly available laboratory equipment. This makes crucial data available to almost any formulation laboratory. The method is most useful for inks composed of metal, ceramic, and metal oxide particles in the range of 100 nm or greater. Because of the density and particle size of these inks, they are readily separable with centrifugation. Further, the spectroscopic signatures of such particles are easy to distinguish from absorbed polymer. The primary limitation of this technique is that the spectroscopy is performed ex-situ on the separated and dried particles as opposed to the particles in dispersion. However, results from attenuated total reflectance spectra of the wet separated particles provide evidence for the validity of the DRIFTS measurement.

Formulation of stable, functional inks is critical to expanding the applications of additive manufacturing. In turn, knowledge of the mechanisms of dispersant/particle bonding is required for effective ink formulation. Diffuse Reflectance Fourier Transform Infrared Spectroscopy (DRIFTS) is presented as a simple, inexpensive way to gain insight into these mechanisms.

In additive manufacturing, or 3D printing, material is deposited drop by drop, to create micron to macroscale layers.  A typical inkjet ink is a colloidal ...

Curtain Flow Column: Optimization of Efficiency and Sensitivity

Active Flow Technology (AFT) is a form of column technology that increases the separation performance of a HPLC column through the use of a specially purpose built multiport end-fitting(s). Curtain Flow (CF) columns belong to the AFT suite of columns, specifically the CF column is designed so that the sample is injected into the radial central region of the bed and a curtain flow of mobile phase surrounding the injection of solute prevents the radial dispersion of the sample to the wall. The column functions as an 'infinite diameter' column. The purpose of the design is to overcome the radial heterogeneity of the column bed, and at the same time maximize the sample load into the radial central region of the column bed, which serves to increase detection sensitivity. The protocol described herein outlines the system and CF column set up and the tuning process for an optimized infinite diameter 'virtual' column.

Here, we present a protocol for the operation and optimization of Active Flow Technology (AFT) column in Curtain flow (CF) mode for enhanced separation performance.

Active Flow Technology (AFT) is a form of column technology that increases the separation performance of a HPLC column through the use of a specially ...

The Preparation and Properties of Thermo-reversibly Cross-linked Rubber Via Diels-Alder Chemistry

A method for using Diels Alder thermo-reversible chemistry as cross-linking tool for rubber products is demonstrated. In this work, a commercial ethylene-propylene rubber, grafted with maleic anhydride, is thermo-reversibly cross-linked in two steps. The pending anhydride moieties are first modified with furfurylamine to graft furan groups to the rubber backbone. These pendant furan groups are then cross-linked with a bis-maleimide via a Diels-Alder coupling reaction. Both reactions can be performed under a broad range of experimental conditions and can easily be applied on a large scale. The material properties of the resulting Diels-Alder cross-linked rubbers are similar to a peroxide-cured ethylene/propylene/diene rubber (EPDM) reference. The cross-links break at elevated temperatures (&#62; 150 &#176;C) via the retro-Diels-Alder reaction and can be reformed by thermal annealing at lower temperatures (50-70 &#176;C). Reversibility of the system was proven with infrared spectroscopy, solubility tests and mechanical properties. Recyclability of the material was also shown in a practical way, i.e., by cutting a cross-linked sample into small parts and compression molding them into new samples displaying comparable mechanical properties, which is not possible for conventionally cross-linked rubbers.

A simple two-step approach involving rubber modification and cross-linking yields fully reworkable, elastic rubber products.

A method for using Diels Alder thermo-reversible chemistry as cross-linking tool for rubber products is demonstrated. In this work, a commercial ...

Synthesis of Platinum-nickel Nanowires and Optimization for Oxygen Reduction Performance

Platinum-nickel (Pt-Ni) nanowires were developed as fuel cell electrocatalysts, and were optimized for the performance and durability in the oxygen reduction reaction. Spontaneous galvanic displacement was used to deposit Pt layers onto Ni nanowire substrates. The synthesis approach produced catalysts with high specific activities and high Pt surface areas. Hydrogen annealing improved Pt and Ni mixing and specific activity. Acid leaching was used to preferentially remove Ni near the nanowire surface, and oxygen annealing was used to stabilize near-surface Ni, improving durability and minimizing Ni dissolution. These protocols detail the optimization of each post-synthesis processing step, including hydrogen annealing to 250 &#176;C, exposure to 0.1 M nitric acid, and oxygen annealing to 175 &#176;C. Through these steps, Pt-Ni nanowires produced increased activities more than an order of magnitude than Pt nanoparticles, while offering significant durability improvements. The presented protocols are based on Pt-Ni systems in the development of fuel cell catalysts. These techniques have also been used for a variety of metal combinations, and can be applied to develop catalysts for a number of electrochemical processes.

&#160; The protocol describes the synthesis and electrochemical testing of platinum-nickel nanowires. Nanowires were synthesized by the galvanic displacement of a nickel nanowire template. Post-synthesis processing, including hydrogen annealing, acid leaching, and oxygen annealing were used to optimize nanowire performance and durability in the oxygen reduction reaction.

Platinum-nickel (Pt-Ni) nanowires were developed as fuel cell electrocatalysts, and were optimized for the performance and durability in the oxygen ...

Color Spot Test As a Presumptive Tool for the Rapid Detection of Synthetic Cathinones

Synthetic cathinones are a large class of new psychoactive substances (NPS) that are increasingly prevalent in drug seizures made by law enforcement and other border protection agencies globally. Color testing is a presumptive identification technique indicating the presence or absence of a particular drug class using rapid and uncomplicated chemical methods. Owing to their relatively recent emergence, a color test for the specific identification of synthetic cathinones is not currently available. In this study, we introduce a protocol for the presumptive identification of synthetic cathinones, employing three aqueous reagent solutions: copper(II) nitrate, 2,9-dimethyl-1,10-phenanthroline (neocuproine) and sodium acetate. Small pin-head sized amounts (approximately 0.1-0.2 mg) of the suspected drugs are added to the wells of a porcelain spot plate, and each reagent is then added dropwise sequentially before heating on a hotplate. A color change from very light blue to yellow-orange after 10 min indicates the likely presence of synthetic cathinones. The highly stable and specific test reagent has the potential for use in the presumptive screening of unknown samples for synthetic cathinones in a forensic laboratory. However, the nuisance of an added heating step for the color change result limits the test to laboratory application and decreases the likelihood of an easy translation to field testing.

Here we present a simple, inexpensive, and selective chemical spot test protocol for the detection of synthetic cathinones, a class of new psychoactive substances. The protocol is suitable for use in various areas of law enforcement that encounter illicit material.

Synthetic cathinones are a large class of new psychoactive substances (NPS) that are increasingly prevalent in drug seizures made by law enforcement and ...

Characterization of Synthetic Polymers via Matrix Assisted Laser Desorption Ionization Time of Flight (MALDI-TOF) Mass Spectrometry

There are many techniques that can be employed in the characterization of synthetic homopolymers, but few provide as useful of information for end group analysis as matrix-assisted laser desorption ionization time of flight mass spectrometry (MALDI-TOF MS). This tutorial demonstrates methods for optimization of the sample preparation, spectral acquisition, and data analysis of synthetic polymers using MALDI-TOF MS. Critical parameters during sample preparation include the selection of the matrix, identification of an appropriate cationization salt, and tuning the relative proportions of the matrix, cation, and analyte. The acquisition parameters, such as mode (linear or reflector), polarization (positive or negative), acceleration voltage, and delay time, are also important. Given some knowledge of the chemistry involved to synthesize the polymer and optimizing both the data acquisition parameters and the sample preparation conditions, spectra should be obtained with sufficient resolution and mass accuracy to enable the unambiguous determination of the end groups of most homopolymers (masses below 10,000) in addition to the repeat unit mass and the overall molecular weight distribution. Though demonstrated on a limited set of polymers, these general techniques are applicable to a much wider range of synthetic polymers for determining mass distributions, though end group determination is only possible for homopolymers with narrow dispersity.

Characterization of Synthetic Polymers via Matrix Assisted Laser Desorption Ionization Time of Flight MALDI-TOF Mass Spectrometry

A protocol for the matrix-assisted laser desorption ionization time of flight mass spectrometry (MALDI-TOF MS) characterization of synthetic polymers is described including the optimization of sample preparation, spectral acquisition, and data analysis.

There are many techniques that can be employed in the characterization of synthetic homopolymers, but few provide as useful of information for end group ...

A Facile Synthetic Method to Obtain Bismuth Oxyiodide Microspheres Highly Functional for the Photocatalytic Processes of Water Depuration

Bismuth oxyhalide (BiOI) is a promising material for sunlight-driven-environmental photocatalysis. Given that the physical structure of this kind of materials is highly related to its photocatalytic performance, it is necessary to standardize the synthetic methods in order to obtain the most functional architectures and, thus, the highest photocatalytic efficiency. Here, we report a reliable route to obtain BiOI microspheres via the solvothermal process, using Bi(NO3)3 and potassium iodide (KI) as precursors, and ethylene glycol as a template. The synthesis is standardized in a 150 mL autoclave, at 126 &#176;C for 18 h. This results in 2-3 &#181;m-sized mesoporous microspheres, with a relevant specific surface area (61.3 m2/g). Shortening the reaction times in the synthesis results in amorphous structures, while higher temperatures lead to a slight increase in the porosity of the microspheres, with no effect in the photocatalytic performance. The materials are photo-active under UV-A/visible light irradiation for the degradation of the antibiotic ciprofloxacin in water. This method has demonstrated to be effective in interlaboratory tests, obtaining similar BiOI microspheres in Mexican and Chilean research groups.

This article describes a synthetic method to obtain bismuth oxyiodide microspheres, which are highly functional to perform the photocatalytic removal of organic pollutants, such as ciprofloxacin, in water under UV-A/visible light irradiation.

Bismuth oxyhalide (BiOI) is a promising material for sunlight-driven-environmental photocatalysis. Given that the physical structure of this kind of ...