We thank Dr. Guido Scarabelli and Hongyang Li for extensive testing throughout development as well as the Bio3D user community and the University of Bergen structural bioinformatics workshop participants for feedback and comments that have improved this application.

NOTE: A typical Bio3D-web session proceeds through five consecutive and dependent steps (see Figure 1 for a schematic representation). Each step is implemented as a consecutive navigation tab of the web application namely SEARCH, ALIGN, FIT, PCA, and eNMA.
1. Structure Search and Selection (SEARCH)
<ol>
	<li>Input structure
	<ol>
		<li>Obtain the PDB ID of adenylate kinase (Adk), e.g. by searching the PDB [http://www.rcsb.org/pdb]. Alternatively, obtain the protein amino acid sequence of interest, e.g. from UniProt [http://uniprot.org].</li>
		<li>Enter the four character long PDB ID for Adk (e.g. 1AKE), or paste a protein sequence, to the text box in the &#34;Input structure or sequence&#34; panel.</li>
	</ol>
	</li>
	<li>Hit selection
	<ol>
		<li>Click the blue &#34;Next&#34; (Hit selection) button in the first panel or simply scroll down to panel B) &#34;Hit selection&#34; for further analysis.</li>
		<li>Make sure the &#34;Limit total number of included structures&#34; slider is set to its maximum value to include all structures above the cutoff.</li>
		<li>Lower the &#34;Adjust the inclusion BitScore cutoff&#34; to include more distantly related hits, or increase it to exclude.</li>
	</ol>
	</li>
	<li>Optional hit filtering
	<ol>
		<li>Click the blue &#34;Next&#34; (Hit selection) button in the first panel or simply scroll down to panel C) &#34;Optional filtering of related structures for further analysis&#34;.</li>
		<li>Make sure the selected hits represent relevant structures by inspecting details of the table, e.g. PDB name, species, and bound ligands.</li>
		<li>Manually refine the selected subset of structures if needed by clicking the rows of the table. 
		&#8203;NOTE: Rows highlighted with a blue color depict PDB IDs selected for further analysis in subsequent tabs.</li>
	</ol>
	</li>
</ol>
2. Multiple Sequence Alignment Analysis (ALIGN)
<ol>
	<li>Click the ALIGN tab to perform sequence alignment of the selected structures from the SEARCH tab.</li>
	<li>Alignment summary
	<ol>
		<li>Review the alignment summary in panel A) &#34;Alignment summary&#34;. Make sure that the regions of interest are aligned and not masked by gaps in one or more structures.</li>
		<li>If needed, toggle the &#34;Display alignment editing options&#34; and remove unwanted PDB IDs, e.g. PDBs with missing residues.</li>
	</ol>
	</li>
	<li>Sequence alignment analysis
	<ol>
		<li>Click the blue &#34;Next&#34; (Analysis) button to perform sequence-based clustering analysis of the collected structures.</li>
		<li>Select the plot option Dendrogram. Adjust the Cluster into K groups slider to partition the structures into k groups.</li>
		<li>Optionally change the clustering method if desired by toggling the More clustering and output options checkbox.</li>
	</ol>
	</li>
	<li>Residue conservation analysis
	<ol>
		<li>Click the blue &#34;Next&#34; (Conservation) button to calculate the column-wise residue conservation.</li>
		<li>Select the Aligned structure sets to generate a plot of the residue conservation at each alignment position.</li>
		<li>Select Structures aligned with PFAM seed alignment to show conservation calculated with respect to the associated PFAM seed alignment containing representative members of the family.</li>
	</ol>
	</li>
	<li>Sequence alignment display
	<ol>
		<li>Click the blue &#34;Next&#34; (Alignment) button to show the full sequence alignment with the in-browser alignment visualization tool.</li>
	</ol>
	</li>
</ol>
3. Structure Fitting and Analysis (FIT)
<ol>
	<li>Perform structure superimposition by entering the FIT tab.</li>
	<li>Structure superposition
	<ol>
		<li>Toggle the &#34;Show PDBs&#34; checkbox to visualize the aligned protein structures in-browser.</li>
		<li>Make sure the protein structures are superimposed to corresponding and relevant regions by visual inspections. Click and drag the mouse over the structures to rotate, and scroll to zoom.</li>
		<li>Adjust the coloring of the structures by clicking on the &#34;Color options&#34;. Coloring options include alignment position, structural variability per position, RMSD cluster groups, sequence cluster groups, aligned regions and secondary structure.</li>
		<li>Download the superposed structures as either conventional PDB files or as a single PyMOL session file for visualization in a specialized molecular viewer program.</li>
	</ol>
	</li>
	<li>Structure analysis
	<ol>
		<li>Click the blue &#34;Next&#34; (Analysis) button to perform structure-based clustering of the collected PDB structures.</li>
		<li>Toggle the RMSD Heatmap in the Plot options dropdown menu.</li>
		<li>Adjust the clustering options, including the clustering method itself, through toggling the &#34;More clustering and output options&#34; checkbox. 
		NOTE: Pairwise RMSD data can also be visualized as a dendrogram, a histogram or a heat map.</li>
	</ol>
	</li>
	<li>Residue fluctuations
	<ol>
		<li>Click the blue &#34;Next&#34; (RMSF) button to view the structural variability of each residue (shown as an RMSF plot) with major secondary structure elements shown in the marginal regions of the x-axis.</li>
		<li>Toggle the Show B-factors checkbox to overlay crystallographic B-factors of the reference structure onto the RMSF plot.</li>
	</ol>
	</li>
