Influence of non-bonded parameters on the quality of NMR structures: a new force field for NMR structure calculation

Deciphering the structure and mechanism of action of computer-designed mastoparan peptides

Mastoparans are cationic peptides with multifunctional pharmacological properties. Mastoparan-R1 and mastoparan-R4 were computationally designed based on native mastoparan-L from wasps and have improved therapeutic potential for the control of bacterial infections. Here, we evaluated whether these peptides maintain their activity against Escherichia coli strains under a range of salt concentrations. We found that mastoparan-R1 and mastoparan-R4 preserved their activity under the conditions tested, including having antibacterial activities at physiological salt concentrations. The overall structure of the peptides was investigated using circular dichroism spectroscopy in a range of solvents. No significant changes in secondary structure were observed (random coil in aqueous solutions and α-helix in hydrophobic and anionic environments). The three-dimensional structures of mastoparan-R1 and mastoparan-R4 were elucidated through nuclear magnetic resonance spectroscopy, revealing amphipathic α-helical segments for Leu3-Ile13 (mastoparan-R1) and Leu3-Ile14 (mastoparan-R4). Possible membrane-association mechanisms for mastoparan-R1 and mastoparan-R4 were investigated through surface plasmon resonance and leakage studies with synthetic POPC and POPC/POPG (4:1) lipid bilayers. Mastoparan-L had the highest affinity for both membrane systems, whereas the two analogs had weaker association, but improved selectivity for lysing anionic membranes. This finding was also supported by molecular dynamics simulations, in which mastoparan-R1 and mastoparan-R4 were found to have greater interactions with bacteria-like membranes compared with model mammalian membranes. Despite having a few differences in their functional and structural profiles, the mastoparan-R1 analog stood out with the highest activity, greater bacteriostatic potential, and selectivity for lysing anionic membranes. This study reinforces the potential of mastoparan-R1 as a drug candidate.

An Electrostatically-steered Conformational Selection Mechanism Promotes SARS-CoV-2 Spike Protein Variation

After two years since the outbreak, the COVID-19 pandemic remains a global public health emergency. SARS-CoV-2 variants with substitutions on the spike (S) protein emerge increasing the risk of immune evasion and cross-species transmission. Here, we analyzed the evolution of the S protein as recorded in 276,712 samples collected before the start of vaccination efforts. Our analysis shows that most variants destabilize the S protein trimer, increase its conformational heterogeneity and improve the odds of the recognition by the host cell receptor. Most frequent substitutions promote overall hydrophobicity by replacing charged amino acids, reducing stabilizing local interactions in the unbound S protein trimer. Moreover, our results identify "forbidden" regions that rarely show any sequence variation, and which are related to conformational changes occurring upon fusion. These results are significant for understanding the structure and function of SARS-CoV-2 related proteins which is a critical step in vaccine development and epidemiological surveillance.

How to assess the structural dynamics of transcription factors by integrating sparse NMR and EPR constraints with molecular dynamics simulations

We review recent advances in modeling structural ensembles of transcription factors from nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopic data, integrated with molecular dynamics (MD) simulations. We focus on approaches that confirm computed conformational ensembles by sparse constraints obtained from magnetic resonance. This combination enables the deduction of functional and structural protein models even if nuclear Overhauser effects (NOEs) are too scarce for conventional structure determination. We highlight recent insights into the folding-upon-DNA binding transitions of intrinsically disordered transcription factors that could be assessed using such integrative approaches.

Adeno-Associated Virus (AAV-DJ)-Cryo-EM Structure at 1.56 Å Resolution

Adeno-associated virus is the leading viral vector for gene therapy. AAV-DJ is a recombinant variant developed for tropism to the liver. The AAV-DJ structure has been determined to 1.56 Å resolution through cryo-electron microscopy (cryo-EM). Only apoferritin is reported in preprints at 1.6 Å or higher resolution, and AAV-DJ nearly matches the highest resolutions ever attained through X-ray diffraction of virus crystals. However, cryo-EM has the advantage that most of the hydrogens are clear, improving the accuracy of atomic refinement, and removing ambiguity in hydrogen bond identification. Outside of secondary structures where hydrogen bonding was predictable a priori, the networks of hydrogen bonds coming from direct observation of hydrogens and acceptor atoms are quite different from those inferred even at 2.8 Å resolution. The implications for understanding viral assembly mean that cryo-EM will likely become the favored approach for high resolution structural virology.

Solution structure of the autophagy-related protein LC3C reveals a polyproline II motif on a mobile tether with phosphorylation site

(Macro-)autophagy is a compartmental degradation pathway conserved from yeast to mammals. The yeast protein Atg8 mediates membrane tethering/hemifusion and cargo recruitment and is essential for autophagy. The human MAP1LC3/GABARAP family proteins show high sequence identity with Atg8, but MAP1LC3C is distinguished by a conspicuous amino-terminal extension with unknown functional significance. We have determined the high-resolution three-dimensional structure and measured the backbone dynamics of MAP1LC3C by NMR spectroscopy. From Ser18 to Ala120, MAP1LC3C forms an α-helix followed by the ubiquitin-like tertiary fold with two hydrophobic binding pockets used by MAP1LC3/GABARAP proteins to recognize targets presenting LC3-interacting regions (LIRs). Unlike other MAP1LC3/GABARAP proteins, the amino-terminal region of MAP1LC3C does not form a stable helix α₁ but a "sticky arm" consisting of a polyproline II motif on a flexible linker. Ser18 at the interface between this linker and the structural core can be phosphorylated in vitro by protein kinase A, which causes additional conformational heterogeneity as monitored by NMR spectroscopy and molecular dynamics simulations, including changes in the LIR-binding interface. Based on these results we propose that the amino-terminal polyproline II motif mediates specific interactions with the microtubule cytoskeleton and that Ser18 phosphorylation modulates the interplay of MAP1LC3C with its various target proteins.

Periplasmic Expression of 4/7 α-Conotoxin TxIA Analogs in E. coli Favors Ribbon Isomer Formation - Suggestion of a Binding Mode at the α7 nAChR

Peptides derived from animal venoms provide important research tools for biochemical and pharmacological characterization of receptors, ion channels, and transporters. Some venom peptides have been developed into drugs (such as the synthetic ω-conotoxin MVIIA, ziconotide) and several are currently undergoing clinical trials for various clinical indications. Challenges in the development of peptides include their usually limited supply from natural sources, cost-intensive chemical synthesis, and potentially complicated stereoselective disulfide-bond formation in the case of disulfide-rich peptides. In particular, if extended structure–function analysis is performed or incorporation of stable isotopes for NMR studies is required, the comparatively low yields and high costs of synthesized peptides might constitute a limiting factor. Here we investigated the expression of the 4/7 α-conotoxin TxIA, a potent blocker at α3β2 and α7 nicotinic acetylcholine receptors (nAChRs), and three analogs in the form of maltose binding protein fusion proteins in Escherichia coli. Upon purification via nickel affinity chromatography and release of the toxins by protease cleavage, HPLC analysis revealed one major peak with the correct mass for all peptides. The final yield was 1–2 mg of recombinant peptide per liter of bacterial culture. Two-electrode voltage clamp analysis on oocyte-expressed nAChR subtypes demonstrated the functionality of these peptides but also revealed a 30 to 100-fold potency decrease of expressed TxIA compared to chemically synthesized TxIA. NMR spectroscopy analysis of TxIA and two of its analogs confirmed that the decreased activity was due to an alternative disulfide linkage rather than the missing C-terminal amidation, a post-translational modification that is common in α-conotoxins. All peptides preferentially formed in the ribbon conformation rather than the native globular conformation. Interestingly, in the case of the α7 nAChR, but not the α3β2 subtype, the loss of potency could be rescued by an R5D substitution. In conclusion, we demonstrate efficient expression of functional but alternatively folded ribbon TxIA variants in E. coli and provide the first structure–function analysis for a ribbon 4/7-α-conotoxin at α7 and α3β2 nAChRs. Computational analysis based on these data provide evidence for a ribbon α-conotoxin binding mode that might be exploited to design ligands with optimized selectivity.

