Luba Tchertanov

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme that converts vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents supplied by its redox partner through thiol–disulphide exchange reactions. The functionally related molecular complexes assembled during this process have never been described, except for a proposed de novo model of a ‘precursor’ complex of hVKORC1 associated with protein disulphide isomerase (PDI). Using numerical approaches (in silico modelling and molecular dynamics simulation), we generated alternative 3D models for each molecular complex bonded either covalently or non-covalently. These models differ in the orientation of the PDI relative to hVKORC1 and in the cysteine residue involved in forming protein–protein disulphide bonds. Based on a comparative analysis of these models’ shape, folding, and conformational dynamics, the most probable putative complexes, mimicking the ‘precursor’, ‘intermediate’, and ‘successor’ states, were suggested. In addition, we propose using these complexes to develop the ‘allo-network drugs’ necessary for treating blood diseases.

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme transforming vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents delivered by its redox partner through thiol-disulfide exchange reactions. The luminal loop (L-loop) is the principal mediator of hVKORC1 activation, and it is a region frequently harbouring numerous missense mutations. Four L-loop hVKORC1 mutants, suggested in vitro as either resistant (A41S, H68Y) or completely inactive (S52W, W59R), were studied in the oxidised state by numerical approaches (in silico). The DYNASOME and POCKETOME of each mutant were characterised and compared to the native protein, recently described as a modular protein composed of the structurally stable transmembrane domain (TMD) and the intrinsically disordered L-loop, exhibiting quasi-independent dynamics. The DYNASOME of mutants revealed that L-loop missense point mutations impact not only its folding and dynamics, but also those of the TMD, highlighting a strong mutation-specific interdependence between these domains. Another consequence of the mutation-induced effects manifests in the global changes (geometric, topological, and probabilistic) of the newly detected cryptic pockets and the alternation of the recognition properties of the L-loop with its redox protein. Based on our results, we postulate that (i) intra-protein allosteric regulation and (ii) the inherent allosteric regulation and cryptic pockets of each mutant depend on its DYNASOME; and (iii) the recognition of the redox protein by hVKORC1 (INTERACTOME) depend on their DYNASOME. This multifaceted description of proteins produces “omics” data sets, crucial for understanding the physiological processes of proteins and the pathologies caused by alteration of the protein properties at various “omics” levels. Additionally, such characterisation opens novel perspectives for the development of “allo-network drugs” essential for the treatment of blood disorders.

Receptor tyrosine kinase (RTK) KIT is key regulator of cellular signalling, and its deregulation contributes to the development and progression of many serious diseases. Several mutations lead to the constitutive activation of the cytoplasmic domain of KIT, causing the aberrant intracellular signalling observed in malignant tumours. Elucidating the molecular basis of mutation-induced effects at the atomistic level is absolutely required. We report the first dynamic 3D model (DYNASOME) of the full-length cytoplasmic domain of the oncogenic mutant KIT^D816V generated through unbiased long-timescale MD simulations under conditions mimicking the natural environment of KIT. The comparison of the structural and dynamical properties of multidomain KIT^D816V with those of wild type KIT (KIT^WT) allowed us to evaluate the impact of the D816V mutation on each protein domain, including multifunctional well-ordered and intrinsically disordered (ID) regions. The two proteins were compared in terms of free energy landscape and intramolecular coupling. The increased intrinsic disorder and gain of coupling within each domain and between distant domains in KIT^D816V demonstrate its inherent self-regulated constitutive activation. The search for pockets revealed novel allosteric pockets (POCKETOME) in each protein, KIT^D816V and KIT^WT. These pockets open an avenue for the development of new highly selective allosteric modulators specific to KIT^D816V.

The kinase insert domain (KID) of RTK KIT is a key recruitment region for downstream signalling proteins (DSPs). KID, as a multisite phosphorylation region, provides alternative recognition sites for DSPs and activates them by binding a phosphotyrosine (pY) to their SH2 domains. Significant steric, biochemical, and biophysical requirements must be fulfilled by each pair of interacting proteins as the adaptation of their configurations is mandatory for the selective activation of DSPs. The accurate 3D atomistic models obtained by modelling and molecular dynamics (MD) simulations of phosphorylated KID (p-KID) have been delivered to describe KID INTERACTOME. By taking phosphorylated KID^pY721 and the N-terminal SH2 domain of phosphatidylinositol 3-kinase (PI3K), a physiological partner of KID, we showed the two proteins are intrinsically disordered. Using 3D models of both proteins, we probe alternative orientations of KID^pY721 relative to the SH2 binding pocket using automatic docking (HADDOCK) and intuitive user-guided docking. This modelling yields to two possible models of the functionally related non-covalent complex KID^pY721/SH2, where one can be regarded as the first precursor to probe PI3K activation via KIT KID. We suggest that such generation of a KID/SH2 complex is best suited for future studies of the post-transduction effects of RTK KIT.

Receptor tyrosine kinases (RTKs) are modular membrane proteins possessing both well-folded and disordered domains acting together in ligand-induced activation and regulation of post-transduction processes that tightly couple extracellular and cytoplasmic events. They ensure the fine-turning control of signal transmission by signal transduction. Deregulation of RTK KIT, including overexpression and gain of function mutations, has been detected in several human cancers. In this paper, we analysed by in silico techniques the Kinase Insert Domain (KID), a key platform of KIT transduction processes, as a generic macrocycle (KID^GC), a cleaved isolated polypeptide (KID^C), and a natively fused TKD domain (KID^D). We assumed that these KID species have similar structural and dynamic characteristics indicating the intrinsically disordered nature of this domain. This finding means that both polypeptides, cyclic KID^GC and linear KID^C, are valid models of KID integrated into the RTK KIT and will be helpful for further computational and empirical studies of post-transduction KIT events.

