ZDOCK server: interactive docking prediction of protein-protein complexes and symmetric multimers

Target-Driven Positive Selection at Hot Spots of Scorpion Toxins Uncovers Their Potential in Design of Insecticides

Positive selection sites (PSSs), a class of amino acid sites with an excess of nonsynonymous to synonymous substitutions, are indicators of adaptive molecular evolution and have been detected in many protein families involved in a diversity of biological processes by statistical approaches. However, few studies are conducted to evaluate their functional significance and the driving force behind the evolution (i.e., agent of selection). Scorpion α-toxins are a class of multigene family of peptide neurotoxins affecting voltage-gated Na(+ )(Nav) channels, whose members exhibit differential potency and preference for insect and mammalian Nav channels. In this study, we undertook a systematical molecular dissection of nearly all the PSSs newly characterized in the Mesobuthus α-toxin family and a two-residue insertion ((19)AlaPhe(20)) located within a positively selected loop via mutational analysis of α-like MeuNaTxα-5, one member affecting both insect and mammalian Nav channels. This allows to identify hot-spot residues on its functional face involved in interaction with the receptor site of Nav channels, which comprises two PSSs (Ile(40) and Leu(41)) and the small insertion, both located on two spatially separated functional loops. Mutations at these hot-spots resulted in a remarkably decreased anti-mammalian activity in MeuNaTxα-5 with partially impaired or enhanced insecticide activity, suggesting the potential of PSSs in designing promising candidate insecticides from scorpion α-like toxins. Based on an experiment-guided toxin-channel complex model and high evolutionary variability in the receptor site of predators and prey of scorpions, we provide new evidence for target-driven adaptive evolution of scorpion toxins to deal with their targets' diversity.

Predicted structure of MIF/CD74 and RTL1000/CD74 complexes

Macrophage migration inhibitory factor (MIF) is a key cytokine in autoimmune and inflammatory diseases that attracts and then retains activated immune cells from the periphery to the tissues. MIF exists as a homotrimer and its effects are mediated through its primary receptor, CD74 (the class II invariant chain that exhibits a highly structured trimerization domain), present on class II expressing cells. Although a number of binding residues have been identified between MIF and CD74 trimers, their spatial orientation has not been established. Using a docking program in silico, we have modeled binding interactions between CD74 and MIF as well as CD74 and a competitive MIF inhibitor, RTL1000, a partial MHC class II construct that is currently in clinical trials for multiple sclerosis. These analyses revealed 3 binding sites on the MIF trimer that each were predicted to bind one CD74 trimer through interactions with two distinct 5 amino acid determinants. Surprisingly, predicted binding of one CD74 trimer to a single RTL1000 antagonist utilized the same two 5 residue determinants, providing strong suggestive evidence in support of the MIF binding regions on CD74. Taken together, our structural modeling predicts a new MIF(CD74)3 dodecamer that may provide the basis for increased MIF potency and the requirement for ~3-fold excess RTL1000 to achieve full antagonism.

Modeling of the Binding of Peptide Blockers to Voltage-Gated Potassium Channels: Approaches and Evidence

Modeling of the structure of voltage-gated potassium (KV) channels bound to peptide blockers aims to identify the key amino acid residues dictating affinity and provide insights into the toxin-channel interface. Computational approaches open up possibilities for in silico rational design of selective blockers, new molecular tools to study the cellular distribution and functional roles of potassium channels. It is anticipated that optimized blockers will advance the development of drugs that reduce over activation of potassium channels and attenuate the associated malfunction. Starting with an overview of the recent advances in computational simulation strategies to predict the bound state orientations of peptide pore blockers relative to KV-channels, we go on to review algorithms for the analysis of intermolecular interactions, and then take a look at the results of their application.

Homology modeling, docking, and molecular dynamics simulation of the receptor GALR2 and its interactions with galanin and a positive allosteric modulator

Galanin receptor type 2 (GALR2) is a class A G-protein-coupled receptor (GPCR), and it has been reported that orthosteric ligands and positive allosteric modulators (PAMs) of GALR2 could potentially be used to treat epilepsy. So far, the X-ray structure of this receptor has not been resolved, and knowledge of the 3D structure of GALR2 may prove informative in attempts to design novel ligands and to explore the mechanism for the allosteric modulation of this receptor. In this study, homology modeling was used to obtain several GALR2 models using known templates. ProSA-web Z-scores and Ramachandran plots as well as pre-screening against a test dataset of known compounds were all utilized to select the best model of GALR2. Molecular dockings of galanin (a peptide) and a nonpeptide ligand were carried out to choose the (GALR2 model)-galanin complex that showed the closest agreement with the corresponding experimental data. Finally, a 50-ns MD simulation was performed to study the interactions between the GALR2 model and the synthetic and endogenous ligands. The results from docking and MD simulation showed that, besides the reported residues, Tyr160(4.60), Ile105(3.32), Ala274(7.35), and Tyr163(ECL2) also appear to play important roles in the binding of galanin. The potential allosteric binding pockets in the GALR2 model were then investigated via MD simulation. The results indicated that the mechanism for the allosteric modulation caused by PAMs is the binding of the PAM at pocket III, which is formed by galanin, ECL2, TM2, TM3, and ECL1; this results in the disruption of the Na(+)-binding site and/or the Na(+) ion pathway, leading to GALR2 agonism.

FRODOCK 2.0: fast protein–protein docking server

Summary: The prediction of protein–protein complexes from the structures of unbound components is a challenging and powerful strategy to decipher the mechanism of many essential biological processes. We present a user-friendly protein–protein docking server based on an improved version of FRODOCK that includes a complementary knowledge-based potential. The web interface provides a very effective tool to explore and select protein–protein models and interactively screen them against experimental distance constraints. The competitive success rates and efficiency achieved allow the retrieval of reliable potential protein–protein binding conformations that can be further refined with more computationally demanding strategies.

Availability and Implementation: The server is free and open to all users with no login requirement at http://frodock.chaconlab.org

Contact:  pablo@chaconlab.org

Supplementary information:  Supplementary data are available at Bioinformatics online.

Inhibition of PRL-2·CNNM3 Protein Complex Formation Decreases Breast Cancer Proliferation and Tumor Growth

