Distance Restraints Research Articles

A Simple and Efficient Protein Structure Refinement Method.

Improving the quality of a given protein structure can serve as the ultimate solution for accurate protein structure prediction, and seeking such a method is currently a challenge in computational structural biology. In order to promote and encourage much needed such efforts, CASP (Critical Assessment of Structure Prediction) has been providing an ideal computational experimental platform, where it was reported only recently (since CASP10) that systematic protein structure refinement is possible by carrying out extensive (approximately millisecond) MD simulations with proper restraints generated from the given structure. Using an explicit solvent model and much reduced positional and distance restraints than previously exercised, we propose a refinement protocol that combines a series of short (5 ns) MD simulations with energy minimization procedures. Testing and benchmarking on 54 CASP8-10 refinement targets and 34 CASP11 refinement targets shows quite promising results. Using only a small fraction of MD simulation steps (nanosecond versus millisecond), systematic protein structure refinement was demonstrated in this work, indicating that refinement of a given model can be achieved using a few hours of desktop computing.

Folding Membrane Proteins by Deep Transfer Learning.

Computational elucidation of membrane protein (MP) structures is challenging partially due to lack of sufficient solved structures for homology modeling. Here, we describe a high-throughput deep transfer learning method that first predicts MP contacts by learning from non-MPs and then predicts 3D structure models using the predicted contacts as distance restraints. Tested on 510 non-redundant MPs, our method has contact prediction accuracy at least 0.18 better than existing methods, predicts correct folds for 218 MPs, and generates 3D models with root-mean-square deviation (RMSD) less than 4 and 5Å for 57 and 108 MPs, respectively. A rigorous blind test in the continuous automated model evaluation project shows that our method predicted high-resolution 3D models for two recent test MPs of 210 residues with RMSD ∼2Å. We estimated that our method could predict correct folds for 1,345-1,871 reviewed human multi-pass MPs including a few hundred new folds, which shall facilitate the discovery of drugs targeting at MPs.

Exploring the binding pathways of the 14-3-3ζ protein: Structural and free-energy profiles revealed by Hamiltonian replica exchange molecular dynamics with distancefield distance restraints.

The 14-3-3 protein family performs regulatory functions in eukaryotic organisms by binding to a large number of phosphorylated protein partners. Whilst the binding mode of the phosphopeptides within the primary 14-3-3 binding site is well established based on the crystal structures of their complexes, little is known about the binding process itself. We present a computational study of the process by which phosphopeptides bind to the 14-3-3ζ protein. Applying a novel scheme combining Hamiltonian replica exchange molecular dynamics and distancefield restraints allowed us to map and compare the most likely phosphopeptide-binding pathways to the 14-3-3ζ protein. The most important structural changes to the protein and peptides involved in the binding process were identified. In order to bind phosphopeptides to the primary interaction site, the 14-3-3ζ adopted a newly found wide-opened conformation. Based on our findings we additionally propose a secondary interaction site on the inner surface of the 14-3-3ζ dimer, and a direct interference on the binding process by the flexible C-terminal tail. A minimalistic model was designed to allow for the efficient calculation of absolute binding affinities. Binding affinities calculated from the potential of mean force along the binding pathway are in line with the available experimental estimates for two of the studied systems.

The Exact Nuclear Overhauser Enhancement: Recent Advances.

Although often depicted as rigid structures, proteins are highly dynamic systems, whose motions are essential to their functions. Despite this, it is difficult to investigate protein dynamics due to the rapid timescale at which they sample their conformational space, leading most NMR-determined structures to represent only an averaged snapshot of the dynamic picture. While NMR relaxation measurements can help to determine local dynamics, it is difficult to detect translational or concerted motion, and only recently have significant advances been made to make it possible to acquire a more holistic representation of the dynamics and structural landscapes of proteins. Here, we briefly revisit our most recent progress in the theory and use of exact nuclear Overhauser enhancements (eNOEs) for the calculation of structural ensembles that describe their conformational space. New developments are primarily targeted at increasing the number and improving the quality of extracted eNOE distance restraints, such that the multi-state structure calculation can be applied to proteins of higher molecular weights. We then review the implications of the exact NOE to the protein dynamics and function of cyclophilin A and the WW domain of Pin1, and finally discuss our current research and future directions.

