Iron inhibits Escherichia coli topoisomerase I activity by targeting the first two zinc‐binding sites in the C‐terminal domain

Escherichia coli DNA topoisomerase I (TopA) contains a 67 kDa N‐terminal catalytic domain and a 30 kDa C‐terminal zinc‐binding region (ZD domain) which has three adjacent tetra‐cysteine zinc‐binding motifs. Previous studies have shown that E. coli TopA can bind both iron and zinc, and that iron binding in TopA results in failure to unwind the negatively supercoiled DNA. Here, we report that each E. coli TopA monomer binds one atom of iron via the first two zinc‐binding motifs in ZD domain and both the first and second zinc‐binding motifs are required for iron binding in TopA. The site‐directed mutagenesis studies further reveal that while the mutation of the third zinc‐binding motif has very little effect on TopA's activity, mutation of the first two zinc‐binding motifs in TopA greatly diminishes the topoisomerase activity in vitro and in vivo, indicating that the first two zinc‐binding motifs in TopA are crucial for its function. The DNA‐binding activity assay and intrinsic tryptophan fluorescence measurements show that iron binding in TopA may decrease the single‐stranded (ss) DNA‐binding activity of ZD domain and also change the protein structure of TopA, which subsequently modulate topoisomerase activity.     

A distance‐dependent atomic knowledge‐based potential and force for discrimination of native structures from decoys

The purpose of this article is to introduce a novel model for discriminating correctly folded proteins from well designed decoy structures using mechanical interatomic forces. In our model, we consider a protein as a collection of springs and the force imposed to each atom is calculated. A potential function is obtained from statistical contact preferences within known protein structures. Combining this function with the spring equation, the interatomic forces are calculated. Finally, we consider a structure and define a score function on the 3D structure of a protein. We compare the force imposed to each atom of a protein with the corresponding atom in the other structures. We then assign larger scores to those atoms with lower forces. The total score is the sum of partial scores of atoms. The optimal structure is assumed to be the one with the highest score in the data set. To evaluate the performance of our model, we apply it on several decoy sets. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

All‐atom and coarse‐grained simulations of the forced unfolding pathways of the SNARE complex

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP‐25), is thought to drive membrane fusion by assembling into a four‐helix bundle through a zippering process. In support of the above zippering model, a recent single‐molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C‐terminal domain → the N‐terminal domain. To offer detailed structural insights to this unzipping process, we have performed all‐atom and coarse‐grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all‐atom force field (with an implicit solvent model) or a coarse‐grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse‐grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair‐wise residue–residue distance fluctuations collected from all‐atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino‐acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins. Proteins 2014; 82:1376–1386. © 2014 Wiley Periodicals, Inc.     

Analysis of proton wires in the enzyme active site suggests a mechanism of c‐di‐GMP hydrolysis by the EAL domain phosphodiesterases

We report for the first time a hydrolysis mechanism of the cyclic dimeric guanosine monophosphate (c‐di‐GMP) by the EAL domain phosphodiesterases as revealed by molecular simulations. A model system for the enzyme‐substrate complex was prepared on the base of the crystal structure of the EAL domain from the BlrP1 protein complexed with c‐di‐GMP. The nucleophilic hydroxide generated from the bridging water molecule appeared in a favorable position for attack on the phosphorus atom of c‐di‐GMP. The most difficult task was to find a pathway for a proton transfer to the O3' atom of c‐di‐GMP to promote the O3'? P bond cleavage. We show that the hydrogen bond network extended over the chain of water molecules in the enzyme active site and the Glu359 and Asp303 side chains provides the relevant proton wires. The suggested mechanism is consistent with the structural, mutagenesis, and kinetic experimental studies on the EAL domain phosphodiesterases. Proteins 2016; 84:1670–1680. © 2016 Wiley Periodicals, Inc.     

Single‐molecule pair studies of the interactions of the α‐GalNAc (Tn‐antigen) form of porcine submaxillary mucin with soybean agglutinin

