Comprehensive Clinical Diagnostic Pipelines Reveal New Variants in Alpha-1 Antitrypsin Deficiency.

Alpha-1 antitrypsin deficiency (AATD) is an underdiagnosed disorder associated with mutations in the SERPINA1 gene encoding alpha-1 antitrypsin (AAT). Severe AATD can manifest as pulmonary emphysema and progressive liver disease. Besides the most common pathogenic variants S (E264V) and Z (E342K), many rarer genetic variants of AAT have been found in patients and in the general population. Here we report a panel of new SERPINA1 variants, including 4 null and 16 missense alleles, identified among a cohort of individuals with suspected AATD whose phenotypic follow-up showed inconclusive or atypical results. Because the pathogenic significance of the missense variants was unclear purely on the basis of clinical data, the integration of computational, biochemical, and cellular studies was used to define the associated risk of disease. Established pathogenicity predictors and structural analysis identified a panel of candidate damaging mutations that were characterized by expression in mammalian cell models. Polymer formation, intracellular accumulation, and secretory efficiency were evaluated experimentally. Our results identified two AAT mutants with a Z-like polymerogenic severe deficiency profile (Smilano and Mcampolongo) and three milder variants (Xsarezzo, Pdublin, and Ctiberias). Overall, the experimentally determined behavior of the variants was in agreement with the pathogenicity scores of the REVEL (an ensemble method for predicting the pathogenicity of rare missense variants) predictor, supporting the utility of this bioinformatic tool in the initial assessment of newly identified amino acid substitutions of AAT. Our study, in addition to describing 20 new SERPINA1 variants, provides a model for a multidisciplinary approach to classification of rare AAT variants and their clinical impact on individuals with rare AATD genotypes.

Study of ferritin self-assembly and heteropolymer formation by the use of Fluorescence Resonance Energy Transfer (FRET) technology.

The high stability and strong self-assembly properties made ferritins the most used proteins for nanotechnological applications. Human ferritins are made of 24 subunits of the H- and L-type that coassemble in an almost spherical nanocage 12nm across, delimiting a large cavity. The mechanism and kinetics of ferritin self-assembly and why H/L heteropolymers formation is favored over the homopolymers remain unclarified. In order to study this, we used the Fluorescence Resonance Energy Transfer (FRET) tool by binding multiple donor or acceptor Alexa Fluor fluorophores on the outer surface of human H and L ferritins and then denaturing and reassembling them in different proportions and conditions. The FRET efficiency increase from <0.3 of the disassembled to >0.7 in the assembled allowed to study the assembly kinetics. We found that their assembly was complete in about one hour, and that the initial rate of self-assembly of H/L heteropolymers was slightly faster than that of the H/H homopolymers. Then, by adding various proportions of unlabeled H or L-chains to the FRET system we found that the presence of the L-chains displaced the formation of H-H dimers more efficiently than that of the H-chains. This favored formation of H/L heterodimers, which is the initial step in ferritin self-assembly, contributes to explain the preferred formation of H/L heteropolymers over the H or L homopolymers. Moreover, we found that the H-chains arrange at distant positions on the heteropolymeric shell until they reach a number above eight, when they start to co-localize.

Looking at human cytosolic sialidase NEU2 structural features with an interdisciplinary approach.

Circular dichroism (CD) spectra at variable temperatures have been recorded for human cytosolic sialidase NEU2 in buffered water solutions and in the presence of divalent cations. The results show the prevalence of ß-strands together with a considerable amount of α-helical structure, while in the solid state, from both previous X-ray diffraction analysis and our CD data on film samples, the content of ß-strands is higher. In solution, a significant change in CD spectra occurs with an increase in temperature, related to a decrease in α-helix content and a slight increase in ß-strand content. In the same range of temperatures, the enzymatic activity decreases. Although the overall structure of the protein appears to be particularly stable, molecular dynamics simulations performed at various temperatures evidence local conformational changes possibly relevant for explaining the relative lability of enzymatic activity.

Disease-associated mutations in the coil 2B domain of human lamin A/C affect structural properties that mediate dimerization and intermediate filament formation.

