Assessing nanobody interaction with SARS-CoV-2 Nsp9.

The interaction between SARS-CoV-2 non-structural protein Nsp9 and the nanobody 2NSP90 was investigated by NMR spectroscopy using the paramagnetic perturbation methodology PENELOP (Paramagnetic Equilibrium vs Nonequilibrium magnetization Enhancement or LOss Perturbation). The Nsp9 monomer is an essential component of the replication and transcription complex (RTC) that reproduces the viral gRNA for subsequent propagation. Therefore preventing Nsp9 recruitment in RTC would represent an efficient antiviral strategy that could be applied to different coronaviruses, given the Nsp9 relative invariance. The NMR results were consistent with a previous characterization suggesting a 4:4 Nsp9-to-nanobody stoichiometry with the occurrence of two epitope pairs on each of the Nsp9 units that establish the inter-dimer contacts of Nsp9 tetramer. The oligomerization state of Nsp9 was also analyzed by molecular dynamics simulations and both dimers and tetramers resulted plausible. A different distribution of the mapped epitopes on the tetramer surface with respect to the former 4:4 complex could also be possible, as well as different stoichiometries of the Nsp9-nanobody assemblies such as the 2:2 stoichiometry suggested by the recent crystal structure of the Nsp9 complex with 2NSP23 (PDB ID: 8dqu), a nanobody exhibiting essentially the same affinity as 2NSP90. The experimental NMR evidence, however, ruled out the occurrence in liquid state of the relevant Nsp9 conformational change observed in the same crystal structure.

Exploring the role of cyclodextrins as a cholesterol scavenger: a molecular dynamics investigation of conformational changes and thermodynamics.

This study presents a comprehensive analysis of the cholesterol binding mechanism and conformational changes in cyclodextrin (CD) carriers, namely ßCD, 2HPßCD, and MßCD. The results revealed that the binding of cholesterol to CDs was spontaneous and thermodynamically favorable, with van der Waals interactions playing a dominant role, while Coulombic interactions have a negligible contribution. The solubility of cholesterol/ßCD and cholesterol/MßCD complexes was lower compared to cholesterol/2HPßCD complex due to stronger vdW and Coulombic repulsion between water and CDs. Hydrogen bonding was found to have a minor role in the binding process. The investigation of mechanisms and kinetics of binding demonstrated that cholesterol permeates into the CD cavities completely. Replicas consideration indicated that while the binding to 2HPßCD occurred perpendicularly and solely through positioning cholesterol's oxygen toward the primary hydroxyl rim (PHR), the mechanism of cholesterol binding to ßCD and MßCD could take place with the orientation of oxygen towards both rims. Functionalization resulted in decreased cavity polarity, increased constriction tendency, and altered solubility and configuration of the carrier. Upon cholesterol binding, the CDs expanded, increasing the cavity volume in cholesterol-containing systems. The effects of cholesterol on the relative shape anisotropy (κ2) and asphericity parameter (b) in cyclodextrins were investigated. ßCD exhibited a spherical structure regardless of cholesterol presence, while 2HPßCD and MßCD displayed more pronounced non-sphericity in the absence of cholesterol. Loading cholesterol transformed 2HPßCD and MßCD into more spherical shapes, with increased probabilities of higher κ2. MßCD showed a higher maximum peak of κ2 compared to 2HPßCD after cholesterol loading, while 2HPßCD maintained a significant maximum peak at 0.2 for b.

Novel IGFALS mutations with predicted pathogenetic effects by the analysis of AlphaFold structure.

