Biomedical subjects

Martin Zacharias

Publications and source records attributed to Martin Zacharias.

32 records · Page 2Linked to original sources

PubMed2005

ATTRACT: protein-protein docking in CAPRI using a reduced protein model.

Protein-protein complex structures have been predicted for CAPRI Rounds 3 and 5 using a reduced protein model. Proteins are represented by up to 3 pseudoatoms per amino acid. The docking approach termed ATTRACT is based on energy minimization in translational and rotational degrees of freedom of one protein with respect to another protein. The reduced protein model allows one to perform systematic docking minimization of many thousand start structures in reasonable computer time. Flexibility of critical surface side-chains can be accounted for by a multiple conformational copy approach. The multicopy approach allows simultaneous adjustment of side-chain conformations and optimization of translational and rotational degrees of freedom of one protein with respect to the partner during docking. For 3 (Targets 8, 14, and 19) out of 5 CAPRI targets, the approach resulted in predictions in close agreement with experiment [root-mean-square deviation (RMSD) of backbone atoms within 10 A of the protein-protein interface < 1.8 A]. The comparison of predicted and experimental structures of the CAPRI targets indicates that besides local conformational changes (e.g., changes in side-chain conformations), global conformational changes of the protein backbone can be critical for complex formation. These conformational changes not accounted for during docking are a likely reason for the unrealistic predictions in 2 cases (Targets 9 and 18).

Algorithms↗

PubMed2005

A novel cGUUAg tetraloop structure with a conserved yYNMGg-type backbone conformation from cloverleaf 1 of bovine enterovirus 1 RNA.

The 5'-terminal cloverleaf (CL)-like RNA structures are essential for the initiation of positive- and negative-strand RNA synthesis of entero- and rhinoviruses. SLD is the cognate RNA ligand of the viral proteinase 3C (3C(pro)), which is an indispensable component of the viral replication initiation complex. The structure of an 18mer RNA representing the apical stem and the cGUUAg D-loop of SLD from the first 5'-CL of BEV1 was determined in solution to a root-mean-square deviation (r.m.s.d.) (all heavy atoms) of 0.59 A (PDB 1Z30). The first (antiG) and last (synA) nucleotide of the D-loop forms a novel 'pseudo base pair' without direct hydrogen bonds. The backbone conformation and the base-stacking pattern of the cGUUAg-loop, however, are highly similar to that of the coxsackieviral uCACGg D-loop (PDB 1RFR) and of the stable cUUCGg tetraloop (PDB 1F7Y) but surprisingly dissimilar to the structure of a cGUAAg stable tetraloop (PDB 1MSY), even though the cGUUAg BEV D-loop and the cGUAAg tetraloop differ by 1 nt only. Together with the presented binding data, these findings provide independent experimental evidence for our model [O. Ohlenschlager, J. Wohnert, E. Bucci, S. Seitz, S. Hafner, R. Ramachandran, R. Zell and M. Gorlach (2004) Structure, 12, 237-248] that the proteinase 3C(pro) recognizes structure rather than sequence.

3C Viral Proteases↗

PubMed2005

Design and biological evaluation of linear and cyclic phosphopeptide ligands of the N-terminal SH2 domain of protein tyrosine phosphatase SHP-1.

In an effort to gain further insight into the conformational and topographical requirements for recognition by the N-terminal SH2 domain of protein tyrosine phosphatase SHP-1, we synthesized a series of linear and cyclic peptides derived from the sequence surrounding phosphotyrosine 2267 in the receptor tyrosine kinase Ros (EGLNpYMVL). A molecular modeling approach was used to suggest peptide modifications sterically compatible with the N-SH2-peptide binding groove and possibly enhanced binding affinities compared to the parent peptide. The potencies of the synthesized compounds were evaluated by assaying their ability to stimulate phosphatase activity as well as by their binding affinities to the GST-fused N-SH2 domain of SHP-1. In the series of linear peptides, structural modifications of Ros pY2267 in positions pY + 1 to pY + 3 by amino acid residues structurally related to Phe, for example l-erythro/threo-Abu(betaPh) (5a, 5b), yielded ligands with increased binding affinity. The incorporation of d-amino acid residues at pY + 1 and pY + 3 led to inactive peptides. The replacement of Phe in both pY + 1 and pY + 3 by Tic (1,2,3,4-tetrahydroisoquinoline-3-carboxylic acid) was also not tolerated due to steric hindrance. Cyclic peptides (13, 14) that were linked via residues in positions pY - 1 (Lys) and pY + 2 (Asp/Glu) and contained a Gly residue in the bridging unit displayed much lower potencies for the stimulation of SHP-1 activity but increased binding affinities compared to Ros pY2267. They partially competed with Ros pY2267 in the activation assay. Such cyclic structures may serve as scaffolds for competitive SHP-1 inhibitor design targeting N-SH2 domain-protein interactions that block SHP-1 activation.