Mucins form a group of heavily O‐glycosylated biologically important glycoproteins that are involved in a variety of biological functions, including modulating immune response, inflammation, and adhesion. Mucins are also involved in cancer and metastasis and often express diagnostic cancer antigens. Recently, a modified porcine submaxillary mucin (Tn‐PSM) containing GalNAcα1‐O‐Ser/Thr residues was shown to bind to soybean agglutinin (SBA) with ～10⁶‐fold enhanced affinity relative to GalNAcα1‐O‐Ser, the pancarcinoma carbohydrate antigen. In this study, dynamic force spectroscopy is used to investigate molecular pairs of SBA and Tn‐PSM. A number of force jumps that demonstrate unbinding or rebinding events were observed up to a distance equal to 2.0 μm, consistent with the length of the mucin chain. The unbinding force increased from 103 to 402 pN with increasing force loading rate. The position of the activation barrier in the energy landscape of the interaction was 0.1 nm. The lifetime of the SBA–TnPSM complex in the absence of applied force was determined to be in the range 1.3–1.9 s. Kinetic parameters describing the rate of dissociation of other sugar lectin interactions are in the range 3.3 × 10^?3–2.5 × 10^?3 s. The long lifetime of the SBA‐TnPSM complex is compatible with a binding model in which lectin molecules “bind and jump” from α‐GalNAc residue to α‐GalNAc residue along the polypeptide chain of Tn‐PSM before dissociating. These findings have important implications for the molecular recognition properties of mucins. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 719–728, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

REMO: A new protocol to refine full atomic protein models from C‐alpha traces by optimizing hydrogen‐bonding networks

Protein structure prediction approaches usually perform modeling simulations based on reduced representation of protein structures. For biological utilizations, it is an important step to construct full atomic models from the reduced structure decoys. Most of the current full atomic model reconstruction procedures have defects which either could not completely remove the steric clashes among backbone atoms or generate final atomic models with worse topology similarity relative to the native structures than the reduced models. In this work, we develop a new protocol, called REMO, to generate full atomic protein models by optimizing the hydrogen‐bonding network with basic fragments matched from a newly constructed backbone isomer library of solved protein structures. The algorithm is benchmarked on 230 nonhomologous proteins with reduced structure decoys generated by I‐TASSER simulations. The results show that REMO has a significant ability to remove steric clashes, and meanwhile retains good topology of the reduced model. The hydrogen‐bonding network of the final models is dramatically improved during the procedure. The REMO algorithm has been exploited in the recent CASP8 experiment which demonstrated significant improvements of the I‐TASSER models in both atomic‐level structural refinement and hydrogen‐bonding network construction. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Folding of Top7 in unbiased all‐atom Monte Carlo simulations

For computational studies of protein folding, proteins with both helical and β‐sheet secondary structure elements are very challenging, as they expose subtle biases of the physical models. Here, we present reproducible folding of a 92 residue α/β protein (residues 3–94 of Top7, PDB ID: 1QYS) in computer simulations starting from random initial conformations using a transferable physical model which has been previously shown to describe the folding and thermodynamic properties of about 20 other smaller proteins of different folds. Top7 is a de novo designed protein with two α‐helices and a five stranded β‐sheet. Experimentally, it is known to be unusually stable for its size, and its folding transition distinctly deviates from the two‐state behavior commonly seen in natural single domain proteins. In our all‐atom implicit solvent parallel tempering Monte Carlo simulations, Top7 shows a rapid transition to a group of states with high native‐like secondary structure, and a much slower subsequent transition to the native state with a root mean square deviation of about 3.5 Å from the experimentally determined structure. Consistent with experiments, we find Top7 to be thermally extremely stable, although the simulations also find a large number of very stable non‐native states with high native‐like secondary structure. Proteins 2013; 81:1446–1456. © 2013 Wiley Periodicals, Inc.     

Molecular modeling of the binding modes of the iron‐sulfur protein to the Jac1 co‐chaperone from Saccharomyces cerevisiae by all‐atom and coarse‐grained approaches

