An overview of tools for the validation of protein NMR structures

Restraint validation of biomolecular structures determined by NMR in the Protein Data Bank

Biomolecular structure analysis from experimental NMR studies generally relies on restraints derived from a combination of experimental and knowledge-based data. A challenge for the structural biology community has been a lack of standards for representing these restraints, preventing the establishment of uniform methods of model-vs-data structure validation against restraints and limiting interoperability between restraint-based structure modeling programs. The NEF and NMR-STAR formats provide a standardized approach for representing commonly used NMR restraints. Using these restraint formats, a standardized validation system for assessing structural models of biopolymers against restraints has been developed and implemented in the wwPDB OneDep data deposition-validation-biocuration system. The resulting wwPDB restraint violation report provides a model vs. data assessment of biomolecule structures determined using distance and dihedral restraints, with extensions to other restraint types currently being implemented. These tools are useful for assessing NMR models, as well as for assessing biomolecular structure predictions based on distance restraints.

Restraint Validation of Biomolecular Structures Determined by NMR in the Protein Data Bank

Biomolecular structure analysis from experimental NMR studies generally relies on restraints derived from a combination of experimental and knowledge-based data. A challenge for the structural biology community has been a lack of standards for representing these restraints, preventing the establishment of uniform methods of model-vs-data structure validation against restraints and limiting interoperability between restraint-based structure modeling programs. The NMR exchange (NEF) and NMR-STAR formats provide a standardized approach for representing commonly used NMR restraints. Using these restraint formats, a standardized validation system for assessing structural models of biopolymers against restraints has been developed and implemented in the wwPDB OneDep data deposition-validation-biocuration system. The resulting wwPDB Restraint Violation Report provides a model vs. data assessment of biomolecule structures determined using distance and dihedral restraints, with extensions to other restraint types currently being implemented. These tools are useful for assessing NMR models, as well as for assessing biomolecular structure predictions based on distance restraints.

Methods for Structure Determination of SH2 Domain-Phosphopeptide Complexes by NMR

Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique to solve the structure of biomolecular complexes at atomic resolution in solution. Small proteins such as Src-homology 2 (SH2) domains have fast tumbling rates and long-lived NMR signals, making them particularly suited to be studied by standard NMR methods. SH2 domains are modular proteins whose function is the recognition of sequences containing phosphotyrosines. In this chapter, we describe the application of NMR to assess the interaction between SH2 domains and phosphopeptides and determine the structure of the resulting complexes.

Core packing of well‐defined X‐ray and NMR structures is the same

Numerous studies have investigated the differences and similarities between protein structures determined by solution NMR spectroscopy and those determined by X-ray crystallography. A fundamental question is whether any observed differences are due to differing methodologies or to differences in the behavior of proteins in solution versus in the crystalline state. Here, we compare the properties of the hydrophobic cores of high-resolution protein crystal structures and those in NMR structures, determined using increasing numbers and types of restraints. Prior studies have reported that many NMR structures have denser cores compared with those of high-resolution X-ray crystal structures. Our current work investigates this result in more detail and finds that these NMR structures tend to violate basic features of protein stereochemistry, such as small non-bonded atomic overlaps and few Ramachandran and sidechain dihedral angle outliers. We find that NMR structures solved with more restraints, and which do not significantly violate stereochemistry, have hydrophobic cores that have a similar size and packing fraction as their counterparts determined by X-ray crystallography at high resolution. These results lead us to conclude that, at least regarding the core packing properties, high-quality structures determined by NMR and X-ray crystallography are the same, and the differences reported earlier are most likely a consequence of methodology, rather than fundamental differences between the protein in the two different environments.

In silico designing of putative peptides for targeting pathological protein Htt in Huntington's disease

