This study explored the association of deleterious variants in pharmacodynamics (PD) genes with statin response and adverse effects in patients with familial hypercholesterolemia (FH) and analyzed their potential effects on protein structure and stability. Clinical and laboratory data were obtained from 144 adult FH patients treated with statins. A panel of 32 PD genes was analyzed by exon-targeted gene sequencing. Deleterious variants were identified using prediction algorithms and their structural effects were analyzed by molecular modeling studies. A total of 102 variants were predicted as deleterious (83 missense, 8 stop-gain, 4 frameshift, 1 indel, 6 splicing). The variants ABCA1 rs769705621 (indel), LPA rs41267807 (p.Tyr2023Cys) and KIF6 rs20455 (p.Trp719Arg) were associated with reduced low-density lipoprotein cholesterol (LDLc) response to statins, and the LPL rs1801177 (p.Asp36Asn) with increased LDLc response (P < 0.05). LPA rs3124784 (p.Arg2016Cys) was predicted to increase statin response (P = 0.022), and ABCA1 rs769705621 to increase the risk of statin-related adverse events (SRAE) (P = 0.027). LPA p.Arg2016Cys and LPL p.Asn36Asp maintained interactions with solvent, LPA p.Tyr2023Cys reduced intramolecular interaction with Gln1987, and KIF6 p.Trp719Arg did not affect intramolecular interactions. DDMut analysis showed that LPA p.Arg2016Cys and p.Tyr2023Cys and LPL p.Asp36Asn caused energetically favorable changes, and KIF6 p.Trp719Arg resulted in unfavorable energetic changes, affecting protein stability. Deleterious variants in ABCA1, LPA, LPL and KIF6 are associated with variability in LDLc response to statins, and ABCA1 rs769705621 is associated with SRAE risk in FH patients. Molecular modeling studies suggest that LPA p.Tyr2023Cys and KIF6 p.Trp719Arg disturb protein conformational structure and stability.