Structural Architecture of the Nucleosome Remodeler ISWI Determined from Cross-Linking, Mass Spectrometry, SAXS, and Modeling

Chromatin remodeling factors assume critical roles by regulating access to nucleosomal DNA. To determine the architecture of the Drosophila ISWI remodeling enzyme, we developed an integrative structural approach that combines protein cross-linking, mass spectrometry, small-angle X-ray scattering, and computational modeling. The resulting structural model shows the ATPase module in a resting state with both ATPase lobes twisted against each other, providing support for a conformation that was recently trapped by crystallography. The autoinhibiting NegC region does not protrude from the ATPase module as suggested previously. The regulatory NTR domain is located near both ATPase lobes. The full-length enzyme is flexible and can adopt a compact structure in solution with the C-terminal HSS domain packing against the ATPase module. Our data imply a series of conformational changes upon activation of the enzyme and illustrate how the NTR, NegC, and HSS domains contribute to regulation of the ATPase module.

Structural dissection of an interaction between transcription initiation and termination factors implicated in promoter-terminator cross-talk

Functional cross-talk between the promoter and terminator of a gene has long been noted. Promoters and terminators are juxtaposed to form gene loops in several organisms, and gene looping is thought to be involved in transcriptional regulation. The general transcription factor IIB (TFIIB) and the C-terminal domain phosphatase Ssu72, essential factors of the transcription preinitiation complex and the mRNA processing and polyadenylation complex, respectively, are important for gene loop formation. TFIIB and Ssu72 interact both genetically and physically, but the molecular basis of this interaction is not known. Here we present a crystal structure of the core domain of TFIIB in two new conformations that differ in the relative distance and orientation of the two cyclin-like domains. The observed extraordinary conformational plasticity may underlie the binding of TFIIB to multiple transcription factors and promoter DNAs that occurs in distinct stages of transcription, including initiation, reinitiation, and gene looping. We mapped the binding interface of the TFIIB-Ssu72 complex using a series of systematic, structure-guided in vitro binding and site-specific photocross-linking assays. Our results indicate that Ssu72 competes with acidic activators for TFIIB binding and that Ssu72 disrupts an intramolecular TFIIB complex known to impede transcription initiation. We also show that the TFIIB-binding site on Ssu72 overlaps with the binding site of symplekin, a component of the mRNA processing and polyadenylation complex. We propose a hand-off model in which Ssu72 mediates a conformational transition in TFIIB, accounting for the role of Ssu72 in transcription reinitiation, gene looping, and promoter-terminator cross-talk.

Identification of a Ligand-Binding Site on the Staphylococcus aureus DnaG Primase C-Terminal Domain

The interface between the DnaG primase C-terminal domain (CTD) and the N-terminal domain of DnaB helicase is essential for bacterial DNA replication because it allows coordinated priming of DNA synthesis at the replication fork while the DNA is being unwound. Because these two proteins are conserved in all bacteria and distinct from those in eukaryotes, their interface is an attractive antibiotic target. To learn more about this interface, we determined the solution structure and dynamics of the DnaG primase CTD from Staphylococcus aureus, a medically important bacterial species. Comparison with the known primase CTD structures shows there are two biologically relevant conformations, an open conformation that likely binds to DnaB helicase and a closed conformation that does not. The S. aureus primase CTD is in the closed conformation, but nuclear magnetic resonance (NMR) dynamic studies indicate there is considerable movement in the linker between the two subdomains and that N564 is the most dynamic residue within the linker. A high-throughput NMR ligand affinity screen identified potential binding compounds, among which were acycloguanosine and myricetin. Although the affinity for these compounds and adenosine was in the millimolar range, all three bind to a common pocket that is present only on the closed conformation of the CTD. This binding pocket is at the opposite end of helices 6 and 7 from N564, the key hinge residue. The identification of this binding pocket should allow the development of stronger-binding ligands that can prevent formation of the CTD open conformation that binds to DnaB helicase.

High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Structure determination of proteins by NMR is unique in its ability to measure restraints, very accurately, in environments and under conditions that closely mimic those encountered in vivo. For example, advances in solid-state NMR methods enable structure determination of membrane proteins in detergent-free lipid bilayers, and of large soluble proteins prepared by sedimentation, while parallel advances in solution NMR methods and optimization of detergent-free lipid nanodiscs are rapidly pushing the envelope of the size limit for both soluble and membrane proteins. These experimental advantages, however, are partially squandered during structure calculation, because the commonly used force fields are purely repulsive and neglect solvation, Van der Waals forces and electrostatic energy. Here we describe a new force field, and updated energy functions, for protein structure calculations with EEFx implicit solvation, electrostatics, and Van der Waals Lennard-Jones forces, in the widely used program Xplor-NIH. The new force field is based primarily on CHARMM22, facilitating calculations with a wider range of biomolecules. The new EEFx energy function has been rewritten to enable OpenMP parallelism, and optimized to enhance computation efficiency. It implements solvation, electrostatics, and Van der Waals energy terms together, thus ensuring more consistent and efficient computation of the complete nonbonded energy lists. Updates in the related python module allow detailed analysis of the interaction energies and associated parameters. The new force field and energy function work with both soluble proteins and membrane proteins, including those with cofactors or engineered tags, and are very effective in situations where there are sparse experimental restraints. Results obtained for NMR-restrained calculations with a set of five soluble proteins and five membrane proteins show that structures calculated with EEFx have significant improvements in accuracy, precision, and conformation, and that structure refinement can be obtained by short relaxation with EEFx to obtain improvements in these key metrics. These developments broaden the range of biomolecular structures that can be calculated with high fidelity from NMR restraints.

Elucidating Mechanisms of Molecular Recognition Between Human Argonaute and miRNA Using Computational Approaches

MicroRNA (miRNA) and Argonaute (AGO) protein together form the RNA-induced silencing complex (RISC) that plays an essential role in the regulation of gene expression. Elucidating the underlying mechanism of AGO-miRNA recognition is thus of great importance not only for the in-depth understanding of miRNA function but also for inspiring new drugs targeting miRNAs. In this chapter we introduce a combined computational approach of molecular dynamics (MD) simulations, Markov state models (MSMs), and protein-RNA docking to investigate AGO-miRNA recognition. Constructed from MD simulations, MSMs can elucidate the conformational dynamics of AGO at biologically relevant timescales. Protein-RNA docking can then efficiently identify the AGO conformations that are geometrically accessible to miRNA. Using our recent work on human AGO2 as an example, we explain the rationale and the workflow of our method in details. This combined approach holds great promise to complement experiments in unraveling the mechanisms of molecular recognition between large, flexible, and complex biomolecules.

Protein-protein and peptide-protein docking and refinement using ATTRACT in CAPRI

