Identification and mechanistic basis of non-ACE2 blocking neutralizing antibodies from COVID-19 patients with deep RNA sequencing and molecular dynamics simulations.

Variants of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continue to cause disease and impair the effectiveness of treatments. The therapeutic potential of convergent neutralizing antibodies (NAbs) from fully recovered patients has been explored in several early stages of novel drugs. Here, we identified initially elicited NAbs (Ig Heavy, Ig lambda, Ig kappa) in response to COVID-19 infection in patients admitted to the intensive care unit at a single center with deep RNA sequencing (>100 million reads) of peripheral blood as a diagnostic tool for predicting the severity of the disease and as a means to pinpoint specific compensatory NAb treatments. Clinical data were prospectively collected at multiple time points during ICU admission, and amino acid sequences for the NAb CDR3 segments were identified. Patients who survived severe COVID-19 had significantly more of a Class 3 antibody (C135) to SARS-CoV-2 compared to non-survivors (15059.4 vs. 1412.7, p = 0.016). In addition to highlighting the utility of RNA sequencing in revealing unique NAb profiles in COVID-19 patients with different outcomes, we provided a physical basis for our findings via atomistic modeling combined with molecular dynamics simulations. We established the interactions of the Class 3 NAb C135 with the SARS-CoV-2 spike protein, proposing a mechanistic basis for inhibition via multiple conformations that can effectively prevent ACE2 from binding to the spike protein, despite C135 not directly blocking the ACE2 binding motif. Overall, we demonstrate that deep RNA sequencing combined with structural modeling offers the new potential to identify and understand novel therapeutic(s) NAbs in individuals lacking certain immune responses due to their poor endogenous production. Our results suggest a possible window of opportunity for administration of such NAbs when their full sequence becomes available. A method involving rapid deep RNA sequencing of patients infected with SARS-CoV-2 or its variants at the earliest infection time could help to develop personalized treatments using the identified specific NAbs.

Principles of target DNA cleavage and the role of Mg2+ in the catalysis of CRISPR-Cas9.

At the core of the CRISPR-Cas9 genome-editing technology, the endonuclease Cas9 introduces site-specific breaks in DNA. However, precise mechanistic information to ameliorating Cas9 function is still missing. Here, multi-microsecond molecular dynamics, free-energy and multiscale simulations are combined with solution NMR and DNA cleavage experiments to resolve the catalytic mechanism of target DNA cleavage. We show that the conformation of an active HNH nuclease is tightly dependent on the catalytic Mg2+, unveiling its cardinal structural role. This activated Mg2+-bound HNH is consistently described through molecular simulations, solution NMR and DNA cleavage assays, revealing also that the protonation state of the catalytic H840 is strongly affected by active site mutations. Finally, ab-initio QM(DFT)/MM simulations and metadynamics establish the catalytic mechanism, showing that the catalysis is activated by H840 and completed by K866, rationalising DNA cleavage experiments. This information is critical to enhance the enzymatic function of CRISPR-Cas9 toward improved genome-editing.

Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9.

CRISPR-Cas9 (clustered regularly interspaced short palindromic repeat and associated Cas9 protein) is a molecular tool with transformative genome editing capabilities. At the molecular level, an intricate allosteric signaling is critical for DNA cleavage, but its role in the specificity enhancement of the Cas9 endonuclease is poorly understood. Here, multi-microsecond molecular dynamics is combined with solution NMR and graph theory-derived models to probe the allosteric role of key specificity-enhancing mutations. We show that mutations responsible for increasing the specificity of Cas9 alter the allosteric structure of the catalytic HNH domain, impacting the signal transmission from the DNA recognition region to the catalytic sites for cleavage. Specifically, the K855A mutation strongly disrupts the allosteric connectivity of the HNH domain, exerting the highest perturbation on the signaling transfer, while K810A and K848A result in more moderate effects on the allosteric communication. This differential perturbation of the allosteric signal correlates to the order of specificity enhancement (K855A > K848A ~ K810A) observed in biochemical studies, with the mutation achieving the highest specificity most strongly perturbing the signaling transfer. These findings suggest that alterations of the allosteric communication from DNA recognition to cleavage are critical to increasing the specificity of Cas9 and that allosteric hotspots can be targeted through mutational studies for improving the system's function.