Drug Design↗

PubMed2005

Hereditary lattice corneal dystrophy is associated with corneal amyloid deposits enclosing C-terminal fragments of keratoepithelin.

PURPOSE: To investigate the molecular basis of hereditary lattice corneal dystrophy (LCD) type IIIA associated with corneal amyloid deposits afflicting several members of a four-generation family. METHODS: Histologic, immunohistochemical and biochemical studies were performed on corneal tissue samples obtained after perforating keratoplasty. DNA was extracted from peripheral blood leukocytes. All exons of the keratoepithelin-encoding TGFBI gene were amplified and sequenced. The presence of a mutation was confirmed by digestion of the isolated PCR product with the restriction enzyme AlwNI. RESULTS: The cornea of the index patient (II-1) contained large patchy deposits of amyloid, which were immunoreactive for the C terminus of keratoepithelin. Western blot analysis of the polypeptide chains extracted from the amyloid deposits of paraffin-embedded tissue revealed that these represented mainly fragments of the full-length protein. The smallest fragments were 6.5 and 6.9 kDa. DNA analyses of the TGFBI gene revealed a heterozygous T-->C transition at the second position of codon 540 in exon 12, indicating that replacement of phenylalanine by serine (Phe540Ser) leads to dominant disease. The mutation creates a new restriction site for the enzyme AlwNI. Five of the examined family members carried this mutation. Three of them (aged >/=41 years) had the disease, two family members (aged <20 years) do not yet show any clinical symptoms. An additional inconsequential single-nucleotide polymorphism (T1667C) was found at the third position of the same codon (Phe540Phe) in three unaffected family members. CONCLUSIONS: This is the first report of a single-nucleotide mutation at codon 540 of TGFBI leading to LCD, and the first to demonstrate that the amyloid deposits in LCD contain proteolytic fragments of keratoepithelin.

Adult↗

PubMed2004

Influence of a fluorobenzene nucleobase analogue on the conformational flexibility of RNA studied by molecular dynamics simulations.

Chemically modified bases are frequently used to stabilize nucleic acids, to study the driving forces for nucleic acid structure formation and to tune DNA and RNA hybridization conditions. In particular, fluorobenzene and fluorobenzimidazole base analogues can act as universal bases able to pair with any natural base and to stabilize RNA duplex formation. Although these base analogues are compatible with an A-form RNA geometry, little is known about the influence on the fine structure and conformational dynamics of RNA. In the present study, nano-second molecular dynamics (MD) simulations have been performed to characterize the dynamics of RNA duplexes containing a central 1'-deoxy-1'-(2,4-difluorophenyl)-beta-D-ribofuranose base pair or opposite to an adenine base. For comparison, RNA with a central uridine:adenine pair and a 1'-deoxy-1'-(phenyl)-beta-D-ribofuranose opposite to an adenine was also investigated. The MD simulations indicate a stable overall A-form geometry for the RNAs with base analogues. However, the presence of the base analogues caused a locally enhanced mobility of the central bases inducing mainly base pair shear and opening motions. No stable 'base-paired' geometry was found for the base analogue pair or the base analogue:adenine pairs, which explains in part the universal base character of these analogues. Instead, the conformational fluctuations of the base analogues lead to an enhanced accessibility of the bases in the major and minor grooves of the helix compared with a regular base pair.

Adenine↗

PubMed2004

The human T-cell lymphotropic virus type-I dimerization initiation site forms a hairpin loop, unlike previously characterized retroviral dimerization motifs.