Progressive Stereo Locking (PSL): A Residual Dipolar Coupling Based Force Field Method for Determining the Relative Configuration of Natural Products and Other Small Molecules.

Establishing the relative configuration of a bioactive natural product represents the most challenging part in determining its structure. Residual dipolar couplings (RDCs) are sensitive probes of the relative spatial orientation of internuclear vectors. We adapted a force field structure calculation methodology to allow free sampling of both R and S configurations of the stereocenters of interest. The algorithm uses a floating alignment tensor in a simulated annealing protocol to identify the conformations and configurations that best fit experimental RDC and distance restraints (from NOE and J-coupling data). A unique configuration (for rigid molecules) or a very small number of configurations (for less rigid molecules) of the structural models having the lowest chiral angle energies and reasonable magnitudes of the alignment tensor are provided as the best predictions of the unknown configuration. For highly flexible molecules, the progressive locking of their stereocenters into their statistically dominant R or S state dramatically reduces the number of possible relative configurations. The result is verified by checking that the same configuration is obtained by initiating the locking from different regions of the molecule. For all molecules tested having known configurations (with conformations ranging from mostly rigid to highly flexible), the method accurately determined the correct configuration.

Integrating NOE and RDC using sum-of-squares relaxation for protein structure determination.

We revisit the problem of protein structure determination from geometrical restraints from NMR, using convex optimization. It is well-known that the NP-hard distance geometry problem of determining atomic positions from pairwise distance restraints can be relaxed into a convex semidefinite program (SDP). However, often the NOE distance restraints are too imprecise and sparse for accurate structure determination. Residual dipolar coupling (RDC) measurements provide additional geometric information on the angles between atom-pair directions and axes of the principal-axis-frame. The optimization problem involving RDC is highly non-convex and requires a good initialization even within the simulated annealing framework. In this paper, we model the protein backbone as an articulated structure composed of rigid units. Determining the rotation of each rigid unit gives the full protein structure. We propose solving the non-convex optimization problems using the sum-of-squares (SOS) hierarchy, a hierarchy of convex relaxations with increasing complexity and approximation power. Unlike classical global optimization approaches, SOS optimization returns a certificate of optimality if the global optimum is found. Based on the SOS method, we proposed two algorithms-RDC-SOS and RDC-NOE-SOS, that have polynomial time complexity in the number of amino-acid residues and run efficiently on a standard desktop. In many instances, the proposed methods exactly recover the solution to the original non-convex optimization problem. To the best of our knowledge this is the first time SOS relaxation is introduced to solve non-convex optimization problems in structural biology. We further introduce a statistical tool, the Cramér-Rao bound (CRB), to provide an information theoretic bound on the highest resolution one can hope to achieve when determining protein structure from noisy measurements using any unbiased estimator. Our simulation results show that when the RDC measurements are corrupted by Gaussian noise of realistic variance, both SOS based algorithms attain the CRB. We successfully apply our method in a divide-and-conquer fashion to determine the structure of ubiquitin from experimental NOE and RDC measurements obtained in two alignment media, achieving more accurate and faster reconstructions compared to the current state of the art.

A Peptoid Square Helix via Synergistic Control of Backbone Dihedral Angles.