A simple approach for reconstruction of non-uniformly sampled pseudo-3D NMR data for accurate measurement of spin relaxation parameters.

We explain how to conduct a pseudo-3D relaxation series NUS measurement so that it can be reconstructed by existing 3D NUS reconstruction methods to give accurate relaxation values. We demonstrate using reconstruction algorithms IST and SMILE that this 3D approach allows lower sampling densities than for independent 2D reconstructions. This is in keeping with the common finding that higher dimensionality increases signal sparsity, enabling lower sampling density. The approach treats the relaxation series as ordinary 3D time-domain data whose imaginary part in the pseudo-dimension is zero, and applies any suitably linear 3D NUS reconstruction method accordingly. Best results on measured and simulated data were achieved using acquisitions with 9 to 12 planes and exponential spacing in the pseudo-dimension out to ~ 2 times the inverse decay time. Given these criteria, in typical cases where 2D reconstructions require 50% sampling, the new 3D approach generates spectra reliably at sampling densities of 25%.

1H, 13C, 15 N backbone resonance assignment of the recognition lobe subdomain 3 (Rec3) from Streptococcus pyogenes CRISPR-Cas9.

Rec3 is a subdomain of the recognition (Rec) lobe within CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-associated protein Cas9 that is involved in nucleic acid binding and is critical to HNH endonuclease activation. Here, we report the backbone resonance assignments of an engineered construct of the Rec3 subdomain from Streptococcus pyogenes Cas9. We also analyze backbone chemical shift data to predict secondary structure and an overall fold that is consistent with that of Rec3 from the full-length S. pyogenes Cas9 protein.

Simultaneous Bilateral Fractures of the Tali.

A 22-year-old male US Military Academy cadet fell while sidestepping across the 8-ft-high bar portion of the indoor obstacle course. The cadet, in immense pain, was unable to bear weight immediately after the fall. Following examination by physical therapists within the fitness center, emergency medical service personnel transported the cadet to the emergency department for definitive care. Radiographs demonstrated impaction fractures of the tali, while computed tomography imaging of his ankles demonstrated Hawkins type 2 fractures bilaterally, indicating talar neck displacement and dislocation of the subtalar joints. J Orthop Sports Phys Ther 2020;50(7):409. doi:10.2519/jospt.2020.9098.

Acute Fracture of the Third Metacarpal in a Handball Player.

A 20-year-old right hand-dominant male military cadet presented to the direct-access physical therapy clinic complaining of pain and swelling of his right hand, which was injured while competing in a team handball match the day before. Due to suspicion of a third metacarpal fracture, fluoroscopy was performed on the cadet's hand in the physical therapy clinic, and an apparent oblique fracture was noted. Confirmatory radiographs were ordered and the cadet was referred for orthopaedic consultation. J Orthop Sports Phys Ther 2020;50(1):44. doi:10.2519/jospt.2020.9034.

NMR and computational methods for molecular resolution of allosteric pathways in enzyme complexes.

Allostery is a ubiquitous biological mechanism in which a distant binding site is coupled to and drastically alters the function of a catalytic site in a protein. Allostery provides a high level of spatial and temporal control of the integrity and activity of biomolecular assembles composed of proteins, nucleic acids, or small molecules. Understanding the physical forces that drive allosteric coupling is critical to harnessing this process for use in bioengineering, de novo protein design, and drug discovery. Current microscopic models of allostery highlight the importance of energetics, structural rearrangements, and conformational fluctuations, and in this review, we discuss the synergistic use of solution NMR spectroscopy and computational methods to probe these phenomena in allosteric systems, particularly protein-nucleic acid complexes. This combination of experimental and theoretical techniques facilitates an unparalleled detection of subtle changes to structural and dynamic equilibria in biomolecules with atomic resolution, and we provide a detailed discussion of specialized NMR experiments as well as the complementary methods that provide valuable insight into allosteric pathways in silico. Lastly, we highlight two case studies to demonstrate the adaptability of this approach to enzymes of varying size and mechanistic complexity.

Allosteric Motions of the CRISPR-Cas9 HNH Nuclease Probed by NMR and Molecular Dynamics.

