Biomedical subjects

Martin Zacharias

Publications and source records attributed to Martin Zacharias.

At least 19 recordsLinked to original sources

PubMed2026

Beyond benchmarking: an expert-guided consensus approach to spatially aware clustering.

Spatial omics technologies have revolutionized the study of tissue architecture and cellular heterogeneity by integrating molecular profiles with spatial localization. In spatially resolved transcriptomics, delineating higher-order anatomical structures is critical for understanding how cellular organization affects function. However, the reliability of current benchmarks of spatially aware clustering (SAC) methods is undermined by their narrow focus on Visium and brain tissue datasets and the incorrect interpretation of manual annotation as ground truth. Here we present SACCELERATOR, a community-driven, extensible framework that standardizes data formatting, method integration and metric evaluation, enabling rapid inclusion of new methods and datasets. Our analysis revealed substantial limitations in the generalizability and reproducibility of SAC methods and shows that anatomical labels commonly used as ground truths are often biased, error prone and unsuitable for benchmarking. Rather than ranking methods, we propose a consensus-guided workflow where descriptive spatial metrics highlight high-entropy regions of method disagreement, enabling targeted feedback for tissue experts. Applied to brain and cancer datasets, this approach uncovered biologically meaningful patterns overlooked by individual SAC methods and manual annotations, highlighting the need for iterative, expert-in-the-loop evaluation.

Benchmarking↗

PubMed2007

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential.

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.

Alanine↗

PubMed2006

Toward understanding the inactivation mechanism of monooxygenase P450 BM-3 by organic cosolvents: a molecular dynamics simulation study.

Cytochrome P450 BM-3 from Bacillus megaterium is an extensively studied enzyme for industrial applications. A major focus of current protein engineering research is directed to improving the catalytic performance of P450 BM-3 toward nonnatural substrates of industrial importance in the presence of organic solvents or cosolvents. For the latter reason, it is important to study the effect of organic cosolvent molecules on the structure and dynamics of the enzyme, in particular, the effect of cosolvent molecules on the active site's structure and dynamics. In this paper, we have studied, using molecular dynamics (MD) simulations, the F87A mutant of P450 BM-3 in the presence of DMSO as cosolvent, to understand the role of the F87A substitution for its catalytic activity. This mutant exhibits an altered regioselectivity and substrate specificity compared with wild-type; however, it has lower tolerance toward DMSO. The simulation results offer an explanation for the DMSO sensitivity of the F87A mutant. Our simulation results show that the F87 side chain prevents the disturbance of the water molecule bound to the heme iron by DMSO molecules. The absence of the phenyl ring in F87A mutant promotes interactions of the DMSO molecule with the heme iron resulting in water displacement by DMSO at the catalytic heme center.

Bacillus megaterium↗

PubMed2006

B-DNA under stress: over- and untwisting of DNA during molecular dynamics simulations.

The twist flexibility of DNA is central to its many biological functions. Explicit solvent molecular dynamics simulations in combination with an umbrella sampling restraining potential have been employed to study induced twist deformations in DNA. Simulations allowed us to extract free energy profiles for twist deformations and were performed on six DNA dodecamer duplexes to cover all 10 possible DNA basepair steps. The shape of the free energy curves was similar for all duplexes. The calculated twist deformability was in good agreement with experiment and showed only modest variation for the complete duplexes. However, the response of the various basepair steps on twist stress was highly nonuniform. In particular, pyrimidine/purine steps were much more flexible than purine/purine steps followed by purine/pyrimidine steps. It was also possible to extract correlations of twist changes and other helical as well as global parameters of the DNA molecules. Twist deformations were found to significantly alter the local as well as global shape of the DNA modulating the accessibility for proteins and other ligands. Severe untwisting of DNA below an average of 25 degrees per basepair step resulted in the onset of a global structural transition with a significantly smaller twist at one end of the DNA compared to the other.

Base Pairing↗

PubMed2006