The continued expansion of the fields of macromolecular chemistry and nanoscience has motivated the development of new secondary structures that can serve as architectural elements of innovative materials, molecular machines, biological probes, and even commercial medicines. Synthetic foldamers are particularly attractive systems for developing such elements because they are specifically designed to facilitate synthetic manipulation and functional diversity. However, relatively few predictive design principles exist that permit both rational and modular control of foldamer secondary structure, while maintaining the capacity for facile diversification of displayed functionality. We demonstrate here that the synergistic application of two such principles in the design of peptoid foldamers yields a new and unique secondary structure that we term an "η-helix" due to its repeating turns, which are highly reminiscent of peptide β-turns. Solution-phase structures of η-helices were obtained by simulated annealing using NOE-derived distance restraints, and the NMR spectra of a series of designed η-helices were altogether consistent with the primary adoption of this structure. The structure is resilient to solvent and temperature changes, and accommodates diversification without requiring postsynthetic manipulation. The unique shape, broad structural stability, and synthetic accessibility of η-helices could facilitate their utilization in a wide range of applications.

Selective 1H-1H Distance Restraints in Fully Protonated Proteins by Very Fast Magic-Angle Spinning Solid-State NMR.

Very fast magic-angle spinning (MAS > 80 kHz) NMR combined with high-field magnets has enabled the acquisition of proton-detected spectra in fully protonated solid samples with sufficient resolution and sensitivity. One of the primary challenges in structure determination of protein is observing long-range 1H-1H contacts. Here we use band-selective spin-lock pulses to obtain selective 1H-1H contacts (e.g., HN-HN) on the order of 5-6 Å in fully protonated proteins at 111 kHz MAS. This approach is a major advancement in structural characterization of proteins given that magnetization can be selectively transferred between protons that are 5-6 Å apart despite the presence of other protons at shorter distance. The observed contacts are similar to those previously observed only in perdeuterated proteins with selective protonation. Simulations and experiments show the proposed method has performance that is superior to that of the currently used methods. The method is demonstrated on GB1 and a β-barrel membrane protein, AlkL.

High Precision FRET at Single-molecule Level for Biomolecule Structure Determination.

A protocol on how to perform high-precision interdye distance measurements using Förster resonance energy transfer (FRET) at the single-molecule level in multiparameter fluorescence detection (MFD) mode is presented here. MFD maximizes the usage of all "dimensions" of fluorescence to reduce photophysical and experimental artifacts and allows for the measurement of interdye distance with an accuracy up to ~1 Å in rigid biomolecules. This method was used to identify three conformational states of the ligand-binding domain of the N-methyl-D-aspartate (NMDA) receptor to explain the activation of the receptor upon ligand binding. When comparing the known crystallographic structures with experimental measurements, they agreed within less than 3 Å for more dynamic biomolecules. Gathering a set of distance restraints that covers the entire dimensionality of the biomolecules would make it possible to provide a structural model of dynamic biomolecules.

CompleXView: a server for the interpretation of protein abundance and connectivity information to identify protein complexes.

The molecular understanding of cellular processes requires the identification and characterization of the involved protein complexes. Affinity-purification and mass spectrometric analysis (AP–MS) are performed on a routine basis to detect proteins assembled in complexes. In particular, protein abundances obtained by quantitative mass spectrometry and direct protein contacts detected by crosslinking and mass spectrometry (XL–MS) provide complementary datasets for revealing the composition, topology and interactions of modules in a protein network. Here, we aim to combine quantitative and connectivity information by a webserver tool in order to infer protein complexes. In a first step, modeling protein abundances and functional annotations from Gene Ontology (GO) results in a network which, in a second step, is integrated with connectivity data from XL–MS analysis in order to complement and validate the protein complexes in the network. The output of our integrative approach is a quantitative protein interaction map which is supplemented with topological information of the detected protein complexes. compleXView is built up by two independent modules which are dedicated to the analysis of label-free AP–MS data and to the visualization of the detected complexes in a network together with crosslink-derived distance restraints. compleXView is available to all users without login requirements at http://xvis.genzentrum.lmu.de/compleXView.

Oligomeric Structure of Anabaena Sensory Rhodopsin in a Lipid Bilayer Environment by Combining Solid-State NMR and Long-range DEER Constraints

