Novel insights into congenital surfactant dysfunction disorders by in silico analysis of ABCA3 proteins

Guo-Liang Xiao·Yuan Gao·Hu Hao·Tao Wei·Chun Hong·Yue Wang·Ying-Yi Lin·Xiu-Fang Chi·Ying Liu·Hong-Yi Gao·Chuan Nie

Surfactants produced by type II alveolar epithelial cells(AT2 cells) are usually present in inclusion organelles called lamellar bodies (LBs).The ATP-binding cassette subfamily A member 3 (ABCA3) transporter primarily exists in AT2 cells and is generally considered to be one of the critical regulators of biogenesis of LBs and surfactant metabolism in the lungs [1—4].ABCA3mutations are the most common cause of congenital surfactant dysfunction disorders(CSDDs),resulting in fatal neonatal respiratory distress and pediatric or adult interstitial lung disease [3,5—7].More than 200 disease-associatedABCA3variants have been identified in symptomatic infants and children [8].The estimated prevalence of deleteriousABCA3mutations in the population is 1/70—1/33,with a predicted disease incidence of 1/20,000—1/4400 [7].

CSDD caused byABCA3mutations lacks specific therapies,and patients are mostly subjected to ineffective conventional treatments [9].A combination of hydroxychloroquine,corticosteroids,and azithromycin may improve symptoms and delay disease progression in some cases [10—14];however,evidence-based guidelines are currently lacking.Lung transplantation appears to be a promising way to treat respiratory distress syndrome due to ABCA3 defects [8,15].Nevertheless,finding a suitable lung donor in the neonatal period is a considerable challenge for patients.Thus,gene therapy may be a new method for treating these diseases [3].Therefore,it is crucial to understand the genetic characteristics and pathogenic mechanisms of ABCA3.

The mechanisms underlying compound heterozygosity remain unclear.Since compound heterozygosity mutations play an essential role in the currently reported cases,additional biophysical analyses are needed to predict the effect of gene mutations on the structure and function of the ABCA3 protein,such as homology modeling or 3D model construction for visualization of the target gene sequence.

In general,comprehensive functional analysis,by studying how mutations affect protein variant functions,is essential for understanding the molecular pathogenesis of CSDD and interpreting their clinical significance.This study reports a novel compound heterozygous mutation caused by two missense mutations.In silico analysis elucidated the effects of these novel compound heterozygous variants on protein function using 3D structures.Furthermore,computational analysis of the wild-type and variant proteins revealed the harmful nature of these mutations.The findings of this study can help increase the understanding of compound heterozygous mutations inABCA3.In addition,through impact and domain function analysis before and after amino acid mutation,we can better understand how amino acid substitutions contribute to disease occurrence.

The male neonate was born to a 30-year-old mother(gravida 1,para 1) by vaginal delivery at the 41st week of gestation with 3300 g birth weight.The Apgar scores were 9,10,and 10 at 1,5,and 10 minutes,respectively.Two hours after birth,the infant presented with a weak cry,shallow breathing,hypotonia,and cyanosis,and oxygen saturation decreased to 65%—72%.After mechanical ventilation with 60% oxygen and a mean airway pressure of 15 mmHg,oxygen saturation increased to 90%.Cardiac ultrasonography showed an atrial septal defect (4.0 mm),patent ductus arteriosus (4.2 mm),and persistent pulmonary hypertension (75 mmHg).Arterial blood gas showed pH 7.281,pO226.2 mmHg,pCO267.4 mmHg,and HCO3-25.8 mmol/L.Nitric oxide (NO) inhalation improved persistent pulmonary hypertension,and the oxygen requirement decreased to 40%with SpO292%.Consecutive chest radiographs revealed diffuse lung disease similar to neonatal hyaline membrane disease (Supplementary Fig.1).In addition,chest X-ray examination did not provide valuable information to assist in the diagnosis.Computed tomography revealed blurred bronchovascular bundles and decreased transparency in both lungs.In addition,we detected a large patchy highdensity shadow in the anterior segment of the left upper lobe of the right lung and on the lower lobe of both lungs(Supplementary Fig.2).Pulmonary hypertension improved to 35 mmHg after NO inhalation,according to the cardiac ultrasonography findings on day 3.However,hypoxemia and carbon dioxide retention persisted after ventilation with high-frequency oscillatory ventilation or volume-guaranteed pressure controlled (tide volume 6 mL/kg).High-throughput genetic sequencing technology detectedUreaplasmainfection in the blood and sputum,which was cleared one week after azithromycin treatment.