Distance distributions of short polypeptides recovered by fluorescence resonance energy transfer in the 10 A domain.

Fluorescence resonance energy transfer (FRET) between tryptophan (Trp) as donor and 2,3-diazabicyclo[2.2.2]oct-2-ene (Dbo) as acceptor was studied by steady-state and time-resolved fluorescence spectroscopy. The unique feature of this FRET pair is its exceptionally short Förster radius (10 A), which allows one to recover distance distributions in very short structureless peptides. The technique was applied to Trp-(GlySer)n-Dbo-NH2 peptides with n = 0-10, for which the average probe/quencher distance ranged between 8.7 and 13.7 A experimentally (in propylene glycol, analysis according to wormlike chain model) and 8.6-10.2 A theoretically (for n = 0-6, GROMOS96 molecular dynamics simulations). The larger FRET efficiency in steady-state compared to time-resolved fluorescence experiments was attributed to a static quenching component, suggesting that a small but significant part (ca. 10%) of the conformations are already in van der Waals contact when excitation occurs.

Bridged Bicyclo Compounds, Heterocyclic↗

PubMed2006

Sequence specificity of SHP-1 and SHP-2 Src homology 2 domains. Critical roles of residues beyond the pY+3 position.

A combinatorial phosphotyrosyl (pY) peptide library was screened to determine the amino acid preferences at the pY+4 to pY+6 positions for the four SH2 domains of protein-tyrosine phosphatases SHP-1 and SHP-2. Individual binding sequences selected from the library were resynthesized and their binding affinities and specificities to various SH2 domains were further evaluated by SPR studies, stimulation of SHP-1 and SHP-2 phosphatase activity, and in vitro pulldown assays. These studies reveal that binding of a pY peptide to the N-SH2 domain of SHP-2 is greatly enhanced by a large hydrophobic residue (Trp, Tyr, Met, or Phe) at the pY+4 and/or pY+5 positions, whereas binding to SHP-1 N-SH2 domain is enhanced by either hydrophobic or positively charged residues (Arg, Lys, or His) at these positions. Similar residues at the pY+4 to pY+6 positions are also preferred by SHP-1 and SHP-2 C-SH2 domains, although their influence on the overall binding affinities is much smaller compared with the N-SH2 domains. A structural model was generated to qualitatively interpret the contribution of the pY+4 and pY+5 residues to the overall binding affinity. Examination of pY motifs from known SHP-1 and SHP-2-binding proteins shows that many of the pY motifs contain a hydrophobic or positively charged residue(s) at the pY+4 and pY+5 positions.

Amino Acid Sequence↗

PubMed2006

Minor groove deformability of DNA: a molecular dynamics free energy simulation study.

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.

Algorithms↗

PubMed2006

Accounting for loop flexibility during protein-protein docking.

Although reliable docking can now be achieved for systems that do not undergo important induced conformational change upon association, the presence of flexible surface loops, which must adapt to the steric and electrostatic properties of a partner, generally presents a major obstacle. We report here the first docking method that allows large loop movements during a systematic exploration of the possible arrangements of the two partners in terms of position and rotation. Our strategy consists in taking into account an ensemble of possible loop conformations by a multi-copy representation within a reduced protein model. The docking process starts from regularly distributed positions and orientations of the ligand around the whole receptor. Each starting configuration is submitted to energy minimization during which the best-fitting loop conformation is selected based on the mean-field theory. Trials were carried out on proteins with significant differences in the main-chain conformation of the binding loop between isolated form and complexed form, which were docked to their partner considered in their bound form. The method is able to predict complexes very close to the crystal complex both in terms of relative position of the two partners and of the geometry of the flexible loop. We also show that introducing loop flexibility on the isolated protein form during systematic docking largely improves the predictions of relative position of the partners in comparison with rigid-body docking.

Binding Sites↗

PubMed2006

Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study.

Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.

Adenine↗

PubMed2006

Enzymatic degradation of beta- and mixed alpha,beta-oligopeptides.