Huntington's disease is a neurodegenerative disease caused by CAG repeat in the first exon of HTT (Huntingtin) gene, leading to abnormal form of Htt protein containing enlarged polyglutamine strands of variable length that stick together to form aggregates and is toxic to brain causing brain damage. Complete reversal of brain damage is not possible till date but recovery may be possible by peptide therapy. The peptide-based therapy for Huntington's disease includes both poly Q peptide as well as non poly Q peptides like (QBP1)2, p42, Exendin 4, ED11, CaM, BiP, Leuprorelin peptide. The novel approach that is currently being tested in this article is the peptide-based therapy to target the mutated protein. This approach is based on the principle of preventing the aggregation of mutant Htt by blocking the potential sites responsible for protein aggregation and thereby ameliorating the disease symptoms. Herein, we have screened a variety of potential peptides that were known to prevent the protein aggregation, comparatively analyzed their binding affinity with homology modeled Htt protein, designed novel peptides based upon conservation analysis among screened potential peptides as a therapeutic agent, comparatively analyzed the therapeutic potential of novel peptides against modeled Htt protein for investigating the therapeutic prospects of Huntington's disease. We have designed a peptide for the therapy of Huntington's disease by comparing several peptides, which are already in use for Huntington's disease.

NMR-Based Methods for Protein Analysis

Nuclear magnetic resonance (NMR) spectroscopy is a well-established method for analyzing protein structure, interaction, and dynamics at atomic resolution and in various sample states including solution state, solid state, and membranous environment. Thanks to rapid NMR methodology development, the past decade has witnessed a growing number of protein NMR studies in complex systems ranging from membrane mimetics to living cells, which pushes the research frontier further toward physiological environments and offers unique insights in elucidating protein functional mechanisms. In particular, in-cell NMR has become a method of choice for bridging the huge gap between structural biology and cell biology. Herein, we review the recent developments and applications of NMR methods for protein analysis in close-to-physiological environments, with special emphasis on in-cell protein structural determination and the analysis of protein dynamics, both difficult to be accessed by traditional methods.

A method for validating the accuracy of NMR protein structures

We present a method that measures the accuracy of NMR protein structures. It compares random coil index [RCI] against local rigidity predicted by mathematical rigidity theory, calculated from NMR structures [FIRST], using a correlation score (which assesses secondary structure), and an RMSD score (which measures overall rigidity). We test its performance using: structures refined in explicit solvent, which are much better than unrefined structures; decoy structures generated for 89 NMR structures; and conventional predictors of accuracy such as number of restraints per residue, restraint violations, energy of structure, ensemble RMSD, Ramachandran distribution, and clashscore. Restraint violations and RMSD are poor measures of accuracy. Comparisons of NMR to crystal structures show that secondary structure is equally accurate, but crystal structures are typically too rigid in loops, whereas NMR structures are typically too floppy overall. We show that the method is a useful addition to existing measures of accuracy.

Molecular Mechanism by which Cobra Venom Cardiotoxins Interact with the Outer Mitochondrial Membrane

Cardiotoxin CTII from Najaoxiana cobra venom translocates to the intermembrane space (IMS) of mitochondria to disrupt the structure and function of the inner mitochondrial membrane. At low concentrations, CTII facilitates ATP-synthase activity, presumably via the formation of non-bilayer, immobilized phospholipids that are critical in modulating ATP-synthase activity. In this study, we investigated the effects of another cardiotoxin CTI from Najaoxiana cobra venom on the structure of mitochondrial membranes and on mitochondrial-derived ATP synthesis. By employing robust biophysical methods including 31P-NMR and 1H-NMR spectroscopy, we analyzed the effects of CTI and CTII on phospholipid packing and dynamics in model phosphatidylcholine (PC) membranes enriched with 2.5 and 5.0 mol% of cardiolipin (CL), a phospholipid composition that mimics that in the outer mitochondrial membrane (OMM). These experiments revealed that CTII converted a higher percentage of bilayer phospholipids to a non-bilayer and immobilized state and both cardiotoxins utilized CL and PC molecules to form non-bilayer structures. Furthermore, in order to gain further understanding on how cardiotoxins bind to mitochondrial membranes, we employed molecular dynamics (MD) and molecular docking simulations to investigate the molecular mechanisms by which CTII and CTI interactively bind with an in silico phospholipid membrane that models the composition similar to the OMM. In brief, MD studies suggest that CTII utilized the N-terminal region to embed the phospholipid bilayer more avidly in a horizontal orientation with respect to the lipid bilayer and thereby penetrate at a faster rate compared with CTI. Molecular dynamics along with the Autodock studies identified critical amino acid residues on the molecular surfaces of CTII and CTI that facilitated the long-range and short-range interactions of cardiotoxins with CL and PC. Based on our compiled data and our published findings, we provide a conceptual model that explains a molecular mechanism by which snake venom cardiotoxins, including CTI and CTII, interact with mitochondrial membranes to alter the mitochondrial membrane structure to either upregulate ATP-synthase activity or disrupt mitochondrial function.