Oligomerization of membrane proteins is common in nature. Here, we combine spin-labeling double electron–electron resonance (DEER) and solid-state NMR (ssNMR) spectroscopy to refine the structure of an oligomeric integral membrane protein, Anabaena sensory rhodopsin (ASR), reconstituted in a lipid environment. An essential feature of such a combined approach is that it provides structural distance restraints spanning a range of ca 3–60Å while using the same sample preparation (i.e., mutations, paramagnetic labeling, and reconstitution in lipid bilayers) for both ssNMR and DEER. Direct modeling of the multispin effects on DEER signal allowed for the determination of the oligomeric order and for obtaining long-range DEER distance restraints between the ASR trimer subunits that were used to refine the ssNMR structure of ASR. The improved structure of the ASR trimer revealed a more compact packing of helices and side chains at the intermonomer interface, compared to the structure determined using the ssNMR data alone. The extent of the refinement is significant when compared with typical helix movements observed for the active states of homologous proteins. Our combined approach of using complementary DEER and NMR measurements for the determination of oligomeric structures would be widely applicable to membrane proteins where paramagnetic tags can be introduced. Such a method could be used to study the effects of the lipid membrane composition on protein oligomerization and to observe structural changes in protein oligomers upon drug, substrate, and co-factor binding.

C-1311 (Symadex), a potential anti-cancer drug, intercalates into DNA between A and G moieties. NMR-derived and MD-refined stereostructure of the d(GAGGCCTC)2:C-1311 complex

Imidazoacridinone C-1311 (Symadex®) is an antitumor agent which has been recommended for Phase II clinical trials a few years ago. Previously, it was shown experimentally that during the initial stage of its action C-1311 forms stable intercalation complexes with DNA duplexes. Herein, a NMR-derived stereostructure of d(GAGGCCTC)2:C-1311 complex was reported. The ligand was found locating itself between A and G moieties, forming symmetrical DNA:drug 1:2 mol/mol complex. Intercalation site was located upon the DNA-ligand proton/proton dipolar couplings observed in the NOESY spectrum and the performed MD simulations. NMR-derived stereostructure was hence refined by restrained MD using distance restraints obtained from the NOESY data and the result was compared with MD-derived structure of the proposed complex, obtained from the calculations performed with distance restraints applied only for hydrogen bonds in the terminal GC base pairs. The results of both simulations were coherent. Basing on the observed C-1311's intercalation sites and on our previous results concerning the d(CGATCG)2:C-1311 complex, we stated that AG/GA sequences are the preferred binding sites of imidazoacridinone C-1311.

NMR-based method of small changes reveals how DNA mutator APOBEC3A interacts with its single-stranded DNA substrate.

APOBEC3 proteins are double-edged swords. They deaminate cytosine to uracil in single-stranded DNA and provide protection, as part of our innate immune system, against viruses and retrotransposons, but they are also involved in cancer evolution and development of drug resistance. We report a solution-state model of APOBEC3A interaction with its single-stranded DNA substrate obtained with the ‘method of small changes’. This method compares pairwise the 2D 15N–1H NMR spectra of APOBEC3A bearing a deactivating mutation E72A in the presence of 36 slightly different DNA substrates. From changes in chemical shifts of peptide N–H moieties, the positions of each nucleotide relative to the protein can be identified. This provided distance restraints for molecular-dynamic simulations to derive a 3-D molecular model of the APOBEC3A–ssDNA complex. The model reveals that loops 1 and 7 of APOBEC3A move to accommodate substrate binding, indicating an important role for protein-DNA dynamics. Overall, our method may prove useful to study other DNA–protein complexes where crystallographic techniques or full NMR structure calculations are hindered by weak binding or other problems. Subsequent to submission, an APOBEC3A structure with a bound DNA oligomer was published and coordinates released, which has provided an unbiased validation of the ‘method of small changes’.

Exploring RNA conformational space under sparse distance restraints

