Overexpression and RNAi-mediated Knockdown of Two 3ß-hydroxy-Δ5-steroid dehydrogenase Genes in Digitalis lanata Shoot Cultures Reveal Their Role in Cardenolide Biosynthesis.

3ß-hydroxy-Δ5-steroid dehydrogenases (3ßHSDs) are supposed to be involved in 5ß-cardenolide biosynthesis. Here, a novel 3ßHSD (Dl3ßHSD2) was isolated from Digitalis lanata shoot cultures and expressed in E. coli. Recombinant Dl3ßHSD1 and Dl3ßHSD2 shared 70% amino acid identity, reduced various 3-oxopregnanes and oxidised 3-hydroxypregnanes, but only rDl3ßHSD2 converted small ketones and secondary alcohols efficiently. To explain these differences in substrate specificity, we established homology models using borneol dehydrogenase of Salvia rosmarinus (6zyz) as the template. Hydrophobicity and amino acid residues in the binding pocket may explain the difference in enzyme activities and substrate preferences. Compared to Dl3ßHSD1, Dl3ßHSD2 is weakly expressed in D. lanata shoots. High constitutive expression of Dl3ßHSDs was realised by Agrobacterium-mediated transfer of Dl3ßHSD genes fused to the CaMV-35S promotor into the genome of D. lanata wild type shoot cultures. Transformed shoots (35S:Dl3ßHSD1 and 35S:Dl3ßHSD2) accumulated less cardenolides than controls. The levels of reduced glutathione (GSH), which is known to inhibit cardenolide formation, were higher in the 35S:Dl3ßHSD1 lines than in the controls. In the 35S:Dl3ßHSD1 lines cardenolide levels were restored after adding of the substrate pregnane-3,20-dione in combination with buthionine-sulfoximine (BSO), an inhibitor of GSH formation. RNAi-mediated knockdown of the Dl3ßHSD1 yielded several shoot culture lines with strongly reduced cardenolide levels. In these lines, cardenolide biosynthesis was fully restored after addition of the downstream precursor pregnan-3ß-ol-20-one, whereas upstream precursors such as progesterone had no effect, indicating that no shunt pathway could overcome the Dl3ßHSD1 knockdown. These results can be taken as the first direct proof that Dl3ßHSD1 is indeed involved in 5ß-cardenolide biosynthesis.

Key interactions in the trimolecular complex consisting of the rheumatoid arthritis-associated DRB1*04:01 molecule, the major glycosylated collagen II peptide and the T-cell receptor.

OBJECTIVES: Rheumatoid arthritis (RA) is an autoimmune disease strongly associated with the major histocompatibility complex (MHC) class II allele DRB1*04:01, which encodes a protein that binds self-peptides for presentation to T cells. This study characterises the autoantigen-presenting function of DRB1*04:01 (HLA-DRA*01:01/HLA-DRB1*04:01) at a molecular level for prototypic T-cell determinants, focusing on a post-translationally modified collagen type II (Col2)-derived peptide. METHODS: The crystal structures of DRB1*04:01 molecules in complex with the peptides HSP70289-306, citrullinated CILP982-996 and galactosylated Col2259-273 were determined on cocrystallisation. T cells specific for Col2259-273 were investigated in peripheral blood mononuclear cells from patients with DRB1*04:01-positive RA by cytofluorometric detection of the activation marker CD154 on peptide stimulation and binding of fluorescent DRB1*0401/Col2259-273 tetramer complexes. The cDNAs encoding the T-cell receptor (TCR) α-chains and ß-chains were cloned from single-cell sorted tetramer-positive T cells and transferred via a lentiviral vector into TCR-deficient Jurkat 76 cells. RESULTS: The crystal structures identified peptide binding to DRB1*04:01 and potential side chain exposure to T cells. The main TCR recognition sites in Col2259-273 were lysine residues that can be galactosylated. RA T-cell responses to DRB1*04:01-presented Col2259-273 were dependent on peptide galactosylation at lysine 264. Dynamic molecular modelling of a functionally characterised Col2259-273-specific TCR complexed with DRB1*04:01/Col2259-273 provided evidence for differential allosteric T-cell recognition of glycosylated lysine 264. CONCLUSIONS: The MHC-peptide-TCR interactions elucidated in our study provide new molecular insights into recognition of a post-translationally modified RA T-cell determinant with a known dominant role in arthritogenic and tolerogenic responses in murine Col2-induced arthritis.