</ol>
4. Principal Component Analysis (PCA)
<ol>
	<li>Perform principal component analysis by entering the &#34;PCA&#34; tab.</li>
	<li>Visualization of the principal components
	<ol>
		<li>Toggle the &#34;Show PC Trajectory&#34; checkbox to visualize motions described by the PCs with the in-browser visualization tool.</li>
		<li>Make sure &#34;Principal Component 1&#34; is chosen from the first drop down menu.</li>
		<li>To visualize the motions described by other PCs, choose the desired PC from the &#34;Choose Principal Component&#34; drop down menu.</li>
		<li>Change the coloring of the trajectory from the &#34;Color options&#34; drop down menu.</li>
		<li>Choose &#34;Variability per position&#34; from the &#34;Color options&#34; to color by displacement magnitude.</li>
		<li>Click the &#34;Download PDB trajectory&#34; button in the &#34;Principal Component Visualization&#34; panel to obtain a trajectory view of the motion described by the PCs.</li>
		<li>Click the button &#34;Download PyMOL&#34; session file to generate a PyMOL session file giving the motions as a vector field.</li>
	</ol>
	</li>
	<li>Conformer analysis
	<ol>
		<li>Project the individual structures onto two selected PCs by clicking on the blue &#34;Next&#34; (Plot) button.</li>
		<li>Make sure &#34;PC on X-axis&#34; is set to 1, and &#34;PC on Y-axis&#34; to 2. To project the structures onto other PCs, adjust the PC numbering accordingly.</li>
		<li>Choose &#34;Cluster by PC Subspace&#34; to color the structures in the plot by PC-based clustering; &#34;RMSD&#34; to color by &#34;RMSD-based&#34; clustering; and &#34;Sequence&#34; to color by sequence based clustering.</li>
		<li>Click on any individual points in the plot to label the structures. Alternatively, highlight one or more structures in the table &#34;PCA conformer plot annotation&#34; below the plot.</li>
		<li>Slide the PCs in subspace slider to including more/less PCs for the clustering algorithm.</li>
	</ol>
	</li>
	<li>Residue contributions
	<ol>
		<li>Calculate the residue contributions to the individual PCs by clicking the blue &#34;Next&#34; (Residue contributions) button.</li>
		<li>Plot the contributions for additional PCs by including the PC number in the &#34;Choose Principal Component&#34; text box.</li>
		<li>Toggle the &#34;Spread lines&#34; checkbox avoid plotting the residue contributions on top of each other.</li>
		<li>Toggle off the &#34;Multiline plot&#34; checkbox to plot the residue contributions in separate plots.</li>
		<li>Toggle the &#34;Show RMSF&#34; to include the RMSF values (from the FIT tab).</li>
	</ol>
	</li>
</ol>
5. Ensemble Normal Mode Analysis (eNMA)
<ol>
	<li>Click the eNMA tab to initiate normal modes (NMs) calculation.</li>
	<li>Filter structure
	<ol>
		<li>Adjust the number of structures by lowering or increasing the &#34;Cutoff&#34; for structure inclusion/exclusion.</li>
		<li>Click the green &#34;Run Ensemble NMA&#34; to start the NMA calculation.</li>
	</ol>
	</li>
	<li>Normal modes visualization
	<ol>
		<li>Scroll down to the second panel of the eNMA tab (Normal Modes Visualization) for visualization of the NMs. 
		&#8203;NOTE: By default, the NM with the highest overlap (similarity) to PC-1 is displayed in the visualization window.</li>
		<li>To visualize the motions described by other NMs or other PDB structures, choose the desired NM and structure from the &#34;Choose Mode&#34; and &#34;Show NMs for structure&#34; drop down menus, respectively.</li>
	</ol>
	</li>
	<li>Residue fluctuations
	<ol>
		<li>Click the blue &#34;Next&#34; (Fluctuations) button to calculate the residue-wise fluctuations of structures selected for eNMA.</li>
		<li>Toggle the &#34;Cluster by RMSD&#34; to color the fluctuation profiles by RMSD-based clustering.</li>
		<li>Toggle the &#34;Cluster by RMSIP&#34; to color the fluctuation profiles by RMSIP-based clustering.</li>
		<li>Toggle the &#34;Spread lines&#34; checkbox to plot the grouped fluctuation profiles apart from each other.</li>
	</ol>
	</li>
	<li>Comparing NMA and PCA
	<ol>
		<li>Click the blue &#34;Next&#34; (PCA-vs-NMA) button to calculate the similarity between the individual NMs and PCs.</li>
		<li>Select a PDB ID from the &#34;Compare NMs of structure&#34; drop down menu to calculate the similarity between the NMs of this structure to the PCs calculated in the PCA tab.</li>
	</ol>
	</li>
	<li>Overlap analysis
	<ol>
		<li>Click the blue &#34;Next&#34; (Overlap analysis) button to calculate the overlap between calculated NMs and the structure difference vector between two selected structures.</li>
		<li>Select a &#39;reference&#39; PDB ID from the &#34;Compare NMs of structure&#34; drop down menu and or one or more PDB IDs in the structure table for the pairwise comparison with the reference PDB.</li>
	</ol>
	</li>
	<li>Clustering analysis
	<ol>
		<li>Click the blue &#34;Next&#34; (Clustering) button to perform structure clustering based on pair-wise NM similarity (RMSIP).</li>
	</ol>
	</li>
</ol>

