Search PubMed⌕ Search

Biomedical subjects

M Nilges

Publications and source records attributed to M Nilges.

At least 19 recordsLinked to original sources
PubMed2000

Efficient sampling in collective coordinate space.

Collective motions in biological macromolecules have been shown to be important for function. The most important collective motions occur on slow time scales, which poses a sampling problem in dynamic simulation of biomolecules. We present a novel method for efficient conformational sampling. The method combines the simulation of an ensemble of concurrent trajectories with restraints acting on the ensemble of structures as a whole. Two properties of the ensemble may be restrained: (i) the variance of the ensemble and (ii) the average position of the ensemble. Both properties are defined in a subspace of collective coordinate space spanned by an arbitrary number of modes. We show that weak restraints on the ensemble variance suffice for an increase in sampling efficiency along soft modes by two orders of magnitudes. The resulting trajectories exhibit virtually the same structural quality as trajectories generated by restraint-free-molecular dynamics simulation, as judged by standard structure validation tools. The method is used to probe the resistance of a structure against conformational changes along collective modes and clearly distinguishes soft from stiff modes. Further applications are discussed. Proteins 2000;39:82-88.

Computer Simulation↗
PubMed2000

Unraveling the symmetry ambiguity in a hexamer: calculation of the R6 human insulin structure.

Crystallographic and NMR studies of insulin have revealed a highly flexible molecule with a range of different aggregation and structural states; the importance of these states for the function of the hormone is still unclear. To address this question, we have studied the solution structure of the insulin R6 symmetric hexamer using NMR spectroscopy. Structure determination of symmetric oligomers by NMR is complicated due to 'symmetry ambiguity' between intra- and intermonomer NOEs, and between different classes of intermonomer NOEs. Hence, to date, only two symmetric tetramers and one symmetric pentamer (VTB, B subunit of verotoxin) have been solved by NMR: there has been no other symmetric hexamer or higher-order oligomer. Recently, we reported a solution structure for R6 insulin hexamer. However, in that study, a crystal structure was used as a reference to resolve ambiguities caused by the threefold symmetry; the same method was used in solving VTB. Here, we have successfully recalculated R6 insulin using the symmetry-ADR method, a computational strategy in which ambiguities are resolved using the NMR data alone. Thus the obtained structure is a refinement of the previous R6 solution structure. Correlated motions in the final structural ensemble were analysed using a recently developed principal component method; this suggests the presence of two major conformational substates. The study demonstrates that the solution structure of higher-order symmetric oligomers can be determined unambiguously from NMR data alone, using the symmetry-ADR method. This success bodes well for future NMR studies of higher-order symmetric oligomers. The correlated motions observed in the structural ensemble suggest a new insight into the mechanism of phenol exchange and the T6 <--> R6 transition of insulin in solution.

Binding Sites↗
PubMed2000

Refinement of the protein backbone angle psi in NMR structure calculations.

Cross-correlated relaxation rates involving the Calpha-Halpha dipolar interaction and the carbonyl (C') chemical shift anisotropy (CSA) have been measured using two complementary 3D experiments. We show that the protein backbone angle psi can be directly refined against such cross-correlated relaxation rates (gammaHalphaCalpha,C') and the three-bond H/D isotope effect on the Calpha chemical shifts (3 deltaCalpha(ND)). By simultaneously using both experimental parameters as restraints during NMR structure calculations, a unique value for the backbone angle psi is defined. We have applied the new refinement method to the alpha-Spectrin SH3 domain (a beta-sheet protein) and to the Sgs1p HRDC domain (an alpha-helical protein) and show that the quality of the NMR structures is substantially improved, judging from the atomic coordinate precision and the Ramachandran map. In addition, the psi-refined NMR structures of the SH3 domain deviate less from the 1.8 A crystal structure, suggesting an improved accuracy. The proposed refinement method can be used to significantly improve the quality of NMR structures and will be applicable to larger proteins.

Carbon↗
PubMed1999

The three-dimensional structure of the HRDC domain and implications for the Werner and Bloom syndrome proteins.