The oncogenic phosphatase of regenerating liver 2 (PRL-2) has been shown to regulate intracellular magnesium levels by forming a complex through an extended amino acid loop present in the Bateman module of the CNNM3 magnesium transporter. Here we identified highly conserved residues located on this amino acid loop critical for the binding with PRL-2. A single point mutation (D426A) of one of those critical amino acids was found to completely disrupt PRL-2·human Cyclin M 3 (CNNM3) complex formation. Whole-cell voltage clamping revealed that expression of CNNM3 influenced the surface current, whereas overexpression of the binding mutant had no effect, indicating that the binding of PRL-2 to CNNM3 is important for the activity of the complex. Interestingly, overexpression of the CNNM3 D426A-binding mutant in cancer cells decreased their ability to proliferate under magnesium-deprived situations and under anchorage-independent growth conditions, demonstrating a PRL-2·CNNM3 complex-dependent oncogenic advantage in a more stringent environment. We further confirmed the importance of this complex in vivo using an orthotopic xenograft breast cancer model. Finally, because molecular modeling showed that the Asp-426 side chain in CNNM3 buries into the catalytic cavity of PRL-2, we showed that a PRL inhibitor could abrogate complex formation, resulting in a decrease in proliferation of human breast cancer cells. In summary, we provide evidence that this fundamental regulatory aspect of PRL-2 in cancer cells could potentially lead to broadly applicable and innovative therapeutic avenues.

iBET: Immersive visualization of biological electron-transfer dynamics

Recently, we presented a computational framework named VizBET to simulate and visualize biological electron-transfer (ET) dynamics. The visualization process was encapsulated as a plugin to the Visual Molecular Dynamics (VMD) software. However, the user's ability to understand complex, multidimensional ET pathways was severely limited when visualized in 2D on traditional computer monitors. To provide a more accurate representation with enhanced depth perception, we here present an extension of VizBET named iBET to render the VMD model of ET dynamics in a commodity virtual reality (VR) platform. The paper describes detailed procedures to export VMD models into the Unity game engine and render it in an Oculus Rift head mounted display. With the increasing availability of low-cost VR systems like the Rift and rich programmability of game engines, the iBET framework provides a powerful means to explore and understand not only biological ET processes but also a unique experiential tool for broad scientific communities.

cDNA sequence and protein bioinformatics analyses of MSTN in African catfish (Clarias gariepinus)

Myostatin, also known as growth differentiation factor 8, has been identified as a potent negative regulator of skeletal muscle growth. The purpose of this study was to characterize and predict function of the myostatin gene of the African catfish (Cg-MSTN). Expression of Cg-MSTN was determined at three growth stages to establish the relationship between the levels of MSTN transcript and skeletal muscle growth. The partial cDNA sequence of Cg-MSTN was cloned by using published information from its congener walking catfish (Cm-MSTN). The Cg-MSTN was 1194 bp in length encoding a protein of 397 amino acids. The deduced MSTN sequence exhibited key functional sites similar to those of other members of the TGF-β superfamily, especially, the proteolytic processing site (RXXR motif) and nine conserved cysteines at the C-terminal. Expression of MSTN appeared to be correlated with muscle development and growth of African catfish. Protein bioinformatics revealed that the primary sequence of Cg-MSTN shared 98 % sequence identity with that of walking catfish Cm-MSTN with only two different residues, [Formula: see text]. and [Formula: see text]. The proposed model of Cg-MSTN revealed the key point mutation [Formula: see text] causing a 7.35 Å shorter distance between the N- and C-lobes and an approximately 11° narrow angle than those of Cm-MSTN. The substitution of a proline residue near the proteolytic processing site which altered the structure of myostatin may play a critical role in reducing proteolytic activity of this protein in African catfish.

A potential therapeutic peptide-based neutralizer that potently inhibits Shiga toxin 2 in vitro and in vivo

Shiga toxin 2 (Stx2) is a major virulence factor in infections with Stx-producing Escherichia coli (STEC), which can cause serious clinical complications in humans, such as hemolytic uremic syndrome (HUS). Recently, we screened and identified two peptide-based Stx2 neutralizers, TF-1 and WA-8, which specifically and directly bind to Stx2. Computer simulations suggested that the majority of TF-1 or WA-8 binds tightly at the receptor-binding site 3 of Stx2. The two peptides also effectively inhibited the cytotoxic activity of Stx2 by blocking the binding of Stx2 to target cells. TF-1 exhibits remarkable therapeutic potency in both mice and rat toxicity models. In mice toxicity models, TF-1 provided full protection when mice were injected with 5 LD50 of Stx2. In rat toxicity models, TF-1 reduced fatal tissue damage and completely protected rats from the lethal challenges of Stx2. In these rats, TF-1 significantly decreased the concentration of Stx2 in blood and diminished tissue distribution levels of Stx2. Furthermore, TF-1 effectively protected rats from the pathological effects caused by Stx2, especially in the kidney, thymus, adrenal gland, and lung. Taken together, these results indicate that TF-1 is a promising therapeutic agent against the pathogenicity of Stx2.

Computational analysis of the CB1 carboxyl-terminus in the receptor-G protein complex

Despite the important role of the carboxyl-terminus (Ct) of the activated brain cannabinoid receptor one (CB1) in the regulation of G protein signaling, a structural understanding of interactions with G proteins is lacking. This is largely due to the highly flexible nature of the CB1 Ct that dynamically adapts its conformation to the presence of G proteins. In the present study, we explored how the CB1 Ct can interact with the G protein by building on our prior modeling of the CB1-Gi complex (Shim, Ahn, and Kendall, The Journal of Biological Chemistry 2013;288:32449-32465) to incorporate a complete CB1 Ct (Glu416(Ct)-Leu472(Ct)). Based on the structural constraints from NMR studies, we employed ROSETTA to predict tertiary folds, ZDOCK to predict docking orientation, and molecular dynamics (MD) simulations to obtain two distinct plausible models of CB1 Ct in the CB1-Gi complex. The resulting models were consistent with the NMR-determined helical structure (H9) in the middle region of the CB1 Ct. The CB1 Ct directly interacted with both Gα and Gβ and stabilized the receptor at the Gi interface. The results of site-directed mutagenesis studies of Glu416(Ct), Asp423(Ct), Asp428(Ct), and Arg444(Ct) of CB1 Ct suggested that the CB1 Ct can influence receptor-G protein coupling by stabilizing the receptor at the Gi interface. This research provided, for the first time, models of the CB1 Ct in contact with the G protein.

Screening Platform toward New Anti-HIV Aptamers Set on Molecular Docking and Fluorescence Quenching Techniques

By using a new rapid screening platform set on molecular docking simulations and fluorescence quenching techniques, three new anti-HIV aptamers targeting the viral surface glycoprotein 120 (gp120) were selected, synthesized, and assayed. The use of the short synthetic fluorescent peptide V35-Fluo mimicking the V3 loop of gp120, as the molecular target for fluorescence-quenching binding affinity studies, allowed one to measure the binding affinities of the new aptamers for the HIV-1 gp120 without the need to obtain and purify the full recombinant gp120 protein. The almost perfect correspondence between the calculated Kd and the experimental EC50 on HIV-infected cells confirmed the reliability of the platform as an alternative to the existing methods for aptamer selection and measuring of aptamer-protein equilibria.