Receptors Tyrosine Kinases (RTKs) act as sensors for extracellular ligands, the binding of which triggers dimerization, activation, and autophosphorylation of specific tyrosine (Y) residues in the Cytoplasmic Domain (CD). This leads to the recruitment and activation of multiple downstream signaling proteins, which regulate various aspects of cellular physiology.

Human vitamin K epoxide reductase (hVKORC1) enzymatic activity requires an initial activation by a specific redox protein, a less studied step in the hVKORC1 vital cycle. Significant steric conditions must be met by enzymes, being that to adapt their configurations is mandatory for hVKORC1 activation. We studied, by molecular dynamics (MD) simulations, the folding and conformational plasticity of hVKORC1 in its inactive (fully oxidised) state using available structures, crystallographic and from de novo modelling. According to the obtained results, hVKORC1 is a modular protein composed of the stable transmembrane domain (TMD) and intrinsically disordered luminal (L) loop, possessing the great plasticity/adaptability required to perform various steps of the activation process. The docking (HADDOCK) of Protein Disulfide Isomerase (PDI) onto different hVKORC1 conformations clearly indicated that the most interpretable solutions were found on the target closed L-loop form, a prevalent conformation of hVKORC1’s oxidised state. We also suggest that the cleaved L-loop is an appropriate entity to study hVKORC1 recognition/activation by its redox protein. Additionally, the application of hVKORC1 (membrane protein) in aqueous solution is likely to prove to be very useful in practice in either in silico studies or in vitro experiments.

RTK KIT regulates a variety of crucial cellular processes via its cytoplasmic domain (CD), which is composed of the tyrosine kinase domain, crowned by the highly flexible domains—the juxtamembrane region, kinase insertion domain, and C-tail, which are key recruitment regions for downstream signalling proteins. To prepare a structural basis for the characterization of the interactions of KIT with its signalling proteins (KIT INTERACTOME), we generated the 3D model of the full-length CD attached to the transmembrane helix. This generic model of KIT in inactive state was studied by molecular dynamics simulation under conditions mimicking the natural environment of KIT. With the accurate atomistic description of the multidomain KIT dynamics, we explained its intrinsic (intra-domain) and extrinsic (inter-domain) disorder and represented the conformational assemble of KIT through free energy landscapes. Strongly coupled movements within each domain and between distant domains of KIT prove the functional interdependence of these regions, described as allosteric regulation, a phenomenon widely observed in many proteins. We suggested that KIT, in its inactive state, encodes all properties of the active protein and its post-transduction events.

The kinase insert domain (KID) of RTK KIT is the key recruitment region for downstream signalling proteins. KID, studied by molecular dynamics simulations as a cleaved polypeptide and as a native domain fused to KIT, showed intrinsic disorder represented by a set of heterogeneous conformations. The accurate atomistic models showed that the helical fold of KID is mainly sequence dependent. However, the reduced fold of the native KID suggests that its folding is allosterically controlled by the kinase domain. The tertiary structure of KID represents a compact array of highly variable α- and 3₁₀-helices linked by flexible loops playing a principal role in the conformational diversity. The helically folded KID retains a collapsed globule-like shape due to non-covalent interactions associated in a ternary hydrophobic core. The free energy landscapes constructed from first principles—the size, the measure of the average distance between the conformations, the amount of helices and the solvent-accessible surface area—describe the KID disorder through a collection of minima (wells), providing a direct evaluation of conformational ensembles. We found that the cleaved KID simulated with restricted N- and C-ends better reproduces the native KID than the isolated polypeptide. We suggest that a cyclic, generic KID would be best suited for future studies of KID f post-transduction effects.

Redox (reduction–oxidation) reactions control many important biological processes in all organisms, both prokaryotes and eukaryotes. This reaction is usually accomplished by canonical disulphide-based pathways involving a donor enzyme that reduces the oxidised cysteine residues of a target protein, resulting in the cleavage of its disulphide bonds. Focusing on human vitamin K epoxide reductase (hVKORC1) as a target and on four redoxins (protein disulphide isomerase (PDI), endoplasmic reticulum oxidoreductase (ERp18), thioredoxin-related transmembrane protein 1 (Tmx1) and thioredoxin-related transmembrane protein 4 (Tmx4)) as the most probable reducers of VKORC1, a comparative in-silico analysis that concentrates on the similarity and divergence of redoxins in their sequence, secondary and tertiary structure, dynamics, intraprotein interactions and composition of the surface exposed to the target is provided. Similarly, hVKORC1 is analysed in its native state, where two pairs of cysteine residues are covalently linked, forming two disulphide bridges, as a target for Trx-fold proteins. Such analysis is used to derive the putative recognition/binding sites on each isolated protein, and PDI is suggested as the most probable hVKORC1 partner. By probing the alternative orientation of PDI with respect to hVKORC1, the functionally related noncovalent complex formed by hVKORC1 and PDI was found, which is proposed to be a first precursor to probe thiol–disulphide exchange reactions between PDI and hVKORC1.