BACKGROUND: The HRDC (helicase and RNaseD C-terminal) domain is found at the C terminus of many RecQ helicases, including the human Werner and Bloom syndrome proteins. RecQ helicases have been shown to unwind DNA in an ATP-dependent manner. However, the specific functional roles of these proteins in DNA recombination and replication are not known. An HRDC domain exists in both of the human RecQ homologues that are implicated in human disease and may have an important role in their function. RESULTS: We have determined the three-dimensional structure of the HRDC domain in the Saccharomyces cerevisiae RecQ helicase Sgs1p by nuclear magnetic resonance (NMR) spectroscopy. The structure resembles auxiliary domains in bacterial DNA helicases and other proteins that interact with nucleic acids. We show that a positively charged region on the surface of the Sgs1p HRDC domain can interact with DNA. Structural similarities to bacterial DNA helicases suggest that the HRDC domain functions as an auxiliary domain in RecQ helicases. Homology models of the Werner and Bloom HRDC domains show different surface properties when compared with Sgs1p. CONCLUSIONS: The HRDC domain represents a structural scaffold that resembles auxiliary domains in proteins that are involved in nucleic acid metabolism. In Sgs1p, the HRDC domain could modulate the helicase function via auxiliary contacts to DNA. However, in the Werner and Bloom syndrome helicases the HRDC domain may have a role in their functional differences by mediating diverse molecular interactions.

Adenosine Triphosphatases↗
PubMed1999

Classification of protein sequences by homology modeling and quantitative analysis of electrostatic similarity.

Protein electrostatics plays a key role in ligand binding and protein-protein interactions. Therefore, similarities or dissimilarities in electrostatic potentials can be used as indicators of similarities or dissimilarities in protein function. We here describe a method to compare the electrostatic properties within protein families objectively and quantitatively. Three-dimensional structures are built from database sequences by comparative modeling. Molecular potentials are then computed for these with a continuum solvation model by finite difference solution of the Poisson-Boltzmann equation or analytically as a multipole expansion that permits rapid comparison of very large datasets. This approach is applied to 104 members of the Pleckstrin homology (PH) domain family. The deviation of the potentials of the homology models from those of the corresponding experimental structures is comparable to the variation of the potential in an ensemble of structures from nuclear magnetic resonance data or between snapshots from a molecular dynamics simulation. For this dataset, the results for analysis of the full electrostatic potential and the analysis using only monopole and dipole terms are very similar. The electrostatic properties of the PH domains are generally conserved despite the extreme sequence divergence in this family. Notable exceptions from this conservation are seen for PH domains linked to a Db1 homology (DH) domain and in proteins with internal PH domain repeats.

Data Interpretation, Statistical↗
PubMed1999

Functionally important correlated motions in the single-stranded DNA-binding protein encoded by filamentous phage Pf3.

To elucidate the interplay between different parts of dimeric single-stranded DNA-binding proteins we have studied the correlated motions in the protein encoded by filamentous phage Pf3 via the combined use of 15N-NMR relaxation experiments, molecular dynamics simulations and essential dynamics calculations. These studies provide insight into the mechanism underlying the protein-DNA binding reaction. The most important motions can be described by a few essential modes. Most outstanding is the correlated symmetric motion of the DNA-binding wings, which are far apart in the structure. This motion determines the access of DNA to the DNA-binding domain. A correlation between the motion of the DNA-binding wing and the complex loop is indicated to play a role in the cooperative binding of the protein to DNA. These motions are in the nanosecond regime in correspondence with the 15N-NMR relaxation experiments.

DNA, Single-Stranded↗
PubMed1999

Exploring protein interiors: the role of a buried histidine in the KH module fold.

The K-homology (KH) module is a novel RNA-binding motif. The structures of a representative KH motif from vigilin (vig-KH6) and of the first KH domain of fmr1 have been recently solved by nuclear magnetic resonance (NMR) and automated assignment-refinement techniques (ARIA). While a hydrophobic residue is found at position 21 in most of the KH modules, a buried His is conserved in all the 15 KH repeats of vigilin. This position must therefore have a key structural role in stabilizing the hydrophobic core. In the present work, we have addressed the following questions in order to obtain a detailed description of the role of His 21: i) what is the exact role of the histidine in the hydrophobic core of vig-KH6? ii) can we define the interactions that allow a conserved buried position to be occupied by a histidine both in vig-KH6 and in the whole vigilin KH sub-family? iii) how is the structure and stability of vig-KH6 influenced by the state of protonation of this histidine? To answer these questions, we have carried out an extensive refinement of the vig-KH6 structure using both an improved ARIA protocol starting from different initial structures and successively running restrained and unrestrained trajectories in water. An analysis of the stability of secondary structural elements, solvent accessibility, and hydrogen bonding patterns allows hypothesis on the structural role of residue His 21 and on the interactions that this residue forms with the environment. The importance of the protonation state of His 21 on the stability of the KH fold was addressed and validated by experimental results.