Protein unfolding as a switch from self-recognition to high-affinity client binding

Stress-specific activation of the chaperone Hsp33 requires the unfolding of a central linker region. This activation mechanism suggests an intriguing functional relationship between the chaperone's own partial unfolding and its ability to bind other partially folded client proteins. However, identifying where Hsp33 binds its clients has remained a major gap in our understanding of Hsp33's working mechanism. By using site-specific Fluorine-19 nuclear magnetic resonance experiments guided by in vivo crosslinking studies, we now reveal that the partial unfolding of Hsp33's linker region facilitates client binding to an amphipathic docking surface on Hsp33. Furthermore, our results provide experimental evidence for the direct involvement of conditionally disordered regions in unfolded protein binding. The observed structural similarities between Hsp33's own metastable linker region and client proteins present a possible model for how Hsp33 uses protein unfolding as a switch from self-recognition to high-affinity client binding.

Under stress conditions the molecular chaperone Hsp33 is activated to process unfolded proteins. Here, the authors use in vivo and in vitro crosslinking and ¹⁹F-NMR to elucidate the binding site for misfolded proteins and are able to propose a model for its mechanism of action.

Interactions of a class IIb bacteriocin with a model lipid bilayer, investigated through molecular dynamics simulations

The emergence of antibiotic resistant microorganisms poses an alarming threat to global health. Antimicrobial peptides (AMPs) are considered a possible effective alternative to conventional antibiotic therapies. An understanding of the mechanism of action of AMPs is needed in order to better control and optimize their bactericidal activity. Plantaricin EF is a heterodimeric AMP, consisting of two peptides Plantaricin E (PlnE) and Plantaricin F (PlnF). We studied the behavior of these peptides on the surface of a model lipid bilayer. We identified the residues that facilitate peptide-peptide interactions. We also identified residues that mediate interactions of the dimer with the membrane. PlnE interacts with the membrane through amino acids at both its termini, while only the N terminus of PlnF approaches the membrane. By comparing the activity of single-site mutants of the two-peptide bacteriocin and the simulations of the bacteriocin on the surface of a model lipid bilayer, structure activity relationships are proposed. These studies allow us to generate hypotheses that relate biophysical interactions observed in simulations with the experimentally measured activity. We find that single-site amino acid substitutions result in markedly stronger antimicrobial activity when they strengthen the interactions between the two peptides, while, concomitantly, they weaken peptide-membrane association. This effect is more pronounced in the case of the PlnE mutant (G20A), which interacts the strongest with PlnF and the weakest with the membrane while displaying the highest activity.

Principles of assembly reveal a periodic table of protein complexes

Structural insights into protein complexes have had a broad impact on our understanding of biological function and evolution. In this work, we sought a comprehensive understanding of the general principles underlying quaternary structure organization in protein complexes. We first examined the fundamental steps by which protein complexes can assemble, using experimental and structure-based characterization of assembly pathways. Most assembly transitions can be classified into three basic types, which can then be used to exhaustively enumerate a large set of possible quaternary structure topologies. These topologies, which include the vast majority of observed protein complex structures, enable a natural organization of protein complexes into a periodic table. On the basis of this table, we can accurately predict the expected frequencies of quaternary structure topologies, including those not yet observed. These results have important implications for quaternary structure prediction, modeling, and engineering.

A human monoclonal antibody against HPV16 recognizes an immunodominant and neutralizing epitope partially overlapping with that of H16.V5

The presence of neutralizing epitopes in human papillomavirus (HPV) L1 virus-like particles (VLPs) is the structural basis of prophylactic vaccines. An anti-HPV16 neutralizing monoclonal antibody (N-mAb) 26D1 was isolated from a memory B cell of a human vaccinee. The pre-binding of heparan sulfate to VLPs inhibited the binding of both N-mAbs to the antigen, indicating that the epitopes are critical for viral cell attachment/entry. Hybrid VLP binding with surface loop swapping between types indicated the essential roles of the DE and FG loops for both 26D1 (DEa in particular) and H16.V5 binding. Specifically, Tyr(135) and Val(141) on the DEa loop were shown to be critical residues for 26D1 binding via site-directed mutagenesis. Partially overlap between the epitopes between 26D1 and H16.V5 was shown using pairwise epitope mapping, and their binding difference is demonstrated to be predominantly in DE loop region. In addition, 26D1 epitope is immunodominant epitope recognized by both antibodies elicited by the authentic virus from infected individuals and polyclonal antibodies from vaccinees. Overall, a partially overlapping but distinct neutralizing epitope from that of H16.V5 was identified using a human N-mAb, shedding lights to the antibody arrays as part of human immune response to vaccination and infection.

Molecular insight into the activity of LasR protein from Pseudomonas aeruginosa in the regulation of virulence gene expression by this organism

Pseudomonas aeruginosa is an opportunistic human pathogen. This organism attacks human patients suffering from diseases like AIDS, cancer, cystic fibrosis, etc. One of the important virulent factors produced by this organism is Hydrogen Cyanide. This is expressed from the genes encoded by the hcnABC operon. The expressions of the genes encoded by hcnABC operon are mediated mainly by the interactions of LasR protein with the corresponding promoter region of the hcnABC operon. The LasR protein acts as a dimer and binds to the promoter DNA with the help of an autoinducer ligand. However, till date the detailed molecular mechanism of how the LasR protein interacts with the promoter DNA is not clearly known. Therefore, in this work, an attempt has been made to analyze the mode of interactions of the LasR protein with the promoter DNA region of the hcnABC operon. We analyzed the three dimensional structure of the LasR protein from Pseudomonas aeruginosa and docked the protein with the autoinducer ligand. We then docked the ligand-bound-LasR-protein as well the LasR-protein-without-the-autoinducer-ligand on to the promoter DNA region of hcnABC operon. We analyzed the details of the interaction profiles of LasR protein with the autoinducer ligand. We also deciphered the details of the LasR promoter-DNA interactions. We compared the modes of DNA bindings by the LasR protein in presence and absence of the autoinducer ligand and tried to analyze the molecular details of the binding of LasR protein with the promoter DNA region of hcnABC operon during hcnABC gene expression. This study may therefore pave the pathway for future experiments to determine the relative effects of the amino acid residues of LasR protein in DNA binding during the transcription of hcnABC operon.

Rigid-Docking Approaches to Explore Protein-Protein Interaction Space

Protein-protein interactions play core roles in living cells, especially in the regulatory systems. As information on proteins has rapidly accumulated on publicly available databases, much effort has been made to obtain a better picture of protein-protein interaction networks using protein tertiary structure data. Predicting relevant interacting partners from their tertiary structure is a challenging task and computer science methods have the potential to assist with this. Protein-protein rigid docking has been utilized by several projects, docking-based approaches having the advantages that they can suggest binding poses of predicted binding partners which would help in understanding the interaction mechanisms and that comparing docking results of both non-binders and binders can lead to understanding the specificity of protein-protein interactions from structural viewpoints. In this review we focus on explaining current computational prediction methods to predict pairwise direct protein-protein interactions that form protein complexes.