PURPOSE: According to the American College of Medical Genetics (ACMG) classification, variants of uncertain significance (VUS) are gene variations whose impact on the disease risk is not yet known. VUS, therefore, represent an unmet need for genetic counselling. Aim of the study is the use the AlphaFold artificial intelligence algorithm to predict the impact of novel mutations of the IGFALS gene, detected in a subject with short stature and initially classified as VUS according to the ACMG classification. METHODS: A short-stature girl and her parents have been investigated. IGFALS mutations have been detected through clinical exome and confirmed by Sanger sequencing. The potential presence of co-occurring gene alterations was investigated in the proband by whole exome and CGH array. Structure of the ALS protein (encoded by the IGFALS gene) was evaluated through the AlphaFold artificial intelligence algorithm. RESULTS: Two IGFALS variants were found in the proband: c.1349T > C (p.Leu450Pro) and c.1363_1365delCTC (p.Leu455del), both classified as VUS, according to ACMG. Parents' analysis highlighted the in trans position of the two variants. AlphaFold showed that the mutated positions were found the concave side a horseshoe structure of the ALS protein, likely interfering with protein-protein interactions. According to a loss of function (LoF) effect of the two variants, reduced levels of the IGF1 and IGFBP-3 proteins, as well as a growth hormone (GH) excess were detected in the proband's serum. CONCLUSIONS: By using the AlphaFold structure we were able to predict two IGFALS gene mutations initially classified as VUS, as potentially pathogenetic. Our proof-of-concept showed a potential application of AlphaFold as tool to a better inform VUS interpretation of genetic tests.

Biallelic ATOH1 Gene Variant in Siblings With Pontocerebellar Hypoplasia, Developmental Delay, and Hearing Loss.

Background and Objectives: To report on the novel association of biallelic variant in atonal basic helix-loop-helix transcription factor 1 (ATOH1) gene and pontocerebellar hypoplasia (PCH), severe global developmental delay, intellectual disability, and hearing loss in a family with 2 affected siblings. Methods: A detailed clinical assessment and exome sequencing of peripheral blood sample were performed. Segregation analysis with Sanger sequencing and structural modeling of the variant was performed to support the pathogenicity of the variant. Results: A homozygous missense variant (NM_005172.1:c.481C>G) in the ATOH1 gene was identified in the proband and his affected sister. The segregation analysis subsequently confirmed its segregation with an apparently recessive PCH in this family. ATOH1 encodes for the atonal basic helix-loop-helix (bHLH) transcription factor 1, a core transcription factor in the developing cerebellum, brainstem, and dorsal spinal cord, and in the ear. The identified variant results in the p.(Arg161Gly) amino acid substitution in the evolutionarily conserved DNA-binding bHLH domain of the ATOH1 protein. Biallelic missense variants in this domain were previously reported to result in disordered cerebellar development and hearing loss in animal models. In silico homology modeling revealed that p.Arg161Gly in ATOH1 protein probably disrupts a salt bridge with DNA backbone phosphate and increases the flexibility of the bHLH helix-both of which together affect the binding capability of the bHLH domain to the DNA. Discussion: Based on the sequencing results and evidence from structural modeling of the identified variant, as well as with previous reports of ATOH1 gene disruption, we conclude that ATOH1 may represent a novel candidate gene associated with the phenotype of PCH, global developmental delay, and hearing loss in humans.

The corona of protein-gold nanoparticle systems: the role of ionic strength.

The nature of the nanoparticle-protein corona is emerging as a key aspect in determining the impact of nanomaterials on proteins and in general on the biological response. We previously demonstrated that citrate-stabilized gold nanoparticles (Cit-AuNPs) interact with ß2-microglobulin (ß2m) preserving the protein native structure. Moreover, Cit-AuNPs are able to hinder in vitro fibrillogenesis of a ß2m pathologic variant, namely D76N, by reducing the oligomeric association of the protein in solution. Here, we clarify the characteristics of the interaction between ß2m and Cit-AuNPs by means of different techniques, i.e. surface enhanced Raman spectroscopy, NMR and quartz crystal microbalance with dissipation monitoring. All the results obtained clearly show that by simply changing the ionic strength of the medium it is possible to switch from a labile and transient nature of the protein-NP adduct featuring the so-called soft corona, to a more "hard" interaction with a layer of proteins having a longer residence time on the NP surface. This confirms that the interaction between ß2m and Cit-AuNPs is dominated by electrostatic forces which can be tuned by modifying the ionic strength.

NMR-Based Analysis of Nanobodies to SARS-CoV-2 Nsp9 Reveals a Possible Antiviral Strategy Against COVID-19.