Carrier Proteins↗
PubMed1999

Molecular dynamics and accuracy of NMR structures: effects of error bounds and data removal.

The effect of internal dynamics on the accuracy of nuclear magnetic resonance (NMR) structures was studied in detail using model distance restraint sets (DRS) generated from a 6.6 nanosecond molecular dynamics trajectory of bovine pancreatic trypsin inhibitor. The model data included the effects of internal dynamics in a very realistic way. Structure calculations using different error estimates were performed with iterative removal of systematically violated restraints. The accuracy of each calculated structure was measured as the atomic root mean square (RMS) difference to the optimized average structure derived from the trajectory by structure factors refinement. Many of the distance restraints were derived from NOEs that were significantly affected by internal dynamics. Depending on the error bounds used, these distance restraints seriously distorted the structure, leading to deviations from the coordinate average of the dynamics trajectory even in rigid regions. Increasing error bounds uniformly for all distance restraints relieved the strain on the structures. However, the accuracy did not improve. Significant improvement of accuracy was obtained by identifying inconsistent restraints with violation analysis, and excluding them from the calculation. The highest accuracy was obtained by setting bounds rather tightly, and removing about a third of the restraints. The limiting accuracy for all backbone atoms was between 0.6 and 0.7 A. Also, the precision of the structures increased with removal of inconsistent restraints, indicating that a high precision is not simply the consequence of tight error bounds but of the consistency of the DRS. The precision consistently overestimated the accuracy.

Amino Acids↗
PubMed1999

Influence of internal dynamics on accuracy of protein NMR structures: derivation of realistic model distance data from a long molecular dynamics trajectory.

In order to study the effect of internal dynamics on the accuracy of NMR structures in detail, we generated NOE distance data from a long molecular dynamics trajectory of BPTI. Cross-relaxation rates were calculated from the trajectory by analysis of the appropriate proton-proton vector autocorrelation functions. A criterion for the convergence of correlation functions was developed, and the analysis was restricted to those correlation functions that had converged within the simulation time. Effective distances were determined from the calculated cross-relaxation rates. Internal dynamics affected the derived distances in a realistic way, since they were subject both to radial averaging (which increases the cross-relaxation rate) and angular averaging (which decreases the cross-relaxation rate). The comparison of the effective distances with average distance between the protons during the trajectory showed that for most the effects of angular and distance averaging essentially cancel out. For these distances, the effective distance derived from an NOE is therefore a very good estimate of the average distance, or the distance in the average structure. However, for about 10% of the distances, the effective distance was more than 10% larger than the average distance, while for about 5%, it was more than 10% smaller, in some cases by more than 2 A. Little correlation is observed between the effects on cross-relaxation rates to different protons of the same residue. The results of this analysis have implications for the way structures are calculated from NOE distance data. For many distances, the assumption of a rigid structure is valid, and large error bounds would result in the loss of too much information content. On the other hand, the error bounds very often employed are not wide enough for some of the effects seen in our study.

Animals↗
PubMed1999

The PH superfold: a structural scaffold for multiple functions.

Pleckstrin homology (PH) domains form a structurally conserved family that is associated with many regulatory pathways within the cell. Domains with a nearly identical fold are found in other families that share no sequence similarity, suggesting the existence of a stable PH superfold. The PH domains generally function as regulated membrane-binding modules that bind to inositol lipids and respond to upstream signals by targeting the host proteins to the correct cellular sites. The other domains with a similar fold, such as the phosphotyrosine binding domains, recognize protein ligands.

Amino Acid Motifs↗
PubMed1999