One of the main and most astonishing characteristics of peptides comprised of beta-amino acids with proteinogenic side chains is their extraordinarily high stability towards enzymatic degradation. So far, only certain microbial enzymes have been shown to cleave N-terminal beta(3)-homoamino acid residues from peptides. In this work, the L-aminopeptidase-D-amidase/esterase (DmpA) from Ochrobactrum anthropi LMG7991 is compared to two closely related beta-peptidyl aminopeptidases (BapA), which originate from Sphingosinicella strains, and to microsomal leucine aminopeptidase (LAP) as a reference. All four enzymes are aminopeptidases cleaving N-terminal amino acids from small peptides. Degradation experiments reveal that DmpA and both BapA enzymes exhibit unique, but clearly distinct substrate specificities and preferences. DmpA also cleaves beta- and mixed alpha,beta-peptides and amides, but a short side chain of the N-terminal beta-amino acid residue seems to be a prerequisite, since only peptides carrying N-terminal betahGly and beta(3)hAla are hydrolyzed with good efficiencies. Both beta-peptidyl aminopeptidases cleave beta-amino acids from a variety of beta-peptides and mixed alpha,beta-peptides, but they do not accept alpha-amino acids in the N-terminal position. Astonishingly, DmpA exhibited much higher catalytical rates for the mixed dipeptide carnosine (H-betahGly-His-OH) than for any other substrate described until now.

Aminopeptidases↗

PubMed2006

Phosphopeptide ligands of the SHP-1 N-SH2 domain: effects on binding and stimulation of phosphatase activity.

Src homology 2 (SH2)-domain-mediated interactions with phosphotyrosine (pY)-containing ligands are critical for the regulation of SHP-1 phosphatase activity. Peptides based on a binding site from receptor tyrosine kinase Ros (EGLN-pY2267-MVL, 1) have recently been shown to bind to the SHP-1 N-terminal SH2 domain (N-SH2) with considerably high affinity. In addition, two peptides cyclized between positions -1 and +2 relative to pY (EGLc[K(COCH(2)NH)pYMX]L-NH(2), 2: X=D, 3: X=E) bound to the N-SH2 domain, but did not activate the enzyme and even partially prevented stimulation of SHP-1 activity by the physiological ligand 1. These findings prompted us to further examine the determinants for optimal binding to the N-SH2 domain and for the stimulation and inhibition of SHP-1 activity. Herein we demonstrate that combining the preferred residues in both pY+1 (such as Phe or norleucine, Nle) and pY+3 (such as homophenylalanine, Hfe) leads to highly efficient activating ligands of SHP-1. Particularly in the context of the cyclic peptides 7 (EGLc[K(COCH(2)NH)pYFD]Hfe-NH(2)) and 8 (EGLc[K(COCH(2)NH)pYNleD]HfeL-NH(2)), the incorporation of these residues resulted in high-affinity ligands with a significantly increased ability to stimulate SHP-1 activity. We suggest that different binding modes (according to consensus sequences class I and II) are responsible for obtaining either activating (7 and 8) or nonactivating (2 and 3) ligands. Peptides such as 7 and 8 that bind in the extended fashion of the type II mode activate the phosphatase through complete filling of the cavity for pY+3. In contrast, peptides such as 2 and 3 that bind in the class I mode do not activate the enzyme because they allow more conformational space at pY+3. Therefore, their binding does not force the conformational transition necessary to trigger the dissociation of N-SH2 and the catalytic domain.

Amino Acid Sequence↗

PubMed2006

RNA kink-turns as molecular elbows: hydration, cation binding, and large-scale dynamics.

