Computational Characterization of Pathogenic LMNA Missense Variants: Structural Instability, Altered Binding, and Conformational Dynamics.

BACKGROUND

Mutations in the LMNA gene underlie a broad spectrum of laminopathies, including muscular dystrophies, cardiomyopathies, and premature aging syndromes; however, the molecular mechanisms by which missense variants disrupt Lamin A structural integrity remain incompletely characterized. Systematic computational approaches for prioritizing pathogenic variants and elucidating their structural consequences are critically needed.

METHODS

An integrated multistep in silico framework was employed to investigate the structural and functional consequences of LMNA missense variants. Variant prioritization was performed using the Evo2 nucleotide language model via delta log-likelihood scoring, followed by bioinformatic annotation using SIFT, PANTHER-PSEP, PhD-SNP, and E-SNPs&GO. Protein stability assessment was conducted with DynaMut, INPS-MD, I-Mutant2.0, and MUpro. Variants localized within globular domains-N456D, N456T, and G465D-together with the known pathogenic variant M540T as a positive control, were selected for three-dimensional structural modeling using PyMOL and AlphaFold2, molecular docking with lonafarnib as a reference ligand via AutoDock Vina, and 100 ns molecular dynamics simulations using GROMACS with the Amber ff14SB force field. Conformational dynamics were characterized through principal component analysis and free-energy surface construction.

RESULTS

Evo2-based screening of the full LMNA coding sequence identified 50 high-priority loss-of-function variants, of which N456D, N456T, and G465D were retained for structural investigation based on their globular domain localization and multitool pathogenicity predictions. All three variants were consistently predicted to alter physicochemical properties and reduce structural stability relative to wild-type Lamin A. Molecular docking revealed mutation-dependent changes in lonafarnib binding profiles. The known pathogenic control M540T exhibited comparable structural and dynamic behavior, supporting the reliability of the prioritization workflow. Molecular dynamics analyses demonstrated altered RMSD trajectories, increased residue-level flexibility, and modified hydrogen bonding patterns in mutant systems. Free-energy landscape analyses revealed expanded conformational basins, particularly pronounced in the G465D variant, indicating increased structural plasticity.

CONCLUSION

This integrated computational framework provides a systematic strategy for prioritizing pathogenic LMNA variants and characterizing their structural consequences at the atomic level. The identified variants-N456D, N456T, and G465D-represent structurally disruptive substitutions consistent with the established role of globular domain destabilization in other laminopathy-associated variants, offering testable hypotheses for experimental validation in cellular and animal models.