The ATTRACT coarse-grained docking approach in combination with various types of atomistic, flexible refinement methods has been applied to predict protein-protein and peptide-protein complexes in CAPRI rounds 28-36. For a large fraction of CAPRI targets (12 out of 18), at least one model of acceptable or better quality was generated, corresponding to a success rate of 67%. In particular, for several peptide-protein complexes excellent predictions were achieved. In several cases, a combination of template-based modeling and extensive molecular dynamics-based refinement yielded medium and even high quality solutions. In one particularly challenging case, the structure of an ubiquitylation enzyme bound to the nucleosome was correctly predicted as a set of acceptable quality solutions. Based on the experience with the CAPRI targets, new interface refinement approaches and methods for ab-initio peptide-protein docking have been developed. Failures and possible improvements of the docking method with respect to scoring and protein flexibility will also be discussed. Proteins 2017; 85:391-398. © 2016 Wiley Periodicals, Inc.

A protocol for the refinement of NMR structures using simultaneously pseudocontact shift restraints from multiple lanthanide ions

The binding of paramagnetic metal ions to proteins produces a number of different effects on the NMR spectra of the system. In particular, when the magnetic susceptibility of the metal ion is anisotropic, pseudocontact shifts (PCSs) arise and can be easily measured. They constitute very useful restraints for the solution structure determination of metal-binding proteins. In this context, there has been great interest in the use of lanthanide(III) ions to induce PCSs in diamagnetic proteins, e.g. through the replacement native calcium(II) ions. By preparing multiple samples in each of which a different ion of the lanthanide series is introduced, it is possible to obtain multiple independent PCS datasets that can be used synergistically to generate protein structure ensembles (typically called bundles). For typical NMR-based determination of protein structure, it is necessary to perform an energetic refinement of such initial bundles to obtain final structures whose geometric quality is suitable for deposition in the PDB. This can be conveniently done by using restrained molecular dynamics simulations (rMD) in explicit solvent. However, there are no available protocols for rMD using multiple PCS datasets as part of the restraints. In this work, we extended the PCS module of the AMBER MD package to handle multiple datasets and tuned a previously developed protocol for NMR structure refinement to achieve consistent convergence with PCS restraints. Test calculations with real experimental data show that this new implementation delivers the expected improvement of protein geometry, resulting in final structures that are of suitable quality for deposition. Furthermore, we observe that also initial structures generated only with traditional restraints can be successfully refined using traditional and PCS restraints simultaneously.

Development of a μO-Conotoxin Analogue with Improved Lipid Membrane Interactions and Potency for the Analgesic Sodium Channel NaV1.8

The μO-conotoxins MrVIA, MrVIB, and MfVIA inhibit the voltage-gated sodium channel NaV1.8, a well described target for the treatment of pain; however, little is known about the residues or structural elements that define this activity. In this study, we determined the three-dimensional structure of MfVIA, examined its membrane binding properties, performed alanine-scanning mutagenesis, and identified residues important for its activity at human NaV1.8. A second round of mutations resulted in (E5K,E8K)MfVIA, a double mutant with greater positive surface charge and greater affinity for lipid membranes compared with MfVIA. This analogue had increased potency at NaV1.8 and was analgesic in the mouse formalin assay.

Structure determination of helical filaments by solid-state NMR spectroscopy

The controlled formation of filamentous protein complexes plays a crucial role in many biological systems and represents an emerging paradigm in signal transduction. The mitochondrial antiviral signaling protein (MAVS) is a central signal transduction hub in innate immunity that is activated by a receptor-induced conversion into helical superstructures (filaments) assembled from its globular caspase activation and recruitment domain. Solid-state NMR (ssNMR) spectroscopy has become one of the most powerful techniques for atomic resolution structures of protein fibrils. However, for helical filaments, the determination of the correct symmetry parameters has remained a significant hurdle for any structural technique and could thus far not be precisely derived from ssNMR data. Here, we solved the atomic resolution structure of helical MAVS(CARD) filaments exclusively from ssNMR data. We present a generally applicable approach that systematically explores the helical symmetry space by efficient modeling of the helical structure restrained by interprotomer ssNMR distance restraints. Together with classical automated NMR structure calculation, this allowed us to faithfully determine the symmetry that defines the entire assembly. To validate our structure, we probed the protomer arrangement by solvent paramagnetic resonance enhancement, analysis of chemical shift differences relative to the solution NMR structure of the monomer, and mutagenesis. We provide detailed information on the atomic contacts that determine filament stability and describe mechanistic details on the formation of signaling-competent MAVS filaments from inactive monomers.

Structural investigations of the p53/p73 homologs from the tunicate species Ciona intestinalis reveal the sequence requirements for the formation of a tetramerization domain

Most members of the p53 family of transcription factors form tetramers. Responsible for determining the oligomeric state is a short oligomerization domain consisting of one β-strand and one α-helix. With the exception of human p53 all other family members investigated so far contain a second α-helix as part of their tetramerization domain. Here we have used nuclear magnetic resonance spectroscopy to characterize the oligomerization domains of the two p53-like proteins from the tunicate Ciona intestinalis, representing the closest living relative of vertebrates. Structure determination reveals for one of the two proteins a new type of packing of this second α-helix on the core domain that was not predicted based on the sequence, while the other protein does not form a second helix despite the presence of crucial residues that are conserved in all other family members that form a second helix. By mutational analysis, we identify a proline as well as large hydrophobic residues in the hinge region between both helices as the crucial determinant for the formation of a second helix.

Solution NMR characterization of chemokine CXCL8/IL-8 monomer and dimer binding to glycosaminoglycans: structural plasticity mediates differential binding interactions

Structural plasticity plays a major role in determining differential binding of CXCL8 monomer and dimer to glycosaminoglycans (GAGs) and that dimer is the high-affinity GAG ligand. We propose that these properties play important roles in orchestrating in vivo chemokine-mediated neutrophil function.

Chemokine CXCL8/interleukin-8 (IL-8) plays a crucial role in directing neutrophils and oligodendrocytes to combat infection/injury and tumour cells in metastasis development. CXCL8 exists as monomers and dimers and interaction of both forms with glycosaminoglycans (GAGs) mediate these diverse cellular processes. However, very little is known regarding the structural basis underlying CXCL8–GAG interactions. There are conflicting reports on the affinities, geometry and whether the monomer or dimer is the high-affinity GAG ligand. To resolve these issues, we characterized the binding of a series of heparin-derived oligosaccharides [heparin disaccharide (dp2), heparin tetrasaccharide (dp4), heparin octasaccharide (dp8) and heparin 14-mer (dp14)] to the wild-type (WT) dimer and a designed monomer using solution NMR spectroscopy. The pattern and extent of binding-induced chemical shift perturbation (CSP) varied between dimer and monomer and between longer and shorter oligosaccharides. NMR-based structural models show that different interaction modes coexist and that the nature of interactions varied between monomer and dimer and oligosaccharide length. MD simulations indicate that the binding interface is structurally plastic and provided residue-specific details of the dynamic nature of the binding interface. Binding studies carried out under conditions at which WT CXCL8 exists as monomers and dimers provide unambiguous evidence that the dimer is the high-affinity GAG ligand. Together, our data indicate that a set of core residues function as the major recognition/binding site, a set of peripheral residues define the various binding geometries and that the structural plasticity of the binding interface allows multiplicity of binding interactions. We conclude that structural plasticity most probably regulates in vivo CXCL8 monomer/dimer–GAG interactions and function.

A Practical Implicit Membrane Potential for NMR Structure Calculations of Membrane Proteins