The iron‐sulfur protein 1 (Isu1) and the J‐type co‐chaperone Jac1 from yeast are part of a huge ATP‐dependent system, and both interact with Hsp70 chaperones. Interaction of Isu1 and Jac1 is a part of the iron‐sulfur cluster biogenesis system in mitochondria. In this study, the structure and dynamics of the yeast Isu1–Jac1 complex has been modeled. First, the complete structure of Isu1 was obtained by homology modeling using the I‐TASSER server and YASARA software and thereafter tested for stability in the all‐atom force field AMBER. Then, the known experimental structure of Jac1 was adopted to obtain initial models of the Isu1–Jac1 complex by using the ZDOCK server for global and local docking and the AutoDock software for local docking. Three most probable models were subsequently subjected to the coarse‐grained molecular dynamics simulations with the UNRES force field to obtain the final structures of the complex. In the most probable model, Isu1 binds to the left face of the Γ‐shaped Jac1 molecule by the β‐sheet section of Isu1. Residues L₁₀₅, L₁₀₉, and Y₁₆₃ of Jac1 have been assessed by mutation studies to be essential for binding (Ciesielski et al., J Mol Biol 2012; 417:1–12). These residues were also found, by UNRES/molecular dynamics simulations, to be involved in strong interactions between Isu1 and Jac1 in the complex. Moreover, N₉₅, T₉₈, P₁₀₂, H₁₁₂, V₁₅₉, L₁₆₇, and A₁₇₀ of Jac1, not yet tested experimentally, were also found to be important in binding. Proteins 2015; 83:1414–1426. © 2015 Wiley Periodicals, Inc.     

Directed self‐immobilization of alkaline phosphatase on micro‐patterned substrates via genetically fused metal‐binding peptide

Current biotechnological applications such as biosensors, protein arrays, and microchips require oriented immobilization of enzymes. The characteristics of recognition, self‐assembly and ease of genetic manipulation make inorganic binding peptides an ideal molecular tool for site‐specific enzyme immobilization. Herein, we demonstrate the utilization of gold binding peptide (GBP1) as a molecular linker genetically fused to alkaline phosphatase (AP) and immobilized on gold substrate. Multiple tandem repeats (n = 5, 6, 7, 9) of gold binding peptide were fused to N‐terminus of AP (nGBP1‐AP) and the enzymes were expressed in E. coli cells. The binding and enzymatic activities of the bi‐functional fusion constructs were analyzed using quartz crystal microbalance spectroscopy and biochemical assays. Among the multiple‐repeat constructs, 5GBP1‐AP displayed the best bi‐functional activity and, therefore, was chosen for self‐immobilization studies. Adsorption and assembly properties of the fusion enzyme, 5GBP1‐AP, were studied via surface plasmon resonance spectroscopy and atomic force microscopy. We demonstrated self‐immobilization of the bi‐functional enzyme on micro‐patterned substrates where genetically linked 5GBP1‐AP displayed higher enzymatic activity per area compared to that of AP. Our results demonstrate the promising use of inorganic binding peptides as site‐specific molecular linkers for oriented enzyme immobilization with retained activity. Directed assembly of proteins on solids using genetically fused specific inorganic‐binding peptides has a potential utility in a wide range of biosensing and bioconversion processes. Biotechnol. Bioeng. 2009;103: 696–705. © 2009 Wiley Periodicals, Inc.     

Self‐guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self‐guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrC^r). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson–Boltzmann calculations of the pK_a values of the active site residues suggest an increase in the pK_a of His‐84 on phosphorylation of Asp‐54. In SGLD simulations of the phosphorylated active form with charged His‐84, the average position of the regulatory helix α4 is found closer to the starting structure than in simulations with the neutral His‐84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix α4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix α4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Human topoisomerase I poisoning by protoberberines: potential roles for both drug-DNA and drug-enzyme interactions

