@article{lapin_nevzorov_2022, title={Validation of protein backbone structures calculated from NMR angular restraints using Rosetta (vol 73, pg 229, 2019)}, ISSN={["1573-5001"]}, DOI={10.1007/s10858-022-00398-w}, journal={JOURNAL OF BIOMOLECULAR NMR}, author={Lapin, Joel and Nevzorov, Alexander A.}, year={2022}, month={Aug} }
 @article{lapin_awosanya_esteves_nevzorov_2021, title={H-1/C-13/N-15 triple-resonance experiments for structure determinaton of membrane proteins by oriented-sample NMR}, volume={111}, ISSN={["1527-3326"]}, DOI={10.1016/j.ssnmr.2020.101701}, abstractNote={The benefits of triple-resonance experiments for structure determination of macroscopically oriented membrane proteins by solid-state NMR are discussed. While double-resonance 1H/15N experiments are effective for structure elucidation of alpha-helical domains, extension of the method of oriented samples to more complex topologies and assessing side-chain conformations necessitates further development of triple-resonance (1H/13C/15N) NMR pulse sequences. Incorporating additional spectroscopic dimensions involving 13C spin-bearing nuclei, however, introduces essential complications arising from the wide frequency range of the 1H-13C dipolar couplings and 13C CSA (>20 ​kHz), and the presence of the 13C-13C homonuclear dipole-dipole interactions. The recently reported ROULETTE-CAHA pulse sequence, in combination with the selective z-filtering, can be used to evolve the structurally informative 1H-13C dipolar coupling arising from the aliphatic carbons while suppressing the signals from the carbonyl and methyl regions. Proton-mediated magnetization transfer under mismatched Hartman-Hahn conditions (MMHH) can be used to correlate 13C and 15N nuclei in such triple-resonance experiments for the subsequent 15N detection. The recently developed pulse sequences are illustrated for n-acetyl Leucine (NAL) single crystal and doubly labeled Pf1 coat protein reconstituted in magnetically aligned bicelles. An interesting observation is that in the case of 15N-labeled NAL measured at 13C natural abundance, the triple (1H/13C/15N) MMHH scheme predominantly gives rise to long-range intermolecular magnetization transfers from 13C to 15N spins; whereas direct Hartmann-Hahn 13C/15N transfer is entirely intramolecular. The presented developments advance NMR of oriented samples for structure determination of membrane proteins and liquid crystals.}, journal={SOLID STATE NUCLEAR MAGNETIC RESONANCE}, author={Lapin, Joel and Awosanya, Emmanuel O. and Esteves, Richard J. A. and Nevzorov, Alexander A.}, year={2021}, month={Feb} }
 @article{lapin_nevzorov_2020, title={Computer-generated pulse sequences for H-1-N-15 and H-1(alpha)-C-13(alpha) separated local-field experiments}, volume={317}, ISSN={["1096-0856"]}, DOI={10.1016/j.jmr.2020.106794}, abstractNote={High-resolution separated local field (SLF) experiments are employed in oriented-sample solid state NMR to measure angular-dependent heteronuclear dipolar couplings for structure determination. While traditionally these experiments have been designed analytically by determining cycles of pulses with specific phases and durations to achieve cancellation of the homonuclear dipolar terms in the average Hamiltonian, recent work has introduced a computational approach to optimizing linewidths of the 1H-15N dipolar resonances. Accelerated by GPU processors, a computer algorithm searches for the optimal parameters by simulating numerous 1H-15N NMR spectra. This approach, termed ROULETTE, showed promising results by developing a new pulse sequence (ROULETTE-1.0) exhibiting 18% sharper mean linewidths than SAMPI4 for an N-acetyl Leucine (NAL) crystal. Herein, we expand on this previous work to improve the performance of the 1H-15N SLF experiment and extend the work beyond the original approach to new SLF experiments. The new algorithm, in addition to finding pulse durations and phases, now searches for the optimal on/off application scheme of radio frequency irradiation on each channel. This constitutes true de novo optimization, effectively optimizing every aspect of a pulse sequence instead of just phases and durations. With an improved ROULETTE algorithm, we have found a new 1H-15N pulse sequence, termed ROULETTE-2.0, yielding 32% sharper mean linewidths than SAMPI4 for NAL crystal at 500 MHz 1H frequency. Whereas both SAMPI4 and ROULETTE-1.0 have a window where the rf power on the I-channel is turned off, the new pulse sequence is entirely windowless. Furthermore, the reliability of the algorithm has been greatly improved in terms of avoiding false positives, i.e. well-performing pulse sequences in silica that fail to render narrow resonances in experiment. The program has been extended to the 13Cα-1Hα SLF experiments, using a 6 subdwell architecture similar to the 1H-15N optimization. Compared to the PISEMA pulse sequence, the mean 13Cα-1Hα linewidth is 17% sharper for the new pulse sequence, termed ROULETTE-CAHA. In addition to superior performance, the work demonstrates the broad applicability of the algorithm and its adaptability to different NMR experiments and spin systems.}, journal={JOURNAL OF MAGNETIC RESONANCE}, author={Lapin, Joel and Nevzorov, Alexander A.}, year={2020}, month={Aug} }
 @article{lapin_nevzorov_2020, title={De novo NMR pulse sequence design using Monte-Carlo optimization techniques}, volume={310}, ISBN={1096-0856}, DOI={10.1016/j.jmr.2019.106641}, abstractNote={Separated Local Field (SLF) experiments have been routinely used for measuring 1H-15N heteronuclear dipolar couplings in oriented-sample solid-state NMR for structure determination of proteins. In the on-going pursuit of designing better-performing SLF pulse sequences (e.g. by increasing the number of subdwells, and varying the rf amplitudes and phases), analytical treatment of the relevant average Hamiltonian terms may become cumbersome and/or nearly impossible. Numerical simulations of NMR experiments using GPU processors can be employed to rapidly calculate spectra for moderately sized spin systems, which permit an efficient numeric optimization of pulse sequences by the Monte Carlo Simulated Annealing protocol. In this work, a computational strategy was developed to find the optimal phases and timings that substantially improve the 1H-15N dipolar linewidths over a broad range of dipolar couplings as compared to SAMPI4. More than 100 pulse sequences were developed de novo and tested on an N-acetyl Leucine crystal. Seventeen distinct pulse sequences were shown to produce sharper mean linewidths than SAMPI4. Overall, these pulse sequences have more variable parameters (involving non-quadrature phases) and do not involve symmetry between the odd and even dwells, which would likely preclude their rigorous analytical treatment. The top performing pulse sequence, termed ROULETTE-1, has 18% sharper mean linewidths than SAMPI4 when run on an N-acetyl Leucine crystal. This sequence was also shown to be robust over a broad range of 1H carrier frequencies and various crystal orientations. The performance of such an optimized pulse sequence was also illustrated on 15N Leucine-labeled Pf1 coat protein reconstituted in magnetically aligned bicelles. For the optimized pulse sequence the mean peak width was 14% sharper than SAMPI4, which in turn yielded a better signal to noise ratio, 20:1 vs. 17:1. This method is potentially extendable to de novo development of a variety of NMR experiments.}, journal={JOURNAL OF MAGNETIC RESONANCE}, author={Lapin, Joel and Nevzorov, Alexander A.}, year={2020}, month={Jan} }
 @article{awosanya_lapin_nevzorov_2020, title={NMR "Crystallography" for Uniformly (C-13, N-15)-Labeled Oriented Membrane Proteins}, volume={59}, ISSN={["1521-3773"]}, DOI={10.1002/anie.201915110}, abstractNote={In oriented-sample (OS) solid-state NMR of membrane proteins, the angular-dependent dipolar couplings and chemical shifts provide direct input for structure calculations. However, so far only the  1 H- 15 N dipolar couplings and  15 N chemical shifts have been routinely assessed in oriented  15 N-labeled samples. The main stumbling block for extending this technique to membrane proteins of arbitrary topology has remained in the lack of additional experimental restraints. We have developed a new experimental triple-resonance NMR technique, which was applied to uniformly doubly ( 15 N,  13 C) labeled Pf1 coat protein in magnetically aligned DMPC/DHPC bicelles. The previously inaccessible  1 H α - 13 C α  dipolar have been measured, which make it possible to determine the torsion angles between the peptide planes without assuming α-helical structure  a priori .  The fitting of three angular restraints per peptide plane and filtering by Rosetta scoring functions has yielded a consensus α-helical transmembrane structure for Pf1 protein.}, number={9}, journal={ANGEWANDTE CHEMIE-INTERNATIONAL EDITION}, author={Awosanya, Emmanuel O. and Lapin, Joel and Nevzorov, Alexander A.}, year={2020}, month={Feb}, pages={3554–3557} }
 @article{lapin_nevzorov_2019, title={Validation of protein backbone structures calculated from NMR angular restraints using Rosetta}, volume={73}, ISSN={["1573-5001"]}, DOI={10.1007/s10858-019-00251-7}, abstractNote={Multidimensional solid-state NMR spectra of oriented membrane proteins can be used to infer the backbone torsion angles and hence the overall protein fold by measuring dipolar couplings and chemical shift anisotropies, which depend on the orientation of each peptide plane with respect to the external magnetic field. However, multiple peptide plane orientations can be consistent with a given set of angular restraints. This ambiguity is further exacerbated by experimental uncertainty in obtaining and interpreting such restraints. The previously developed algorithms for structure calculations using angular restraints typically involve a sequential walkthrough along the backbone to find the torsion angles between the consecutive peptide plane orientations that are consistent with the experimental data. This method is sensitive to experimental uncertainty in interpreting the peak positions of as low as ± 10 Hz, often yielding high structural RMSDs for the calculated structures. Here we present a significantly improved version of the algorithm which includes the fitting of several peptide planes at once in order to prevent propagation of error along the backbone. In addition, a protocol has been devised for filtering the structural solutions using Rosetta scoring functions in order to find the structures that both fit the spectrum and satisfy bioinformatics restraints. The robustness of the new algorithm has been tested using synthetic angular restraints generated from the known structures for two proteins: a soluble protein 2gb1 (56 residues), chosen for its diverse secondary structure elements, i.e. an alpha-helix and two beta-sheets, and a membrane protein 4a2n, from which the first two transmembrane helices (having a total of 64 residues) have been used. Extensive simulations have been performed by varying the number of fitted planes, experimental error, and the number of NMR dimensions. It has been found that simultaneously fitting two peptide planes always shifted the distribution of the calculated structures toward lower structural RMSD values as compared to fitting a single torsion-angle pair. For each protein, irrespective of the simulation parameters, Rosetta was able to distinguish the most plausible structures, often having structural RMSDs lower than 2 Å with respect to the original structure. This study establishes a framework for de-novo protein structure prediction using a combination of solid-state NMR angular restraints and bioinformatics.}, number={5}, journal={JOURNAL OF BIOMOLECULAR NMR}, author={Lapin, Joel and Nevzorov, Alexander A.}, year={2019}, month={May}, pages={229–244} }