The lamin proteins are essential components of the nuclear lamina of eukaryotic cells, that are involved in a complex association mechanism to attain a functional supermolecular structure. Mutations of the lamin A/C gene are associated with several different neuromuscular diseases, and the detailed effect of disease-associated amino acid substitutions on the structure and stability of human lamin dimers is yet unknown. Here we present a structural and thermodynamic characterization by means of molecular dynamics simulations of the effect of pathological mutations (S326T, R331P, R331Q, E347K, E358K, M371K, and R377H) on the association of the coil 2B domains that mediate lamin A/C oligomerization. The structures attained during the simulations, along with the quantification of the contribution of each residue to the dimerization energies, support a lamin association mechanism mediated by homophilic intermolecular interactions promoted by dissociative conformational changes at distinct positions in the coiled coil. The pathogenic mutations can both increase or decrease the stability of lamin A/C dimers, and a possible correlation between the effect of the amino acid substitutions and disease onset and severity is presented.

Electrical and mechanical anharmonicities from NIR-VCD spectra of compounds exhibiting axial and planar chirality: the cases of (S)-2,3-pentadiene and methyl-d(3) (R)- and (S)-[2.2]paracyclophane-4-carboxylate.

The IR and Near infrared (NIR) vibrational circular dichroism (VCD) spectra of molecules endowed with noncentral chirality have been investigated. Data for fundamental, first, and second overtone regions of (S)-2,3-pentadiene, exhibiting axial chirality, and methyl-d(3) (R)- and (S)-[2.2]paracyclophane-4-carboxylate, exhibiting planar chirality have been measured and analyzed. The analysis of NIR and IR VCD spectra was based on the local-mode model and the use of density functional theory (DFT), providing mechanical and electrical anharmonic terms for all CH-bonds. The comparison of experimental and calculated spectra is satisfactory and allows one to monitor fine details in the asymmetric charge distribution in the molecules: these details consist in the harmonic frequencies, in the principal anharmonicity constants, in both the atomic polar and axial tensors and in their first and second derivatives with respect to the CH-stretching coordinates.

Estrogen-Like Effect of Mitotane Explained by Its Agonist Activity on Estrogen Receptor-α.

Mitotane is the cornerstone of medical treatment of adrenocortical carcinoma. Estrogenic-like side effects frequently occur in patients, and previous studies explored the chemical nature of the interaction between estrogen receptor-α (ER-α) and toxic compounds, including the DDD derivatives. We used molecular docking and molecular dynamics (MD) simulations to explore the possible interaction between mitotane and the ER-α receptor and the induced conformational changes. The ER-α expressing MCF-7 cells were exposed to mitotane with/without tamoxifen, and the cell viability/proliferation was evaluated by MTT assay and direct count. The transient ER-α silencing was performed using two ER-α siRNA (50 nM) and verified by Western blot. MDA-MB-231 cells were used as a negative control. Mitotane showed a similar docking configuration to 17ß-estradiol and bisphenol A (BPA) and a significant binding affinity to ER-α. MD simulations showed that mitotane preserves the active conformation of ER-α more than both BPA and Bisphenol C, classifying it as an agonist. Exposure of MCF-7 cells to mitotane led to the concentration-dependent increase of cell viability and proliferation, which was reduced in the presence of tamoxifen and nullified by the transient ER-α knock-down. Integrating bioinformatics approaches with cell biology and pharmacological methods, we demonstrated that mitotane directly binds and activates ER-α.

Sodium bis(2-ethylhexyl)sulfosuccinate self-aggregation in vacuo: molecular dynamics simulation.

Molecular dynamics (MD) simulations were conducted for systems in vacuo consisting of n AOT(-) anions (bis(2-ethylhexyl)sulfosuccinate ions) and n+/- 1 or n Na(+) ions up to n = 20. For n = 15, positively charged systems with Li(+), K(+), and Cs(+) cations were also considered. All systems were observed to form reverse micelle-like aggregates whose centre is occupied by cations and polar heads in a very compact solid-like way, while globally the aggregate has the form of an elongated and rather flat ellipsoid. Various types of statistical analyses were carried out on the systems to enlighten structural and dynamical properties including gyration radius, atomic pair correlation functions, atomic B-factor and moment of inertia tensor. For completeness and comparison the stability of reverse micelle is tested in the case of neutral n = 20 system in CCl(4) solution.