21-Hydroxypregnane 21-O-malonylation, a crucial step in cardenolide biosynthesis, can be achieved by substrate-promiscuous BAHD-type phenolic glucoside malonyltransferases from Arabidopsis thaliana and homolog proteins from Digitalis lanata.

Three putative 21-hydroxypregnane 21-O-malonyltransferases (21MaT) from Digitalis lanata were partially purified. Two of them were supposed to be BAHD-type enzymes. We were unable to purify them in quantities necessary for reliable sequencing. We identified two genes in A. thaliana coding for substrate-promiscuous BAHD-type phenolic glucoside malonyltransferases (AtPMaT1, AtPMaT2) and docked various 21-hydroxypregnanes into the substrate-binding site of a homology model built on the BAHD template 2XR7 (NtMaT1 from N. tabacum). Recombinant forms of Atpmat1 and Atpmat2 were expressed in E. coli and the recombinant enzymes characterized with regard to their substrate preferences. They were shown to malonylate various 21-hydroxypregnanes. The Atpmat1 sequence was used to identify candidate genes in Digitalis lanata (Dlmat1 to Dlmat4). Dlmat1 and Dlmat2 were also expressed in E. coli and shown to possess 21-hydroxypregnane 21-O-malonyltransferase activity.

Synthesis, Radiosynthesis and Biological Evaluation of Buprenorphine-Derived Phenylazocarboxamides as Novel µ-Opioid Receptor Ligands.

Targeted structural modifications have led to a novel type of buprenorphine-derived opioid receptor ligand displaying an improved selectivity profile for the µ-OR subtype. On this basis, it is shown that phenylazocarboxamides may serve as useful bioisosteric replacements for the widely occurring cinnamide units, without loss of OR binding affinity or subtype selectivity. This study further includes functional experiments pointing to weak partial agonist properties of the novel µ-OR ligands, as well as docking and metabolism experiments. Finally, the unique bifunctional character of phenylazocarboxylates, herein serving as precursors for the azocarboxamide subunit, was exploited to demonstrate the accessibility of an 18 F-fluorinated analogue.

Induction of the tetracycline repressor: characterization by molecular-dynamics simulations.

Extensive molecular-dynamics simulations show that the distance between the centers of gravity of the two equivalent helices 3 in the DNA-binding heads of the dimer of the tetracycline-repressor protein (TetR) can be used as a reliable diagnostic of induction. This is not, however, true for X-ray structures, but only for molecular-dynamics simulations. This is suggested to be because TetR is inherently flexible along the coordinate of the allosteric change (as is always likely to be the case for allosteric proteins), so that crystal-packing forces can determine the conformation of the protein. However, the time scale of the allosteric rearrangement in the absence of DNA-complexation is found to be of the order of tens of nanoseconds, so that rearrangements can be observed reproducibly in 100 ns simulations. Metastable (pre-equilibrium) conformations of TetR have been observed for up to 60 ns. The likely equilibrium processes and key features of the TetR system are discussed.

Diarylpropane-1,3-dione derivatives as TetR-inducing tetracycline mimetics: Synthesis and biological investigations.

Synthesis, biological investigations and molecular docking studies of nonantibiotic and nontetracyclic inducers that feature a minimal key motif of the natural lead tetracycline are presented. The diarylpropane-1,3-dione motif was identified as the minimal substructure responsible for TetR induction by tetracyclines. The first nontetracyclic surrogates of the natural tetracyclines displayed significant inducing effects for TetR(BD)S135L, whereby the chlorohydroxyphenyl-substituted beta-diketone 31 displayed the highest activity. Interestingly, antibiotic activity could not be detected for 31. Homology modeling based on the X-ray structure of 7-chlorotetracycline bound to TetR indicated analogous binding modes for the natural inducer and the synthetic diarylpropane-1,3-dione derivatives.