The formation of genomic RNA dimers during the retroviral life cycle is essential for optimal viral replication and infectivity. The sequences and RNA structures responsible for this interaction are located in the untranslated 5' leader RNA, along with other cis-acting signals. Dimer formation occurs by specific interaction between identical structural motifs. It is believed that an initial kissing hairpin forms following self-recognition by autocomplementary RNA loops, leading to formation of an extended stable duplex. The dimerization initiation site (DIS) of the deltaretrovirus human T-cell lymphotropic virus type-I (HTLV-I) has been previously localized to a 14-nucleotide sequence predicted to contain an RNA stem loop. Biochemical probing of the monomeric RNA structure using RNAse T1, RNAse V1, RNAse U2, lead acetate, and dimethyl sulfate has led to the generation of the first structural map of the HTLV-I DIS. A comprehensive data set of individual nucleotide modifications reveals that the structural motif responsible for HTLV-I RNA dimerization forms a trinucleotide RNA loop, unlike any previously characterized retroviral dimerization motif. Molecular modeling demonstrates that this can be formed by an unusual C:synG base pair closing the loop. Comparative phylogeny indicates that such a motif may also exist in other deltaretroviruses.

5' Untranslated Regions↗

PubMed2004

Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP.

Most current docking methods to identify possible ligands and putative binding sites on a receptor molecule assume a rigid receptor structure to allow virtual screening of large ligand databases. However, binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a bound ligand. An approach is presented that allows relaxation of the protein conformation in precalculated soft flexible degrees of freedom during ligand-receptor docking. For the immunosuppressant FK506-binding protein FKBP, the soft flexible modes are extracted as principal components of motion from a molecular dynamics simulation. A simple penalty function for deformations in the soft flexible mode is used to limit receptor protein deformations during docking that avoids a costly recalculation of the receptor energy by summing over all receptor atom pairs at each step. Rigid docking of the FK506 ligand binding to an unbound FKBP conformation failed to identify a geometry close to experiment as favorable binding site. In contrast, inclusion of the flexible soft modes during systematic docking runs selected a binding geometry close to experiment as lowest energy conformation. This has been achieved at a modest increase of computational cost compared to rigid docking. The approach could provide a computationally efficient way to approximately account for receptor flexibility during docking of large numbers of putative ligands and putative docking geometries.

Binding Sites↗

PubMed2004

Molecular determinants for high-affinity block of human EAG potassium channels by antiarrhythmic agents.

Undesired block of human ERG1 potassium channels is the basis for cardiac side effects of many different types of drugs. Therefore, it is important to know exactly why some drugs particularly bind to these channels with high affinity. Upon expression in mammalian cells and Xenopus laevis oocytes, we investigated the inhibition of the closely related hEAG1 and hEAG2 channels by agents that have previously been reported to block hERG1 channels. Clofilium inhibited hEAG1 and hERG1 with the same potency, whereas hEAG2 was about 150-fold less sensitive to this antiarrhythmic agent. The molecular determinants for this difference are residues Ser436 and Val437 in the inner cavity of the pore and Ala453, which is located in S6 (i.e., remote from the inner cavity). A modeling approach that allowed for partial conformational relaxation of hEAG model structures upon ligand docking suggests that high-affinity block of ether à go-go channels is mediated by an anchoring of the clofilium alkane tail between S6 and the pore helices. In qualitative agreement with experiments, the mutations of hEAG1 residues Ser436 and Val437 to the corresponding larger hEAG2 residues (Thr432, Ile433) resulted in reduced sterical fit between the ligand and the binding cavity. The model is further supported by functional assays involving (+)-N-[1'-(6-cyano-1,2,3,4-tetrahydro-2(R)-naphthalenyl)-3,4-dihydro-4(R)-hydroxyspiro(2H-1-benzopyran-2,4'-piperidin)-6-yl]methanesulfonamide monohydrochloride (MK-499), terfenadine, quinidine, and tetrabutylammonium that are differentially affected by mutations in the pore pocket.

Amino Acid Sequence↗

PubMed2004

Tapasin and other chaperones: models of the MHC class I loading complex.