<ol><li>Kornev, A. P., Taylor, S. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamics-Driven+Allostery+in+Protein+Kinases%5BTitle%5D)+AND+%22Trends+Biochem.+Sci%22%5BJournal%5D)">Dynamics-Driven Allostery in Protein Kinases.</a> Trends Biochem. Sci. 40 (11), 628-647 (2015).</li>
<li>Yao, X. -Q., Grant, B. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Domain-opening+and+dynamic+coupling+in+the+%CE%B1-subunit+of+heterotrimeric+G+proteins%5BTitle%5D)+AND+%22Biophys.+J%22%5BJournal%5D)">Domain-opening and dynamic coupling in the α-subunit of heterotrimeric G proteins.</a> Biophys. J. 105 (2), L08-L10 (2013).</li>
<li>Henzler-Wildman, K. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Intrinsic+motions+along+an+enzymatic+reaction+trajectory%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Intrinsic motions along an enzymatic reaction trajectory.</a> Nature. 450 (7171), 838-844 (2007).</li>
<li>Boehr, D., Nussinov, R., Wright, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+role+of+dynamic+conformational+ensembles+in+biomolecular+recognition%5BTitle%5D)+AND+%22Nat.+Chem.+Biol%22%5BJournal%5D)">The role of dynamic conformational ensembles in biomolecular recognition.</a> Nat. Chem. Biol. 5 (11), 789-796 (2009).</li>
<li>Teilum, K., Olsen, J. G., Kragelund, B. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Functional+aspects+of+protein+flexibility%5BTitle%5D)+AND+%22Cell+Mol+Life+Sci%22%5BJournal%5D)">Functional aspects of protein flexibility.</a> Cell Mol Life Sci. 66 (14), 2231-2247 (2009).</li>
<li>Henzler-Wildman, K., Kern, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+personalities+of+proteins%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Dynamic personalities of proteins.</a> Nature. 450 (7172), 964-972 (2007).</li>
<li>Grant, B. J., Gorfe, A. A., McCammon, J. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Large+conformational+changes+in+proteins%3A+signaling+and+other+functions%5BTitle%5D)+AND+%22Curr.+Opin.+Struct.+Biol%22%5BJournal%5D)">Large conformational changes in proteins: signaling and other functions.</a> Curr. Opin. Struct. Biol. 20 (2), 142-147 (2010).</li>
<li>Grant, B. J., Rodrigues, A. P. C., ElSawy, K. M., McCammon, J. A., Caves, L. S. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Bio3d%3A+an+R+package+for+the+comparative+analysis+of+protein+structures%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Bio3d: an R package for the comparative analysis of protein structures.</a> Bioinformatics. 22 (21), 2695-2696 (2006).</li>
<li>Bakan, A., Meireles, L. M., Bahar, I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ProDy%3A+protein+dynamics+inferred+from+theory+and+experiments%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">ProDy: protein dynamics inferred from theory and experiments.</a> Bioinformatics. 27 (11), 1575-1577 (2011).</li>
<li>Zimmermann, M. T., Kloczkowski, A., Jernigan, R. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MAVENs%3A+motion+analysis+and+visualization+of+elastic+networks+and+structural+ensembles%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">MAVENs: motion analysis and visualization of elastic networks and structural ensembles.</a> BMC Bioinformatics. 12 (1), 264(2011).</li>
<li>Yang, L. -W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((oGNM%3A+online+computation+of+structural+dynamics+using+the+Gaussian+Network+Model%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">oGNM: online computation of structural dynamics using the Gaussian Network Model.</a> Nucleic Acids Res. 34, 24-31 (2006).</li>
<li>Suhre, K., Sanejouand, Y. -H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ElNemo%3A+a+normal+mode+web+server+for+protein+movement+analysis+and+the+generation+of+templates+for+molecular+replacement%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">ElNemo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement.</a> Nucleic Acids Res. 32, W610-W614 (2004).</li>
<li>Tiwari, S. P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((WEBnm%40+v2.0%3A+Web+server+and+services+for+comparing+protein+flexibility%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">WEBnm@ v2.0: Web server and services for comparing protein flexibility.</a> BMC Bioinformatics. 15 (1), 427(2014).</li>
<li>Hrabe, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PDBFlex%3A+exploring+flexibility+in+protein+structures%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">PDBFlex: exploring flexibility in protein structures.</a> Nucleic Acids Res. 44, D423-D428 (2016).</li>
<li>Skjærven, L., Jariwala, S., Yao, X. -Q., Grant, B. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Online+interactive+analysis+of+protein+structure+ensembles+with+Bio3D-web%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Online interactive analysis of protein structure ensembles with Bio3D-web.</a> Bioinformatics. , first published online (2016).</li>
<li>Skjærven, L., Yao, X., Scarabelli, G., Grant, B. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrating+protein+structural+dynamics+and+evolutionary+analysis+with+Bio3D%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">Integrating protein structural dynamics and evolutionary analysis with Bio3D.</a> BMC Bioinformatics. 15 (399), 1-11 (2014).</li>
<li>Eddy, S. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Accelerated+Profile+HMM+Searches%5BTitle%5D)+AND+%22PLoS+Comput.+Biol%22%5BJournal%5D)">Accelerated Profile HMM Searches.</a> PLoS Comput. Biol. 7 (10), (2011).</li>
<li>Edgar, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MUSCLE%3A+multiple+sequence+alignment+with+high+accuracy+and+high+throughput%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">MUSCLE: multiple sequence alignment with high accuracy and high throughput.</a> Nucleic Acids Res. 32 (5), 1792-1797 (2004).</li>
<li>Berman, H. M. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/10592235">The Protein Data Bank.</a> Nucleic Acids Res. 28 (1), 235-242 (2000).</li>
<li>Finn, R. D., et al. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3965110/">Pfam: the protein families database.</a> Nucleic Acids Res. 42, D222-D230 (2014).</li>
<li>Kerns, S. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+energy+landscape+of+adenylate+kinase+during+catalysis%5BTitle%5D)+AND+%22Nat.+Struct.+Mol.+Biol%22%5BJournal%5D)">The energy landscape of adenylate kinase during catalysis.</a> Nat. Struct. Mol. Biol. 22 (2), 124-131 (2015).</li>
</ol>

The authors declare that they have no competing financial interests.