The highly anisotropic environment of the lipid bilayer membrane imposes significant constraints on the structures and functions of membrane proteins. However, NMR structure calculations typically use a simple repulsive potential that neglects the effects of solvation and electrostatics, because explicit atomic representation of the solvent and lipid molecules is computationally expensive and impractical for routine NMR-restrained calculations that start from completely extended polypeptide templates. Here, we describe the extension of a previously described implicit solvation potential, eefxPot, to include a membrane model for NMR-restrained calculations of membrane protein structures in XPLOR-NIH. The key components of eefxPot are an energy term for solvation free energy that works together with other nonbonded energy functions, a dedicated force field for conformational and nonbonded protein interaction parameters, and a membrane function that modulates the solvation free energy and dielectric screening as a function of the atomic distance from the membrane center, relative to the membrane thickness. Initial results obtained for membrane proteins with structures determined experimentally in lipid bilayer membranes show that eefxPot affords significant improvements in structural quality, accuracy, and precision. Calculations with eefxPot are straightforward to implement and can be used to both fold and refine structures, as well as to run unrestrained molecular-dynamics simulations. The potential is entirely compatible with the full range of experimental restraints measured by various techniques. Overall, it provides a useful and practical way to calculate membrane protein structures in a physically realistic environment.

Improved reliability, accuracy and quality in automated NMR structure calculation with ARIA

In biological NMR, assignment of NOE cross-peaks and calculation of atomic conformations are critical steps in the determination of reliable high-resolution structures. ARIA is an automated approach that performs NOE assignment and structure calculation in a concomitant manner in an iterative procedure. The log-harmonic shape for distance restraint potential and the Bayesian weighting of distance restraints, recently introduced in ARIA, were shown to significantly improve the quality and the accuracy of determined structures. In this paper, we propose two modifications of the ARIA protocol: (1) the softening of the force field together with adapted hydrogen radii, which is meaningful in the context of the log-harmonic potential with Bayesian weighting, (2) a procedure that automatically adjusts the violation tolerance used in the selection of active restraints, based on the fitting of the structure to the input data sets. The new ARIA protocols were fine-tuned on a set of eight protein targets from the CASD-NMR initiative. As a result, the convergence problems previously observed for some targets was resolved and the obtained structures exhibited better quality. In addition, the new ARIA protocols were applied for the structure calculation of ten new CASD-NMR targets in a blind fashion, i.e. without knowing the actual solution. Even though optimisation of parameters and pre-filtering of unrefined NOE peak lists were necessary for half of the targets, ARIA consistently and reliably determined very precise and highly accurate structures for all cases. In the context of integrative structural biology, an increasing number of experimental methods are used that produce distance data for the determination of 3D structures of macromolecules, stressing the importance of methods that successfully make use of ambiguous and noisy distance data.

Markov State Models Reveal a Two-Step Mechanism of miRNA Loading into the Human Argonaute Protein: Selective Binding followed by Structural Re-arrangement

Argonaute (Ago) proteins and microRNAs (miRNAs) are central components in RNA interference, which is a key cellular mechanism for sequence-specific gene silencing. Despite intensive studies, molecular mechanisms of how Ago recognizes miRNA remain largely elusive. In this study, we propose a two-step mechanism for this molecular recognition: selective binding followed by structural re-arrangement. Our model is based on the results of a combination of Markov State Models (MSMs), large-scale protein-RNA docking, and molecular dynamics (MD) simulations. Using MSMs, we identify an open state of apo human Ago-2 in fast equilibrium with partially open and closed states. Conformations in this open state are distinguished by their largely exposed binding grooves that can geometrically accommodate miRNA as indicated in our protein-RNA docking studies. miRNA may then selectively bind to these open conformations. Upon the initial binding, the complex may perform further structural re-arrangement as shown in our MD simulations and eventually reach the stable binary complex structure. Our results provide novel insights in Ago-miRNA recognition mechanisms and our methodology holds great potential to be widely applied in the studies of other important molecular recognition systems.

In RNA interference, Argonaute proteins and microRNAs together form the functional core that regulates the gene expression with high sequence specificity. Elucidating the detailed mechanism of molecular recognition between Argonaute proteins and microRNAs is thus important not only for the fundamental understanding of RNA interference, but also for the further development of microRNA-based therapeutic application. In this work, we propose a two-step model to understand the mechanism of microRNA loading into human Argonaute-2: selective binding followed by structural re-arrangement. Our model is based on the results from a combined approach of molecular dynamics simulations, Markov State Models and protein-RNA docking. In particular, we identify a metastable open state of apo hAgo2 in rapid equilibrium with other states. Some of conformations in this open state have largely exposed RNA binding groove that can accommodate microRNA. We further show that the initial Argonaute-microRNA binding complex undergoes structural re-arrangement to reach stable binary crystal structure. These results provide novel insights into the underlying mechanism of Argonaute-microRNA recognition. In addition, our method is readily applicable to the investigation of other complex molecular recognition events such as protein-protein interactions and protein-ligand binding.

Solution-state NMR structure of the putative morphogene protein BolA (PFE0790c) from Plasmodium falciparum