Converting the bis-FeIV state of the diheme enzyme MauG to Compound I decreases the reorganization energy for electron transfer

The electron transfer (ET) properties of two types of high-valent hemes were studied within the same protein matrix; the bis-Fe(IV) state of MauG and the Compound I state of Y294H MauG. The latter is formed as a consequence of mutation of the tyrosine which forms the distal axial ligand of the six-coordinate heme that allows it to stabilize Fe(IV) in the absence of an external ligand. The rates of the ET reaction of each high-valent species with the type I copper protein, amicyanin, were determined at different temperatures and analysed by ET theory. The reaction with bis-Fe(IV) wild-type (WT) MauG exhibited a reorganization energy (λ) that was 0.39 eV greater than that for the reaction of Compound I Y295H MauG. It is concluded that the delocalization of charge over the two hemes in the bis-Fe(IV) state is responsible for the larger λ, relative to the Compound I state in which the Fe(V) equivalent is isolated on one heme. Although the increase in λ decreases the rate of ET, the delocalization of charge decreases the ET distance to its natural substrate protein, thus increasing the ET rate. This describes how proteins can balance different ET properties of complex redox cofactors to optimize each system for its particular ET or catalytic reaction.

Refining the Structural Model of a Heterohexameric Protein Complex: Surface Induced Dissociation and Ion Mobility Provide Key Connectivity and Topology Information

Toyocamycin nitrile hydratase (TNH) is a protein hexamer that catalyzes the hydration of toyocamycin to produce sangivamycin. The structure of hexameric TNH and the arrangement of subunits within the complex, however, have not been solved by NMR or X-ray crystallography. Native mass spectrometry (MS) clearly shows that TNH is composed of two copies each of the α, β, and γ subunits. Previous surface induced dissociation (SID) tandem mass spectrometry on a quadrupole time-of-flight (QTOF) platform suggests that the TNH hexamer is a dimer composed of two αβγ trimers; furthermore, the results suggest that α-β interact most strongly (Blackwell et al. Anal. Chem. 2011, 83, 2862-2865). Here, multiple complementary MS based approaches and homology modeling have been applied to refine the structure of TNH. Solution-phase organic solvent disruption coupled with native MS agrees with the previous SID results. By coupling surface induced dissociation with ion mobility mass spectrometry (SID/IM), further information on the intersubunit contacts and relative interfacial strengths are obtained. The results show that TNH is a dimer of αβγ trimers, that within the trimer the α, β subunits bind most strongly, and that the primary contact between the two trimers is through a γ-γ interface. Collisional cross sections (CCSs) measured from IM experiments are used as constraints for postulating the arrangement of the subunits represented by coarse-grained spheres. Covalent labeling (surface mapping) together with protein complex homology modeling and docking of trimers to form hexamer are utilized with all the above information to propose the likely quaternary structure of TNH, with chemical cross-linking providing cross-links consistent with the proposed structure. The novel feature of this approach is the use of SID-MS with ion mobility to define complete connectivity and relative interfacial areas of a heterohexameric protein complex, providing much more information than is available from solution disruption. That information, when combined with CCS-guided coarse-grained modeling and covalent labeling restraints for homology modeling and trimer-trimer docking, provides atomic models of a previously uncharacterized heterohexameric protein complex.

Structural and functional interactions between six-transmembrane μ-opioid receptors and β2-adrenoreceptors modulate opioid signaling

The primary molecular target for clinically used opioids is the μ-opioid receptor (MOR). Besides the major seven-transmembrane (7TM) receptors, the MOR gene codes for alternatively spliced six-transmembrane (6TM) isoforms, the biological and clinical significance of which remains unclear. Here, we show that the otherwise exclusively intracellular localized 6TM-MOR translocates to the plasma membrane upon coexpression with β₂-adrenergic receptors (β₂-ARs) through an interaction with the fifth and sixth helices of β₂-AR. Coexpression of the two receptors in BE(2)-C neuroblastoma cells potentiates calcium responses to a 6TM-MOR ligand, and this calcium response is completely blocked by a selective β₂-antagonist in BE(2)-C cells, and in trigeminal and dorsal root ganglia. Co-administration of 6TM-MOR and β₂-AR ligands leads to substantial analgesic synergy and completely reverses opioid-induced hyperalgesia in rodent behavioral models. Together, our results provide evidence that the heterodimerization of 6TM-MOR with β₂-AR underlies a molecular mechanism for 6TM cellular signaling, presenting a unique functional responses to opioids. This signaling pathway may contribute to the hyperalgesic effects of opioids that can be efficiently blocked by β₂-AR antagonists, providing a new avenue for opioid therapy.

Binding of the Covalent Flavin Assembly Factor to the Flavoprotein Subunit of Complex II

Escherichia coli harbors two highly conserved homologs of the essential mitochondrial respiratory complex II (succinate:ubiquinone oxidoreductase). Aerobically the bacterium synthesizes succinate:quinone reductase as part of its respiratory chain, whereas under microaerophilic conditions, the quinol:fumarate reductase can be utilized. All complex II enzymes harbor a covalently bound FAD co-factor that is essential for their ability to oxidize succinate. In eukaryotes and many bacteria, assembly of the covalent flavin linkage is facilitated by a small protein assembly factor, termed SdhE in E. coli. How SdhE assists with formation of the covalent flavin bond and how it binds the flavoprotein subunit of complex II remain unknown. Using photo-cross-linking, we report the interaction site between the flavoprotein of complex II and the SdhE assembly factor. These data indicate that SdhE binds to the flavoprotein between two independently folded domains and that this binding mode likely influences the interdomain orientation. In so doing, SdhE likely orients amino acid residues near the dicarboxylate and FAD binding site, which facilitates formation of the covalent flavin linkage. These studies identify how the conserved SdhE assembly factor and its homologs participate in complex II maturation.

Modification of a deoxynivalenol-antigen-mimicking nanobody to improve immunoassay sensitivity by site-saturation mutagenesis