Molecular dynamics simulations on human fibulin-4 mutants D203A and E126K reveal conformational changes in EGF domains potentially responsible for enhanced protease lability and impaired extracellular matrix assembly.

Fibulin-4 is a 50 kDa glycoprotein of elastic fibers and plays an important role in development and function of elastic tissues. Fibulin-4 consists of a tandem array of five calcium-binding epidermal growth factor-like modules flanked by N- and C-terminal domains. Mutations in the human fibulin-4 gene EFEMP2 have been identified in patients affected with various arteriopathies including aneurysm, arterial tortuosity, or stenosis, but the molecular basis of most genotype-phenotype correlations is unknown. Here we present biochemical and computer modelling approaches designed to gain further insight into changes in structure and function of two fibulin-4 mutations (E126K and D203A), which are potentially involved in Ca2+ binding in the EGF2 and EGF4 domain, respectively. Using recombinantly produced fibulin-4 mutant and wild type proteins we show that both mutations introduced additional protease cleavage sites, impaired extracellular assembly into fibers, and affected binding to to fibrillin-1, latent TGF-ß-binding proteins, and the lysyl oxidase LOXL2. Molecular dynamics studies indicated that the E126K and D203A mutations do not necessarily result in a direct loss of the complexed Ca2+ ion after 500 ns simulation time, but in significantly enhanced fluctuations within the connecting loop between EGF3 and EGF4 domains and other conformational changes. In contrast, intentionally removing Ca2+ from EGF4 (D203A ΔCa) predicted dramatic changes in the protein structure. These results may explain the changes in protease cleavage sites, reduced secretion and impaired extracellular assembly of the E126K and D203A fibulin-4 mutants and provide further insight into understanding the molecular basis of the associated clinical phenotypes.

Distinct acidic clusters and hydrophobic residues in the alternative splice domains X1 and X2 of alpha7 integrins define specificity for laminin isoforms.

The binding specificity of alpha7beta1 integrins for different laminin isoforms is defined by the X1 and X2 splice domains located in the beta-propeller domain of the alpha7 subunit. In order to gain insight into the mechanism of specific laminin-integrin interactions, we defined laminin-binding epitopes of the alpha7X1 and -X2 domains by single amino acid substitutions and domain swapping between X1 and X2. The interaction of mutated, recombinantly prepared alpha7X1beta1 and alpha7X2beta1 heterodimers with various laminin isoforms was studied by surface plasmon resonance and solid phase binding assays. The data show that distinct clusters of surface-exposed acidic residues located in different positions of the X1 and the X2 loops are responsible for the specific recognition of laminins. These residues are conserved between the respective X1 or X2 splice domains of the alpha7 chains of different species, some also in the corresponding X1/X2 splice domains of alpha6 integrin. Interestingly, ligand binding was also modulated by mutating surface-exposed hydrophobic residues (alpha7X1L205, alpha7X2Y208) at positions corresponding to the fibronectin binding synergy site in alpha5beta1 integrin. Mutations in X1 that affected binding to laminin-1 also affected binding to laminin-8 and -10, but not to the same extent, thus allowing conclusions on the specific role of individual surface epitopes in the selective recognition of laminin-1 versus laminins -8 and -10. The role of the identified epitopes was confirmed by molecular dynamics simulations of wild-type integrins and several inactivating mutations. The analysis of laminin isoform interactions with various X1/X2 chimaera lend further support to the key role of negative surface charges and pointed to an essential contribution of the N-terminal TARVEL sequence of the X1 domain for recognition of laminin-8 and -10. In conclusion, specific surface epitopes containing charged and hydrophobic residues are essential for ligand binding and define specific interactions with laminin isoforms.

Monitoring biological membrane-potential changes: a CI QM/MM study.