Following the entry into the host cell, SARS-CoV-2 replication is mediated by the replication transcription complex (RTC) assembled through a number of nonstructural proteins (Nsps). A monomeric form of Nsp9 is particularly important for RTC assembly and function. In the present study, 136 unique nanobodies targeting Nsp9 are generated. Several nanobodies belonging to different B-cell lineages are expressed, purified, and characterized. Results from immunoassays applied to purified Nsp9 and neat saliva from coronavirus disease (COVID-19) patients show that these nanobodies effectively and specifically recognize both recombinant and endogenous Nsp9. Nuclear magnetic resonance analyses supported by molecular dynamics reveal a composite Nsp9 oligomerization pattern and demonstrate that both nanobodies stabilize the tetrameric form of wild-type Nsp9 also identifying the epitopes on the tetrameric assembly. These results can have important implications in the potential use of these nanobodies to combat viral replication.

A novel de novo HDAC8 missense mutation causing Cornelia de Lange syndrome.

BACKGROUND: Cornelia de Lange syndrome (CdLS) is a rare and clinically variable syndrome characterized by growth impairment, multi-organ anomalies, and a typical set of facial dysmorphisms. Here we describe a 2-year-old female child harboring a novel de novo missense variant in HDAC8, whose phenotypical score, according to the recent consensus on CdLS clinical diagnostic criteria, allowed the diagnosis of a non-classic CdLS. METHODS: Clinical exome sequencing was performed on the trio, identifying a de novo heterozygous variant in HDAC8 (NM_018486; c. 356C>G p.Thr119Arg). Molecular modeling was performed to evaluate putative functional consequence of the HDAC8 protein. RESULTS: The variant HDAC8 c.356C>G is classified as pathogenic following the ACMG (American College of Medical Genetics and Genomics)/AMP (Association for Molecular Pathology) guidelines. By molecular modeling, we confirmed the deleterious effect of this variant, since the amino acid change compromises the conformational flexibility of the HDAC8 loop required for optimal catalytic function. CONCLUSION: We described a novel Thr119Arg mutation in HDAC8 in a patient displaying the major phenotypic traits of the CdLS. Our results suggest that a modest change outside an active site is capable of triggering global structural changes that propagate through the protein scaffold to the catalytic site, creating de facto haploinsufficiency.

Structure of Nanobody Nb23.

BACKGROUND: Nanobodies, or VHHs, are derived from heavy chain-only antibodies (hcAbs) found in camelids. They overcome some of the inherent limitations of monoclonal antibodies (mAbs) and derivatives thereof, due to their smaller molecular size and higher stability, and thus present an alternative to mAbs for therapeutic use. Two nanobodies, Nb23 and Nb24, have been shown to similarly inhibit the self-aggregation of very amyloidogenic variants of ß2-microglobulin. Here, the structure of Nb23 was modeled with the Chemical-Shift (CS)-Rosetta server using chemical shift assignments from nuclear magnetic resonance (NMR) spectroscopy experiments, and used as prior knowledge in PONDEROSA restrained modeling based on experimentally assessed internuclear distances. Further validation was comparatively obtained with the results of molecular dynamics trajectories calculated from the resulting best energy-minimized Nb23 conformers. METHODS: 2D and 3D NMR spectroscopy experiments were carried out to determine the assignment of the backbone and side chain hydrogen, nitrogen and carbon resonances to extract chemical shifts and interproton separations for restrained modeling. RESULTS: The solution structure of isolated Nb23 nanobody was determined. CONCLUSIONS: The structural analysis indicated that isolated Nb23 has a dynamic CDR3 loop distributed over different orientations with respect to Nb24, which could determine differences in target antigen affinity or complex lability.

A novel de novo NIPA1 missense mutation associated to hereditary spastic paraplegia.

SPG6 accounts for 1% of autosomal dominant Hereditary Spastic Paraplegia (HSP) and is caused by pathogenic variants in NIPA1, which encodes a magnesium transporter located in plasma membrane and early endosomes, implicated in neuronal development and maintenance. Here we report a 39-year-old woman affected by progressive gait disturbance associated to absence seizures episodes within childhood. Clinical exome sequencing identified a likely pathogenic de novo heterozygous variant in NIPA1 (NM_144599.5 c.249 C > G; p.Asn83Lys). Molecular modelling was performed to evaluate putative functional consequence of the NIPA1 protein. Indeed, the Asn83Lys modification is predicted to induce a significant perturbation of the protein structure, altering signal transduction or small-molecule transport by modulating the length of the second transmembrane domain. This is the first study reporting a SPG6-affected patient harbouring the NIPA1 p.Asn83Lys mutation.