The presence of Kink-turns (Kt) at key functional sites in the ribosome (e.g., A-site finger and L7/L12 stalk) suggests that some Kink-turns can confer flexibility on RNA protuberances that regulate the traversal of tRNAs during translocation. Explicit solvent molecular dynamics demonstrates that Kink-turns can act as flexible molecular elbows. Kink-turns are associated with a unique network of long-residency static and dynamical hydration sites that is intimately involved in modulating their conformational dynamics. An implicit solvent conformational search confirms the flexibility of Kink-turns around their X-ray geometries and identifies a second low-energy region with open structures that could correspond to Kink-turn geometries seen in solution experiments. An extended simulation of Kt-42 with the factor binding site (helices 43 and 44) shows that the local Kt-42 elbow-like motion fully propagates beyond the Kink-turn, and that there is no other comparably flexible site in this rRNA region. Kink-turns could mediate large-scale adjustments of distant RNA segments.

Binding Sites↗

PubMed2006

Serum-stable RNA aptamers to an invariant surface domain of live African trypanosomes.

African trypanosomes are extracellular blood parasites that cause sleeping sickness in humans and Nagana in cattle. The therapeutics used to control and treat these diseases are very ineffective and thus, the development of new drugs is urgently needed. We have previously suggested to use trypanosome-specific RNA aptamers as tools for the development of novel trypanocidal compounds. Here, we report the selection of a 2'-NH(2)-modified RNA aptamer that binds to live trypanosomes with an affinity of 70 +/- 15 nM. The aptamer adopts a stable G-quartet structure and has a half-life in human serum of > 30 h. RNA binding is restricted to the flagellar attachment zone, located between the cell body and the flagellum of the parasite. We demonstrate that antigen-tagged preparations of the aptamer can bind to live trypanosomes and that they can be used to re-direct immunoglobulins to the parasite surface.

Animals↗

PubMed2005

A statistical analysis of random mutagenesis methods used for directed protein evolution.

We have developed a statistical method named MAP (mutagenesis assistant program) to equip protein engineers with a tool to develop promising directed evolution strategies by comparing 19 mutagenesis methods. Instead of conventional transition/transversion bias indicators as benchmarks for comparison, we propose to use three indicators based on the subset of amino acid substitutions generated on the protein level: (1) protein structure indicator; (2) amino acid diversity indicator with a codon diversity coefficient; and (3) chemical diversity indicator. A MAP analysis for a single nucleotide substitution was performed for four genes: (1) heme domain of cytochrome P450 BM-3 from Bacillus megaterium (EC 1.14.14.1); (2) glucose oxidase from Aspergillus niger (EC 1.1.3.4); (3) arylesterase from Pseudomonas fluorescens (EC 3.1.1.2); and (4) alcohol dehydrogenase from Saccharomyces cerevisiae (EC 1.1.1.1). Based on the MAP analysis of these four genes, 19 mutagenesis methods have been evaluated and criteria for an ideal mutagenesis method have been proposed. The statistical analysis showed that existing gene mutagenesis methods are limited and highly biased. An average amino acid substitution per residue of only 3.15-7.4 can be achieved with current random mutagenesis methods. For the four investigated gene sequences, an average fraction of amino acid substitutions of 0.5-7% results in stop codons and 4.5-23.9% in glycine or proline residues. An average fraction of 16.2-44.2% of the amino acid substitutions are preserved, and 45.6% (epPCR method) are chemically different. The diversity remains low even when applying a non-biased method: an average of seven amino acid substitutions per residue, 2.9-4.7% stop codons, 11.1-16% glycine/proline residues, 21-25.8% preserved amino acids, and 55.5% are amino acids with chemically different side-chains. Statistical information for each mutagenesis method can further be used to investigate the mutational spectra in protein regions regarded as important for the property of interest.

Amino Acid Substitution↗

PubMed2005

Accounting for global protein deformability during protein-protein and protein-ligand docking.