MHC (major histocompatibility complex) class I molecules bind intracellular virus-derived peptides in the endoplasmic reticulum (ER) and present them at the cell surface to cytotoxic T lymphocytes. Peptide-free class I molecules at the cell surface, however, could lead to aberrant T cell killing. Therefore, cells ensure that class I molecules bind high-affinity ligand peptides in the ER, and restrict the export of empty class I molecules to the Golgi apparatus. For both of these safeguard mechanisms, the MHC class I loading complex (which consists of the peptide transporter TAP, the chaperones tapasin and calreticulin, and the protein disulfide isomerase ERp57) plays a central role. This article reviews the actions of accessory proteins in the biogenesis of class I molecules, specifically the functions of the loading complex in high-affinity peptide binding and localization of class I molecules, and the known connections between these two regulatory mechanisms. It introduces new models for the mode of action of tapasin, the role of the class I loading complex in peptide editing, and the intracellular localization of class I molecules.

Animals↗

PubMed2004

Conformational flexibility of the MHC class I alpha1-alpha2 domain in peptide bound and free states: a molecular dynamics simulation study.

Major histocompatibility complex class I proteins play a key role in the recognition and presentation of peptide antigens to the host immune system. The structure of various major histocompatibility complex class I proteins has been determined experimentally in complex with several antigenic peptides. However, the structure in the unbound (empty) form is not known. To study the conformational dynamics of the empty major histocompatibility complex class I molecule comparative molecular dynamics simulations have been performed starting from the crystal structure of a peptide bound class I peptide-binding domain in the presence and absence of a peptide ligand. Simulations including the bound peptide stayed close to the experimental start structure at both simulation temperatures (300 and 355 K) during the entire simulation of 26 ns. Several independent simulations in the absence of peptide indicate that the empty domain may not adopt a single defined conformation but is conformationally significantly more heterogeneous in particular within the alpha-helices that flank the peptide binding cleft. The calculated conformational dynamics along the protein chain correlate well with available spectroscopic data and with the observed site-specific sensitivity of the empty class I protein to proteolytic digestion. During the simulations at 300 K the binding region for the peptide N-terminus stayed close to the conformation in the bound state, whereas the anchor region for the C-terminus showed significantly larger conformational fluctuations. This included a segment at the beginning of the second alpha-helix in the domain that is likely to be involved in the interaction with the chaperone protein tapasin during the peptide-loading process. The simulation studies further indicate that peptide binding at the C- and N-terminus may follow different mechanisms that involve different degrees of induced conformational changes in the peptide-binding domain. In particular binding of the peptide C-terminus may require conformational stabilization by chaperone proteins during peptide loading.

Binding Sites↗

PubMed2003

Effective dephosphorylation of Src substrates by SHP-1.

The protein-tyrosine phosphatase SHP-1 is a negative regulator of multiple signal transduction pathways. We observed that SHP-1 effectively antagonized Src-dependent phosphorylations in HEK293 cells. This occurred by dephosphorylation of Src substrates, because Src activity was unaffected in the presence of SHP-1. One reason for efficient dephosphorylation was activation of SHP-1 by Src. Recombinant SHP-1 had elevated activity subsequent to phosphorylation by Src in vitro, and SHP-1 variants with mutated phosphorylation sites in the C terminus, SHP-1 Y538F, and SHP-1 Y538F,Y566F were less active toward Src-generated phosphoproteins in intact cells. A second reason for efficient dephosphorylation is the substrate selectivity of SHP-1. Pull-down experiments with different GST-SHP-1 fusion proteins revealed efficient interaction of Src-generated phosphoproteins with the SHP-1 catalytic domain rather than with the SH2 domains. Phosphopeptides that correspond to good Src substrates were efficiently dephosphorylated by SHP-1 in vitro. Phosphorylated "optimal Src substrate" AEEEIpYGEFEA (where pY is phosphotyrosine) and a phosphopeptide corresponding to a recently identified Src phosphorylation site in p120 catenin, DDLDpY(296)GMMSD, were excellent SHP-1 substrates. Docking of these phosphopeptides into the catalytic domain of SHP-1 by molecular modeling was consistent with the biochemical data and explains the efficient interaction. Acidic residues N-terminal of the phosphotyrosine seem to be of major importance for efficient substrate interaction. Residues C-terminal of the phosphotyrosine probably contribute to the substrate selectivity of SHP-1. We propose that activation of SHP-1 by Src and complementary substrate specificities of SHP-1 and Src may lead to very transient Src signals in the presence of SHP-1.