The structure in solution of the b domain of protein disulfide isomerase.

Protein disulfide isomerase (PDI) is a multifunctional protein of the endoplasmic reticulum, which catalyzes the formation, breakage and rearrangement of disulfide bonds during protein folding. It consists of four domains designated a, b, b and a. Both a and a domains contains an active site with the sequence motif -Cys-Gly-His-Cys-involved directly in thiol-disulfide exchange reactions. As expected these domains have structures very similar to the ubiquitous redox protein thioredoxin. A low-resolution NMR structure of the b domain revealed that this domain adopts a fold similar to the PDI a domain and thioredoxin [Kemmink, J., Darby, N.J., Dijkstra, K., Nilges, M. and Creighton, T.E. (1997) Curr. Biol. 7, 239-245]. A refined ensemble of solution structures based on the input of 1865 structural restraints shows that the structure of PDI b is well defined throughout the complete protein except for about 10 residues at the C-terminus of the sequence. 15N relaxation data show that these residues are disordered and not part of this structural domain. Therefore the domain boundaries of PDI can now be fixed with reasonable precision. Structural comparison of the PDI b domain with thioredoxin and PDI a reveals several features important for thiol-disulfide exchange activity.

Amino Acid Sequence↗
PubMed1998

The three-dimensional structure of a type I module from titin: a prototype of intracellular fibronectin type III domains.

BACKGROUND: Titin is a huge protein ( approximately 3 MDa) that is present in the contractile unit (sarcomere) of striated muscle and has a key role in muscle assembly and elasticity. Titin is mainly composed of two types of module (type I and II). Type I modules are found exclusively in the region of titin localised in the A band, where they are arranged in a super-repeat pattern that correlates with the ultrastructure of the thick filament. No structure of a titin type I module has been reported so far. RESULTS: We have determined the structure of a representative type I module, A71, using nuclear magnetic resonance (NMR) spectroscopy. The structure has the predicted fibronectin type III fold. Titin-specific conserved residues are either located at the putative module-module interfaces or along one side of the protein surface. Several proline residues that contribute to two stretches in a polyproline II helix conformation are solvent-exposed and line up as a continuous ribbon extending over more than two-thirds of the module surface. Homology models of the type I module N-terminal to A71 (A70) and the double module A70-A71 were used to discuss possible intermodule interactions and their role in module-module orientation. CONCLUSIONS: As residues at the module-module interfaces are highly conserved, we speculate that similar interactions govern all of the interfaces between type I modules in titin. This conservation would lead to a regular multiple array of similar surface structures. Such an arrangement would allow arrays of contiguous type I modules to expose multiple proline stretches in a highly regular way and these may act as binding sites for other thick filament proteins.

Amino Acid Sequence↗
PubMed1998

Crystallography & NMR system: A new software suite for macromolecular structure determination.

A new software suite, called Crystallography & NMR System (CNS), has been developed for macromolecular structure determination by X-ray crystallography or solution nuclear magnetic resonance (NMR) spectroscopy. In contrast to existing structure-determination programs, the architecture of CNS is highly flexible, allowing for extension to other structure-determination methods, such as electron microscopy and solid-state NMR spectroscopy. CNS has a hierarchical structure: a high-level hypertext markup language (HTML) user interface, task-oriented user input files, module files, a symbolic structure-determination language (CNS language), and low-level source code. Each layer is accessible to the user. The novice user may just use the HTML interface, while the more advanced user may use any of the other layers. The source code will be distributed, thus source-code modification is possible. The CNS language is sufficiently powerful and flexible that many new algorithms can be easily implemented in the CNS language without changes to the source code. The CNS language allows the user to perform operations on data structures, such as structure factors, electron-density maps, and atomic properties. The power of the CNS language has been demonstrated by the implementation of a comprehensive set of crystallographic procedures for phasing, density modification and refinement. User-friendly task-oriented input files are available for nearly all aspects of macromolecular structure determination by X-ray crystallography and solution NMR.

Computer Simulation↗
PubMed1998

Heteronuclear relaxation study of the PH domain of beta-spectrin: restriction of loop motions upon binding inositol trisphosphate.