In recent decades, new less-invasive, nonlinear optical methods have been proposed and optimized for monitoring fast physiological processes in biological cells. One of these methods allows the action potential (AP) in cardiomyocytes or neurons to be monitored by means of second-harmonic generation (SHG). We now present the first, to our knowledge, simulations of the dependency of the intensity of the second harmonic (I(SHG)) on variations of the transmembrane potential (TMP) in a cardiomyocyte during an action potential (AP). For this, an amphiphilic potential-sensitive styryl dye molecule with nonlinear optical properties was embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer, replacing one of the phospholipid molecules. External electrical fields with different strengths were applied across the membrane to simulate the AP of a heart-muscle cell. We used a combined classical/quantum mechanical approach to model the structure and the spectroscopic properties of the embedded chromophore. Two 10 ns molecular dynamics (MD) simulations provided input geometries for semiempirical molecular orbital (QM/MM) single-point configuration interaction (CI) calculations, which were used to calculate the wavelengths and oscillator strengths of electronic transitions in the di-8-ANEPPS dye molecule. The results were then used in a sum-over-states treatment to calculate the second-order hyperpolarizability. The square of the hyperpolarizability scales with the intensity of the second harmonic, which is used to monitor the action potentials of cardiomyocytes experimentally. Thus, we computed changes in the intensity of the second harmonic (DeltaI(SHG)) as function of TMP changes. Our results agree well with experimental measurements.

Molecular dynamics simulations of the tetracycline-repressor protein: the mechanism of induction.

Molecular dynamics simulations on the tetracycline-repressor (TetR) protein, both in the absence of an inducer and complexed with the inducers tetracycline and 5a,6-anhydrotetracycline, show significant differences in the structures and dynamics of the induced and non-induced forms of the protein. Calpha-density-difference plots, low-frequency normal vibrations and inter-residue interaction energies all point to a common mechanism of induction. The inducer displaces Asp156 from the magnesium ion in the binding pocket, leading to a short cascade of rearrangements of salt bridges that results in the allosteric change. The increased flexibility of the induced form of the protein is suggested to contribute to the decrease in binding affinity to DNA on induction.

Molecular dynamics characterization of the structures and induction mechanisms of a reverse phenotype of the tetracycline receptor.

Molecular-dynamics simulations have been used to investigate the mechanism of induction of a mutant (revTetR) of the tetracycline repressor protein (TetR) that shows the reverse phenotype (i.e., it is induced in the absence of tetracyclines and not in their presence). Low-frequency, normal-mode analyses demonstrate that the reverse phenotype is reproduced by the simulations on the basis of criteria established for wild-type TetR. The reverse phenotype is caused by the fact that the DNA-binding heads in revTetR are closer than the ideal distance needed for DNA-binding when no inducer is present. This distance increases on binding an inducer. Whereas this distance increase makes the interhead distance too large in wild-type TetR, it increases to the ideal value in revTetR. Thus, the mechanism of induction is the same for the two proteins, but the consequences are reversed because of the smaller interhead distance in revTetR when no inducer is present.

Structural changes and binding characteristics of the tetracycline-repressor binding site on induction.

The binding motif (pharmacophore) for induction and the changes in the structure of the binding site that accompany induction have been determined from molecular-dynamics simulations on the tetracycline-repressor signal-transduction protein. The changes and the induction mechanism are discussed and compared with conclusions drawn from earlier X-ray structures. The differences in inducer strength of tetracycline and 5a,6-anhydrotetracycline are discussed with respect to their interaction in the MD simulations.

SCRF-DFT and NMR comparison of tetracycline and 5a,6-anhydrotetracycline in solution.

A combination of structures, energies, and spectral data calculated using density functional theory (DFT) with experimental NMR data has been used to assign conformational equilibria for tetracycline and 5a,6-anhydrotetracycline in water at pH 1, 7, and 10 and in chloroform (5a,6-anhydrotetracycline) and methanol (tetracycline). The results suggest that tetracycline always prefers the extended conformation but that 5a,6-anhydrotetracycline exists in water as a mixture of the two conformers and in chloroform exclusively in the twisted conformation. The conformational equilibria are also shown to be pH dependent.

Progesterone 5ß-reductases/iridoid synthases (PRISE): gatekeeper role of highly conserved phenylalanines in substrate preference and trapping is supported by molecular dynamics simulations.