Amino Acid Sequence↗

PubMed2003

Sequence context dependence of tandem guanine:adenine mismatch conformations in RNA: a continuum solvent analysis.

Guanine:adenine (G:A) mismatches and in particular tandem G:A (tG:A) mismatches are frequently observed in biological RNA molecules and can serve as sites for tertiary interaction, metal binding and protein recognition. Depending on the surrounding sequence tG:A mismatches can adopt different basepairing topologies. In the sequence context (5'-) GGAC (tandem G:A in bold) a face-to-face (imino or Watson-Crick-like) pairing is preferred whereas in the CGAG context, G and A adopt a sheared arrangement. Systematic conformational searches with a generalized Born continuum model and molecular dynamics simulations including explicit water molecules and ions have been used to generate face-to-face and sheared tG:A mismatches in both CGAG and GGAC sequence contexts. Conformations from both approaches were evaluated using the same force field and a Poisson-Boltzmann continuum solvent model. Although the substate analysis predicted the sheared arrangement to be energetically preferred in both sequence contexts, a significantly greater preference of the sheared form was found for the CGAG context. In agreement with the experimental observation, the analysis of molecular dynamics trajectories indicated a preference of the sheared form in the case of the CGAG-context and a favorization of the face-to-face form in the case of the GGAC context. The computational studies allowed to identify energetic contributions that stabilize or destabilize the face-to-face and sheared tandem mismatch topologies. The calculated nonpolar solvation and Lennard-Jones packing interaction were found to stabilize the sheared topology independent of the sequence context. Electrostatic contributions are predicted to make the most significant contribution to the sequence context dependence on the structural preference of tG:A mismatches.

Adenine↗

PubMed2003

Protein-protein docking with a reduced protein model accounting for side-chain flexibility.

A protein-protein docking approach has been developed based on a reduced protein representation with up to three pseudo atoms per amino acid residue. Docking is performed by energy minimization in rotational and translational degrees of freedom. The reduced protein representation allows an efficient search for docking minima on the protein surfaces within. During docking, an effective energy function between pseudo atoms has been used based on amino acid size and physico-chemical character. Energy minimization of protein test complexes in the reduced representation results in geometries close to experiment with backbone root mean square deviations (RMSDs) of approximately 1 to 3 A for the mobile protein partner from the experimental geometry. For most test cases, the energy-minimized experimental structure scores among the top five energy minima in systematic docking studies when using both partners in their bound conformations. To account for side-chain conformational changes in case of using unbound protein conformations, a multicopy approach has been used to select the most favorable side-chain conformation during the docking process. The multicopy approach significantly improves the docking performance, using unbound (apo) binding partners without a significant increase in computer time. For most docking test systems using unbound partners, and without accounting for any information about the known binding geometry, a solution within approximately 2 to 3.5 A RMSD of the full mobile partner from the experimental geometry was found among the 40 top-scoring complexes. The approach could be extended to include protein loop flexibility, and might also be useful for docking of modeled protein structures.

Algorithms↗

PubMed2002

Molecular dynamics simulations on the free and complexed N-terminal SH2 domain of SHP-2.

Signal transduction events are often mediated by small protein domains such as SH2 (Src homology 2) domains that recognize phosphotyrosines (pY) and flanking sequences. In case of the SHP-2 receptor tyrosine phosphatase an N-terminal SH2 domain binds and inactivates the phosphatase (PTP) domain. The pY-peptide-binding site on the N-terminal SH2 domain does not overlap with the PTP binding region. Nevertheless, pY-peptide binding causes domain dissociation and phosphatase activation. Comparative multi-nanosecond molecular dynamics simulations on the N-SH2 domain in ligand-bound and free states have been performed to study the allosteric mechanism that leads to domain dissociation upon pY-peptide binding. Significant ligand-dependent differences in the conformational flexibility of regions that are involved in SH2-PTP domain association have been observed. The results support a mechanism of signal transduction where SH2-peptide binding modulates the domain flexibility and reduces its capacity to fit into the entrance of the PTP catalytic domain of SHP-2.

Allosteric Regulation↗