Case Report: Novel Compound Heterozygous RNASEH2B Mutations Cause Aicardi-Goutières Syndrome.

Aicardi-Goutières Syndrome (AGS) is a rare disorder characterized by neurological and immunological signs. In this study we have described a child with a phenotype consistent with AGS carrying a novel compound heterozygous mutation in RNASEH2B gene. Next Generation Sequencing revealed two heterozygous variants in RNASEH2B gene. We also highlighted a reduction of RNase H2B transcript and protein levels in all the family members. Lower protein levels of RNase H2A have been observed in all the members of the family as well, whereas a deep depletion of RNase H2C has only been identified in the affected child. The structural analysis showed that both mutations remove many intramolecular contacts, possibly introducing conformational rearrangements with a decrease of the stability of RNase H2B and strongly destabilizing the RNase H2 complex. Taken together, these results highlight the importance of an integrated diagnostic approach which takes into consideration clinical, genetic, and molecular analyses.

Entropy of Two-Molecule Correlated Translational-Rotational Motions Using the kth Nearest Neighbor Method.

The entropy associated with rotations, translations, and their coupled motions provides an important contribution to the free energy of many physicochemical processes such as association and solvation. The kth nearest neighbor method, which offers a convenient way to estimate the entropy in high-dimensional spaces, has been previously applied for translational-rotational entropy estimation. Here, we explore the possibility of extending the kth nearest neighbor method to the computation of the entropy of correlated translation-rotations of two molecules, i.e., in the product space of two translation-rotations, both referred to the same independent reference system, which is relevant for all cases in which the correlated translational-rotational motion of more than two molecules is involved. Numerical tests show that, albeit the relatively high dimensionality (12) of the space, the kth nearest neighbor approach provides an accurate estimate for the entropy of two correlated translational-rotational motions, even when computed from a limited number of samples.

Amplifying the spectrum of SPAST gene mutations.

Hereditary spastic paraplegias (HSPs) include a group of neurodegenerative disorders characterized by slowly progressive spasticity and weakness of the lower extremities, caused by axon degeneration of corticospinal tracts. Spastic paraplegia type 4 (SPG4) is the most common autosomal dominant form of HSP and is caused by mutations in the SPAST gene. SPAST gene encodes for the protein spastin, a member of the ATPases Associated with a variety of cellular Activity (AAA) family.We describe a newly variant in SPAST gene, within an Italian family affected by pure HSP. In particular, we found a heterozygous intragenic microdeletion of 3T in exon 13 of SPG4 gene. The 3T deletion results in a mutated protein with a unique leucine residues deletion at the protein position 508, in the AAA ATPase domain. This variant is not registered in any public database either as rare normal variant nor as mutation in SPAST gene and the importance of this aminoacid is confirmed by the absolute conservation in multiple alignments with diverse species. We conclude that the novel SPAST gene variant identified is probably pathogenic and destabilizes the precise arrangement of the nucleotide binding domain, with a consequent loss-of-function of the mutated spastin protein.

Insights into a Protein-Nanoparticle System by Paramagnetic Perturbation NMR Spectroscopy.

BACKGROUND: The interaction between proteins and nanoparticles is a very relevant subject because of the potential applications in medicine and material science in general. Further interest derives from the amyloidogenic character of the considered protein, ß2-microglobulin (ß2m), which may be regarded as a paradigmatic system for possible therapeutic strategies. Previous evidence showed in fact that gold nanoparticles (AuNPs) are able to inhibit ß2m fibril formation in vitro. METHODS: NMR (Nuclear Magnetic Resonance) and ESR (Electron Spin Resonance) spectroscopy are employed to characterize the paramagnetic perturbation of the extrinsic nitroxide probe Tempol on ß2m in the absence and presence of AuNPs to determine the surface accessibility properties and the occurrence of chemical or conformational exchange, based on measurements conducted under magnetization equilibrium and non-equilibrium conditions. RESULTS: The nitroxide perturbation analysis successfully identifies the protein regions where protein-protein or protein-AuNPs interactions hinder accessibility or/and establish exchange contacts. These information give interesting clues to recognize the fibrillation interface of ß2m and hypothesize a mechanism for AuNPs fibrillogenesis inhibition. CONCLUSIONS: The presented approach can be advantageously applied to the characterization of the interface in protein-protein and protein-nanoparticles interactions.