A nanobody (N-28) which can act as a deoxynivalenol (DON) antigen has been generated, and its residues Thr102-Ser106 were identified to bind with anti-DON monoclonal antibody by alanine-scanning mutagenesis. Site-saturation mutagenesis was used to analyze the plasticity of five residues and to improve the sensitivity of the N-28-based immunoassay. After mutagenesis, three mutants were selected by phage immunoassay and were sequenced. The half-maximal inhibitory concentrations of the immunoassay based on mutants N-28-T102Y, N-28-V103L, and N-28-Y105F were 24.49 ± 1.0, 51.83 ± 2.5, and 35.65 ± 1.6 ng/mL, respectively, showing the assay was, respectively, 3.2, 1.5, and 2.2 times more sensitive than the wild-type-based assay. The best mutant, N-28-T102Y, was used to develop a competitive phage ELISA to detect DON in cereals with high specificity and accuracy. In addition, the structural properties of N-28-T102Y and N-28 were investigated, revealing that the affinity of N-28-T102Y decreased because of increased steric hindrance with the large side chain. The lower-binding-affinity antigen mimetic may contribute to the improvement of the sensitivity of competitive immunoassays. These results demonstrate that nanobodies would be a favorable tool for engineering. Moreover, our results have laid a solid foundation for site-saturation mutagenesis of antigen-mimicking nanobodies to improve immunoassay sensitivity for small molecules.

Identification of the HIV-1 Vif and Human APOBEC3G Protein Interface

Human cells express natural antiviral proteins, such as APOBEC3G (A3G), that potently restrict HIV replication. As a counter-defense, HIV encodes the accessory protein Vif, which binds A3G and mediates its proteasomal degradation. Our structural knowledge on how Vif and A3G interact is limited, because a co-structure is not available. We identified specific points of contact between Vif and A3G by using functional assays with full-length A3G, patient-derived Vif variants, and HIV forced evolution. These anchor points were used to model and validate the Vif-A3G interface. The resultant co-structure model shows that the negatively charged β4-α4 A3G loop, which contains primate-specific variation, is the core Vif binding site and forms extensive interactions with a positively charged pocket in HIV Vif. Our data present a functional map of this viral-host interface and open avenues for targeted approaches to block HIV replication by obstructing the Vif-A3G interaction.

Mechanistic Scrutiny Identifies a Kinetic Role for Cytochrome b5 Regulation of Human Cytochrome P450c17 (CYP17A1, P450 17A1)

Cytochrome P450c17 (P450 17A1, CYP17A1) is a critical enzyme in the synthesis of androgens and is now a target enzyme for the treatment of prostate cancer. Cytochrome P450c17 can exhibit either one or two physiological enzymatic activities differentially regulated by cytochrome b5. How this is achieved remains unknown. Here, comprehensive in silico, in vivo and in vitro analyses were undertaken. Fluorescence Resonance Energy Transfer analysis showed close interactions within living cells between cytochrome P450c17 and cytochrome b5. In silico modeling identified the sites of interaction and confirmed that E48 and E49 residues in cytochrome b5 are essential for activity. Quartz crystal microbalance studies identified specific protein-protein interactions in a lipid membrane. Voltammetric analysis revealed that the wild type cytochrome b5, but not a mutated, E48G/E49G cyt b5, altered the kinetics of electron transfer between the electrode and the P450c17. We conclude that cytochrome b5 can influence the electronic conductivity of cytochrome P450c17 via allosteric, protein-protein interactions.

Peptide Functionalized Gold Nanorods for the Sensitive Detection of a Cardiac Biomarker Using Plasmonic Paper Devices

The sensitivity of localized surface plasmon resonance (LSPR) of metal nanostructures to adsorbates lends itself to a powerful class of label-free biosensors. Optical properties of plasmonic nanostructures are dependent on the geometrical features and the local dielectric environment. The exponential decay of the sensitivity from the surface of the plasmonic nanotransducer calls for the careful consideration in its design with particular attention to the size of the recognition and analyte layers. In this study, we demonstrate that short peptides as biorecognition elements (BRE) compared to larger antibodies as target capture agents offer several advantages. Using a bioplasmonic paper device (BPD), we demonstrate the selective and sensitive detection of the cardiac biomarker troponin I (cTnI). The smaller sized peptide provides higher sensitivity and a lower detection limit using a BPD. Furthermore, the excellent shelf-life and thermal stability of peptide-based LSPR sensors, which precludes the need for special storage conditions, makes it ideal for use in resource-limited settings.

Modelling the molecular mechanism of protein-protein interactions and their inhibition: CypD-p53 case study

Cyclophilin D (CypD) is an important regulatory protein involved in mitochondrial membrane permeability transition and cell death. Further, the mitochondrial CypD-p53 axis is an important contributor to necroptosis, a form of programmed necrosis, involved in various cardiovascular and neurological disorders. The CypD ligand, Cyclosporin A (CsA), was identified as an inhibitor of this interaction. In this study, using computational methods, we have attempted to model the CypD-p53 interaction in order to delineate their mode of binding and also to disclose the molecular mechanism, by means of which CsA interferes with this interaction. It was observed that p53 binds at the CsA-binding site of CypD. The knowledge obtained from this modelling was employed to identify novel CypD inhibitors through structure-based methods. Further, the identified compounds were tested by a similar strategy, adopted during the modelling process. This strategy could be applied to study the mechanism of protein-protein interaction (PPI) inhibition and to identify novel PPI inhibitors.

Screening and Identifying a Novel ssDNA Aptamer against Alpha-fetoprotein Using CE-SELEX

Alpha-fetoprotein (AFP) is a liver cancer associated protein and has long been utilized as a serum tumor biomarker of disease progression. AFP is usually detected in HCC patients by an antibody based system. Recently, however, aptamers generated from systematic evolution of ligands by exponential enrichment (SELEX) were reported to have an alternative potential in targeted imaging, diagnosis and therapy. In this study, AFP-bound ssDNA aptamers were screened and identified using capillary electrophoresis (CE) SELEX technology. After cloning, sequencing and motif analysis, we successfully confirmed an aptamer, named AP273, specifically targeting AFP. The aptamer could be used as a probe in AFP immunofluorescence imaging in HepG2, one AFP positive cancer cell line, but not in A549, an AFP negative cancer cell line. More interesting, the aptamer efficiently inhibited the migration and invasion of HCC cells after in vivo transfection. Motif analysis revealed that AP273 had several stable secondary motifs in its structure. Our results indicate that CE-SELEX technology is an efficient method to screen specific protein-bound ssDNA, and AP273 could be used as an agent in AFP-based staining, diagnosis and therapy, although more works are still needed.

Signaling and Adaptation Modulate the Dynamics of the Photosensoric Complex of Natronomonas pharaonis