Vein Patterning 1 (VEP1)-encoded progesterone 5ß-reductases/iridoid synthases (PRISE) belong to the short-chain dehydrogenase/reductase superfamily of proteins. They are characterized by a set of highly conserved amino acids in the substrate-binding pocket. All PRISEs are capable of reducing the activated C=C double bond of various enones enantioselectively and therefore have a potential as biocatalysts in bioorganic synthesis. Here, recombinant forms of PRISEs of Arabidopsis thaliana and Digitalis lanata were modified using site-directed mutagenesis (SDM). In rDlP5ßR, a set of highly conserved amino acids in the vicinity of the catalytic center was individually substituted for alanine resulting in considerable to complete loss of enone reductase activity. F153 and F343, which can be found in most PRISEs known, are located at the outer rim of the catalytic cavity and seem to be involved in substrate binding and their role was addressed in a series of SDM experiments. The wild-type PRISE accepted progesterone (large hydrophobic 1,4-enone) as well as 2-cyclohexen-1-one (small hydrophilic 1,4-enone), whereas the double mutant rAtP5ßR_F153A_F343A converted progesterone much better than the wild-type enzyme but almost lost its capability of reducing 2-cyclohexen-1-one. Recombinant Draba aizoides P5ßR (rDaP5ßR) has a second pair of phenylalanines at position 156 and 345 at the rim of the binding site. These two phenylalanines were introduced into rAtP5ßR_F153A_F343A and the resulting quadruple mutant rAtP5ßR_F153A_F343A_V156F_V345F partly recovered the ability to reduce 2-cyclohexen-1-one. These results can best be explained by assuming a trapping mechanism in which phenylalanines at the rim of the substrate-binding pocket are involved. The dynamic behavior of individual P5ßRs and mutants thereof was investigated by molecular dynamics simulations and all calculations supported the 'gatekeeper' role of phenylalanines at the periphery of the substrate-binding pocket. Our findings provide structural and mechanistic explanations for the different substrate preferences seen among the natural PRISEs and help to explain the large differences in catalytic efficiency found for different types of 1,4-enones.

Functional consequence of fibulin-4 missense mutations associated with vascular and skeletal abnormalities and cutis laxa.

Fibulin-4 is a 60kDa calcium binding glycoprotein that has an important role in development and integrity of extracellular matrices. It interacts with elastin, fibrillin-1 and collagen IV as well as with lysyl oxidases and is involved in elastogenesis and cross-link formation. To date, several mutations in the fibulin-4 gene (FBLN4/EFEMP2) are known in patients whose major symptoms are vascular deformities, aneurysm, cutis laxa, joint laxity, or arachnodactyly. The pathogenetic mechanisms how these mutations translate into the clinical phenotype are, however, poorly understood. In order to elucidate these mechanisms, we expressed fibulin-4 mutants recombinantly in HEK293 cells, purified the proteins in native forms and analyzed alterations in protein synthesis, secretion, matrix assembly, and interaction with other proteins in relation to wild type fibulin-4. Our studies show that different mutations affect these properties in multiple ways, resulting in fibulin-4 deficiency and/or impaired ability to form elastic fibers. The substitutions E126K and C267Y impaired secretion of the protein, but not mRNA synthesis. Furthermore, the E126K mutant showed less resistance to proteases, reduced binding to collagen IV and fibrillin-1, as well as to LTBP1s and LTBP4s. The A397T mutation introduced an extra O-glycosylation site and deleted binding to LTBP1s. We show that fibulin-4 binds stronger than fibulin-3 and -5 to LTBP1s, 3, and 4s, and to the lysyl oxidases LOX and LOXL1; the binding of fibulin-4 to the LOX propeptide was strongly reduced by the mutation E57K. These findings show that different mutations in the fibulin-4 gene result in different molecular defects affecting secretion rates, protein stability, LOX-induced cross-linking, or binding to other ECM components and molecules of the TGF-ß pathway, and thus illustrate the complex role of fibulin-4 in connective tissue assembly.

Modeling the similarity and divergence of dopamine D2-like receptors and identification of validated ligand-receptor complexes.

