In-silico prediction of TGF-β1 non-synonymous variants and their impact on binding affinity to Fresolimumab

TGF-β1 is a potent immunoregulatory cytokine that plays diverse roles in development, bone healing, fibrosis, and cancer. However, characterizing TGF-β1 gene variants is challenging because the structural and functional consequences of these variants are still undetermined. In this study, we aimed to perform an in-silico analysis of TGF-β1 non-synonymous variants and their pathogenic effects on the TGF-β1 protein. A total of 10,252 TGF-β1 SNPs were collected from the NCBI dbSNP database and in-silico tools (SIFT, PROVEAN, Mutation Taster, ClinVar, PolyPhen-2, CScape, MutPred, and ConSurf) were used. The in-silico predicted potential variants were further investigated for their binding to the TGF-β1 targeting drug “Fresolimumab”. Molecular docking was performed using HADDOCK and confirmed by PRODIGY and PDBsum. The in-silico analysis predicted four potential TGF-β1 nsSNPs: E47G in the LAP domain of the propeptide and I22T, L28F, and E35D in the mature TGF-β1 peptide. HADDOCK and molecular dynamics simulations revealed that the I22T and E35D variants have higher binding affinity for Fresolimumab as compared to the wild type and L28F variants. Molecular dynamics simulations (100 ns) and principal component analysis showed that TGF-β1 variants influenced the protein structure and caused variations in the internal dynamics of protein complexes with the antibody. Among them, the E35D variant significantly destabilized the TGF-β1 protein structure, resulting in rearrangement in the binding site and affecting the interactions with the Fresolimumab. This study identified four variants that can affect the TGF-β1 protein structure and result in functional consequences such as impaired response to Fresolimumab.

: TGF-β1 is a potent immunoregulatory cytokine that plays diverse roles in development, bone healing, fibrosis, and cancer. However, characterizing TGF-β1 gene variants is challenging because the structural and functional consequences of these variants are still undetermined. In this study, we aimed to perform an in-silico analysis of TGF-β1 non-synonymous variants and their pathogenic effects on the TGF-β1 protein. A total of 10,252 TGF-β1 SNPs were collected from the NCBI dbSNP database and in-silico tools (SIFT, PROVEAN, Mutation Taster, ClinVar, PolyPhen-2, CScape, MutPred, and ConSurf) were used. The in-silico predicted potential variants were further investigated for their binding to the TGF-β1 targeting drug "Fresolimumab". Molecular docking was performed using HADDOCK and confirmed by PRODIGY and PDBsum. The in-silico analysis predicted four potential TGF-β1 nsSNPs: E47G in the LAP domain of the propeptide and I22T, L28F, and E35D in the mature TGF-β1 peptide. HADDOCK and molecular dynamics simulations revealed that the I22T and E35D variants have higher binding affinity for Fresolimumab as compared to the wild type and L28F variants. Molecular dynamics simulations (100 ns) and principal component analysis showed that TGF-β1 variants influenced the protein structure and caused variations in the internal dynamics of protein complexes with the antibody. Among them, the E35D variant significantly destabilized the TGF-β1 protein structure, resulting in rearrangement in the binding site and affecting the interactions with the Fresolimumab. This study identified four variants that can affect the TGF-β1 protein structure and result in functional consequences such as impaired response to Fresolimumab.Communicated by Ramaswamy H. Sarma.