Motile bacteria and archaea respond to chemical and physical stimuli seeking optimal conditions for survival. To this end transmembrane chemo- and photoreceptors organized in large arrays initiate signaling cascades and ultimately regulate the rotation of flagellar motors. To unravel the molecular mechanism of signaling in an archaeal phototaxis complex we performed coarse-grained molecular dynamics simulations of a trimer of receptor/transducer dimers, namely NpSRII/NpHtrII from Natronomonas pharaonis. Signaling is regulated by a reversible methylation mechanism called adaptation, which also influences the level of basal receptor activation. Mimicking two extreme methylation states in our simulations we found conformational changes for the transmembrane region of NpSRII/NpHtrII which resemble experimentally observed light-induced changes. Further downstream in the cytoplasmic domain of the transducer the signal propagates via distinct changes in the dynamics of HAMP1, HAMP2, the adaptation domain and the binding region for the kinase CheA, where conformational rearrangements were found to be subtle. Overall these observations suggest a signaling mechanism based on dynamic allostery resembling models previously proposed for E. coli chemoreceptors, indicating similar properties of signal transduction for archaeal photoreceptors and bacterial chemoreceptors.

Achaea and bacteria can “see” and “sniffle”, they have photo- and chemosensors that measure the environment. On the cell poles, these sensor proteins form large arrays built of several thousands of different receptors. The receptors comprise extracellular or transmembrane sensory domains and elongated homodimeric coiled-coil bundles, which transduce the signals from the membrane across ~20 nm to a conserved cytoplasmic signaling subdomain in an unknown manner. In our study we performed coarse-grained molecular dynamics simulations of the phototactic receptor/transducer complex from Natronomonas pharaonis. Comparing fully methylated and demethylated complexes reveals an interconversion between states of different dynamics along the coiled-coil bundle, which might represent the essential characteristics of the signal transfer from the membrane to the binding sites of the downstream kinase CheA.

HIV-1 Recruits UPF1 but Excludes UPF2 to Promote Nucleocytoplasmic Export of the Genomic RNA

Unspliced, genomic HIV-1 RNA (vRNA) is a component of several ribonucleoprotein complexes (RNP) during the viral replication cycle. In earlier work, we demonstrated that the host upframeshift protein 1 (UPF1), a key factor in nonsense-mediated mRNA decay (NMD), colocalized and associated to the viral structural protein Gag during viral egress. In this work, we demonstrate a new function for UPF1 in the regulation of vRNA nuclear export. OPEN ACCESS Biomolecules 2015, 5 2809 We establish that the nucleocytoplasmic shuttling of UPF1 is required for this function and demonstrate that UPF1 exists in two essential viral RNPs during the late phase of HIV-1 replication: the first, in a nuclear export RNP that contains Rev, CRM1, DDX3 and the nucleoporin p62, and the second, which excludes these nuclear export markers but contains Gag in the cytoplasm. Interestingly, we observed that both UPF2 and the long isoform of UPF3a, UPF3aL, but not the shorter isoforms UPF3aS and UPF3b, are excluded from the UPF1-Rev-CRM1-DDX3 complex as they are negative regulators of vRNA nuclear export. In silico protein-protein docking analyses suggest that Rev binds UPF1 in a region that overlaps the UPF2 binding site, thus explaining the exclusion of this negative regulatory factor by HIV-1 that is necessary for vRNA trafficking. This work uncovers a novel and unique regulatory circuit involving several UPF proteins that ultimately regulate vRNA nuclear export and trafficking.

Structural Insights into the MMACHC-MMADHC Protein Complex Involved in Vitamin B12 Trafficking

Conversion of vitamin B12 (cobalamin, Cbl) into the cofactor forms methyl-Cbl (MeCbl) and adenosyl-Cbl (AdoCbl) is required for the function of two crucial enzymes, mitochondrial methylmalonyl-CoA mutase and cytosolic methionine synthase, respectively. The intracellular proteins MMACHC and MMADHC play important roles in processing and targeting the Cbl cofactor to its destination enzymes, and recent evidence suggests that they may interact while performing these essential trafficking functions. To better understand the molecular basis of this interaction, we have mapped the crucial protein regions required, indicate that Cbl is likely processed by MMACHC prior to interaction with MMADHC, and identify patient mutations on both proteins that interfere with complex formation, via different mechanisms. We further report the crystal structure of the MMADHC C-terminal region at 2.2 Å resolution, revealing a modified nitroreductase fold with surprising homology to MMACHC despite their poor sequence conservation. Because MMADHC demonstrates no known enzymatic activity, we propose it as the first protein known to repurpose the nitroreductase fold solely for protein-protein interaction. Using small angle x-ray scattering, we reveal the MMACHC-MMADHC complex as a 1:1 heterodimer and provide a structural model of this interaction, where the interaction region overlaps with the MMACHC-Cbl binding site. Together, our findings provide novel structural evidence and mechanistic insight into an essential biological process, whereby an intracellular "trafficking chaperone" highly specific for a trace element cofactor functions via protein-protein interaction, which is disrupted by inherited disease mutations.

Structural analysis of inter-genus complexes of V-antigen and its regulator and their stabilization by divalent metal ions

Gram-negative bacteria like Yersinia, Pseudomonas, and Aeromonas need type III secretion system (T3SS) for their pathogenicity. V-antigen and its regulator are essential for functioning of T3SS. There is significant functional conservation amongst V-antigen and its regulator belonging to the Ysc family. In this study, we have structurally characterized the inter-genus complexes of V-antigen and its regulator. ConSurf analysis demonstrates that V-antigens belonging to the Ysc family show high structural identity predominantly confined to the two long helical regions. The regulator of V-antigen shows high conservation in its first intramolecular coiled-coil domain, responsible for interaction with V-antigen. ∆LcrG(1-70) localizes within the groove formed by long helices of LcrV, as observed in PcrV-∆PcrG(13-72) interaction. Inter-genus complexes of LcrV-PcrG and PcrV-LcrG exhibited elongated conformation and 1:1 heterodimeric state like the native complex of PcrV-PcrG and LcrV-LcrG. Both native and inter-genus complexes showed rigid tertiary structure, solvent-exposed hydrophobic patches, and cooperative melting behavior with high melting temperature. LcrV-PcrG and PcrV-LcrG showed nanomolar affinity of interaction, identical to PcrV-PcrG interaction, but stronger than LcrV-LcrG interaction. Calcium (a secretion blocker of T3SS) propels all the complexes towards a highly monodisperse form. Calcium and magnesium increase the helicity of the native and inter-genus complexes, and causes helix-helix stabilization. Stabilization of helices leads to a slight increase in the melting temperature by 1.5-2.0 °C. However, calcium does not alter the affinity of interaction of V-antigen and its regulator, emphasizing the effect of divalent of cations at the structural level without any regulatory implications. Therefore, the structural conservation of these inter-genus complexes could be the basis for their functional complementation.

Computational elucidation of potential antigenic CTL epitopes in Ebola virus