The structural dynamics of protein ligand binding sites is one factor determining the specificity towards related ligands. In this context, the spectrin PH domain, which binds to a number of phosphatidylinositol lipid head groups, was investigated with respect to the dynamics of the binding loops. The latter were found to be of intermediate flexibility on a picosecond to nanosecond time-scale in the free protein and become more rigid upon ligand binding. Significant 15N and proton chemical shift changes occur in the binding loops. The internal correlation time, determined from 15N heteronuclear relaxation data using the standard model-free approach, decreases upon ligand binding. For several residues a concomitant rise in the generalized order parameter is observed. This is interpreted as a dampening effect of the ligand on a slow loop motion, while a fast component is not affected. Molecular dynamics simulations were performed to further investigate this situation. In fact, two time-scales of loop motions in the free state are observed in a 9 ns molecular dynamics trajectory. Agreement with generalized order parameters obtained from the experiment improves when a subtrajectory is analyzed that excludes rare dihedral transitions.

Binding Sites↗
PubMed1998

Are there non-trivial dynamic cross-correlations in proteins?

The analysis of internal motion in ensembles of flexible molecules in coordinate space requires the removal of overall motion by a least-squares fitting procedure of the Cartesian coordinates. It has been demonstrated that the choice of the atom set used for fitting influences the picture of the internal motion of BPTI. We have performed essential dynamics analyses of a 1 ns molecular dynamics trajectory of the single-stranded DNA-binding protein from the Pf3 phage using either all alpha-carbon atoms or the least mobile ones for fitting the trajectory prior to the analysis. We found that covariances of atoms separated by long distances were significantly reduced in the latter case; the overall overlap of essential spaces was still high. In the second part we present a method that does not introduce and bias caused by overall motion: principal component analysis in distance space. Non-trivial dynamic cross-correlations were preserved in distance space, which answers the question posed in the title in the affirmative. However, cross-correlations were throughout smaller than those detected by standard essential dynamics analyses.

Bacteriophages↗
PubMed1998

Essential spaces defined by NMR structure ensembles and molecular dynamics simulation show significant overlap.

Large concerted motions of proteins which span its "essential space," are an important component of protein dynamics. We investigate to what extent structure ensembles generated with standard structure calculation techniques such as simulated annealing can capture these motions by comparing them to long-time molecular dynamics (MD) trajectories. The motions are analyzed by principal component analysis and compared using inner products of eigenvectors of the respective covariance matrices. Two very different systems are studied, the beta-spectrin PH domain and the single-stranded DNA binding protein (ssDBP) from the filamentous phage Pf3. A comparison of the ensembles from NMR and MD shows significant overlap of the essential spaces, which in the case of ssDBP is extraordinarily high. The influence of variations in the specifications of distance restraints is investigated. We also study the influence of the selection criterion for the final structure ensemble on the definition of mobility. The results suggest a modified criterion that improves conformational sampling in terms of amplitudes of correlated motion.

Computer Simulation↗
PubMed1997

Solution structure of the spectrin repeat: a left-handed antiparallel triple-helical coiled-coil.

Cytoskeletal proteins belonging to the spectrin family have an elongated structure composed of repetitive units. The three-dimensional solution structure of the 16th repeat from chicken brain alpha-spectrin (R16) has been determined by NMR spectroscopy and distance geometry-simulated annealing calculations. We used a total of 1035 distance restraints, which included 719 NOE-based values obtained by applying the ambiguous restraints for iterative assignment (ARIA) method. In addition, we performed a direct refinement against 1H-chemical shifts. The final ensemble of 20 structures shows an average RMSD of 1.52 A from the mean for the backbone atoms, excluding loops and N and C termini. R16 is made up of three antiparallel alpha-helices separated by two loops, and folds into a left-handed coiled-coil. The basic unit of spectrin is an antiparallel heterodimer composed of two homologous chains, beta and alpha. These assemble a tetramer via a mechanism that relies on the completion of a single repeat by association of the partial repeats located at the C terminus of the beta-chain (two helices) and at the N terminus of the alpha-chain (one helix). This tetramer is the assemblage able to cross-link actin filaments. Model building by homology of the "tetramerization" repeat from human erythrocyte spectrin illuminates the possible role of point mutations which cause hemolytic anemias.

Amino Acid Sequence↗