Protozoa of the genus Plasmodium are responsible for malaria, which is perhaps the most important parasitic disease to infect mankind. The emergence of Plasmodium strains resistant to current therapeutics and prophylactics makes the development of new treatment strategies urgent. Among the potential targets for new antimalarial drugs is the BolA-like protein PFE0790c from Plasmodium falciparum (Pf-BolA). While the function of BolA is unknown, it has been linked to cell morphology by regulating transcription in response to stress. Using an NMR-based method, an ensemble of 20 structures of Pf-BolA was determined and deposited in the PDB (PDB entry 2kdn). The overall topology of the Pf-BolA structure, α1-β1-β2-η1-α2/η2-β3-α3, with the β-strands forming a mixed β-sheet, is similar to the fold observed in other BolA structures. A helix-turn-helix motif similar to the class II KH fold associated with nucleic acid-binding proteins is present, but contains an FXGXXXL signature sequence that differs from the GXXG signature sequence present in class II KH folds, suggesting that the BolA family of proteins may use a novel protein-nucleic acid interface. A well conserved arginine residue, Arg50, hypothesized to play a role in governing the formation of the C-terminal α-helix in the BolA family of proteins, is too distant to form polar contacts with any side chains in this α-helix in Pf-BolA, suggesting that this conserved arginine may only serve a role in guiding the orientation of this C-terminal helix in some BolA proteins. A survey of BolA structures suggests that the C-terminal helix may not have a functional role and that the third helix (α2/η2) has a `kink' that appears to be conserved among the BolA protein structures. Circular dichroism spectroscopy shows that Pf-BolA is fairly robust, partially unfolding when heated to 353 K and refolding upon cooling to 298 K.

iATTRACT: simultaneous global and local interface optimization for protein-protein docking refinement

Protein-protein interactions are abundant in the cell but to date structural data for a large number of complexes is lacking. Computational docking methods can complement experiments by providing structural models of complexes based on structures of the individual partners. A major caveat for docking success is accounting for protein flexibility. Especially, interface residues undergo significant conformational changes upon binding. This limits the performance of docking methods that keep partner structures rigid or allow limited flexibility. A new docking refinement approach, iATTRACT, has been developed which combines simultaneous full interface flexibility and rigid body optimizations during docking energy minimization. It employs an atomistic molecular mechanics force field for intermolecular interface interactions and a structure-based force field for intramolecular contributions. The approach was systematically evaluated on a large protein-protein docking benchmark, starting from an enriched decoy set of rigidly docked protein-protein complexes deviating by up to 15 Å from the native structure at the interface. Large improvements in sampling and slight but significant improvements in scoring/discrimination of near native docking solutions were observed. Complexes with initial deviations at the interface of up to 5.5 Å were refined to significantly better agreement with the native structure. Improvements in the fraction of native contacts were especially favorable, yielding increases of up to 70%.

Bayesian weighting of statistical potentials in NMR structure calculation

The use of statistical potentials in NMR structure calculation improves the accuracy of the final structure but also raises issues of double counting and possible bias. Because statistical potentials are averaged over a large set of structures, they may not reflect the preferences of a particular structure or data set. We propose a Bayesian method to incorporate a knowledge-based backbone dihedral angle potential into an NMR structure calculation. To avoid bias exerted through the backbone potential, we adjust its weight by inferring it from the experimental data. We demonstrate that an optimally weighted potential leads to an improvement in the accuracy and quality of the final structure, especially with sparse and noisy data. Our findings suggest that no universally optimal weight exists, and that the weight should be determined based on the experimental data. Other knowledge-based potentials can be incorporated using the same approach.

A practical implicit solvent potential for NMR structure calculation

The benefits of protein structure refinement in water are well documented. However, performing structure refinement with explicit atomic representation of the solvent molecules is computationally expensive and impractical for NMR-restrained structure calculations that start from completely extended polypeptide templates. Here we describe a new implicit solvation potential, EEFx (Effective Energy Function for XPLOR-NIH), for NMR-restrained structure calculations of proteins in XPLOR-NIH. The key components of EEFx are an energy term for solvation energy that works together with other nonbonded energy functions, and a dedicated force field for conformational and nonbonded protein interaction parameters. The initial results obtained with EEFx show that significant improvements in structural quality can be obtained. EEFx is computationally efficient and can be used both to fold and refine structures. Overall, EEFx improves the quality of protein conformation and nonbonded atomic interactions. Moreover, such benefits are accompanied by enhanced structural precision and enhanced structural accuracy, reflected in improved agreement with the cross-validated dipolar coupling data. Finally, implementation of EEFx calculations is straightforward and computationally efficient. Overall, EEFx provides a useful method for the practical calculation of experimental protein structures in a physically realistic environment.

Semienzymatic cyclization of disulfide-rich peptides using Sortase A

Disulfide-rich cyclic peptides have generated great interest in the development of peptide-based therapeutics due to their exceptional stability toward chemical, enzymatic, or thermal attack. In particular, they have been used as scaffolds onto which bioactive epitopes can be grafted to take advantage of the favorable biophysical properties of disulfide-rich cyclic peptides. To date, the most commonly used method for the head-to-tail cyclization of peptides has been native chemical ligation. In recent years, however, enzyme-mediated cyclization has become a promising new technology due to its efficiency, safety, and cost-effectiveness. Sortase A (SrtA) is a bacterial enzyme with transpeptidase activity. It recognizes a C-terminal penta-amino acid motif, LPXTG, and cleaves the amide bond between Thr and Gly to form a thioacyl-linked intermediate. This intermediate undergoes nucleophilic attack by an N-terminal poly-Gly sequence to form an amide bond between the Thr and N-terminal Gly. Here, we demonstrate that sortase A can successfully be used to cyclize a variety of small disulfide-rich peptides, including the cyclotide kalata B1, α-conotoxin Vc1.1, and sunflower trypsin inhibitor 1. These peptides range in size from 14 to 29 amino acids and contain three, two, or one disulfide bond, respectively, within their head-to-tail cyclic backbones. Our findings provide proof of concept for the potential broad applicability of enzymatic cyclization of disulfide-rich peptides with therapeutic potential.

The structure of Plasmodium yoelii merozoite surface protein 119, antibody specificity and implications for malaria vaccine design

Merozoite surface protein 1 (MSP1) has been identified as a target antigen for protective immune responses against asexual blood stage malaria, but effective vaccines based on MSP1 have not been developed so far. We have modified the sequence of Plasmodium yoelii MSP119 (the C-terminal region of the molecule) and examined the ability of the variant proteins to bind protective monoclonal antibodies and to induce protection by immunization. In parallel, we examined the structure of the protein and the consequences of the amino acid changes. Naturally occurring sequence polymorphisms reduced the binding of individual protective antibodies, indicating that they contribute to immune evasion, but immunization with these variant proteins still provided protective immunity. One variant that resulted in the localized distortion of a loop close to the N-terminus of MSP119 almost completely ablated protection by immunization, indicating the importance of this region of MSP119 as a target for protective immunity and in vaccine development.

A novel α4/7-conotoxin LvIA from Conus lividus that selectively blocks α3β2 vs. α6/α3β2β3 nicotinic acetylcholine receptors

This study was performed to discover and characterize the first potent α3β2-subtype-selective nicotinic acetylcholine receptor (nAChR) ligand. A novel α4/7-conotoxin, α-CTxLvIA, was cloned from Conus lividus. Its pharmacological profile at Xenopus laevis oocyte-expressed rat nAChR subtypes was determined by 2-electrode voltage-clamp electrophysiology, and its 3-dimensional (3D) structure was determined by NMR spectroscopy. α-CTx LvIA is a 16-aa C-terminally-amidated peptide with 2-disulfide bridges. Using rat subunits expressed in Xenopus oocytes, we found the highest affinity of α-CTxLvIA was for α3β2 nAChRs (IC50 8.7 nM), where blockade was reversible within 2 min. IC50 values were >100 nM at α6/α3β2β3, α6/α3β4, and α3β4 nAChRs, and ≥3 μM at all other subtypes tested. α3β2 vs. α6β2 subtype selectivity was confirmed for human-subunit nAChRs with much greater preference (300-fold) for α3β2 over α6β2 nAChRs. This is the first α-CTx reported to show high selectivity for human α3β2 vs. α6β2 nAChRs. α-CTxLvIA adopts two similarly populated conformations water: one (assumed to be bioactive) is highly structured, whereas the other is mostly random coil in nature. Selectivity differences with the similarly potent, but less selective, α3β2 nAChR antagonist α-CTx PeIA probably reside within the three residues, which differ in loop 2, given their otherwise similar 3D structures. α4/7-CTx LvIA is a new, potent, selective α3β2 nAChR antagonist, which will enable detailed studies of α3β2 nAChR structure, function, and physiological roles.

Structural and functional studies of a phosphatidic acid-binding antifungal plant defensin MtDef4: identification of an RGFRRR motif governing fungal cell entry

MtDef4 is a 47-amino acid cysteine-rich evolutionary conserved defensin from a model legume Medicago truncatula. It is an apoplast-localized plant defense protein that inhibits the growth of the ascomycetous fungal pathogen Fusarium graminearum in vitro at micromolar concentrations. Little is known about the mechanisms by which MtDef4 mediates its antifungal activity. In this study, we show that MtDef4 rapidly permeabilizes fungal plasma membrane and is internalized by the fungal cells where it accumulates in the cytoplasm. Furthermore, analysis of the structure of MtDef4 reveals the presence of a positively charged γ-core motif composed of β₂ and β₃ strands connected by a positively charged RGFRRR loop. Replacement of the RGFRRR sequence with AAAARR or RGFRAA abolishes the ability of MtDef4 to enter fungal cells, suggesting that the RGFRRR loop is a translocation signal required for the internalization of the protein. MtDef4 binds to phosphatidic acid (PA), a precursor for the biosynthesis of membrane phospholipids and a signaling lipid known to recruit cytosolic proteins to membranes. Amino acid substitutions in the RGFRRR sequence which abolish the ability of MtDef4 to enter fungal cells also impair its ability to bind PA. These findings suggest that MtDef4 is a novel antifungal plant defensin capable of entering into fungal cells and affecting intracellular targets and that these processes are mediated by the highly conserved cationic RGFRRR loop via its interaction with PA.

Ionic strength-dependent conformations of a ubiquitin-like small archaeal modifier protein (SAMP2) from Haloferax volcanii

Ubiquitin-like proteins play important roles in diverse biological processes. In this study, we present an unexpected finding that a ubiquitin-like small archaeal modifier protein (SAMP2) from Haloferax volcanii adopts two distinct states under low ionic condition. One of these is similar to the β-grasp structure conserved in ubiquitin-like proteins from eukaryotes; the other is disordered, like prokaryotic ubiquitin-like protein, Pup. Furthermore, our study reveals that the conformation of SAMP2 is dependent on ionic strength. With the increase of ion concentration, SAMP2 undergoes a conformational conversion from disorder to order, indicating that the ordered conformation is the functional form of SAMP2 under the physiological condition of H. volcanii.

A winged helix domain in human MUS81 binds DNA and modulates the endonuclease activity of MUS81 complexes

The MUS81-EME1 endonuclease maintains metazoan genomic integrity by cleaving branched DNA structures that arise during the resolution of recombination intermediates. In humans, MUS81 also forms a poorly characterized complex with EME2. Here, we identify and determine the structure of a winged helix (WH) domain from human MUS81, which binds DNA. WH domain mutations greatly reduce binding of the isolated domain to DNA and impact on incision activity of MUS81-EME1/EME2 complexes. Deletion of the WH domain reduces the endonuclease activity of both MUS81-EME1 and MUS81-EME2 complexes, and incisions made by MUS81-EME2 are made closer to the junction on substrates containing a downstream duplex, such as fork structures and nicked Holliday junctions. WH domain mutation or deletion in Schizosaccharomyces pombe phenocopies the DNA-damage sensitivity of strains deleted for mus81. Our results indicate an important role for the WH domain in both yeast and human MUS81 complexes.

Dicarba α-conotoxin Vc1.1 analogues with differential selectivity for nicotinic acetylcholine and GABAB receptors

Conotoxins have emerged as useful leads for the development of novel therapeutic analgesics. These peptides, isolated from marine molluscs of the genus Conus, have evolved exquisite selectivity for receptors and ion channels of excitable tissue. One such peptide, α-conotoxin Vc1.1, is a 16-mer possessing an interlocked disulfide framework. Despite its emergence as a potent analgesic lead, the molecular target and mechanism of action of Vc1.1 have not been elucidated to date. In this paper we describe the regioselective synthesis of dicarba analogues of Vc1.1 using olefin metathesis. The ability of these peptides to inhibit acetylcholine-evoked current at rat α9α10 and α3β4 nicotinic acetylcholine receptors (nAChR) expressed in Xenopus oocytes has been assessed in addition to their ability to inhibit high voltage-activated (HVA) calcium channel current in isolated rat DRG neurons. Their solution structures were determined by NMR spectroscopy. Significantly, we have found that regioselective replacement of the native cystine framework with a dicarba bridge can be used to selectively tune the cyclic peptide's innate biological activity for one receptor over another. The 2,8-dicarba Vc1.1 isomer retains activity at γ-aminobutyric acid (GABAB) G protein-coupled receptors, whereas the isomeric 3,16-dicarba Vc1.1 peptide retains activity at the α9α10 nAChR subtype. These singularly acting analogues will enable the elucidation of the biological target responsible for the peptide's potent analgesic activity.

Ionic strength-dependent conformations of a ubiquitin-like small archaeal modifier protein (SAMP1) from Haloferax volcanii

Eukaryotic ubiquitin and ubiquitin-like systems play crucial roles in various cellular biological processes. In this work, we determined the solution structure of SAMP1 from Haloferax volcanii by NMR spectroscopy. Under low ionic conditions, SAMP1 presented two distinct conformations, one folded β-grasp and the other disordered. Interestingly, SAMP1 underwent a conformational conversion from disorder to order with ion concentration increasing, indicating that the ordered conformation is the functional form of SAMP1 under the physiological condition of H. volcanii. Furthermore, SAMP1 could interact with proteasome-activating nucleotidase B, supposing a potential role of SAMP1 in the protein degradation pathway mediated by proteasome.

Conservation of dark recovery kinetic parameters and structural features in the pseudomonadaceae "short" light, oxygen, voltage (LOV) protein family: implications for the design of LOV-based optogenetic tools

In bacteria and fungi, various light, oxygen, voltage (LOV) sensory systems that lack a fused effector domain but instead contain only short N- and C-terminal extensions flanking the LOV core exist. In the prokaryotic kingdom, this so-called "short" LOV protein family represents the third largest LOV photoreceptor family. This observation prompted us to study their distribution and phylogeny as well as their photochemical and structural properties in more detail. We recently described the slow and fast reverting "short" LOV proteins PpSB1-LOV and PpSB2-LOV from Pseudomonas putida KT2440 whose adduct state lifetimes varied by 3 orders of magnitude [Jentzsch, K., Wirtz, A., Circolone, F., Drepper, T., Losi, A., Gärtner, W., Jaeger, K. E., and Krauss, U. (2009) Biochemistry 48, 10321-10333]. We now present evidence of the conservation of similar fast and slow-reverting "short" LOV proteins in different Pseudomonas species. Truncation studies conducted with PpSB1-LOV and PpSB2-LOV suggested that the short N- and C-terminal extensions outside of the LOV core domain are essential for the structural integrity and folding of the two proteins. While circular dichroism and solution nuclear magnetic resonance experiments verify that the two short C-terminal extensions of PpSB1-LOV and PpSB2-LOV form independently folding helical structures in solution, bioinformatic analyses imply the formation of coiled coils of the respective structural elements in the context of the dimeric full-length proteins. Given their prototypic architecture, conserved in most more complex LOV photoreceptor systems, "short" LOV proteins could represent ideally suited building blocks for the design of genetically encoded photoswitches (i.e., LOV-based optogenetic tools).

The structure of latherin, a surfactant allergen protein from horse sweat and saliva

Latherin is a highly surface-active allergen protein found in the sweat and saliva of horses and other equids. Its surfactant activity is intrinsic to the protein in its native form, and is manifest without associated lipids or glycosylation. Latherin probably functions as a wetting agent in evaporative cooling in horses, but it may also assist in mastication of fibrous food as well as inhibition of microbial biofilms. It is a member of the PLUNC family of proteins abundant in the oral cavity and saliva of mammals, one of which has also been shown to be a surfactant and capable of disrupting microbial biofilms. How these proteins work as surfactants while remaining soluble and cell membrane-compatible is not known. Nor have their structures previously been reported. We have used protein nuclear magnetic resonance spectroscopy to determine the conformation and dynamics of latherin in aqueous solution. The protein is a monomer in solution with a slightly curved cylindrical structure exhibiting a ‘super-roll’ motif comprising a four-stranded anti-parallel β-sheet and two opposing α-helices which twist along the long axis of the cylinder. One end of the molecule has prominent, flexible loops that contain a number of apolar amino acid side chains. This, together with previous biophysical observations, leads us to a plausible mechanism for surfactant activity in which the molecule is first localized to the non-polar interface via these loops, and then unfolds and flattens to expose its hydrophobic interior to the air or non-polar surface. Intrinsically surface-active proteins are relatively rare in nature, and this is the first structure of such a protein from mammals to be reported. Both its conformation and proposed method of action are different from other, non-mammalian surfactant proteins investigated so far.

Structure and dynamics of human Nedd4-1 WW3 in complex with the αENaC PY motif

Nedd4-1 (neuronal precursor cell expressed developmentally downregulated gene 4-1) is an E3 ubiquitin ligase that interacts with and negatively regulates the epithelial Na(+) channel (ENaC). The WW domains of Nedd4-1 bind to the ENaC subunits via recognition of PY motifs. Human Nedd4-1 (hNedd4-1) contains four WW domains with the third domain (WW3*) showing the strongest affinity to the PY motif. To understand the mechanism underlying this binding affinity, we have carried out NMR structural and dynamics analyses of the hNedd4-1 WW3* domain in complex with a peptide comprising the C-terminal tail of the human ENaC α-subunit. The structure reveals that the peptide interacts in a similar manner to other WW domain-ENaC peptide structures. Crucial interactions that likely provide binding affinity are the broad XP groove facilitating additional contacts between the WW3* domain and the peptide, compared to similar complexes, and the large surface area buried (83Å(2)) between R430 (WW3*) and L647' (αENaC). This corroborates the model-free analysis of the (15)N backbone relaxation data, which showed that R430 is the most rigid residue in the domain (S(2)=0.90±0.01). Carr-Purcell-Meiboom-Gill relaxation dispersion analysis identified two different conformational exchange processes on the μs-ms time-scale. One of these processes involves residues located at the peptide binding interface, suggesting conformational exchange may play a role in peptide recognition. Thus, both structural and dynamic features of the complex appear to define the high binding affinity. The results should aid interpretation of biochemical data and modeling interfaces between Nedd4-1 and other interacting proteins.

Structural and functional analysis of the DEAF-1 and BS69 MYND domains

DEAF-1 is an important transcriptional regulator that is required for embryonic development and is linked to clinical depression and suicidal behavior in humans. It comprises various structural domains, including a SAND domain that mediates DNA binding and a MYND domain, a cysteine-rich module organized in a Cys(4)-Cys(2)-His-Cys (C4-C2HC) tandem zinc binding motif. DEAF-1 transcription regulation activity is mediated through interactions with cofactors such as NCoR and SMRT. Despite the important biological role of the DEAF-1 protein, little is known regarding the structure and binding properties of its MYND domain.Here, we report the solution structure, dynamics and ligand binding of the human DEAF-1 MYND domain encompassing residues 501-544 determined by NMR spectroscopy. The structure adopts a ββα fold that exhibits tandem zinc-binding sites with a cross-brace topology, similar to the MYND domains in AML1/ETO and other proteins. We show that the DEAF-1 MYND domain binds to peptides derived from SMRT and NCoR corepressors. The binding surface mapped by NMR titrations is similar to the one previously reported for AML1/ETO. The ligand binding and molecular functions of the related BS69 MYND domain were studied based on a homology model and mutational analysis. Interestingly, the interaction between BS69 and its binding partners (viral and cellular proteins) seems to require distinct charged residues flanking the predicted MYND domain fold, suggesting a different binding mode. Our findings demonstrate that the MYND domain is a conserved zinc binding fold that plays important roles in transcriptional regulation by mediating distinct molecular interactions with viral and cellular proteins.

Calibration of Boltzmann distribution priors in Bayesian data analysis

The Boltzmann distribution is commonly used as a prior probability in Bayesian data analysis. Examples include the Ising model in statistical image analysis and the canonical ensemble based on molecular dynamics force fields in protein structure calculation. These models involve a temperature or weighting factor that needs to be inferred from the data. Bayesian inference stipulates to determine the temperature based on the model evidence. This is challenging because the model evidence, a ratio of two high-dimensional normalization integrals, cannot be calculated analytically. We outline a replica-exchange Monte Carlo scheme that allows us to estimate the model evidence by use of multiple histogram reweighting. The method is illustrated for an Ising model and examples in protein structure determination.

Structural characterization of the cyclic cystine ladder motif of θ-defensins

The θ-defensins are, to date, the only known ribosomally synthesized cyclic peptides in mammals, and they have promising antimicrobial bioactivities. The characteristic structural motif of the θ-defensins is the cyclic cystine ladder, comprising a cyclic peptide backbone and three parallel disulfide bonds. In contrast to the cyclic cystine knot, which characterizes the plant cyclotides, the cyclic cystine ladder has not been as well described as a structural motif. Here we report the solution structures and nuclear magnetic resonance relaxation properties in aqueous solution of three representative θ-defensins from different species. Our data suggest that the θ-defensins are more rigid and structurally defined than previously thought. In addition, all three θ-defensins were found to self-associate in aqueous solution in a concentration-dependent and reversible manner, a property that might have a role in their mechanism of action. The structural definition of the θ-defensins and the cyclic cystine ladder will help to guide exploitation of these molecules as structural frameworks for the design of peptide drugs.

Phf19 links methylated Lys36 of histone H3 to regulation of Polycomb activity

Polycomb-group proteins are transcriptional repressors with essential roles in embryonic development. Polycomb repressive complex 2 (PRC2) contains the methyltransferase activity for Lys27. However, the role of other histone modifications in regulating PRC2 activity is just beginning to be understood. Here we show that direct recognition of methylated histone H3 Lys36 (H3K36me), a mark associated with activation, by the PRC2 subunit Phf19 is required for the full enzymatic activity of the PRC2 complex. Using NMR spectroscopy, we provide structural evidence for this interaction. Furthermore, we show that Phf19 binds to a subset of PRC2 targets in mouse embryonic stem cells and that this is required for their repression and for H3K27me3 deposition. These findings show that the interaction of Phf19 with H3K36me2 and H3K36me3 is essential for PRC2 complex activity and for proper regulation of gene repression in embryonic stem cells.

The -Cys-X1-X2-Cys- motif of reduced glutaredoxins adopts a consensus structure that explains the low pK(a) of its catalytic cysteine

The -Cys-X1-X2-Cys- active site motif is central to the function of enzymes of the thioredoxin superfamily, including glutaredoxins. Their chemistry depends on the lowered pK(a) of the N-terminal thiolate cysteine of the -Cys-X1-X2-Cys- sequence; therefore its structure, dynamics, and electrostatics matter. Much information about the glutaredoxin structures was obtained by nuclear magnetic resonance (NMR), yet these various NMR structures produced heterogeneous and discordant views of the -Cys-X1-X2-Cys- motifs. This study addresses these inconsistencies by a computational and experimental investigation of three diverse reduced -Cys-X1-X2-Cys- motifs, from human glutaredoxin 1 (hGrx1), Escherichia coli glutaredoxin 2 (EcGrx2), and T4 virus glutaredoxin (T4Grx). The NMR models do not account for the low pK(a) of the N-terminal cysteine. However, extensive investigations of the NMR conformers by simulations yielded consensus structures for the -Cys-X1-X2-Cys- motif, with well-defined orientations for the cysteines. pK(a) calculations indicated that the consensus structure stabilizes the thiolate by local hydrogen bonds. The consensus structures of EcGrx2 and T4Grx formed the basis for predicting low pK(a) values for their N-terminal cysteines. Subsequent experimental titrations showed that these pK(a) values are <5, supporting the validity of the consensus structure. The simulations also revisited the conformational dynamics of side chains around the -Cys-X1-X2-Cys- motif, which allowed reconciliation of calculated and measured pK(a) values for important hGrx1 mutants. The conformational spread of these side chains, which differs between NMR and molecular dynamics models, is likely to be relevant to substrate recognition. The new structural models determined in this work should prove to be valuable in future molecular studies of the glutaredoxins.

Structural basis for 5'-end-specific recognition of single-stranded DNA by the R3H domain from human Sμbp-2

The R3H domain is a conserved sequence motif in nucleic acid binding proteins. Previously, we reported the solution structure of the R3H domain and identified a putative nucleic acid binding site composed of three conserved basic residues [Liepinsh, E., Leonchiks, A., Sharipo, A., Guignard, L. & Otting, G. (2003). Solution structure of the R3H domain from human Sμbp-2. J. Mol. Biol.326, 217-223]. Here, we determine the binding affinities of mononucleotides and dinucleotides for the R3H domain from human Sμbp-2 (Sμbp2-R3H) and map their binding sites on the protein's surface. Although the binding affinities show up to 260-fold selectivity between different nucleotides, their binding sites and conformations seem very similar. Further, we report the NMR structure of the Sμbp2-R3H in complex with deoxyguanosine 5'-monophosphate (dGMP) mimicking the 5'-end of single-stranded DNA. Pseudocontact shifts from a paramagnetic lanthanide tag attached to residue 731 in the mutant A731C confirmed that binding of dGMP brings a loop of the protein into closer proximity. The structure provides the first structural insight into single-stranded nucleic acid recognition by the R3H domain and shows that the R3H domain specifically binds the phosphorylated 5'-end through electrostatic interactions with the two conserved arginines and stacking interactions with the highly conserved histidine.

Molecular dynamics re-refinement of two different small RNA loop structures using the original NMR data suggest a common structure

Restrained molecular dynamics simulations are a robust, though perhaps underused, tool for the end-stage refinement of biomolecular structures. We demonstrate their utility-using modern simulation protocols, optimized force fields, and inclusion of explicit solvent and mobile counterions-by re-investigating the solution structures of two RNA hairpins that had previously been refined using conventional techniques. The structures, both domain 5 group II intron ribozymes from yeast ai5γ and Pylaiella littoralis, share a nearly identical primary sequence yet the published 3D structures appear quite different. Relatively long restrained MD simulations using the original NMR restraint data identified the presence of a small set of violated distance restraints in one structure and a possibly incorrect trapped bulge nucleotide conformation in the other structure. The removal of problematic distance restraints and the addition of a heating step yielded representative ensembles with very similar 3D structures and much lower pairwise RMSD values. Analysis of ion density during the restrained simulations helped to explain chemical shift perturbation data published previously. These results suggest that restrained MD simulations, with proper caution, can be used to "update" older structures or aid in the refinement of new structures that lack sufficient experimental data to produce a high quality result. Notable cautions include the need for sufficient sampling, awareness of potential force field bias (such as small angle deviations with the current AMBER force fields), and a proper balance between the various restraint weights.

Resolution-by-proxy: a simple measure for assessing and comparing the overall quality of NMR protein structures

In protein X-ray crystallography, resolution is often used as a good indicator of structural quality. Diffraction resolution of protein crystals correlates well with the number of X-ray observables that are used in structure generation and, therefore, with protein coordinate errors. In protein NMR, there is no parameter identical to X-ray resolution. Instead, resolution is often used as a synonym of NMR model quality. Resolution of NMR structures is often deduced from ensemble precision, torsion angle normality and number of distance restraints per residue. The lack of common techniques to assess the resolution of X-ray and NMR structures complicates the comparison of structures solved by these two methods. This problem is sometimes approached by calculating "equivalent resolution" from structure quality metrics. However, existing protocols do not offer a comprehensive assessment of protein structure as they calculate equivalent resolution from a relatively small number (<5) of protein parameters. Here, we report a development of a protocol that calculates equivalent resolution from 25 measurable protein features. This new method offers better performance (correlation coefficient of 0.92, mean absolute error of 0.28 Å) than existing predictors of equivalent resolution. Because the method uses coordinate data as a proxy for X-ray diffraction data, we call this measure "Resolution-by-Proxy" or ResProx. We demonstrate that ResProx can be used to identify under-restrained, poorly refined or inaccurate NMR structures, and can discover structural defects that the other equivalent resolution methods cannot detect. The ResProx web server is available at http://www.resprox.ca.

Structural and mechanistic implications of metal binding in the small heat-shock protein αB-crystallin

The human small heat-shock protein αB-crystallin (αB) rescues misfolded proteins from irreversible aggregation during cellular stress. Binding of Cu(II) was shown to modulate the oligomeric architecture and the chaperone activity of αB. However, the mechanistic basis of this stimulation is so far not understood. We provide here first structural insights into this Cu(II)-mediated modulation of chaperone function using NMR spectroscopy and other biophysical approaches. We show that the α-crystallin domain is the elementary Cu(II)-binding unit specifically coordinating one Cu(II) ion with picomolar binding affinity. Putative Cu(II) ligands are His(83), His(104), His(111), and Asp(109) at the dimer interface. These loop residues are conserved among different metazoans, but also for human αA-crystallin, HSP20, and HSP27. The involvement of Asp(109) has direct implications for dimer stability, because this residue forms a salt bridge with the disease-related Arg(120) of the neighboring monomer. Furthermore, we observe structural reorganization of strands β2-β3 triggered by Cu(II) binding. This N-terminal region is known to mediate both the intermolecular arrangement in αB oligomers and the binding of client proteins. In the presence of Cu(II), the size and the heterogeneity of αB multimers are increased. At the same time, Cu(II) increases the chaperone activity of αB toward the lens-specific protein β(L)-crystallin. We therefore suggest that Cu(II) binding unblocks potential client binding sites and alters quaternary dynamics of both the dimeric building block as well as the higher order assemblies of αB.

Protein structure validation by generalized linear model root-mean-square deviation prediction

Large-scale initiatives for obtaining spatial protein structures by experimental or computational means have accentuated the need for the critical assessment of protein structure determination and prediction methods. These include blind test projects such as the critical assessment of protein structure prediction (CASP) and the critical assessment of protein structure determination by nuclear magnetic resonance (CASD-NMR). An important aim is to establish structure validation criteria that can reliably assess the accuracy of a new protein structure. Various quality measures derived from the coordinates have been proposed. A universal structural quality assessment method should combine multiple individual scores in a meaningful way, which is challenging because of their different measurement units. Here, we present a method based on a generalized linear model (GLM) that combines diverse protein structure quality scores into a single quantity with intuitive meaning, namely the predicted coordinate root-mean-square deviation (RMSD) value between the present structure and the (unavailable) "true" structure (GLM-RMSD). For two sets of structural models from the CASD-NMR and CASP projects, this GLM-RMSD value was compared with the actual accuracy given by the RMSD value to the corresponding, experimentally determined reference structure from the Protein Data Bank (PDB). The correlation coefficients between actual (model vs. reference from PDB) and predicted (model vs. "true") heavy-atom RMSDs were 0.69 and 0.76, for the two datasets from CASD-NMR and CASP, respectively, which is considerably higher than those for the individual scores (-0.24 to 0.68). The GLM-RMSD can thus predict the accuracy of protein structures more reliably than individual coordinate-based quality scores.