We show that the application of a small number of restraints predicted by coevolution analysis can provide a powerful restriction on the conformational freedom of an RNA molecule. The greatest degree of restriction occurs when a contact is predicted between the distal ends of a pair of adjacent stemloops but even with this location additional flexibilities in the molecule can mask the contribution. Multiple cross-links, especially those including a pseudoknot provided the strongest restraint on conformational freedom with the effect being most apparent in topologically simple folds and less so if the fold is more topologically entwined. Little was expected for large structures (over 300 bases) and although a few strong localised restrictions were observed, they contributed little to the restraint of the overall fold. Although contacts predicted using a correlated mutation analysis can provide some powerful restrictions on the conformational freedom of RNA molecules, they are too erratic in their occurrence and distribution to provide a general approach to the problem of RNA 3D structure prediction from sequence.

Conformational Analysis of an Antibacterial Cyclodepsipeptide Active against Mycobacterium tuberculosis by a Combined ROE and RDC Analysis.

Griselimycin (GM) and methylgriselimycin (MGM), naturally produced by microorganisms of the genus Streptomyces, are cyclic depsipeptides composed of ten amino acids. They exhibit antibacterial activity against Mycobacterium species by inhibiting the sliding clamp of prokaryotic DNA polymerase III and are therefore considered as potential anti-tuberculosis drugs. The difference between the peptides is the presence of l-(R)-4-methyl-proline in MGM instead of l-proline in GM at position 8 of the amino acid sequence. Methylation increases both metabolic stability and activity of MGM compared to GM. To get deeper insight into the structure-activity relationship, the solution structure of the cyclic part of MGM was determined using rotating-frame nuclear Overhauser effect (ROE) distance restraints and residual dipolar couplings (RDC). The structure of MGM in solution is compared to the structure of GM in a co-crystal with DNA polymerase III subunit beta. As a result, a highly defined structural model of MGM is obtained, which shows related characteristics to the bound GM.

Mutational analysis of cutinase-like enzyme, Cut190, based on the 3D docking structure with model compounds of polyethylene terephthalate

The cutinase-like enzyme, Cut190, from Saccharomonospora viridis AHK190 can degrade the inner block of polyethylene terephthalate (PET) in the presence of Ca2+, and its mutant, S226P/R228S, exhibited increased activity and higher thermostability. The crystal structures of the Cut190 S226P mutant in the absence and presence of Ca2+ were determined, and revealed the large conformational change induced upon Ca2+ binding. However, the substrate-bound 3D structures of Cut190 remained unknown. In this study, to determine the substrate-binding site and improve the enzyme activity, we first built 3D structures of a PET model compound bound to the crystal structures, using the distance restraints between the scissile carbonyl group of the compound and the catalytic site of the enzyme. We then mutated the putative substrate-binding site predicted from the models, and experimentally determined the enzymatic activities of the mutants for the model substrate poly(butylene succinate-co-adipate). The mutated sites with decreased activity were consistent with the putative binding sites predicted by the 3D model from the Ca2+-bound crystal structure, suggesting that the structure of the Ca2+-bound state represents the active state. Notably, we generated two mutants with significantly increased activities.

Systematic evaluation of CS-Rosetta for membrane protein structure prediction with sparse NOE restraints.

We critically test and validate the CS-Rosetta methodology for de novo structure prediction of α-helical membrane proteins (MPs) from NMR data, such as chemical shifts and NOE distance restraints. By systematically reducing the number and types of NOE restraints, we focus on determining the regime in which MP structures can be reliably predicted and pinpoint the boundaries of the approach. Five MPs of known structure were used as test systems, phototaxis sensory rhodopsin II (pSRII), a subdomain of pSRII, disulfide binding protein B (DsbB), microsomal prostaglandin E2 synthase-1 (mPGES-1), and translocator protein (TSPO). For pSRII and DsbB, where NMR and X-ray structures are available, resolution-adapted structural recombination (RASREC) CS-Rosetta yields structures that are as close to the X-ray structure as the published NMR structures if all available NMR data are used to guide structure prediction. For mPGES-1 and Bacillus cereus TSPO, where only X-ray crystal structures are available, highly accurate structures are obtained using simulated NMR data. One main advantage of RASREC CS-Rosetta is its robustness with respect to even a drastic reduction of the number of NOEs. Close-to-native structures were obtained with one randomly picked long-range NOEs for every 14, 31, 38, and 8 residues for full-length pSRII, the pSRII subdomain, TSPO, and DsbB, respectively, in addition to using chemical shifts. For mPGES-1, atomically accurate structures could be predicted even from chemical shifts alone. Our results show that atomic level accuracy for helical membrane proteins is achievable with CS-Rosetta using very sparse NOE restraint sets to guide structure prediction. Proteins 2017; 85:812-826. © 2016 Wiley Periodicals, Inc.