Bio3D-web can be used to interactively explore and map the structural, dynamic and functional states of proteins from available crystallographic structures. Furthermore, the NMA and PCA based clustering results, together with the annotations and sequence based analysis, can be particularly useful for selecting representative structures for more time-consuming analysis such as ensemble small molecule docking or molecular dynamics simulations. Bio3D-web thus facilitates advanced structural bioinformatics analysis for a broader range of researchers by reducing the required level of technical expertise. The current design of Bio3D-web emphasizes simplicity over exhaustive inclusion of the many analysis methods available in the full standalone Bio3D package. In many cases it is envisaged that researchers will use Bio3D-web to understand general trends in their protein family or superfamily of interest, which may then inform more specialized analyses. Bio3D-web is therefore designed to quickly explore biomolecular structure datasets and to act as a hypothesis-generating tool. We encourage users to further explore their data by providing example Bio3D code in the reproducible report that also stores all query details and analysis results.
In the representative example protocol above, we show the capability of Bio3D-web to reveal the structural features of functional conformational transitions of Adk. Additional applications of Bio3D-web include structural and dynamics analysis of user-uploaded PDB structures. For example, the user can upload new structures or indeed protein sequences for analysis. The analysis steps mentioned earlier, especially the eNMA step, can reveal both local and global trends in protein motions, with collective motions being of functional significance. Comparison with apo structures can also reveal characteristics of unbound to bound conformational transitions. Additional examples of application to a range of different protein families are provided online.
Although all proteins are flexible and dynamic entities, not all proteins have atomic resolution structures available in a range of different states (e.g. active and inactive states). Our view of protein structure space is thus a limited one and hence the insight obtained from tools such as Bio3D-web is necessarily also limited for certain proteins. However, with current technological advances and new initiatives for structural genomics the protocol presented here will increasingly become an important route for gaining insight into important structure-function relationships. A critical step, which is particularly important when analyzing more distantly related proteins, is the potential emergence of alignment errors in the ALIGN tab. Alignment errors will inevitably occur when sequence similarity drops below 30% and the user must in such cases double check and correct the sequence alignment in the ALIGN tab. Alignment errors will possibly result in incorrect superimposed structures in the FIT tab and mask the most relevant conformational variations for the subsequent PCA. In addition, the user should be aware of missing residues in the selected PDB structures as in the current implementation PCA can only be performed on protein residues in which all structures have their corresponding carbon alpha atom resolved. Consequently, if a selected PDB has unresolved residues for a particular region of the protein this region will be omitted from PCA.
Bio3D-web is currently limited to the analysis of single chain PDB structures. Consequently, functional motions occurring at the quaternary level cannot be explored using the current protocol. Although we are currently developing new algorithms to include such analysis in Bio3D-web, the only current option is through conventional Bio3D use.
Bio3D-web is the only online application that makes it possible to query and identify structure sets, interpret their patterns of sequence and structural variability, and extract mechanistic information from both analysis and prediction of their structural plasticity. A wide range of molecular visualization tools and online servers enable researchers to explore and analyse individual biomolecular structures. However, existing tools for analysis of the sequence, structure and dynamics of large heterogeneous protein families often require significant computational expertise and typically remain accessible only to users with relevant programming skills. For example, the Bio3D package requires R 8, ProDy requires python and Maven requires Matlab knowledge 9,10. Bio3D-web in contrast requires no programming knowledge and thus increases the accessibility and decreases the entry barrier to performing advanced comparative sequence, structure and dynamics analysis. Furthermore, the preparation, curation, annotation and clean-up of molecular structures that is frequently necessary for efficient analysis is included with the Bio3D-web service. Additionally, the restriction to performing such analysis on capable computational resources is alleviated by our server instance that enables large-scale analysis of many structures that can be initiated and controlled from any modern web browser.
Open development of Bio3D-web is on-going (see https://bitbucket.org/Grantlab/bio3d). We continue to add new analysis functionality and improve existing methods. Future development will focus on the addition of distance matrix based PCA and torsional PCA, more extensive sequence conservation approaches that include a phylogenetic component, ensemble binding site identification, and new approaches for dynamic network analysis across protein families. In this respect the current web application represents the starting point for many other collaborative structural bioinformatic analysis workflows by enabling reproducible and shareable steps on user defined experimental structure sets. We also plan future support of reconstructed biological unit coordinate sets in addition to individual and multiple chains from the asymmetric unit of PDB structures. Additional features will include enhanced saving and loading of collaborative workspaces together with an undo possibility.
Bio3D-web is an online application for interactive analysis of biomolecular structure data. Bio3D-web runs on any modern Web browser and provides functionality for: (1) The identification of related protein structure sets to user specified thresholds of similarity; (2) Their multiple alignment and structure superposition; (3) Sequence and structure conservation analysis; (4) Inter-conformer relationship mapping with principal component analysis, and (5) comparison of predicted internal dynamics via ensemble normal mode analysis. This integrated functionality provides a complete workflow for the investigation of sequence-structure-dynamic relationships within protein families and superfamilies. In addition to a convenient easy to use dynamic interface for exploring the effects of parameter and method choices, Bio3D-web also records the complete user input and subsequent graphical results of a user&#39;s session. This allows users to easily share and reproduce the sequence of analysis steps that created their results.Bio3D-web is implemented entirely in the R language and is based on the Bio3D and Shiny R packages. It can be run from our online server or installed locally on any computer running R. This includes local server installation to provide a customized multi user instance with access to priority structural datasets such as those common in the pharmaceutical industry. Full source code and extensive documentation is provided under a GPL-3 open-source license from: http://thegrantlab.org/bio3d/webapps

The protein data bank (PDB) now contains more than 120,000 protein structures &#8211; many of which are of the same protein family but resolved under different experimental conditions. These multiple structures represent an invaluable resource for understanding the intricacies of protein form and function. For example, the rigorous comparison of these structure ensembles can reveal important molecular mechanisms 1,2,3 and inform on conformational dynamics involved in processes including ligand binding, enzymatic catalysis and bi-molecular recognition 4,5,6,7. New insights can often be obtained from the detailed large-scale analysis of the sequence, structure and dynamics of protein families. However, this typically requires considerable bioinformatics and computer programming expertise together with familiarity with the protein systems under study. For example, software packages such as Bio3D, ProDy and Maven require programing in R, python and Matlab, respectively 8,9,10. Conversely, online tools for analysis of structural flexibility are generally limited to the investigation of individual structures 11,12. An exception in this regard is the recently developed WebNM@ server, which allows for the comparison of flexibility patterns obtained from normal mode analysis (NMA) of several pre-aligned user specified structures13. However, this server lacks an automated procedure for the identification of structures for comparison, their alignment or further analysis beyond NMA. Another recent contribution is the online PDBFlex database, which presents pre-computed analysis of PDB structures sharing 95% or higher sequence identity14. However, analysis of more diverse structure sets is not currently available.
We have previously presented Bio3D-web &#8211; an easy to use web application for the analysis of protein sequence-structure-dynamic relationships15. Bio3D-web is unique in providing easy to use integrated functionality for the identification, comparison and detailed analysis of large homologous structure sets online. Here we present a detailed protocol for the online investigation of protein sequence-structure-dynamics relationship using Bio3D-web. Bio3D-web provides a variety of functions to support the five major steps of data analysis shown in Figure 1 and discussed in detail below. These steps constitute a workflow that spans from query sequence or structure input, through multiple levels of sequence-structure-dynamic analysis, to summary report generation. Results are available immediately through extensive in-browser visualization and plotting devices, as well as through downloading result files in commonly used formats. In addition to a convenient easy to use dynamic interface for exploring the effects of parameter and method choices, Bio3D-web also records the complete user input and subsequent graphical results of a user&#39;s session as a sharable reproducible report in PDF, DOC and HTML formats. User sessions may be saved and reloaded at future times and complete results downloaded and further interpreted by the Bio3D R package on a user&#39;s local machine.
Bio3D-web is powered by the Bio3D R package for analysis of biomolecular structure, sequence and molecular simulation data 8,16. In particular, Bio3D algorithms for rigid-core identification 8, superposition, principal component analysis (PCA) 8, and ensemble normal mode analysis(eNMA) 16 form the basis of the application. We also utilize Bio3D protocols that depend on pHMMER 17 for the identification of related protein structures, and MUSCLE 18 for multiple sequence alignment. Structure and sequence annotations are derived via Bio3D utilities from the RCSB PDB 19 and PFAM databases 20. Bio3D-web can be run from our online server or installed locally on any computer running R. Bio3D-web is open to all users and is provided free of charge under a GPL-3 open-source license from: http://thegrantlab.org/bio3d/webapps

<table><tbody><tr><td>Bio3D-web</td><td></td><td></td><td></td></tr><tr><td>Web-site</td><td>http://thegrantlab.org/bio3d-web/</td><td></td><td></td></tr><tr><td>Requirements</td><td>Web browser</td><td></td><td></td></tr></tbody></table>

investigating protein sequence-structure-dynamics relationships with bio3d-web

We demonstrate the usage of Bio3D-web for the interactive analysis of biomolecular structure data. The Bio3D-web application provides online functionality for: (1) The identification of related protein structure sets to user specified thresholds of similarity; (2) Their multiple alignment and structure superposition; (3) Sequence and structure conservation analysis; (4) Inter-conformer relationship mapping with principal component analysis, and (5) comparison of predicted internal dynamics via ensemble normal mode analysis. This integrated functionality provides a complete online workflow for investigating sequence-structure-dynamic relationships within protein families and superfamilies.

Adenylate kinase (Adk) is a ubiquitous enzyme that functions to maintain the equilibrium between cytoplasmic nucleotides essential for many cellular processes. Adk operates by catalyzing the reversible transfer of a phosphoryl group from ATP to AMP. This reaction is accompanied by well-studied rate limiting conformational transitions 3,21. Here we analyze all currently available Adk structures with Bio3D-web to reveal detailed features and mechanistic principles of these essential transitions.
We can begin our Bio3D-web analysis of Adk by entering the RCSB PDB code of any known Adk structure. For example, entering the PDB ID 1AKE in panel A of the SEARCH tab returns 167 sequence similar structures from which the top 26 are automatically selected for further analysis (see panel B). The annotation presented in panel C indicates that these selected structures are all from E. coli, were solved by x-ray diffraction in a range of space groups; have a resolution range of 1.63 to 2.8 &#197;, and were co-crystalized with a range of different ligands (including no ligands, AMP, ADP, MG and the inhibitor AP5). Note that additional annotation details can be displayed by clicking on &#34;Show/Hide Columns&#34; option in panel C.
Multiple sequence alignment is performed upon entering the ALIGN tab. The first panel of the ALIGN tab displays a summary of the alignment providing details on the number of sequence rows (equivalent to the number of PDB structures), as well as the number of positions (i.e. alignment columns). This includes a specification of the number of gap and non-gap containing columns. The figure on the right hand side of the first row provides a schematic representation of the sequence alignment. Here the grey areas represent non-gap positions, while white areas in the alignment correspond to gaps. A representation of the sequence conservation is shown above the alignment with red areas indicating well-conserved positions, and white indicating less conserved. Note that the sequences in this figure are ordered based on their similarity provided by the clustering dendrogram on the left hand side. The second panel of this tab further facilitates clustering of the selected PDBs based on their pair-wise sequence similarity, which can be visualized either as a dendrogram or a heat map. By default, a dendrogram (or tree diagram) representing the arrangement of clusters is shown. The y- axis of the dendrogram represents the distance (in terms of sequence identity) between the clusters.
Structure superposition is performed automatically upon entering the FIT tab. The superimposed structures, displayed interactively in panel A, indicate the presence of a relatively rigid core region (encompassing residues 1-29, 68-117, and 161-214; see the &#39;optional core and RMSD details&#39; panel at the bottom of the FIT tab for details). Two more variable nucleotide-binding regions (residues 30-67 and 118-167) are also clearly visible (Figure 2). RMSD-based clustering groups these structures into two distinct conformations.
Clicking on the PCA tab more clearly shows the relationship between the structures in terms of the displacements of these regions that effectively close over the bound nucleotide species in related structures (Figure 2B and 2C). The majority of structures are in the &#39;closed&#39; form (blue in Figure 2C) and are associated with a bound ligand or inhibitor. In contrast more &#39;open&#39; conformations are nucleotide and inhibitor free. This is consistent with the extensive body of research on Adk structure and dynamics indicating that an open configuration of these regions is required for nucleotide binding and a closed conformation for efficient phosphoryl transfer and suppression of detrimental hydrolysis events. It is notable that a single PC captures 97% of the total mean square displacement in this Adk structure set and provides a clear and compelling description of the open to closed transition along with the individual residue contributions to this functional displacement (panel C of the app and Figure 2).
Visiting the NMA tab and increasing the number of structures considered for calculation (via decreasing the cutoff for filtering similar structures) indicates that open state structures display enhanced local and global dynamics in comparison to the closed form structures (Figure 2D and panel C of app). Comparing PCA and NMA results for individual structures (panel D) indicates that the first mode of all open form structures displays a relatively high overlap to PC1 (with a mean value of 0.37 &#177; 0.04). In contrast, closed form structures display lower values (with a mean of 0.30 &#177; 0.01). RMSIP values for open form structures (0.62 &#177; 0.003) are also higher than those of closed structures (0.56 &#177; 0.008). In addition, overlap analysis shows that the first modes of the open state are in good agreement with the conformational change that describes the difference of the open and closed states (panel E). Clustering based on RMSIP values again displays a consistent partitioning of open and closed state structures (panel F).
Collectively these results indicate the existence of two major distinct conformational states for Adk. These differ by a collective low frequency displacement of two nucleotide-binding site regions that display distinct flexibilities upon nucleotide binding.
<img alt="Figure 1" src="/files/ftp_upload/55640/55640fig1.jpg" /> 
Figure 1: Bio3D-web overview with screen shots of the PCA and NMA tabs. Bio3D-web takes a user provided protein structure or sequence as input in the SEARCH tab (1). The server provides a list of related structures, which can be selected for further analysis. (2) The ALIGN tab provides sequence alignment and analysis of the structures selected in the SEARCH tab. (3) In the FIT tab all structures are superimposed and visualized in 3D together with the results of conventional pair-wise structure analysis. (4) Principal component analysis of the structure set is performed in the PCA tab to characterize inter-conformer relationships. (5) Normal mode analysis on each structure can be carried out in the eNMA tab to explore dynamic trends for the available structural states. <a href="http://cloudflare.jove.com/files/ftp_upload/55640/55640fig1large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" src="/files/ftp_upload/55640/55640fig2.jpg" /> 
Figure 2: Results of Bio3D-web analysis of adenylate kinase. (A) Available PDB structures of adenylate kinase superimposed on the identified invariant core. Structures are colored according to RMSD-based clustering provided in the FIT tab. (B) Visualization of the principal components is available from the PCA tab to characterize the major conformational variations in the data set. Here, the trajectory corresponding to the first principal component is shown in tube representation showing the large-scale closing motion of the protein. (C) Structures are projected onto their two first principal components in a conformer plot showing a low-dimensional representation of the conformational variability. Each dot (or structure) is colored according to user specified criteria, in this case PCA-based clustering results. (D) Normal mode analysis in the eNMA tab suggests enhanced local and global dynamics for structures in the open state (red) in comparison to the closed form (blue) structures. <a href="http://cloudflare.jove.com/files/ftp_upload/55640/55640fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>

Watch this Scientific Journal Video about Investigating Protein Sequence-structure-dynamics Relationships with Bio3D-web at JoVE.com

Investigating Protein Sequence-structure-dynamics Relationships with Bio3D-web

A protocol for the online investigation of protein sequence-structure-dynamics relationships using Bio3D-web is presented.

We demonstrate the usage of Bio3D-web for the interactive analysis of biomolecular structure data. The Bio3D-web application provides online functionality ...

investigating-protein-sequence-structure-dynamics-relationships-with

Research

JoVE Journal

Biochemistry

15.4K Views.  University of Michigan Medical School. A protocol for the online investigation of protein sequence-structure-dynamics relationships using Bio3D-web is presented. 

Video: Investigating Protein Sequence-structure-dynamics Relationships with Bio3D-web

Assessment of Immunologically Relevant Dynamic Tertiary Structural Features of the HIV-1 V3 Loop Crown R2 Sequence by ab initio Folding

The antigenic diversity of HIV-1 has long been an obstacle to vaccine design, and this variability is especially pronounced in the V3 loop of the virus' surface envelope glycoprotein. We previously proposed that the crown of the V3 loop, although dynamic and sequence variable, is constrained throughout the population of HIV-1 viruses to an immunologically relevant &beta;-hairpin tertiary structure. Importantly, there are thousands of different V3 loop crown sequences in circulating HIV-1 viruses, making 3D structural characterization of trends across the diversity of viruses difficult or impossible by crystallography or NMR. Our previous successful studies with folding of the V3 crown1, 2 used the ab initio algorithm 3 accessible in the ICM-Pro molecular modeling software package (Molsoft LLC, La Jolla, CA) and suggested that the crown of the V3 loop, specifically from positions 10 to 22, benefits sufficiently from the flexibility and length of its flanking stems to behave to a large degree as if it were an unconstrained peptide freely folding in solution. As such, rapid ab initio folding of just this portion of the V3 loop of any individual strain of the 60,000+ circulating HIV-1 strains can be informative. Here, we folded the V3 loop of the R2 strain to gain insight into the structural basis of its unique properties. R2 bears a rare V3 loop sequence thought to be responsible for the exquisite sensitivity of this strain to neutralization by patient sera and monoclonal antibodies4, 5. The strain mediates CD4-independent infection and appears to elicit broadly neutralizing antibodies. We demonstrate how evaluation of the results of the folding can be informative for associating observed structures in the folding with the immunological activities observed for R2.

Assessment of Immunologically Relevant Dynamic Tertiary Structural Features of the HIV-1 V3 Loop Crown R2 Sequence by ab initio Folding

The crown region of different V3 loop sequences of the surface envelope glycoprotein (gp120) of HIV-1 can be structurally characterized in many cases by in silico folding of positions 10 to 22 of the loop using a state-of-the-art ab initio folding algorithm. Here we demonstrate the folding and evaluation of this region of the V3 loop from the R2 strain of HIV-1, a uniquely neutralization sensitive strain with puzzling functional properties.

The antigenic diversity of HIV-1 has long been an obstacle to vaccine design, and this variability is especially pronounced in the V3 loop of the virus' ...

Immunology and Infection

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.

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 ...

Chemistry

Analyzing Dynamic Protein Complexes Assembled On and Released From Biolayer Interferometry Biosensor Using Mass Spectrometry and Electron Microscopy

In vivo, proteins are often part of large macromolecular complexes where binding specificity and dynamics ultimately dictate functional outputs. In this work, the pre-endosomal anthrax toxin is assembled and transitioned into the endosomal complex. First, the N-terminal domain of a cysteine mutant lethal factor (LFN) is attached to a biolayer interferometry (BLI) biosensor through disulfide coupling in an optimal orientation, allowing protective antigen (PA) prepore to bind (Kd 1 nM). The optimally oriented LFN-PAprepore complex then binds to soluble capillary morphogenic gene-2 (CMG2) cell surface receptor (Kd 170 pM), resulting in a representative anthrax pre-endosomal complex, stable at pH 7.5. This assembled complex is then subjected to acidification (pH 5.0) representative of the late endosome environment to transition the PAprepore into the membrane inserted pore state. This PApore state results in a weakened binding between the CMG2 receptor and the LFN-PApore and a substantial dissociation of CMG2 from the transition pore. The thio-attachment of LFN to the biosensor surface is easily reversed by dithiothreitol. Reduction on the BLI biosensor surface releases the LFN-PAprepore-CMG2 ternary complex or the acid transitioned LFN-PApore complexes into microliter volumes. Released complexes are then visualized and identified using electron microscopy and mass spectrometry. These experiments demonstrate how to monitor the kinetic assembly/disassembly of specific protein complexes using label-free BLI methodologies and evaluate the structure and identity of these BLI assembled complexes by electron microscopy and mass spectrometry, respectively, using easy-to-replicate sequential procedures.

Here we present a protocol to monitor the assembly and disassembly of the anthrax toxin using biolayer interferometry (BLI). Following assembly/disassembly on the biosensor surface, the large protein complexes are released from the surface for visualization and identification of components of the complexes using electron microscopy and mass spectrometry, respectively.

In vivo, proteins are often part of large macromolecular complexes where binding specificity and dynamics ultimately dictate functional outputs. In this ...

Analyzing Protein Architectures and Protein-Ligand Complexes by Integrative Structural Mass Spectrometry

Proteins are an important class of biological macromolecules that play many key roles in cellular functions including gene expression, catalyzing metabolic reactions, DNA repair and replication. Therefore, a detailed understanding of these processes provides critical information on how cells function. Integrative structural MS methods offer structural and dynamical information on protein complex assembly, complex connectivity, subunit stoichiometry, protein oligomerization and ligand binding. Recent advances in integrative structural MS have allowed for the characterization of challenging biological systems including large DNA binding proteins and membrane proteins. This protocol describes how to integrate diverse MS data such as native MS and ion mobility-mass spectrometry (IM-MS) with molecular dynamics simulations to gain insights into a helicase-nuclease DNA repair protein complex. The resulting approach provides a framework for detailed studies of ligand binding to other protein complexes involved in important biological processes.

Mass spectrometry (MS) has emerged as an important tool for the investigation of structure and dynamics of macromolecular assemblies. Here, we integrate MS-based approaches to interrogate protein complex formation and ligand binding.

Proteins are an important class of biological macromolecules that play many key roles in cellular functions including gene expression, catalyzing ...

Modeling an Enzyme Active Site using Molecular Visualization Freeware

Biomolecular visualization skills are paramount to understanding key concepts in the biological sciences, such as structure-function relationships and molecular interactions. Various programs allow a learner to manipulate 3D structures, and biomolecular modeling promotes active learning, builds computational skills, and bridges the gap between two dimensional textbook images and the three dimensions of life. A critical skill in this area is to model a protein active site, displaying parts of the macromolecule that can interact with a small molecule, or ligand, in a way that shows binding interactions. In this protocol, we describe this process using four freely available macromolecular modeling programs: iCn3D, Jmol/JSmol, PyMOL, and UCSF ChimeraX. This guide is intended for students seeking to learn the basics of a specific program, as well as instructors incorporating biomolecular modeling into their curriculum. The protocol enables the user to model an active site using a specific visualization program, or to sample several of the free programs available. The model chosen for this protocol is human glucokinase, an isoform of the enzyme hexokinase, which catalyzes the first step of glycolysis. The enzyme is bound to one of its substrates, as well as a non-reactive substrate analog, which allows the user to analyze interactions in the catalytic complex.

A key skill in biomolecular modeling is displaying and annotating active sites in proteins. This technique is demonstrated using four popular free programs for macromolecular visualization: iCn3D, Jmol, PyMOL, and UCSF ChimeraX.

Biomolecular visualization skills are paramount to understanding key concepts in the biological sciences, such as structure-function relationships and ...

Predicting Amino Acid Preferences of Protein-Protein Binding Interfaces

Computational Prediction of Amino Acid Preferences of Potentially Multispecific Peptide-Binding Domains Involved in Protein-Protein Interactions

Many protein-protein interactions involve the binding of short protein segments to peptide-binding domains. Usually, such interactions require the recognition of linear motifs with variable conservation. The combination of highly conserved and more variable regions in the same ligands often contributes to the multispecificity of binding, a common property of enzymes and cell signaling proteins. Characterization of amino acid preferences of peptide-binding domains is important for the design of mediators of protein-protein interactions (PPIs). Computational methods are an efficient alternative to the often costly and cumbersome experimental techniques, enabling the design of potential mediators that can be later validated in downstream experiments. Here, we described a methodology using the Pepspec application of the Rosetta molecular modeling package to predict the amino acid preferences of peptide-binding domains. This methodology is useful when the structure of the receptor protein and the nature of the peptide ligand are both known or can be inferred. The methodology starts with a well-characterized anchor from the ligand, which is extended by randomly adding amino acid residues. The binding affinity of peptides generated this way is then evaluated by flexible-backbone peptide docking in order to select the peptides with the best predicted binding scores. These peptides are then used to calculate amino acid preferences and to optionally compute a position-weight matrix (PWM) that can be used in further studies. To illustrate the application of this methodology, we used the interaction between subunits of human interferon regulatory factor 5 (IRF5), previously known to be multispecific but globally guided by a short conserved motif called pLxIS. The estimated amino acid preferences were consistent with previous knowledge about the IRF5 binding surface. Positions occupied by phosphorylatable serine residues exhibited a high frequency of aspartate and glutamate, likely because their negatively charged side chains are similar to phosphoserine.

Author Spotlight: A Computational Approach to Decipher Amino Acid Preferences in Multispecific Protein-Protein Interactions

We describe a methodology based on sequence diversification to estimate the amino acid preferences of multispecific binding sites in protein-protein interactions (PPIs). In this strategy, thousands of potential peptide ligands are generated and screened in silico, thus overcoming some limitations of available experimental methods.

Many protein-protein interactions involve the binding of short protein segments to peptide-binding domains. Usually, such interactions require the ...

A Protocol for Computer-Based Protein Structure and Function Prediction

Genome sequencing projects have ciphered millions of protein sequence, which require knowledge of their structure and function to improve the understanding of their biological role. Although experimental methods can provide detailed information for a small fraction of these proteins, computational modeling is needed for the majority of protein molecules which are experimentally uncharacterized. The I-TASSER server is an on-line workbench for high-resolution modeling of protein structure and function. Given a protein sequence, a typical output from the I-TASSER server includes secondary structure prediction, predicted solvent accessibility of each residue, homologous template proteins detected by threading and structure alignments, up to five full-length tertiary structural models, and structure-based functional annotations for enzyme classification, Gene Ontology terms and protein-ligand binding sites. All the predictions are tagged with a confidence score which tells how accurate the predictions are without knowing the experimental data. To facilitate the special requests of end users, the server provides channels to accept user-specified inter-residue distance and contact maps to interactively change the I-TASSER modeling; it also allows users to specify any proteins as template, or to exclude any template proteins during the structure assembly simulations. The structural information could be collected by the users based on experimental evidences or biological insights with the purpose of improving the quality of I-TASSER predictions. The server was evaluated as the best programs for protein structure and function predictions in the recent community-wide CASP experiments. There are currently &gt;20,000 registered scientists from over 100 countries who are using the on-line I-TASSER server.

Guidelines for computer based structural and functional characterization of protein using the I-TASSER pipeline is described. Starting from query protein sequence, 3D models are generated using multiple threading alignments and iterative structural assembly simulations. Functional inferences are thereafter drawn based on matches to proteins with known structure and functions.

Genome sequencing projects have ciphered millions of protein sequence, which require knowledge of their structure and function to improve the ...

Biology

Analyzing and Building Nucleic Acid Structures with 3DNA

The 3DNA software package is a popular and versatile bioinformatics tool with capabilities to analyze, construct, and visualize three-dimensional nucleic acid structures. This article presents detailed protocols for a subset of new and popular features available in 3DNA, applicable to both individual structures and ensembles of related structures. Protocol 1 lists the set of instructions needed to download and install the software. This is followed, in Protocol 2, by the analysis of a nucleic acid structure, including the assignment of base pairs and the determination of rigid-body parameters that describe the structure and, in Protocol 3, by a description of the reconstruction of an atomic model of a structure from its rigid-body parameters. The most recent version of 3DNA, version 2.1, has new features for the analysis and manipulation of ensembles of structures, such as those deduced from nuclear magnetic resonance (NMR) measurements and molecular dynamic (MD) simulations; these features are presented in Protocols 4 and 5. In addition to the 3DNA stand-alone software package, the w3DNA web server, located at <a href="http://w3dna.rutgers.edu" target="_blank">http://w3dna.rutgers.edu</a>, provides a user-friendly interface to selected features of the software. Protocol 6 demonstrates a novel feature of the site for building models of long DNA molecules decorated with bound proteins at user-specified locations.

The 3DNA software package is a popular and versatile bioinformatics tool with capabilities to analyze, construct, and visualize three-dimensional nucleic acid structures. This article presents detailed protocols for a subset of new and popular features available in 3DNA, applicable to both individual structures and ensembles of related structures.

The 3DNA software package is a popular and versatile bioinformatics tool with capabilities to analyze, construct, and visualize three-dimensional nucleic ...

Protein WISDOM: A Workbench for In silico De novo Design of BioMolecules

The aim of de novo protein design is to find the amino acid sequences that will fold into a desired 3-dimensional structure with improvements in specific properties, such as binding affinity, agonist or antagonist behavior, or stability, relative to the native sequence. Protein design lies at the center of current advances drug design and discovery. Not only does protein design provide predictions for potentially useful drug targets, but it also enhances our understanding of the protein folding process and protein-protein interactions. Experimental methods such as directed evolution have shown success in protein design. However, such methods are restricted by the limited sequence space that can be searched tractably. In contrast, computational design strategies allow for the screening of a much larger set of sequences covering a wide variety of properties and functionality. We have developed a range of computational de novo protein design methods capable of tackling several important areas of protein design. These include the design of monomeric proteins for increased stability and complexes for increased binding affinity.
To disseminate these methods for broader use we present Protein WISDOM (<a href="http://www.proteinwisdom.org" target="_blank">http://www.proteinwisdom.org</a>), a tool that provides automated methods for a variety of protein design problems. Structural templates are submitted to initialize the design process. The first stage of design is an optimization sequence selection stage that aims at improving stability through minimization of potential energy in the sequence space. Selected sequences are then run through a fold specificity stage and a binding affinity stage. A rank-ordered list of the sequences for each step of the process, along with relevant designed structures, provides the user with a comprehensive quantitative assessment of the design. Here we provide the details of each design method, as well as several notable experimental successes attained through the use of the methods.

Protein WISDOM: A Workbench for In silico De novo Design of BioMolecules

We developed computational de novo protein design methods capable of tackling several important areas of protein design. To disseminate these methods we present Protein WISDOM, an online tool for protein design (<a href="http://www.proteinwisdom.org" target="_blank">http://www.proteinwisdom.org</a>). Starting from a structural template, design of monomeric proteins for increased stability and complexes for increased binding affinity can be performed.

The aim of de novo protein design is to find the amino acid sequences that will fold into a desired 3-dimensional structure with improvements in specific ...

An Integrated Approach for Microprotein Identification and Sequence Analysis

Next-generation sequencing (NGS) has propelled the field of genomics forward and produced whole genome sequences for numerous animal species and model organisms. However, despite this wealth of sequence information, comprehensive gene annotation efforts have proven challenging, especially for small proteins. Notably, conventional protein annotation methods were designed to intentionally exclude putative proteins encoded by short open reading frames (sORFs) less than 300 nucleotides in length to filter out the exponentially higher number of spurious noncoding sORFs throughout the genome. As a result, hundreds of functional small proteins called microproteins (&#60;100 amino acids in length) have been incorrectly classified as noncoding RNAs or overlooked entirely.
Here we provide a detailed protocol to leverage free, publicly available bioinformatic tools to query genomic regions for microprotein-coding potential based on evolutionary conservation. Specifically, we provide step-by-step instructions on how to examine sequence conservation and coding potential using Phylogenetic Codon Substitution Frequencies (PhyloCSF) on the user-friendly University of California Santa Cruz (UCSC) Genome Browser. Additionally, we detail steps to efficiently generate multiple species alignments of identified microprotein sequences to visualize amino acid sequence conservation and recommend resources to analyze microprotein characteristics, including predicted domain structures. These powerful tools can be used to help identify putative microprotein-coding sequences in noncanonical genomic regions or to rule out the presence of a conserved coding sequence with translational potential in a noncoding transcript of interest.

The protocol described here provides detailed instructions on how to analyze genomic regions of interest for microprotein-coding potential using PhyloCSF on the user-friendly UCSC Genome Browser. Additionally, several tools and resources are recommended to further investigate sequence characteristics of identified microproteins to gain insight into their putative functions.

Next-generation sequencing (NGS) has propelled the field of genomics forward and produced whole genome sequences for numerous animal species and model ...

Investigating Protein Sequence-structure-dynamics Relationships with Bio3D-web

Summary

Explore More Videos

PLAYLIST