A novel polyubiquitin chain linkage formed by viral Ubiquitin is resistant to host deubiquitinating enzymes

The Baculoviridae family of viruses encode a viral Ubiquitin (vUb) gene. Though the vUb is homologous to the host eukaryotic Ubiquitin (Ub), its preservation in the viral genome indicates unique functions that are not compensated by the host Ub. We report the structural, biophysical, and biochemical properties of the vUb from Autographa californica multiple nucleo-polyhedrosis virus (AcMNPV). The packing of central helix α1 to the beta-sheet β1-β5 is different between vUb and Ub. Consequently, its stability is lower compared with Ub. However, the surface properties, ubiquitination activity, and the interaction with Ubiquitin-binding domains are similar between vUb and Ub. Interestingly, vUb forms atypical polyubiquitin chain linked by lysine at the 54th position (K54), and the deubiquitinating enzymes are ineffective against the K54-linked polyubiquitin chains. We propose that the modification of host/viral proteins with the K54-linked chains is an effective way selected by the virus to protect the vUb signal from host DeUbiquitinases.

Rare genetic variants in interleukin-37 link this anti-inflammatory cytokine to the pathogenesis and treatment of gout

OBJECTIVE

Gout is characterised by severe interleukin (IL)-1-mediated joint inflammation induced by monosodium urate crystals. Since IL-37 is a pivotal anti-inflammatory cytokine suppressing the activity of IL-1, we conducted genetic and functional studies aimed at elucidating the role of IL-37 in the pathogenesis and treatment of gout.

METHODS

Variant identification was performed by DNA sequencing of all coding bases of IL37 using molecular inversion probe-based resequencing (discovery cohort: gout n=675, controls n=520) and TaqMan genotyping (validation cohort: gout n=2202, controls n=2295). Predictive modelling of the effects of rare variants on protein structure was followed by in vitro experiments evaluating the impact on protein function. Treatment with recombinant IL-37 was evaluated in vitro and in vivo in a mouse model of gout.

RESULTS