Multiple helical conformations of the helix-turn-helix region revealed by NOE-restrained MD simulations of tryptophan aporepressor, TrpR.

The nature of flexibility in the helix-turn-helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature-sensitive mutant showed more helical character in its helix-turn-helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild-type aporepressor are used as time-averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix-turn-helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731-740. © 2016 Wiley Periodicals, Inc.

Mass spectrometry-based cross-linking study shows that the Psb28 protein binds to cytochrome b559 in Photosystem II

Photosystem II (PSII), a large pigment protein complex, undergoes rapid turnover under natural conditions. During assembly of PSII, oxidative damage to vulnerable assembly intermediate complexes must be prevented. Psb28, the only cytoplasmic extrinsic protein in PSII, protects the RC47 assembly intermediate of PSII and assists its efficient conversion into functional PSII. Its role is particularly important under stress conditions when PSII damage occurs frequently. Psb28 is not found, however, in any PSII crystal structure, and its structural location has remained unknown. In this study, we used chemical cross-linking combined with mass spectrometry to capture the transient interaction of Psb28 with PSII. We detected three cross-links between Psb28 and the α- and β-subunits of cytochrome b559, an essential component of the PSII reaction-center complex. These distance restraints enable us to position Psb28 on the cytosolic surface of PSII directly above cytochrome b559, in close proximity to the QB site. Protein-protein docking results also support Psb28 binding in this region. Determination of the Psb28 binding site and other biochemical evidence allow us to propose a mechanism by which Psb28 exerts its protective effect on the RC47 intermediate. This study also shows that isotope-encoded cross-linking with the "mass tags" selection criteria allows confident identification of more cross-linked peptides in PSII than has been previously reported. This approach thus holds promise to identify other transient protein-protein interactions in membrane protein complexes.

The Interval Branch-And-Prune Algorithm for the Protein Structure Determination

A general trend of structural biology is the switch of a rigid description of protein structures, as in the first X-ray structures in the 50's, to a flexible description, in which each protein populate several distinct conformations or even a continuum of conformations. This flexibility was shown in numerous cases to have a crucial importance in the function of proteins. Most of the methods for bio-molecular structure calculations are, up to now, based on a combination of sampling and optimization, which does not allow a systematic exploration of the conformational space. But, in the case of highly flexible bio-molecules, a systematic (or global) exploration of the conformational space would be very welcome.On the other hand, the interval branch-and-prune (iBP) approach has been proposed as a method for allowing a global optimization of molecular structure under distance restraints, the so-called Distance Geometry problem. A recursive implementation of this algorithm has permitted to apply this approach on small structures of proteins in alpha-bundles for which few long-range distance restraints were known. Here, we are going to present the results obtained on a set of protein structures with sizes from 24 to 100 residues, displaying various secondary structures and topology. The distance restraints present on these structures were chosen to contain exclusively short-range information.The results obtained with various sets of distance restraints, including exact values and interval of values, will be presented, in order to experimentally evaluate the complexity of the algorithm on real-case of protein structure determination. Several procedures of acceleration will be used in order to allow a complete exploration of the tree describing the molecular Distance Geometry problem. The efficiency of the exploration will be evaluated using the self-organizing map (SOM) clustering approach, and the quality of obtained conformations with respect to Ramachandran plots and steric clashes will be evaluated.