Author Correction: A novel homozygous mutation in GAD1 gene described in a schizophrenic patient impairs activity and dimerization of GAD67 enzyme.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Comparative analysis of IR and vibrational circular dichroism spectra for a series of camphor-related molecules.

The absorption spectra and vibrational circular dichroism (VCD) spectra in the mid-IR range 1600-950 cm(-1) of 10 camphor-related compounds have been recorded and compared to DFT calculated spectra at the B3PW91/TZ2P level and have been examined together with the corresponding data of the parent molecules. The rigidity of the bridged structure common to all compounds investigated permits (a) identification of three spectroscopic regions in the mid-IR range that can be "used" separately by the interested stereochemist for structural diagnosis and assignment of some major characteristics of the VCD spectra in these regions to what we call "skeletal chiral sense" and (b) recognition of possible conformers for flexible substituent groups, when present. VCD spectra of the 10 molecules have been recorded and analyzed also in the CH-stretching region, 3100-2800 cm(-1). Here, we have been able to identify and characterize features of vibrational excitons by comparison of data within the 10-molecule class. To find a theoretical justification of result (a), we have examined the potential energy distribution of the normal modes in the mid-IR range, the partitioning of the calculated rotational strengths in terms of contributions from all couples of internal coordinates, the angle formed by the two vectors, the electric dipole transition moment and the magnetic dipole transition moment, and finally the overlap of normal modes of different molecules. A discussion is provided as to the usability of the introduced algorithms.

NIR-VCD, vibrational circular dichroism in the near-infrared: experiments, theory and calculations.

The first well documented experiments of Near Infrared Vibrational Circular Dichroism (NIR-VCD) were performed around 1975. We review the thirty year history of NIR-VCD, encompassing both instrumental development and theoretical/computational methods that allow interpretation of experimental spectra, harvesting useful structural information therefrom. We hope to stimulate interest in this still scarcely explored spectroscopy of chiral molecules.

Characterisation of a type II functionally-deficient variant of alpha-1-antitrypsin discovered in the general population.

Lung disease in alpha-1-antitrypsin deficiency (AATD) results from dysregulated proteolytic activity, mainly by neutrophil elastase (HNE), in the lung parenchyma. This is the result of a substantial reduction of circulating alpha-1-antitrypsin (AAT) and the presence in the plasma of inactive polymers of AAT. Moreover, some AAT mutants have reduced intrinsic activity toward HNE, as demonstrated for the common Z mutant, as well as for other rarer variants. Here we report the identification and characterisation of the novel AAT reactive centre loop variant Gly349Arg (p.G373R) present in the ExAC database. This AAT variant is secreted at normal levels in cellular models of AATD but shows a severe reduction in anti-HNE activity. Biochemical and molecular dynamics studies suggest it exhibits unfavourable RCL presentation to cognate proteases and compromised insertion of the RCL into ß-sheet A. Identification of a fully dysfunctional AAT mutant that does not show a secretory defect underlines the importance of accurate genotyping of patients with pulmonary AATD manifestations regardless of the presence of normal levels of AAT in the circulation. This subtype of disease is reminiscent of dysfunctional phenotypes in anti-thrombin and C1-inibitor deficiencies so, accordingly, we classify this variant as the first pure functionally-deficient (type II) AATD mutant.

Molecular dynamics simulation of aqueous solutions of 26-unit segments of p(NIPAAm) and of p(NIPAAm) "doped" with amino acid based comonomers.