Correction: Exploring exchange processes in proteins by paramagnetic perturbation of NMR spectra.

Correction for 'Exploring exchange processes in proteins by paramagnetic perturbation of NMR spectra' by Yamanappa Hunashal et al., Phys. Chem. Chem. Phys., 2020, 22, 6247-6259, DOI: .

Missense NR2F1 variant in monozygotic twins affected with the Bosch-Boonstra-Schaaf optic atrophy syndrome.

BACKGROUND: The Bosch-Boonstra-Schaaf optic atrophy syndrome (BBSOAS) is an autosomal-dominant disorder (OMIM615722) mostly characterized by optic atrophy and/or hypoplasia, mild intellectual disability, hypotonia, seizures/infantile epilepsy. This disorder is caused by loss-of-function alterations of NR2F1 (i.e., either whole gene deletions or single nucleotide variants) and, to date, 40 patients have been identified with deletions or mutations in this gene. Here we describe two monozygotic twins harboring a de novo missense variant in the DNA-binding domain of NR2F1 (c.313G>A, p.Gly105Ser), with well-characterized features associated to BBSOAS. METHODS: Patients' DNA was analyzed by exome sequencing identifying the missense variant c.313G>A in NR2F1 (NM_005654.4). Furthermore, molecular modeling was performed to evaluate putative differences in DNA binding between wild-type and mutated NR2F1. RESULTS: The missense variant is predicted to be likely pathogenetic following the ACMG (American College of Medical Genetics and Genomics)/AMP (Association for Molecular Pathology) guidelines. Indeed, dynamic simulation experiments highlighted that the Gly105Ser substitution let the formation of a hydrogen bond between the S105 side chain and R142 and a base (G5) of the DNA sequence, allowing us to hypothesize that the G105 residue might be evolutionary conserved due to the absence of a side chain, besides glycine conformational features. Therefore, the G105S variation seems to cause a stiffening and a possible deformation in the protein-DNA complex due to the interaction of residues R142-S105 and G5 on the DNA, compared to the wild-type. CONCLUSION: In summary, we described two monozygotic twins harboring a novel Gly105Ser mutation in NR2F1 DNA binding domain, displaying the classical phenotype of BBSOAS-affected patients. Our computational data suggest a dominant negative effect of this newly characterized missense variant. To date, this is the first genetic report analyzing in silico structural consequences of NR2F1 Gly105Ser substitution.

Exploring exchange processes in proteins by paramagnetic perturbation of NMR spectra.

The effect of extrinsic paramagnetic probes on NMR relaxation rates for surface mapping of proteins and other biopolymers is a widely investigated and powerful NMR technique. Here we describe a new application of those probes. It relies on the setting of the relaxation delay to generate magnetization equilibrium and off-equilibrium conditions, in order to tailor the extent of steady state signal recovery with and without the water-soluble nitroxide Tempol. With this approach it is possible to identify signals whose relaxation is affected by exchange processes and, from the relative assignments, to map the protein residues involved in association or conformational interconversion processes on a micro-to-millisecond time scale. This finding is confirmed by the comparison with the results obtained from relaxation dispersion measurements. This simple and convenient method allows preliminary inspection to highlight regions where structural or chemical exchange events are operative, in order to focus on quantitative subsequent determinations by transverse relaxation dispersion experiments or analogous NMR relaxation studies, and/or to gain insights into the predictions of calculations.

Generalized Born radii computation using linear models and neural networks.