Cell-mediated immunity is important for the control of Ebola virus infection. We hypothesized that those HLA A0201 and HLA B40 restricted epitopes derived from Ebola virus proteins, would mount a good antigenic response. Here we employed an immunoinformatics approach to identify specific 9mer amino acid which may be capable of inducing a robust cell-mediated immune response in humans. We identified a set of 28 epitopes that had no homologs in humans. Specifically, the epitopes derived from NP, RdRp, GP and VP40 share population coverage of 93.40%, 84.15%, 74.94% and 77.12%, respectively. Based on the other HLA binding specificity and population coverage, seven novel promiscuous epitopes were identified. These 7 promiscuous epitopes from NP, RdRp and GP were found to have world-wide population coverage of more than 95% indicating their potential significance as useful candidates for vaccine design. Epitope conservancy analysis also suggested that most of the peptides are highly conserved (100%) in other virulent Ebola strain (Mayinga-76, Kikwit-95 and Makona-G3816- 2014) and can therefore be further investigated for their immunological relevance and usefulness as vaccine candidates.

Target-Driven Evolution of Scorpion Toxins

It is long known that peptide neurotoxins derived from a diversity of venomous animals evolve by positive selection following gene duplication, yet a force that drives their adaptive evolution remains a mystery. By using maximum-likelihood models of codon substitution, we analyzed molecular adaptation in scorpion sodium channel toxins from a specific species and found ten positively selected sites, six of which are located at the core-domain of scorpion α-toxins, a region known to interact with two adjacent loops in the voltage-sensor domain (DIV) of sodium channels, as validated by our newly constructed computational model of toxin-channel complex. Despite the lack of positive selection signals in these two loops, they accumulated extensive sequence variations by relaxed purifying selection in prey and predators of scorpions. The evolutionary variability in the toxin-bound regions of sodium channels indicates that accelerated substitutions in the multigene family of scorpion toxins is a consequence of dealing with the target diversity. This work presents an example of atypical co-evolution between animal toxins and their molecular targets, in which toxins suffered from more prominent selective pressure from the channels of their competitors. Our discovery helps explain the evolutionary rationality of gene duplication of toxins in a specific venomous species.

The HGF inhibitory peptide HGP-1 displays promising in vitro and in vivo efficacy for targeted cancer therapy

HGF/MET pathway mediates cancer initiation and development. Thus, inhibition on HGF-initiated MET signaling pathway would provide a new approach to cancer targeted therapeutics. In our study, we identified a targeting peptide candidate binding to HGF which was named HGF binding peptide-1 (HGP-1) via bacterial surface display methods coupled with fluorescence-activated cell sorting (FACS). HGP-1 showed the moderate affinity when determined with surface plasmon resonance (SPR) technique and high specificity in binding to HGF while assessed by fluorescence-based ELISA assay. The results from MTT and in vitro migration assay indicated that HGF-dependent cell proliferation and migration could be inhibited by HGP-1. In vivo administration of HGP-1 led to an effective inhibitory effect on tumor growth in A549 tumor xenograft models. Moreover, findings from Western Blots revealed that HGP-1 could down-regulated the phosphorylation levels of MET and ERK1/2 initiated by HGF, which suggested that HGP-1 could disrupt the activation of HGF/MET signaling to influence the cell activity. All the data highlighted the potential of HGP-1 to be a potent inhibitor for HGF/MET signaling.

MST2-RASSF protein-protein interactions through SARAH domains

The detailed, atomistic-level understanding of molecular signaling along the tumor-suppressive Hippo signaling pathway that controls tissue homeostasis by balancing cell proliferation and death through apoptosis is a promising avenue for the discovery of novel anticancer drug targets. The activation of kinases such as Mammalian STE20-Like Protein Kinases 1 and 2 (MST1 and MST2)-modulated through both homo- and heterodimerization (e.g. interactions with Ras association domain family, RASSF, enzymes)-is a key upstream event in this pathway and remains poorly understood. On the other hand, RASSFs (such as RASSF1A or RASSF5) act as important apoptosis activators and tumor suppressors, although their exact regulatory roles are also unclear. We present recent molecular studies of signaling along the Ras-RASSF-MST pathway, which controls growth and apoptosis in eukaryotic cells, including a variety of modern molecular modeling and simulation techniques. Using recently available structural information, we discuss the complex regulatory scenario according to which RASSFs perform dual signaling functions, either preventing or promoting MST2 activation, and thus control cell apoptosis. Here, we focus on recent studies highlighting the special role being played by the specific interactions between the helical Salvador/RASSF/Hippo (SARAH) domains of MST2 and RASSF1a or RASSF5 enzymes. These studies are crucial for integrating atomistic-level mechanistic information about the structures and conformational dynamics of interacting proteins, with information available on their system-level functions in cellular signaling.

Extension of resolution and oligomerization-state studies of 2,4'-dihydroxyacetophenone dioxygenase from Alcaligenes sp. 4HAP

The enzyme 2,4'-dihydroxyacetophenone dioxygenase (DAD) catalyses the conversion of 2,4'-dihydroxyacetophenone to 4-hydroxybenzoic acid and formic acid. This enzyme is a very unusual dioxygenase in that it cleaves a C-C bond in a substituent of the aromatic ring rather than within the ring itself. Whilst it has been shown that DAD is a tetramer in solution, the recently solved crystal structure of the Alcaligenes sp. 4HAP enzyme was in fact dimeric rather than tetrameric. Since the use of limited chymotrypsinolysis, which apparently results in removal of the first 20 or so N-terminal residues of DAD, was necessary for crystallization of the protein, it was investigated whether this was responsible for the change in its oligomerization state. Gel-filtration and analytical ultracentrifugation studies were conducted, which confirmed that chymotrypsinolysed DAD has an apparent molecular weight of around 40 kDa, corresponding to a dimer. In contrast, the native enzyme has a molecular weight in the 70-80 kDa region, as expected for the tetramer. The structural basis for tetramerization has been investigated by the use of several docking servers, and the results are remarkably consistent with the tetrameric structure of a homologous cupin protein from Ralstonia eutropha (PDB entry 3ebr).

A web interface for easy flexible protein-protein docking with ATTRACT

Protein-protein docking programs can give valuable insights into the structure of protein complexes in the absence of an experimental complex structure. Web interfaces can facilitate the use of docking programs by structural biologists. Here, we present an easy web interface for protein-protein docking with the ATTRACT program. While aimed at nonexpert users, the web interface still covers a considerable range of docking applications. The web interface supports systematic rigid-body protein docking with the ATTRACT coarse-grained force field, as well as various kinds of protein flexibility. The execution of a docking protocol takes up to a few hours on a standard desktop computer.

Updates to the Integrated Protein-Protein Interaction Benchmarks: Docking Benchmark Version 5 and Affinity Benchmark Version 2