Computational docking methods are valuable tools aimed to simplify the costly process of drug development and improvement. Most current approaches assume a rigid receptor structure to allow virtual screening of large numbers of possible ligands and putative binding sites on a receptor molecule. However, inclusion of receptor flexibility can be of critical importance since binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a ligand. Recent approaches to efficiently account for receptor flexibility during docking simulations are reviewed. In particular, accounting efficiently for global conformational changes of the protein backbone during docking is a still challenging unsolved problem. An approximate method has recently been suggested that is based on relaxing the receptor conformation during docking in pre-calculated soft collective degrees of freedom (M. Zacharias, Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP, Proteins: Struct., Funct., Genet. 54 (2004) 759-767). Test applications on protein-protein docking and on docking the inhibitor staurosporine to the apo-form of cAMP-dependent protein kinase A catalytic domain indicate significant improvement of docking results compared to rigid docking at a very modest computational demand. Accounting for receptor conformational changes in pre-calculated global degrees of freedom might offer a promising route to improve systematic docking screening simulations.

Binding Sites↗

PubMed2005

Refinement of protein cores and protein-peptide interfaces using a potential scaling approach.

Refinement of side chain conformations in protein model structures and at the interface of predicted protein-protein or protein-peptide complexes is an important step during protein structural modelling and docking. A common approach for side chain prediction is to assume a rigid protein main chain for both docking partners and search for an optimal set of side chain rotamers to optimize the steric fit. However, depending on the target-template similarity in the case of comparative protein modelling and on the accuracy of an initially docked complex, the main chain template structure is only an approximation of a realistic target main chain. An inaccurate rigid main chain conformation can in turn interfere with the prediction of side chain conformations. In the present study, a potential scaling approach (PS-MD) during a molecular dynamics (MD) simulation that also allows the inclusion of explicit solvent has been used to predict side chain conformations on semi-flexible protein main chains. The PS-MD method converges much faster to realistic protein-peptide interface structures or protein core structures than standard MD simulations. Depending on the accuracy of the protein main chain, it also gives significantly better results compared with the standard rotamer search method.

Amino Acid Sequence↗

PubMed2005

Effect of hexafluoroisopropanol alcohol on the structure of melittin: a molecular dynamics simulation study.

The molecular mechanism by which HFIP stabilizes the alpha-helical structure of peptides is not well understood. In the present study, we use melittin as a model to gain insight into the details of the atomistic interactions of HFIP with the peptide. We have performed extensive comparative molecular dynamics simulations (up to 100 nsec) in the absence and in the presence of HFIP. In agreement with recent NMR experiments, the simulations show rapid loss of tertiary structure in water at pH 2 but much higher helicity in 35% HFIP. The MD simulations also indicate that melittin adopts a highly dynamic global structure in 35% HFIP solution with two alpha-helical segments sampling a wide range of angular orientations. The analysis of the HFIP distribution shows the tendency of HFIP to aggregate around the peptide, increasing the local cosolvent concentration to more than two times that in the bulk concentration. The correlation of local peptide structure with HFIP coating suggests that displacement of water at the peptide surface is the main contribution of HFIP in stabilizing the secondary structure of melittin. Finally, a stabilizing effect promoted by the presence of counter-ions was also observed in the simulations.

Computer Simulation↗

PubMed2005

Structural and dynamic properties of cytochrome P450 BM-3 in pure water and in a dimethylsulfoxide/water mixture.

Solvent molecules play an important role for the structural and dynamical properties of proteins. A major focus of current protein engineering is the development of enzymes that are catalytically active in the presence of organic solvents. The monooxygenase P450 BM-3 is one of the best-studied enzymes and promising for industrial applications but with limited activity in the presence of organic solvents or cosolvents. To gain insights into the structural and dynamical properties of the heme domain of this enzyme in solution, molecular dynamics simulations in pure water and in a 14% DMSO/water mixture were performed. The results of the simulations show overall similar structural fluctuations in both solvent systems, with no indication of partial or global unfolding. In 14% DMSO, the regions comprising the helices E, F, and the EF loop (implicated in controlling the entry to the active site channel) undergo a large shift. Significant changes were also observed near the active site access channel at the residue R47. During the simulation, no DMSO molecule penetrated the active site. However, a significant accumulation of DMSO molecules close to the substrate-binding site and to the Flavin Mononucleotide (FMN) reductase domain interface was observed.

Bacillus megaterium↗