MOTIVATION: Implicit solvent models play an important role in describing the thermodynamics and the dynamics of biomolecular systems. Key to an efficient use of these models is the computation of generalized Born (GB) radii, which is accomplished by algorithms based on the electrostatics of inhomogeneous dielectric media. The speed and accuracy of such computations are still an issue especially for their intensive use in classical molecular dynamics. Here, we propose an alternative approach that encodes the physics of the phenomena and the chemical structure of the molecules in model parameters which are learned from examples. RESULTS: GB radii have been computed using (i) a linear model and (ii) a neural network. The input is the element, the histogram of counts of neighbouring atoms, divided by atom element, within 16 Å. Linear models are ca. 8 times faster than the most widely used reference method and the accuracy is higher with correlation coefficient with the inverse of 'perfect' GB radii of 0.94 versus 0.80 of the reference method. Neural networks further improve the accuracy of the predictions with correlation coefficient with 'perfect' GB radii of 0.97 and ca. 20% smaller root mean square error. AVAILABILITY AND IMPLEMENTATION: We provide a C program implementing the computation using the linear model, including the coefficients appropriate for the set of Bondi radii, as Supplementary Material. We also provide a Python implementation of the neural network model with parameter and example files in the Supplementary Material as well. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

OSCP subunit of mitochondrial ATP synthase: role in regulation of enzyme function and of its transition to a pore.

The permeability transition pore (PTP) is a latent, high-conductance channel of the inner mitochondrial membrane. When activated, it plays a key role in cell death and therefore in several diseases. The investigation of the PTP took an unexpected turn after the discovery that cyclophilin D (the target of the PTP inhibitory effect of cyclosporin A) binds to FO F1 (F)-ATP synthase, thus inhibiting its catalytic activity by about 30%. This observation was followed by the demonstration that binding occurs at a particular subunit of the enzyme, the oligomycin sensitivity conferral protein (OSCP), and that F-ATP synthase can form Ca2+ -activated, high-conductance channels with features matching those of the PTP, suggesting that the latter originates from a conformational change in F-ATP synthase. This review is specifically focused on the OSCP subunit of F-ATP synthase, whose unique features make it a potential pharmacological target both for modulation of F-ATP synthase and its transition to a pore. LINKED ARTICLES: This article is part of a themed section on Mitochondrial Pharmacology: Featured Mechanisms and Approaches for Therapy Translation. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.22/issuetoc.

Dynamics and Thermodynamics of Transthyretin Association from Molecular Dynamics Simulations.

Molecular dynamics simulations are used in this work to probe the structural stability and the dynamics of engineered mutants of transthyretin (TTR), i.e., the double mutant F87M/L110M (MT-TTR) and the triple mutant F87M/L110M/S117E (3M-TTR), in relation to wild-type. Free energy analysis from end-point simulations and statistical effective energy functions are used to analyze trajectories, revealing that mutations do not have major impact on protein structure but rather on protein association, shifting the equilibria towards dissociated species. The result is confirmed by the analysis of 3M-TTR which shows dissociation within the first 10 ns of the simulation, indicating that contacts are lost at the dimer-dimer interface, whereas dimers (formed by monomers which pair to form two extended ß-sheets) appear fairly stable. Overall the simulations provide a detailed view of the dynamics and thermodynamics of wild-type and mutant transthyretins and a rationale of the observed effects.

PDB2ENTROPY and PDB2TRENT: Conformational and Translational-Rotational Entropy from Molecular Ensembles.

Entropy calculation is an important step in the postprocessing of molecular dynamics trajectories or predictive models. In recent years the nearest neighbor method has emerged as a powerful method to deal in a flexible way with the dimensionality of the problem. Here we provide two programs, PBD2ENTROPY and PDB2TRENT that compute the conformational and translational-rotational entropy, respectively, based on the nearest neighbor method. PDB2ENTROPY takes in input two files containing the following: (1) conformational ensembles of the same molecule(s) in PDB format and (2) definitions of torsion angles (a default file is provided where additional user definitions can be easily implemented). PDB2TRENT takes in a file containing samples of the complexed molecules, a string specifying atoms providing the reference framework to superimpose samples, and a string specifying atoms used to compute rotation and translation of one molecule with respect to the other. The C programs and sample demonstration data are available on the GitHub repository (URL: http://github.com/federico-fogolari/pdb2entropy and http://github.com/federico-fogolari/pdb2trent ).