Protoberberines represent a structural class of organic cations that induce topoisomerase I-mediated DNA cleavage, a behavior termed topoisomerase I poisoning. We have employed a broad range of biophysical, biochemical, and computer modeling techniques to characterize and cross-correlate the DNA-binding and topoisomerase poisoning properties of four protoberberine analogues that differ with respect to the substituents on their A- and/or D-rings. Our data reveal the following significant features: (i) The binding of the four protoberberines unwinds duplex DNA by approximately 11 degrees, an observation consistent with an intercalative mode of interaction. (ii) Enthalpically favorable interactions, such as stacking interactions between the intercalated ligand and the neighboring base pairs, provide <50% of the thermodynamic driving force for the complexation of the protoberberines to duplex DNA. Computer modeling studies on protoberberine-DNA complexes suggest that only rings C and D intercalate into the host DNA helix, while rings A and B protrude out of the helix interior into the minor groove. (iii) All four protoberberine analogues are topoisomerase I-specific poisons, exhibiting little or no topoisomerase II poisoning activity. (iv) Modifications of the D-ring influence both DNA binding and topoisomerase I poisoning properties. Specifically, transference of a methoxy substituent from the 11- to the 9-position diminishes both DNA binding affinity and topoisomerase I poisoning activity, an observation suggesting that DNA binding is important in the poisoning of topoisomerase I by protoberberines. (v) Modifications of the A-ring have a negligible impact on DNA binding affinity, while exerting a profound influence on topoisomerase I poisoning activity. Specifically, protoberberine analogues containing either 2,3-dimethoxy; 3,4-dimethoxy; or 3, 4-methylenedioxy substituents all bind DNA with a similar affinity. By contrast, these analogues exhibit markedly different topoisomerase I poisoning activities, with these activities following the hierarchy: 3,4-methylenedioxy > 2,3-dimethoxy > 3, 4-dimethoxy. These differences in topoisomerase I poisoning activity may reflect the differing abilities of the analogues to interact with specific functionalities on the enzyme, thereby stabilizing the enzyme in its cleavable state. In the aggregate, our results are consistent with a mechanistic model in which both ligand-DNA and ligand-enzyme interactions are important for the poisoning of topoisomerase I by protoberberines, with the DNA-directed interactions involving ring D and the enzyme-directed interactions involving ring A. It is reasonable to suggest that the poisoning of topoisomerase I by a broad range of other naturally occurring and synthetic ligands may entail a similar mechanism.     

Structure,dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X‐ray studies

Molecular dynamics simulations have been carried out on all the jacalin–carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X‐ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin–carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate‐binding residues are consistent with the known thermodynamic parameters of jacalin–carbohydrate interactions. The simulations, along with X‐ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin–carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structure prediction of helical transmembrane proteins at two length scales

As the first step toward a multi-scale, hierarchical computational approach for membrane protein structure prediction, the packing of transmembrane helices was modeled at the residue and atom levels, respectively. For predictions at the residue level, the helix-helix and helix-membrane interactions were described by a set of knowledge-based energy functions. For predictions at the atom level, CHARMM19 force field was used. To facilitate the system to overcome energy barriers, the Wang-Landau method was employed, where a random walk is performed in the energy space with a uniform probability. Native-like structures were predicted at both levels for two model systems, each of which consists of two transmembrane helices. Interestingly, consistent results were obtained from simulations at the residue and atom levels for the same system, strongly suggesting the feasibility of a hierarchical approach for membrane protein structure predictions.     

Docking of 4‐oxalocrotonate tautomerase substrates: Implications for the catalytic mechanism

The enzyme 4‐oxalocrotonate tautomerase catalyzes the ketonization of dienols, which after further processing become intermediates in the Krebs cycle. The enzyme uses a general acid–base mechanism for proton transfer: the amino‐terminal proline has been shown to function as the catalytic base and Arg39 has been implicated as the catalytic acid. We report the results of molecular docking simulations of 4‐oxalocrotonate tautomerase with two substrates, 2‐hydroxymuconate and 5‐carboxymethyl‐2‐hydroxymuconate. pK_a calculations are also performed for the free enzyme. The predicted binding mode of 2‐hydroxymuconate is in agreement with experimental data. A model for the binding mode of 5‐carboxymethyl‐2‐hydroxymuconate is proposed which explains the lower catalytic efficiency of the enzyme toward this substrate. The pK_a predictions and docking simulations support residue Arg39 as the general acid for the enzyme catalysis. © 1999 John Wiley & Sons, Inc. Biopoly 50: 319–328, 1999     

Coarse‐grained simulations of proton‐dependent conformational changes in lactose permease