Focusing on the similarity and divergence of GPCR subtypes and their ligand interactions, we generated dopamine D2, D3, and D4 receptor models based on the rhodopsin crystal structure and refined these with an extensive MM/MD protocol. After validation by diagnostic experimental data, subtype-specific relative positions of TM1, 2, 6, and 7 and bending angles of TM7 were found. To sample the conformational space of the complex, we performed simulated-annealing runs of the receptor protein with the sub-nanomolar antagonist spiperone. Docking a representative set of ligands, we were able to identify one superior model for each subtype when excellent correlations between predicted energies of binding and experimental affinities (r2 = 0.72 for D2, 0.91 for D3 and 0.77 for D4) could be observed. Further analysis revealed general ligand interactions with ASP3.32 and aromatic residues in TM6/7 and individual key interactions with TM1 and TM2 residues of the D3 and D4 receptor models, respectively.

Molecular dynamics simulation links conformation of a pore-flanking region to hyperekplexia-related dysfunction of the inhibitory glycine receptor.

Inhibitory glycine receptors mediate rapid synaptic inhibition in mammalian spinal cord and brainstem. The previously identified hyperekplexia mutation GLRA1(P250T), located within the intracellular TM1-2 loop of the GlyR alpha1 subunit, results in altered receptor activation and desensitization. Here, elementary steps of ion channel function of alpha1(250) mutants were resolved and shown to correlate with hydropathy and molar volume of residue alpha1(250). Single-channel recordings and rapid activation kinetic studies using laser pulse photolysis showed reduced conductance but similar open probability of alpha1(P250T) mutant channels. Molecular dynamics simulation of a helix-turn-helix motif representing the intracellular TM1-2 domain revealed alterations in backbone conformation, indicating an increased flexibility in these mutants that paralleled changes in elementary steps of channel function. Thus, the architecture of the TM1-2 loop is a critical determinant of ion channel conductance and receptor desensitization.

Three-dimensional modeling of a pre-B-cell receptor.

Signals delivered by the immunoglobulin (Ig)-like pre-B-cell receptor (pre-BCR) are critical for efficient maturation of early precursor B (pre-B) cells. A pre-BCR contains two immunoglobulin mu-heavy chains (muHC), two surrogate light chains (SLC) consisting of the non-covalently associated polypeptides, VpreB and lambda5, and the heterodimeric signaling transducer Igalpha/beta. Although, it is generally accepted that signals initiated from the pre-BCR are required for efficient expansion and differentiation of pre-B cells, the three-dimensional structure of this receptor has not yet been determined by either NMR or X-ray spectroscopy. Therefore, we used indirect computer-assisted molecular modeling techniques to predict for the first time three-dimensional coordinates of the pre-BCR, the conformation of the SLC components, VpreB and lambda5, and the position and flexibility of the so-called non-Ig-like unique tails at the C-terminus of VpreB and the N-terminus of lambda5. Structure prediction revealed that these unique tails of VpreB and lambda5 protrude from the SLC at the position where the CDR3 of a conventional IgLchain would be located. Thus, the unique tails are accessible for ligand binding, which supports the recent finding that the lambda5 unique tail is required for pre-BCR/stroma cell interaction. Further, the non-covalent interaction of the extra beta-strand of lambda5 (beta8) with VpreB is predicted to result in a stabilization of the tertiary structure of VpreB. In summary, three-dimensional computer modeling suggests that the structure of a pre-BCR resembles that of a conventional B-cell receptor (BCR) and that the lambda5 unique tail could be a major binding site for pre-BCR ligands.

In Silico Adoption of an Orphan Nuclear Receptor NR4A1.

A 4.1 µs molecular dynamics simulation of the NR4A1 (hNur77) apo-protein has been undertaken and a previously undetected druggable pocket has become apparent that is located remotely from the 'traditional' nuclear receptor ligand-binding site. A NR4A1/bis-indole ligand complex at this novel site has been found to be stable over 1 µs of simulation and to result in an interesting conformational transmission to a remote loop that has the capacity to communicate with a NBRE within a RXR-α/NR4A1 heterodimer. Several features of the simulations undertaken indicate how NR4A1 can be affected by alternate-site modulators.