We present an updated and integrated version of our widely used protein-protein docking and binding affinity benchmarks. The benchmarks consist of non-redundant, high-quality structures of protein-protein complexes along with the unbound structures of their components. Fifty-five new complexes were added to the docking benchmark, 35 of which have experimentally measured binding affinities. These updated docking and affinity benchmarks now contain 230 and 179 entries, respectively. In particular, the number of antibody-antigen complexes has increased significantly, by 67% and 74% in the docking and affinity benchmarks, respectively. We tested previously developed docking and affinity prediction algorithms on the new cases. Considering only the top 10 docking predictions per benchmark case, a prediction accuracy of 38% is achieved on all 55 cases and up to 50% for the 32 rigid-body cases only. Predicted affinity scores are found to correlate with experimental binding energies up to r=0.52 overall and r=0.72 for the rigid complexes.

Transcripts involved in hemostasis: Exploring salivary complexes from Haementeria vizottoi leeches through transcriptomics, phylogenetic studies and structural features

Throughout evolution, parasites have adapted in order to successfully intervene in the host defense, producing specific peptides and proteins. Interestingly, these peptides and proteins have been exploited as potential drug candidates against several diseases. Furthermore, biotechnology studies and cDNA libraries have remarkably contributed to identify potentially bioactive molecules. In this regard, herein, a cDNA library of salivary complexes from Haementeria vizottoi leeches was constructed, the transcriptome was characterized and a phylogenetic analysis was performed considering antistasin-like and antiplatelet-like proteins. Hundred twenty three transcripts were identified coding for putative proteins involved in animal feeding (representing about 10% of the expression level). These sequences showed similarities with myohemerythrins, carbonic anhydrases, anticoagulants, antimicrobials, proteases and protease inhibitors. The phylogenetic analysis, regarding antistasin-like and antiplatetlet-like proteins, revealed two main clades in the Rhynchobdellida leeches. As expected, the sequences from H. vizottoi have presented high similarities with those types of proteins. Thus, our findings could be helpful not only to identify new coagulation inhibitors, but also to better understand the biological composition of the salivary complexes.

HK022 Nun Requires Arginine-Rich Motif Residues Distinct from λ N

Bacteriophage λ N protein binds boxB RNA hairpins in the nut (N utilization) sites of immediate early λ transcripts and interacts with host factors to suppress transcriptional termination at downstream terminators. In opposition to λ N, the Nun protein of HK022 binds the boxBs of coinfecting λ transcripts, interacts with a similar or identical set of host factors, and terminates transcription to suppress λ replication. Comparison of N-boxB and Nun-boxB nuclear magnetic resonance (NMR) structural models suggests similar interactions, though limited mutagenesis of Nun is available. Here, libraries of Nun's arginine-rich motif (ARM) were screened for the ability to exclude λ coinfection, and mutants were assayed for Nun termination with a boxB plasmid reporter system. Several Nun ARM residues appear to be immutable: Asp26, Arg28, Arg29, Arg32, Trp33, and Arg36. Asp26 and Trp33 appear to be unable to contact boxB and are not found at equivalent positions in λ N ARM. To understand if the requirement of Asp26, Trp33, and Arg36 indicated differences between HK022 Nun termination and λ N antitermination complexes, the same Nun libraries were fused to the activation domain of λ N and screened for clones able to complement N-deficient λ. Mutants were assayed for N antitermination. Surprisingly, Asp26 and Trp33 were still essential when Nun ARM was fused to N. Docking suggests that Nun ARM contacts a hydrophobic surface of the NusG carboxy-terminal domain containing residues necessary for Nun function. These findings indicate that Nun ARM relies on distinct contacts in its ternary complex and illustrate how protein-RNA recognition can evolve new regulatory functions.

IMPORTANCE

λ N protein interacts with host factors to allow λ nut-containing transcripts to elongate past termination signals. A competing bacteriophage, HK022, expresses Nun protein, which causes termination of λ nut transcripts. λ N and HK022 Nun use similar arginine-rich motifs (ARMs) to bind the same boxB RNAs in nut transcripts. Screening libraries of Nun ARM mutants, both in HK022 Nun and in a λ N fusion, revealed amino acids essential to Nun that could bind one or more host factors. Docking suggests that NusG, which is present in both Nun termination and N antitermination, is a plausible partner. These findings could help understand how transcription elongation is regulated and illustrate how subtle differences allow ARMs to evolve new regulatory functions.

Binding interface between the Salmonella σ(S)/RpoS subunit of RNA polymerase and Crl: hints from bacterial species lacking crl

In many Gram-negative bacteria, including Salmonella enterica serovar Typhimurium (S. Typhimurium), the sigma factor RpoS/σ^S accumulates during stationary phase of growth, and associates with the core RNA polymerase enzyme (E) to promote transcription initiation of genes involved in general stress resistance and starvation survival. Whereas σ factors are usually inactivated upon interaction with anti-σ proteins, σ^S binding to the Crl protein increases σ^S activity by favouring its association to E. Taking advantage of evolution of the σ^S sequence in bacterial species that do not contain a crl gene, like Pseudomonas aeruginosa, we identified and assigned a critical arginine residue in σ^S to the S. Typhimurium σ^S-Crl binding interface. We solved the solution structure of S. Typhimurium Crl by NMR and used it for NMR binding assays with σ^S and to generate in silico models of the σ^S-Crl complex constrained by mutational analysis. The σ^S-Crl models suggest that the identified arginine in σ^S interacts with an aspartate of Crl that is required for σ^S binding and is located inside a cavity enclosed by flexible loops, which also contribute to the interface. This study provides the basis for further structural investigation of the σ^S-Crl complex.

Structural neighboring property for identifying protein-protein binding sites

Background

The protein-protein interaction plays a key role in the control of many biological functions, such as drug design and functional analysis. Determination of binding sites is widely applied in molecular biology research. Therefore, many efficient methods have been developed for identifying binding sites. In this paper, we calculate structural neighboring property through Voronoi diagram. Using 6,438 complexes, we study local biases of structural neighboring property on interface.

Results

We propose a novel statistical method to extract interacting residues, and interacting patches can be clustered as predicted interface residues. In addition, structural neighboring property can be adopted to construct a new energy function, for evaluating docking solutions. It includes new statistical property as well as existing energy items. Comparing to existing methods, our approach improves overall F_nat value by at least 3%. On Benchmark v4.0, our method has average I_rmsd value of 3.31Å and overall F_nat value of 63%, which improves upon I_rmsd of 3.89 Å and F_nat of 49% for ZRANK, and I_rmsd of 3.99Å and F_nat of 46% for ClusPro. On the CAPRI targets, our method has average I_rmsd value of 3.46 Å and overall F_nat value of 45%, which improves upon I_rmsd of 4.18 Å and F_nat of 40% for ZRANK, and I_rmsd of 5.12 Å and F_nat of 32% for ClusPro.

Conclusions

Experiments show that our method achieves better results than some state-of-the-art methods for identifying protein-protein binding sites, with the prediction quality improved in terms of CAPRI evaluation criteria.