Fig.1 ATP-binding cassette A3 (ABCA3) gene mutation analysis. a Family pedigree.A novel ABCA3 compound heterozygosity mutation inherited from the parents. b Gene analysis of ABCA3 in the infant and his parents.Infant: T-to-G substitution in exon 11(c.1142T ＞ G);infant’s mother: A-to-G heterozygous mutation in exon 11(c.1142T ＞ G);infant’s father: no mutation in exon 11. c Gene analysis of ABCA3 in the infant and his parents.Infant: G-to-T substitution in exon 8 (c.731G ＞ T);infant’s father: G-to-T substitution in exon 8 (c.731G ＞ T);infant’s mother: no mutation in exon 8

Fig.2 Modelling of AIP-binding cassette A3 (ABCA3) wild type and mutant-type protein. a,b Three-dimension (3D) structure of the ABCA3 protein.Residue 381 was marked in blue,and the structure of mutant L381R was marked in green. c,d 3D structure showing before and after mutation of residue 381 (L381R) in cartoon format.Wild type was labeled in blue,and the mutant was labeled in green;the mutant L381R protein did not reveal any structural changes on a smaller scale as well as compared to that for the wild type. e,f 3D structure showing before and after mutation of residue 244 (R244M) in cartoon format.Wild type was labeled in blue,and the mutant was labeled in green.The 3D structure in the cartoon format with a smaller design indicated that it was located at the end of a β-sheet inside the protein

However,the infant’s condition did not improve significantly.The peak pressure of respiratory support must consistently be maintained above 22 cm H2O to maintain a tidal volume of 6 mL/kg,and carbon dioxide retention was not markedly improved.Therefore,whole-exome sequencing (WES) and left lung biopsy via video-assisted thoracoscopic surgery (VATS) were performed on day 21,with the parents’ informed consent and ethical review,to collect lung tissue for pathological examination.Uneven inflation and slightly scattered white granules appeared on the surface of the left lung during the operation (Supplementary Fig.3).Evidence-based guidelines for hydroxychloroquine use are currently lacking,and we were unable to obtain parental consent for hydroxychloroquine administration.The infant died of respiratory failure at the age of 40 days,and his parents declined autopsy for him.

Fig.3 Positions of the structural domains on the ABCA3 protein,where the two mutations are located. a Positions of the structural domains on the ABCA3 protein,where the two mutations are located are marked in cartoon format in blue and red for the ECD1 structural domain (blue) and TM5 structural domain (red),respectively.To clearly understand the position of the mutations on the structural domains,we have locally zoomed in on the two mutated amino acids,and the mutated amino acids are indicated in green. b Location of ABCA3 in vivo and its possible distribution across the membrane.The structural domains,where the two mutations are located are shown in cartoon format,where ECD1 is shown in blue,TM5 is shown in red,and the other unmutated structural domains are shown in silver surface format. ABCA3 ATP-binding cassette A3,ECD1 extracellular domain 1,TM transmembrane

The infant underwent VATS according to standard protocols for collecting lung tissues.Specimens from the upper and lower lobes of the left lung were formalin-fixed,paraffin-embedded and investigated by conventional histology and immunohistochemistry.

We collected infant and parent blood samples for WES.Genomic DNA from peripheral blood leukocytes derived from the proband was extracted using the QIAamp DNA blood midi kit (Qiagen,Germany).First,1 μg of each genomic DNA sample was fragmented by sonication and purified to yield 200—300 bp fragments.Second,paired-end adapter oligonucleotides from Illumina (San Diego,California,USA) were ligated into shared genomic DNA.Five hundred nanograms of these tailed fragments were then hybridized to the probe library of SureSelect human all Exon V6(Agilent,Germany).Finally,the enrichment libraries were sequenced on an Illumina NovaSeq 6000 sequencer (Illumina) as 150 bp paired-end reads.