We identified four rare variants in IL37 in six of the discovery gout patients; p.(A144P), p.(G174Dfs*16), p.(C181*) and p.(N182S), whereas none emerged in healthy controls (Fisher's exact p-value=0.043). All variants clustered in the functional domain of IL-37 in exon 5 (p-value=5.71×10^-5). Predictive modelling and functional studies confirmed loss of anti-inflammatory functions and we substantiated the therapeutic potential of recombinant IL-37 in the treatment of gouty inflammation. Furthermore, the carrier status of p.(N182S)(rs752113534) was associated with increased risk (OR=1.81, p-value=0.031) of developing gout in hyperuricaemic individuals of Polynesian ancestry.

CONCLUSION

Here, we provide genetic as well as mechanistic evidence for the role of IL-37 in the pathogenesis of gout, and highlight the therapeutic potential of recombinant IL-37 for the treatment of gouty arthritis.

Structural and functional analysis of SMO-1, the SUMO homolog in Caenorhabditis elegans

SUMO proteins are important post-translational modifiers involved in multiple cellular pathways in eukaryotes, especially during the different developmental stages in multicellular organisms. The nematode C. elegans is a well known model system for studying metazoan development and has a single SUMO homolog, SMO-1. Interestingly, SMO-1 modification is linked to embryogenesis and development in the nematode. However, high-resolution information about SMO-1 and the mechanism of its conjugation is lacking. In this work, we report the high-resolution three dimensional structure of SMO-1 solved by NMR spectroscopy. SMO-1 has flexible N-terminal and C-terminal tails on either side of a rigid beta-grasp folded core. While the sequence of SMO-1 is more similar to SUMO1, the electrostatic surface features of SMO-1 resemble more with SUMO2/3. SMO-1 can bind to typical SUMO Interacting Motifs (SIMs). SMO-1 can also conjugate to a typical SUMOylation consensus site as well as to its natural substrate HMR-1. Poly-SMO-1 chains were observed in-vitro even though SMO-1 lacks any consensus SUMOylation site. Typical deSUMOylation enzymes like Senp2 can cleave the poly-SMO-1 chains. Despite being a single gene, the SMO-1 structure allows it to function in a large repertoire of signaling pathways involving SUMO in C. elegans. Structural and functional features of SMO-1 studies described here will be useful to understand its role in development.

Resolution of ab initio shapes determined from small-angle scattering

Spatial resolution is an important characteristic of structural models, and the authors of structures determined by X-ray crystallography or electron cryo-microscopy always provide the resolution upon publication and deposition. Small-angle scattering of X-rays or neutrons (SAS) has recently become a mainstream structural method providing the overall three-dimensional structures of proteins, nucleic acids and complexes in solution. However, no quantitative resolution measure is available for SAS-derived models, which significantly hampers their validation and further use. Here, a method is derived for resolution assessment for ab initio shape reconstruction from scattering data. The inherent variability of the ab initio shapes is utilized and it is demonstrated how their average Fourier shell correlation function is related to the model resolution. The method is validated against simulated data for proteins with known high-resolution structures and its efficiency is demonstrated in applications to experimental data. It is proposed that henceforth the resolution be reported in publications and depositions of ab initio SAS models.

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.

Observing a late folding intermediate of Ubiquitin at atomic resolution by NMR

The study of intermediates in the protein folding pathway provides a wealth of information about the energy landscape. The intermediates also frequently initiate pathogenic fibril formations. While observing the intermediates is difficult due to their transient nature, extreme conditions can partially unfold the proteins and provide a glimpse of the intermediate states. Here, we observe the high resolution structure of a hydrophobic core mutant of Ubiquitin at an extreme acidic pH by nuclear magnetic resonance (NMR) spectroscopy. In the structure, the native secondary and tertiary structure is conserved for a major part of the protein. However, a long loop between the beta strands β3 and β5 is partially unfolded. The altered structure is supported by fluorescence data and the difference in free energies between the native state and the intermediate is reflected in the denaturant induced melting curves. The unfolded region includes amino acids that are critical for interaction with cofactors as well as for assembly of poly-Ubiquitin chains. The structure at acidic pH resembles a late folding intermediate of Ubiquitin and indicates that upon stabilization of the protein's core, the long loop converges on the core in the final step of the folding process.

Analysis of the interface variability in NMR structure ensembles of protein-protein complexes

NMR structures consist in ensembles of conformers, all satisfying the experimental restraints, which exhibit a certain degree of structural variability. We analyzed here the interface in NMR ensembles of protein-protein heterodimeric complexes and found it to span a wide range of different conservations. The different exhibited conservations do not simply correlate with the size of the systems/interfaces, and are most probably the result of an interplay between different factors, including the quality of experimental data and the intrinsic complex flexibility. In any case, this information is not to be missed when NMR structures of protein-protein complexes are analyzed; especially considering that, as we also show here, the first NMR conformer is usually not the one which best reflects the overall interface. To quantify the interface conservation and to analyze it, we used an approach originally conceived for the analysis and ranking of ensembles of docking models, which has now been extended to directly deal with NMR ensembles. We propose this approach, based on the conservation of the inter-residue contacts at the interface, both for the analysis of the interface in whole ensembles of NMR complexes and for the possible selection of a single conformer as the best representative of the overall interface. In order to make the analyses automatic and fast, we made the protocol available as a web tool at: https://www.molnac.unisa.it/BioTools/consrank/consrank-nmr.html.

Computational approaches for inferring the functions of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are ubiquitously involved in cellular processes and often implicated in human pathological conditions. The critical biological roles of these proteins, despite not adopting a well-defined fold, encouraged structural biologists to revisit their views on the protein structure-function paradigm. Unfortunately, investigating the characteristics and describing the structural behavior of IDPs is far from trivial, and inferring the function(s) of a disordered protein region remains a major challenge. Computational methods have proven particularly relevant for studying IDPs: on the sequence level their dependence on distinct characteristics determined by the local amino acid context makes sequence-based prediction algorithms viable and reliable tools for large scale analyses, while on the structure level the in silico integration of fundamentally different experimental data types is essential to describe the behavior of a flexible protein chain. Here, we offer an overview of the latest developments and computational techniques that aim to uncover how protein function is connected to intrinsic disorder.

Molecular-docking study of malaria drug target enzyme transketolase in Plasmodium falciparum 3D7 portends the novel approach to its treatment

BACKGROUND

Malaria has been a major life threatening mosquito borne disease from long since. Unavailability of any effective vaccine and recent emergence of multi drug resistant strains of malaria pathogen Plasmodium falciparum continues to cause persistent deaths in the tropical and sub-tropical region. As a result, demands for new targets for more effective anti-malarial drugs are escalating. Transketolase is an enzyme of the pentose phosphate pathway; a novel pathway which is involved in energy generation and nucleic acid synthesis. Moreover, significant difference in homology between Plasmodium falciparum transketolase (Pftk) and human (Homo sapiens) transketolase makes it a suitable candidate for drug therapy. Our present study is aimed to predict the 3D structure of Plasmodium falciparum transketolase and design an inhibitor against it.

RESULTS

The primary and secondary structural features of the protein is calculated by ProtParam and SOPMA respectively which revealed the protein is composed of 43.3 % alpha helix and 33.04 % random coils along with 15.62 % extended strands, 8.04 % beta turns. The three dimensional structure of the transketolase is constructed using homology modeling tool MODELLAR utilizing several available transketolase structures as templates. The structure is then subjected to deep optimization and validated by structure validation tools PROCHECK, VERIFY 3D, ERRAT, QMEAN. The predicted model scored 0.74 for global model reliability in PROCHECK analysis, which ensures the quality of the model. According to VERIFY 3D the predicted model scored 0.77 which determines good environmental profile along with ERRAT score of 78.313 which is below 95 % rejection limit. Protein-protein and residue-residue interaction networks are generated by STRING and RING server respectively. CASTp server was used to analyze active sites and His 109, Asn 108 and His 515 are found to be more positive site to dock the substrate, in addition molecular docking simulation with Autodock vina determined the estimated free energy of molecular binding was of -6.6 kcal/mol for most favorable binding of 6'-Methyl-Thiamin Diphosphate.

CONCLUSION

This predicted structure of Pftk will serve first hand in the future development of effective Pftk inhibitors with potential anti-malarial activity. However, this is a preliminary study of designing an inhibitor against Plasmodium falciparum 3D7; the results await justification by in vitro and in vivo experimentations.

NMR-based modeling and refinement of protein 3D structures

NMR is a well-established method to characterize the structure and dynamics of biomolecules in solution. High-quality structures can now be produced thanks to both experimental advances and computational developments that incorporate new NMR parameters and improved protocols and force fields in the structure calculation and refinement process. In this chapter, we give a short overview of the various types of NMR data that can provide structural information, and then focus on the structure calculation methodology itself. We discuss and illustrate with tutorial examples "classical" structure calculation, refinement, and structure validation approaches.

Structure calculation, refinement and validation using CcpNmr Analysis

CcpNmr Analysis provides a streamlined pipeline for both NMR chemical shift assignment and structure determination of biological macromolecules. In addition, it encompasses tools to analyse the many additional experiments that make NMR such a pivotal technique for research into complex biological questions. This report describes how CcpNmr Analysis can seamlessly link together all of the tasks in the NMR structure-determination process. It details each of the stages from generating NMR restraints [distance, dihedral, hydrogen bonds and residual dipolar couplings (RDCs)], exporting these to and subsequently re-importing them from structure-calculation software (such as the programs CYANA or ARIA) and analysing and validating the results obtained from the structure calculation to, ultimately, the streamlined deposition of the completed assignments and the refined ensemble of structures into the PDBe repository. Until recently, such solution-structure determination by NMR has been quite a laborious task, requiring multiple stages and programs. However, with the new enhancements to CcpNmr Analysis described here, this process is now much more intuitive and efficient and less error-prone.

A Structured-based Model for the Decreased Activity of Ala222Val and Glu429Ala Methylenetetrahydrofolate Reductase (MTHFR) Mutants

The structure of human Methylenetetrahydrofolate Reductase (MTHFR) is not known either by NMR or by X-ray methods. Phosphorylation seems to play an important role in the functioning of this flavoprotein. MTHFR catalyzes an irreversible reaction in homocysteine metabolism. Phosphorylation decreases the activity of MTHFR by enhancing the sensitivity of the enzyme to SAdenosylmethione. Two common polymorphisms in MTHFR, Ala222Val and Glu429Ala, can result in a number of vascular diseases. Effects of the Glu429Ala polymorphism on the structure of human MTHFR remain undetermined due to limited structural information. Hence, structural models of the MTHFR mutants were constructed using I-TASSER and assessed by PROCHECK, DFIRE and Verify3D tools. A mechanism is further suggested for the decreased activity of the Ala222Val and Glu429Ala mutants due to a decrease in number of serine phosphorylation sites using information gleaned from the molecular models. This provides insights for the understanding of structure-function relationship for MTHFR.

Quality assessment of protein NMR structures

Biomolecular NMR structures are now routinely used in biology, chemistry, and bioinformatics. Methods and metrics for assessing the accuracy and precision of protein NMR structures are beginning to be standardized across the biological NMR community. These include both knowledge-based assessment metrics, parameterized from the database of protein structures, and model versus data assessment metrics. On line servers are available that provide comprehensive protein structure quality assessment reports, and efforts are in progress by the world-wide Protein Data Bank (wwPDB) to develop a biomolecular NMR structure quality assessment pipeline as part of the structure deposition process. These quality assessment metrics and standards will aid NMR spectroscopists in determining more accurate structures, and increase the value and utility of these structures for the broad scientific community.

Recommendations of the wwPDB NMR Validation Task Force

As methods for analysis of biomolecular structure and dynamics using nuclear magnetic resonance spectroscopy (NMR) continue to advance, the resulting 3D structures, chemical shifts, and other NMR data are broadly impacting biology, chemistry, and medicine. Structure model assessment is a critical area of NMR methods development, and is an essential component of the process of making these structures accessible and useful to the wider scientific community. For these reasons, the Worldwide Protein Data Bank (wwPDB) has convened an NMR Validation Task Force (NMR-VTF) to work with wwPDB partners in developing metrics and policies for biomolecular NMR data harvesting, structure representation, and structure quality assessment. This paper summarizes the recommendations of the NMR-VTF, and lays the groundwork for future work in developing standards and metrics for biomolecular NMR structure quality assessment.