We have performed 75-ns molecular dynamics (MD) simulations of aqueous solutions of a 26-unit NIPAAm oligomer at two temperatures, 302 and 315 K, below and above the experimentally determined lower critical solution temperature (LCST) of p(NIPAAm). We have been able to show that at 315 K the oligomer assumes a compact form, while it keeps a more extended form at 302 K. A similar behavior has been demonstrated for a similar NIPAAm oligomer, where two units had been substituted by methacryloyl-l-valine (MAVA) comonomers, one of them being charged and one neutral. For another analogous oligomer, where the same units had been substituted by methacryloyl-l-leucine (MALEU) comonomers, no transition from the extended to the more compact conformation has been found within the same simulation time. Statistical analysis of the trajectories indicates that this transition is related to the dynamics of the oligomer backbone, and to the formation of intramolecular hydrogen bonds and water-bridges between distant units of the solute. In the MAVA case, we have also evidenced an important role of the neutral MAVA comonomer in stabilizing the compact coiled structure. In the MALEU case, the corresponding comonomer is not equally efficacious and, possibly, is even hindering the readjustment of the oligomer backbone. Finally the self-diffusion coefficient of water molecules surrounding the oligomers at the two temperatures for selected relevant times is observed to characteristically depend on the distance from the solute molecules.

Differential Enzymatic Activity of Rat ADAR2 Splicing Variants Is Due to Altered Capability to Interact with RNA in the Deaminase Domain.

In mammals, adenosine (A) to inosine (I) RNA editing is performed by adenosine deaminases acting on RNA (ADAR), ADAR1 and ADAR2 enzymes, encoded by mRNAs that might undergo splicing process. In rat, two splicing events produce several isoforms of ADAR2, called ADAR2a, ADAR2b, ADAR2e, and ADAR2f, but only ADAR2a and ADAR2b are translated into an active protein. In particular, they differ for ten amino acids located in the catalytic domain of ADAR2b. Here, we focused on these two isoforms, analyzing the splicing pattern and their different function during rat neuronal maturation. We found an increase of editing levels in cortical neurons overexpressing ADAR2a compared to those overexpressing ADAR2b. These results indicate ADAR2a isoform as the most active one, as reported for the homologous human short variant. Furthermore, we showed that the differential editing activity is not due to a different dimerization of the two isoforms; it seems to be linked to the ten amino acids loop of ADAR2b that might interfere with RNA binding, occupying the space volume in which the RNA should be present in case of binding. These data might shed light on the complexity of ADAR2 regulations.

A novel homozygous mutation in GAD1 gene described in a schizophrenic patient impairs activity and dimerization of GAD67 enzyme.

Recently, by whole exome sequencing of schizophrenia (SCZ) patients, we identified a subject that was homozygous for a novel missense substitution (c.391 A > G) in the glutamate acid decarboxylase 1 (GAD1) gene. GAD1 encodes for GAD67 enzyme, catalyzing the production of gamma-aminobutyric acid (GABA) from L-glutamic acid. Here, we studied the impact of this mutation on GAD67 activity, dimerization and subcellular localization. Biochemical assay revealed that c.391 A > G reduces GAD67 enzymatic activity by ~30%, probably due to the impaired homodimerization of homozygous mutants as highlighted by proximity ligation assays. The mutational screening of 120 genes of the "GABAergic system" in a cohort of 4,225 SCZ cases and 5,834 controls (dbGaP: phs000473.v1.p2), did not identify other cases that were homozygous for ultra-rare variants in GAD1, but highlighted an increased frequency of cases that were homozygous for rare variants in genes of the GABA system (SCZ: 0.14% vs. Controls: 0.00%; p-value = 0.0055). In conclusion, this study demonstrates the functional impact of c.391 A > G variant and its biological effect makes it a good candidate as risk variant for SCZ. This study also supports an involvement of ultra-rare variants in GABAergic genes in the etiopathogenesis of SCZ.

The pathological Trento variant of alpha-1-antitrypsin (E75V) shows nonclassical behaviour during polymerization.