During lactose/H⁺ symport, the Escherichia coli lactose permease (LacY) undergoes a series of global conformational transitions between inward‐facing (open to cytoplasmic side) and outward‐facing (open to periplasmic side) states. However, the exact local interactions and molecular mechanisms dictating those large‐scale structural changes are not well understood. All‐atom molecular dynamics simulations have been performed to investigate the molecular interactions involved in conformational transitions of LacY, but the simulations can only explore early or partial global structural changes because of the computational limits (< 100 ns). In this work, we implement a hybrid force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid, to investigate the proton‐dependent dynamics and conformational changes of LacY. The effects of the protonation states on two key glutamate residues (Glu325 and Glu269) have been studied. Our results on the salt‐bridge dynamics agreed with all‐atom simulations at early short time period, validating our simulations. From our microsecond simulations, we were able to observe the complete transition from inward‐facing to outward‐facing conformations of LacY. Our results showed that all helices have participated during the global conformational transitions and helical movements of LacY. The inter‐helical distances measured in our simulations were consistent with the double electron‐electron resonance experiments at both cytoplasmic and periplasmic sides. Our simulations indicated that the deprotonation of Glu325 induced the opening of the periplasmics side and partial closure of the cytoplasmic side of LacY, while protonation of the Glu269 caused a stable cross‐domain salt‐bridge (Glu130‐Arg344) and completely closed the cytoplasmic side. Proteins 2016; 84:1067–1074. © 2016 Wiley Periodicals, Inc.     

The Five-To-Six-Coordination Transition of Ferric Human Serum Heme-Albumin Is Allosterically-Modulated by Ibuprofen and Warfarin: A Combined XAS and MD Study

Human serum albumin (HSA) is involved physiologically in heme scavenging; in turn, heme-albumin (HSA-heme-Fe) displays globin-like properties. Here, the allosteric effect of ibuprofen and warfarin on the local atomic structure around the ferric heme-Fe (heme-Fe(III)) atom of HSA-heme-Fe (HSA-heme-Fe(III)) has been probed by Fe-K edge X-ray absorption spectroscopy (XAS). The quantitative analysis of the Fe-K edge extended X-ray absorption fine structure (EXAFS) signals and modeling of the near edge (XANES) spectral features demonstrated that warfarin and ibuprofen binding modify the local structure of the heme-Fe(III). Combined XAS data analysis and targeted molecular dynamics (MD) simulations provided atomic resolution insights of protein structural rearrangements required to accommodate the heme-Fe(III) upon ibuprofen and warfarin binding. In the absence of drugs, the heme-Fe(III) atom is penta-coordinated having distorted 4+1 configuration made by the nitrogen atoms of the porphyrin ring and the oxygen phenoxy atom of the Tyr161 residue. MD simulations show that ibuprofen and warfarin association to the secondary fatty acid (FA) binding site 2 (FA2) induces a reorientation of domain I of HSA-heme-Fe(III), this leads to the redirection of the His146 residue providing an additional bond to the heme-Fe(III) atom, providing the 5+1 configuration. The comparison of Fe-K edge XANES spectra calculated using MD structures with those obtained experimentally confirms the reliability of the proposed structural model. As a whole, combining XAS and MD simulations it has been possible to provide a reliable model of the heme-Fe(III) atom coordination state and to understand the complex allosteric transition occurring in HSA-heme-Fe(III) upon ibuprofen and warfarin binding.     

Template‐free protein structure prediction and quality assessment with an all‐atom free‐energy model

Biophysical forcefields have contributed less than originally anticipated to recent progress in protein structure prediction. Here, we have investigated the selectivity of a recently developed all‐atom free‐energy forcefield for protein structure prediction and quality assessment (QA). Using a heuristic method, but excluding homology, we generated decoy‐sets for all targets of the CASP7 protein structure prediction assessment with <150 amino acids. The decoys in each set were then ranked by energy in short relaxation simulations and the best low‐energy cluster was submitted as a prediction. For four of nine template‐free targets, this approach generated high‐ranking predictions within the top 10 models submitted in CASP7 for the respective targets. For these targets, our de‐novo predictions had an average GDT_S score of 42.81, significantly above the average of all groups. The refinement protocol has difficulty for oligomeric targets and when no near‐native decoys are generated in the decoy library. For targets with high‐quality decoy sets the refinement approach was highly selective. Motivated by this observation, we rescored all server submissions up to 200 amino acids using a similar refinement protocol, but using no clustering, in a QA exercise. We found an excellent correlation between the best server models and those with the lowest energy in the forcefield. The free‐energy refinement protocol may thus be an efficient tool for relative QA and protein structure prediction. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

A multi‐resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand‐binding site

In multi‐resolution simulations, different system components are simultaneously modeled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of ligand binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modeled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse‐grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically. This multi‐resolution methodology can successfully model stable ligand binding, and we further confirm its validity by exploring the reproduction of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function. Proteins 2016; 84:1902–1913. © 2016 Wiley Periodicals, Inc.