CRISPR-Cas9 is a widely employed genome-editing tool with functionality reliant on the ability of the Cas9 endonuclease to introduce site-specific breaks in double-stranded DNA. In this system, an intriguing allosteric communication has been suggested to control its DNA cleavage activity through flexibility of the catalytic HNH domain. Here, solution NMR experiments and a novel Gaussian-accelerated molecular dynamics (GaMD) simulation method are used to capture the structural and dynamic determinants of allosteric signaling within the HNH domain. We reveal the existence of a millisecond time scale dynamic pathway that spans HNH from the region interfacing the adjacent RuvC nuclease and propagates up to the DNA recognition lobe in full-length CRISPR-Cas9. These findings reveal a potential route of signal transduction within the CRISPR-Cas9 HNH nuclease, advancing our understanding of the allosteric pathway of activation. Further, considering the role of allosteric signaling in the specificity of CRISPR-Cas9, this work poses the mechanistic basis for novel engineering efforts aimed at improving its genome-editing capability.

1H, 13C, 15N backbone and side chain resonance assignment of the HNH nuclease from Streptococcus pyogenes CRISPR-Cas9.

HNH is one of two endonuclease domains of the clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein Cas9 that perform site-specific cleavage of double-stranded DNA. We engineered a novel construct of this critical nuclease from Streptococcus pyogenes Cas9 that not only maintains the wild-type amino acid sequence and fold, but displays enhanced thermostability when compared to the full-length Cas9 enzyme. Here, we report backbone and side chain assignments of the HNH nuclease as a foundational step toward the characterization of protein dynamics and allostery in CRISPR-Cas9.

I260Q DNA polymerase ß highlights precatalytic conformational rearrangements critical for fidelity.

DNA polymerase ß (pol ß) fills single nucleotide gaps in DNA during base excision repair and non-homologous end-joining. Pol ß must select the correct nucleotide from among a pool of four nucleotides with similar structures and properties in order to maintain genomic stability during DNA repair. Here, we use a combination of X-ray crystallography, fluorescence resonance energy transfer and nuclear magnetic resonance to show that pol ß's ability to access the appropriate conformations both before and upon binding to nucleotide substrates is integral to its fidelity. Importantly, we also demonstrate that the inability of the I260Q mutator variant of pol ß to properly navigate this conformational landscape results in error-prone DNA synthesis. Our work reveals that precatalytic conformational rearrangements themselves are an important underlying mechanism of substrate selection by DNA pol ß.

Altering the allosteric pathway in IGPS suppresses millisecond motions and catalytic activity.

Imidazole glycerol phosphate synthase (IGPS) is a V-type allosteric enzyme, meaning that its catalytic rate is critically dependent on activation by its allosteric ligand, N'-[(5'-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (PRFAR). The allosteric mechanism of IGPS is reliant on millisecond conformational motions for efficient catalysis. We engineered four mutants of IGPS designed to disrupt millisecond motions and allosteric coupling to identify regions that are critical to IGPS function. Multiple-quantum Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments and NMR chemical shift titrations reveal diminished enzyme flexibility and a reshaping of the allosteric connectivity in each mutant construct, respectively. The functional relevance of the observed motional quenching is confirmed by significant reductions in glutaminase kinetic activity and allosteric ligand binding affinity. This work presents relevant conclusions toward the control of protein allostery and design of unique allosteric sites for potential enzyme inhibitors with regulatory or therapeutic benefit.

Brief early-life non-specific incorporation of deuterium extends mean life span in Drosophila melanogaster without affecting fecundity.

We have investigated the effects of brief, non-specific deuteration of Drosophila melanogaster by including varying percentages of ²H (D) in the H2O used in the food mix consumed during initial development. Up to 22.5% deuterium oxide (D2O) in H2O was administered, with the result that a low percentage of D2O in the water increased mean life span, whereas the highest percentage used (22.5%) reduced life span. After the one-time treatment period, adult flies were maintained ad libitum with food of normal isotopic distribution. At low deuterium levels, where life span extension was observed, there was no observed change in fecundity. Dead flies were assayed for deuterium incorporation by complete hydrolysis in hot 12 N HCl solution followed by subsequent high-performance liquid chromatography/mass spectrometry (HPLC/MS). Isoleucine and leucine residues showed a small, linear dose-dependent incorporation of deuterium at non-exchangeable sites. Although high levels of D2O itself are toxic for other reasons, higher levels of deuterium incorporation, which can be achieved without toxicity by strategies that avoid direct use of D2O, are clearly worth exploring.