After sequencing,reads were aligned to the human reference genome (GRCh37/hg19) with the Burrows-Wheeler aligner [16],and potential duplicate paired-end reads were removed using the genome analysis toolkit (GATK) v.4.2.0.0(https://github.com/broadinstitute/gatk/releases/tag/4.2.2.0).GATK v.4.2.0.0 was used to obtain base quality score recalibration,indel realignment,single-nucleotide variant,and indel discovery.Genotyping was performed using standard hard filtering parameters [17].

Low-quality variants were f lagged and excluded from subsequent analyses.Bamdst v.1.0.9 (https://github.com/shiquan/bamdst) was used to assess the clean data coverage of each sample using default parameters.All variants identified in the affected individuals were annotated with databases,such as refGene (https://www.refgene.com),Avsnp150 (https://www.ncbi.nlm.nih.gov/snp/),gnomAD211 (http://gnomad-sg.org/),ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/),dbnsfp41a (https://sites.google.com/site/jpopgen/dbNSFP),Intervar (https://wintervar.wglab.org/) by snpeff5.0d ( https://sourceforge.net/projects/snpeff/),and annovar 2020 Jun (https://annovar.openbioinformatics.org/en/latest/).Candidate mutational events were then inspected using an integrative genomics viewer[18].The resulting variants were excluded when their frequency exceeded 1/100 in Genome Aggregation (GnomAD).Variants were correlated with patient phenotypes and the results of clinical investigation.All variants were classified according to the American College of Medical Genetics and Genomics guidelines (ACMG) standards and guidelines.

The structural model of wild-type ABCA3 was obtained from the UniProt database (https://www.unipr ot.org/niport/Q99758).Homologous modeling using SWISS-model(https://swiss model.expasy.org/inter active) and AlphaFold2(https://github.com/deepm ind/alpha fold2) assisted by Google’s Colab (https://colab.research.google.com) arithmetic platform were applied to obtain a relatively accurate model.PyMol (https://pymol.org/2/) was used to label relevant residues.Because of the internal storage size limitation(RAM=13 GB),we built only a partial model.

MutPred2 (http://mutpred.mutdb.org/) is a machine learning-based method and software package [19] that integrates genetic and molecular data to probabilistically infer the pathogenicity of amino acid substitutions.It is a sequence-based model that utilizes a methodology based on recent machine learning advances.It is trained based on unlabeled positive data and combines prior and posterior probabilities.This is achieved by providing general pathogenicity predictions and ranked lists of specif ic molecular alterations that may affect the phenotype.In this study,MutPred2 was used to predict and analyze the structure—function relationship of the ABCA3 protein before and after mutation of the two identified amino acid sites.

Histological examination using hematoxylin—eosin(H&E) staining revealed collapse and enlargement of alveoli,increased fibers in the alveolar basement membrane,serous exudate and tissue reaction in the alveoli cavity,angiodysplasia,and inflammatory cell infiltration (Supplementary Fig.4a—d).Immunohistochemical examination showed strong cytokeratin (CK) expression in alveolar epithelial cells but weak cluster of differentiation (CD) 31 and CD34 expression in pulmonary vascular endothelial cells,which indicated pulmonary vascular dysplasia and widened alveolar septa.In addition,the expression of CD68 in alveolar monocytes revealed inflammatory cell infiltration.We also observed alveolar fibrin deposition by Masson and periodic acid-Schiff (PAS) staining (Supplementary Fig.4e—j).

Two novelABCA3mutations were found in the compound heterozygosis of the infant.Heterozygous variations in theABCA3gene,c.1142T ＞ G and c.731G ＞ T,have not been previously reported and are unavailable in the GnomAD (http://gnomad.broadinstitute.org) and Exome Aggregation Consortium (ExAC) (https://exac.broadinstitute.org) databases.The c.1142T ＞ G variant is located in exon 11 ofABCA3,resulting in p.Leu381Arg.The c.731G ＞ T variant was located in exon 8 ofABCA3,resulting in p.Arg244Met.Next-generation sequencing results revealed compound heterozygous mutations of theABCA3gene,with the c.1142T ＞ G mutation originating from the mother and the c.731G ＞ T mutation inherited from the father (Fig. 1).However,the parents had no signs of respiratory distress,even though they carried a distinct mutation inABCA3.Therefore,according to the ACMG[20],these two mutations were classified as variants of uncertain clinical significance.

Owing to the limitation of the hash rate,the confidence of the calculated model constructed by SWISS modeling is low (global model quality estimate=54%).Even with the help of Google’s Colab,only 67.3% and 70.58% of residues were involved in the structural modeling of p.Leu381Arg (L381R) and p.Arg244Met (R244M) by AlphaFold2 and SWISS modeling,respectively.

The structure of ABCA3 predicted using the online tool AlphaFold2 indicated that residue 381 was located on the surface of the protein (Fig. 2 a,b).The 3D structure in the cartoon format showed that residue 381 was in the middle of an α-helix (Fig. 2 c,d).The mutant L381R did not differ substantially at the structural level compared to the wild type.The structure of ABCA3 showed that two caps covered the head and tail of the protein.Residue 244 was not observed on the surface of ABCA3 (Fig. 2 a,b).The 3D structure in the cartoon format with a smaller design indicated that it was located at the end of a β-sheet inside the protein (Fig. 2 e,f).

By reviewing relevant data [21],based on the structural model and delineation of functional regions,we localizedand functionally analyzed the structural domains,where the two mutations in this study were located (Fig. 3 a).R244M was found in the extracellular domain 1 (ECD1) region,an extracellular structural domain with varying lengths among different members of the ABC family and with associated instability.In contrast,the L381R mutation was primarily located in the transmembrane TM 5 part of transmembrane domain 1 (TMD1),a transmembrane structural domain that is in contact with the cytoplasmic end.At the same time,the TM5 portion of the TMD1 structural domain forms a sizablehydrophobic cavity together with other transmembrane proteins.ABCA3 is primarily distributed on the boundary membrane of the lamellipodia in AT2 cells.To envision its distribution location and transmembrane mode,we used Biorender (https://biorender.com) to draw a 2D distribution schematic of the cell membrane (Fig. 3 b).

We predicted the pathogenicity of the two loci before and after mutation using MutPred2.Relevant analytical results were obtained by introducing the mutation of two loci,R244M and L381R,with the complete amino acid sequence of ABCA3 in the calculation of MutPred2.The results were interpreted in conjunction with the previous delineation of the function and position of the structural domains.

The general pathogenicity score of R244M was 0.751,which is above the pathogenicity threshold of 0.5.This result indicated that the R244M substitution has great potential to produce relevant pathogenic effects on the organism(Table 1).The possible effects were divided into three main areas (in the order of likelihood of occurrence calculated by MutPred2): (1) alteration of the originally ordered interface (Pvalue=8.8×10-3,possibility 33%);(2) gain of loop (Pvalue=4.9×10-3,probability 30%);(3) alteration of transmembrane protein structure and function (Pvalue=1.0×10-3,probability 26%).

Table 1 Prediction results of R244M in MutPred2

Mutations in the L381R locus were also identified in this study.The general pathogenicity score of L381R was 0.808,which was higher than the pathogenicity threshold of 0.5.In terms of molecular mechanism changes,the L381R mutation also has numerous implications (Table 2) divided into two main aspects (in order of likelihood of occurrence).The first was the alteration of transmembrane protein structure and function (Pvalue=4.5×10-3,21% probability).The overall structure of the ABCA3 protein L381R before and after the mutation was generated using a self-optimized prediction method(SOPMA) with alignment,online software for computational analysis of protein secondary structure (Table 3).Before the mutation,the percentages of α-helix,β-turn,and random coil were 46.01%,5.69%,and 32.16%,respectively.After the mutation,the percentage of α-helices decreased by 0.12%,β-turns decreased by 0.06%,and random coils increased by 0.18%,45.89%,5.63%,and 32.34%,respectively.The second effect in order of occurrence likelihood is the sulfation deletion of complex 386 (Tyr) (Pvalue=0.05,probability 1%).

Table 2 Prediction results of L381R in MutPred2

Table 3 Statistics of the percentage of secondary structure of proteins that have changed of L381R in self-optimized prediction method(SOPMA)

ABCA3mutations are the most common cause of genetic surfactant dysfunction disorders,resulting in loss of functionof phospholipid transporters involved in pulmonary surfactant function,leading to fatal neonatal respiratory distress [8,15,22].ABCA3mutations include multiple forms,including missense mutations,splice sites,insertions,and deletions [23].These mutations are autosomal recessive and result in distinct clinical phenotypes [8].Although studies have discovered more than 200 diseases associated withABCA3variants in symptomatic infants and children,approximately ? of the reported cases of pathogenicABCA3mutations are associated with compound heterozygous mutations [3].An overview of some previously reportedABCA3mutations is presented in Table 4.Although the mechanisms underlying this phenomenon have not been well-defined [3,9,23,24,25],whether these mutations are dominant-negative remains to be elucidated.It has been suggested that exposure to certain harmful environmental factors (e.g.,smoking and viral infection) is an important factor in the long-term survival of patients with these compound heterozygous mutants as well as single mutations [3].

The lung tissue histological results of the infant were similar to pathological alterations previously reported regardingABCA3mutations homozygous or compound heterozygous[9].In addition,these characteristic pathological changes reduce the diffusional barrier and ventilation—perfusion mismatch,resulting in persistent hypoxemia and refractory carbon dioxide retention.

Whole-exome analysis revealed a compound heterozygote consisting of twoABCA3variants.These two variants are located in exons 11 (c.1142T ＞ G) and 8 (c.731G ＞ T) and formed a compound heterozygous (trans) relationship.Both were missense mutations in the parents.Transformation is rare and has not been previously reported.The GnomAD,ExAC,human gene mutation database,and ClinVar databases do not contain these two variations.The c.1142T ＞ G variant is located in exon 11 of ABCA3,resulting in p.

Leu381Arg (L381R).The c.731G ＞ T variant is located in exon 8 of ABCA3,resulting in p.Arg244Met (R244M).

In our three-dimention (3D) models,the R244M mutation was located in the ECD1 region,which is an extracellular structural domain.Previous research has shown that missense mutations in the ECD1 region ofABCA3lead to the accumulation of ABCA3 in the endoplasmic reticulum,resulting in defective protein synthesis [29].In contrast,the L381R mutation was located in the TM5 region of the TMD1 domain that,together with other TM proteins,forms a sizable hydrophobic cavity connecting the external environment to the intracellular membrane.Substrates transferred by ABCA3 proteins use these conduits to cross the cell membrane.ECD loops and TMDs are the leading substrate-binding sites that allow the trafficking of lipid molecules [30].

Calculated using the MutPred2 algorithm,the alteration of R244M disrupted the original extracellular membrane orderly interface,where ECD1 was located.The mutation may affect the activity of this ECD1 region in the ABCA3 protein,which contains lipid molecule-binding sites.Moreover,the R244M mutation was speculated to affect the structure of the TM2 transmembrane protein,resulting in alteration of the structural part of the channel occupied by this protein.

The L381R mutation was found in TM5,a partial structural domain of TMD1.TMD1 is in contact with the cytoplasmic end and forms a hydrophobic cavity across the membrane and other structures.Therefore,we hypothesize that the changes in extended strand and β-turn occupancy affected the space of the transmembrane channel structure,and the transformation of the transmembrane channel structure space directly affected the functional ABCA3 protein shift,which may lead to disease.

Based on MutPred2 predictions,we believe that the complex heterozygous mutation located in the structural domains of ECD1 and TMD1 ultimately weakens the transport of substrate lipid molecules,thus reducing the performance of surfactants and eventually leading to CSDD.

This study has some limitations.For example,we did not perform an electron microscopy examination to study LBs and their secretion,and we could not collect samples of the three-generation pedigree and conduct protein functional validation.Furthermore,due to the hash rate limitation,not all residues were involved in the structural modeling of the variants.In addition,the mutation results inferred by the MutPred2 and SOPMA algorithms were not certain to occur but instead represented the likelihood of the effect occurrence.Further research with X-ray diffraction or cryoelectron microscopy is needed.

In conclusion,the novel ABCA3 compound heterozygous mutation resulted in structural changes in the TMD1 and ECD1 regions of ABCA3.They may ultimately weaken the transport of lipid molecules,thus reducing the performance of surfactants and eventually leading to CSDD.The development of genotype—phenotype relationships is often unpredictable in genetic diseases,such as CSDD,caused by compound heterozygous mutations.Therefore,bioinformatics can play an essential role in this class of fields.

Supplementary InformationThe online version contains supplementary material available at https://doi.org/10.1007/s12519-022-00645-y.

AcknowledgementsThe authors gratefully acknowledge Guangzhou Daan Clinical Laboratory Center for their technical assistance.

Author contributionsXGL: conceptualization,methodology,formal analysis,writing—original draft.NC: conceptualization,methodology,formal analysis,writing—review and editing.GY: methodology,data curation,investigation,writing—original draft.HH: methodology.WT and GHY: methodology,data curation,investigation,formal analysis,writing—original draft.HC,WY,LYY,CXF and LY: methodology.

FundingThis work was supported by the Guangzhou Science and Technology Program (201804010090),Natural Science Foundation of Guangdong Province (2018A030313605),and Medical Scientific Research Foundation of Guangdong Province (No.B2020128).

Data availabilityThe data sets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Declarations

Conflict of interestNo financial or non-financial benefits have been received or will be received from any party related directly or indirectly to the subject of this article.

Ethical approvalThe study was approved by the Ethics Committee of Guangdong Women and Children’s Hospital (202201019).

Consent for publicationBefore the genetic studies,informed consent was obtained from the parents.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License,which permits use,sharing,adaptation,distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence,and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence,unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use,you will need to obtain permission directly from the copyright holder.To view a copy of this licence,visit http://creativecommons.org/licenses/by/4.0/.