Severe alpha-1-antitrypsin deficiency (AATD) is most frequently associated with the alpha-1-antitrypsin (AAT) Z variant (E342K). ZZ homozygotes exhibit accumulation of AAT as polymers in the endoplasmic reticulum of hepatocytes. This protein deposition can lead to liver disease, with the resulting low circulating levels of AAT predisposing to early-onset emphysema due to dysregulation of elastinolytic activity in the lungs. An increasing number of rare AAT alleles have been identified in patients with severe AATD, typically in combination with the Z allele. Here we report a new mutation (E75V) in a patient with severe plasma deficiency, which we designate Trento. In contrast to the Z mutant, Trento AAT was secreted efficiently when expressed in cellular models but showed compromised conformational stability. Polyacrylamide gel electrophoresis (PAGE) and ELISA-based analyses of the secreted protein revealed the presence of oligomeric species with electrophoretic and immunorecognition profiles different from those of Z and S (E264V) AAT polymers, including reduced recognition by conformational monoclonal antibodies 2C1 and 4B12. This altered recognition was not due to direct effects on the epitope of the 2C1 monoclonal antibody which we localized between helices E and F. Structural analyses indicate the likely basis for polymer formation is the loss of a highly conserved stabilizing interaction between helix C and the posthelix I loop. These results highlight this region as important for maintaining native state stability and, when compromised, results in the formation of pathological polymers that are different from those produced by Z and S AAT.

Experimental and ab-initio calculated vcd spectra of the first OH-stretching overtone of (1R)-(-) and (1S)-(+)-endo-BORNEOL.

The near infrared (NIR) absorption and NIR-vibrational circular dichroism (NIR-VCD) spectra of dilute solutions of the two enantiomers of endo-borneol have been measured in the first OH-stretching overtone region (1600-1300 nm). By density functional theory (DFT) we calculate mechanical parameters, i.e. the harmonic mechanical frequency and the anharmonicity constant for the OH stretching, and anharmonic electrical parameters; i.e. the dependence on OH-bond length of atomic polar tensors and atomic axial tensors. We evaluate transition integrals for the calculations of rotational and dipole strengths by Morse anharmonic wavefunctions depending on mechanical harmonic frequencies and mechanical anharmonicity parameters that are calculated ab initio. Experimental and calculated spectra compare quite well and this fact allows us to associate differently signed NIR-VCD features with different conformational states of the OH-bond. Absorption features for the fundamental and for the second overtone of the OH stretching are also compared with experiment.

Harmonic and anharmonic features of IR and NIR absorption and VCD spectra of chiral 4-X-[2.2]paracyclophanes.

The vibrational absorption spectra and vibrational circular dichroism (VCD) spectra of both enantiomers of 4-X-[2.2]paracyclophanes (X = COOCD3, Cl, I) have been recorded for a few regions in the range of 900-12000 cm(-1). The analysis of the VCD spectra for the two IR regions, 900-1600 cm(-1) and 2800-3200 cm(-1), is conducted by comparing with DFT calculations of the corresponding spectra; the latter region reveals common motifs of vibrational modes for the three molecules for aliphatic CH stretching fundamentals, whereas in the mid-IR region, one is able to identify specific signatures arising from the substituent groups X. In the CH stretching region between 2900 and 2800 cm(-1), we identify and interpret a group of three IR VCD bands due to HCH bending overtone transitions in Fermi resonance with CH stretching fundamental transitions. The analysis of the NIR region between approximately 8000 and approximately 9000 cm(-1) for X = COOCD3 reveals important features of the aromatic CH stretching overtones that are of value since the aromatic CH stretching fundamentals are almost silent. The intensifying of such overtones is attributed to electrical anharmonicity terms, which are evaluated here by ab initio methods and compared with literature data.

Study of conformational properties of a biologically active peptide of fibronectin by circular dichroism, NMR and molecular dynamics simulation.

Circular dichroism (CD), and NMR spectra have been recorded and molecular dynamics (MD) simulations have been performed in water and water-trifluoroethanol (TFE) mixed solvent for a synthetic biologically active 13-amino-acid fragment of human fibronectin and two related peptides. The CD results are interpreted on the basis of statistical analyses of MD trajectories and of ensuing calculations of CD spectra based on Schellman's matrix method. It is observed that the peptide conformation is quite variable in water and loses its mobility with the addition of TFE. (1)H-NOE data were found to be consistent with the most abundant calculated conformation.