Molecular modeling of biomolecules by paramagnetic NMR and computational hybrid methods

Protein structure prediction using residue-resolved protection factors from hydrogen-deuterium exchange NMR

Hydrogen-deuterium exchange (HDX) measured by nuclear magnetic resonance (NMR) provides structural information for proteins relating to solvent accessibility and flexibility. While this structural information is beneficial, the data cannot be used exclusively to elucidate structures. However, the structural information provided by the HDX-NMR data can be supplemented by computational methods. In previous work, we developed an algorithm in Rosetta to predict structures using qualitative HDX-NMR data (categories of exchange rate). Here we expand on the effort, and utilize quantitative protection factors (PFs) from HDX-NMR for structure prediction. From observed correlations between PFs and solvent accessibility/flexibility measures, we present a scoring function to quantify the agreement with HDX data. Using a benchmark set of 10 proteins, an average improvement of 5.13 Å in root-mean-square deviation (RMSD) is observed for cases of inaccurate Rosetta predictions. Ultimately, seven out of 10 predictions are accurate without including HDX data, and nine out of 10 are accurate when using our PF-based HDX score.

An Integrative Approach to Determine 3D Protein Structures Using Sparse Paramagnetic NMR Data and Physical Modeling

Paramagnetic nuclear magnetic resonance (NMR) methods have emerged as powerful tools for structure determination of large, sparsely protonated proteins. However traditional applications face several challenges, including a need for large datasets to offset the sparsity of restraints, the difficulty in accounting for the conformational heterogeneity of the spin-label, and noisy experimental data. Here we propose an integrative approach to structure determination combining sparse paramagnetic NMR with physical modelling to infer approximate protein structural ensembles. We use calmodulin in complex with the smooth muscle myosin light chain kinase peptide as a model system. Despite acquiring data from samples labeled only at the backbone amide positions, we are able to produce an ensemble with an average RMSD of ∼2.8 Å from a reference X-ray crystal structure. Our approach requires only backbone chemical shifts and measurements of the paramagnetic relaxation enhancement and residual dipolar couplings that can be obtained from sparsely labeled samples.

Protein Structure Prediction from NMR Hydrogen-Deuterium Exchange Data

Amide hydrogen-deuterium exchange (HDX) has long been used to determine regional flexibility and binding sites in proteins; however, the data are too sparse for full structural characterization. Experiments that measure HDX rates, such as HDX-NMR, have far higher throughput compared to structure determination via X-ray crystallography, cryo-EM, or a full suite of NMR experiments. Data from HDX-NMR experiments encode information on the protein structure, making HDX a prime candidate to be supplemented by computational algorithms for protein structure prediction. We have developed a methodology to incorporate HDX-NMR data into ab initio protein structure prediction using the Rosetta software framework to predict structures based on experimental agreement. To demonstrate the efficacy of our algorithm, we examined 38 proteins with HDX-NMR data available, comparing the predicted model with and without the incorporation of HDX data into scoring. The root-mean-square deviation (rmsd, a measure of the average atomic distance between superimposed models) of the predicted model improved by 1.42 Å on average after incorporating the HDX-NMR data into scoring. The average rmsd improvement for the proteins where the selected model rmsd changed after incorporating HDX data was 3.63 Å, including one improvement of more than 11 Å and seven proteins improving by greater than 4 Å, with 12/15 proteins improving overall. Additionally, for independent verification, two proteins that were not part of the original benchmark were scored including HDX data, with a dramatic improvement of the selected model rmsd of nearly 9 Å for one of the proteins. Moreover, we have developed a confidence metric allowing us to successfully identify near-native models in the absence of a native structure. Improvement in model selection with a strong confidence measure demonstrates that protein structure prediction with HDX-NMR is a powerful tool which can be performed with minimal additional computational strain and expense.

Applications of Molecular Dynamics Simulation in Structure Prediction of Peptides and Proteins

Compared with rapid accumulation of protein sequences from high-throughput DNA sequencing, obtaining experimental 3D structures of proteins is still much more difficult, making protein structure prediction (PSP) potentially very useful. Currently, a vast majority of PSP efforts are based on data mining of known sequences, structures and their relationships (informatics-based). However, if closely related template is not available, these methods are usually much less reliable than experiments. They may also be problematic in predicting the structures of naturally occurring or designed peptides. On the other hand, physics-based methods including molecular dynamics (MD) can utilize our understanding of detailed atomic interactions determining biomolecular structures. In this mini-review, we show that all-atom MD can predict structures of cyclic peptides and other peptide foldamers with accuracy similar to experiments. Then, some notable successes in reproducing experimental 3D structures of small proteins through MD simulations (some with replica-exchange) of the folding were summarized. We also describe advancements of MD-based refinement of structure models, and the integration of limited experimental or bioinformatics data into MD-based structure modeling.

Predicting Protein Dimer Structures Using MELD × MD

It is challenging to predict the docked conformations of two proteins. Current methods are susceptible to errors from treating proteins as rigid bodies and from an inability to compute relative Boltzmann populations of different docked conformations. Here, we show that by using the ClusPro server as a front end to generate possible protein-protein contacts, and using Modeling Employing Limited Data (MELD) accelerated molecular dynamics (MELD × MD) as a back end for atomistic simulations, we can find 16/20 native dimer structures of small proteins as those having the lowest free energy, starting from good-bound-backbone structures. We show that atomistic MD free energies can be used to identify native protein dimer structures.

Metal-Ligand Recognition Index Determination by NMR Proton Relaxation Study

In this study, we developed and validated a new proposed parameter quantifying the interaction strength between natural and/or synthetic molecules with paramagnetic metal ions. The Metal ion Recognition Index, Miri, is a quantitative parameter to describe the proton environment and to define their involvement in the inner and/or outer sphere of the paramagnetic metal ion. The method is based on the analysis of NMR proton spin-lattice relaxation rates of a specific ligand in both the diamagnetic and paramagnetic conditions. The proposed procedure is also useful to calculate the ligand proton spin-lattice relaxation rate in the paramagnetic bound conditions, which is typically very difficult to determine experimentally. Miri was used to compare the ligand proton involvement toward different paramagnetic species, in particular the Copper(II)-Piroxicam system. Copper(II)-Piroxicam complex is one of the most active anti-inflammatory and anti-arthritic species. Miri provides an opportunity to improve our knowledge of metal-ligand complexes that play a fundamental role in bioinorganic interactions.

MELD-Path Efficiently Computes Conformational Transitions, Including Multiple and Diverse Paths

The molecular actions of proteins occur along reaction coordinates. Current computer methods have limited ability to explore them. We describe a fast protocol called MELD-path that (1) efficiently samples relevant conformational states via MELD, an accelerator of Molecular Dynamics (MD), (2) seeds multiple short MD trajectories from MELD states, and then (3) constructs Markov State Models (MSM) that give the routes and kinetics. We tested the method against extensive (multi μs) MD simulations of the right-handed- to left-handed-helix transition of a 9-mer peptide of AIB, the symmetry of which allows us to establish convergence. MELD-path finds all the metastable states, their correct relative populations, and the full ensemble of routes, not just a single assumed route. For this transition, we find a very broad route structure. MELD-path is highly parallelizable and efficient, yielding the full route map in a few days of computation. We believe MELD-path could be a general and rapid way to explore mechanistic processes in biomolecules on the computer.