The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Spatial arrangement of the three alpha helices in the solution conformation of E. coli lac repressor DNA-binding domain

The relative orientations of the 3 helices in the DNA-binding domain ('headpiece') of lac repressor have been determined using distance constraints obtained from 2-dimensional 1H nuclear Overhauser enhancement spectra. The relative orientations of its helices is similar to that of the central 3 helices in the DNA-binding domain of the lambda repressor of the bacteriophage lambda.     

Analysis of myosin light and heavy chain types in single human skeletal muscle fibers

In this study, myosin types in human skeletal muscle fibers were investigated with electrophoretic techniques. Single fibers were dissected out of lyophilized surgical biopsies and typed by staining for myofibrillar ATPase after preincubation in acid or alkaline buffers. After 14C-labelling of the fiber proteins in vitro by reductive methylation, the myosin light chain pattern was analysed on two-dimensional gels and the myosin heavy chains were investigated by one-dimensional peptide mapping. Surprisingly, human type I fibers, which contained only the slow heavy chain, were found to contain variable amounts of fast myosin light chains in addition to the two slow light chains LC1s and LC2s. The majority of the type I fibers in normal human muscle showed the pattern LC1s, LC2s and LC1f. Further evidence for the existence in human muscle of a hybrid myosin composed of a slow heavy chain with fast and slow light chains comes from the analysis of purified human myosin in the native state by pyrophosphate gel electrophoresis. With this method, a single band corresponding to slow myosin was obtained; this slow myosin had the light chain composition LC1s, LC2s and LC1f. Type IIA and IIB fibers, on the other hand, revealed identical light chain patterns consisting of only the fast light chains LC1f, LC2f and LC3f but were found to have different myosin havy chains. On the basis of the results presented, we suggest that the histochemical ATPase normally used for fibre typing is determined by the myosin heavy chain type (and not by the light chains). Thus, in normal human muscle a number of 'hybrid' myosins were found to occur, namely two extreme forms of fast myosins which have the same light chains but different heavy chains (IIA and IIB) and a continuum of slow forms consisting of the same heavy chain and slow light chains with a variable fast light chain composition. This is consistent with the different physiological roles these fibers are thought to have in muscle contraction.     

Assignment of 1H, 13C, and 15N resonances of YgiT, a putative DNA interacting protein from E. coli, containing one HTH and two CxxC motifs

The 131 residues protein encoded by the open reading frame ygiT of E. coli contains two characteristic domains: a zinc finger protein-like structure with two CxxC motives at its N-terminus and a helix-turn-helix (HTH) motif at its C-terminus. We report the backbone and side chain ¹H, ¹³C, and ¹⁵N resonances assignment of YgiT.     

Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.     

Fully automated sequence-specific resonance assignments of hetero- nuclear protein spectra

Full automation of the analysis of spectra is a prerequisite for high-throughput NMR studies in structural or functional genomics. Sequence-specific assignments often form the major bottleneck. Here, we present a procedure that yields nearly complete backbone and side chain resonance assignments starting from a set of heteronuclear three-dimensional spectra. Neither manual intervention, e.g., to correct lists obtained from peak picking before feeding these to an assignment program, nor protein-specific information, e.g., structures of homologous proteins, were required. By combining two earlier published procedures, AUTOPSY [Koradi et al. (1998) J. Magn. Reson., 135, 288–297] and GARANT [Bartels et al. (1996) J. Biomol. NMR, 7, 207–213], with a new program, PICS, all necessary steps from spectra analyses to sequence-specific assignments were performed fully automatically. Characteristic features of the present approach are a flexible design allowing as input almost any combination of NMR spectra, applicability to side chains, robustness with respect to parameter choices (such as noise levels) and reproducibility. In this study, automated resonance assignments were obtained for the 14 kD blue copper protein azurin from P. aeruginosa using five spectra: HNCACB, HNHA, HCCH-TOCSY, ¹⁵N-NOESY-HSQC and ¹³C-NOESY-HSQC. Peaks from these three-dimensional spectra were filtered and calibrated with the help of two two-dimensional spectra: ¹⁵N-HSQC and ¹³C-HSQC. The rate of incorrect assignments is less than 1.5% for backbone nuclei and about 3.5% when side chain protons are also considered.     

Optimizing resolution in multidimensional NMR by three-way decomposition

Resolution depends on the number of points sampled in a FID; in indirectly detected dimensions it is an important determinant of the total experiment time. Based on the high redundancy present in NMR data, we propose the following timesaving scheme for three-dimensional spectra. An extensive grid of discrete t1- and t2-values is used, which increases resolution while preserving the spectral width. Total experiment time is reduced by avoiding the recording of t3-FIDs for selected pairs of t1 and t2; typically the recording is omitted for about 75% of the (t1,t2) combinations. These data sets are referred to as sparse, and post-experimental processing making optimal use of spectral redundancy provides the missing, non-recorded data. We have previously shown that three-way decomposition (TWD) within the MUNIN approach provides a practical way to process dense NMR data sets. Here, a novel TWD algorithm [Ibraghimov, (2002) Numer. Linear Algebra Appl. 9, 551–565] is used to complement a sparselyrecorded time-domain data set by providing the missing FIDs for all (t1,t2) combinations omitted in the experiment. A necessary condition is that for each t1-value at least a few FIDs are recorded, and similar for each t2-value. The method is demonstrated on non-uniformly sampled ¹⁵N-NOESY-HSQC data sets recorded for the 14 kD protein azurin. The spectra obtained by TWD, reconstruction and ordinary transform to frequency-domain are, in spite of the large number of signals and the high dynamic range typical for NOESYs, highly similar to a corresponding reference spectrum, for which all (t1,t2) combinations were recorded.     

Measles virus matrix protein specifies apical virus release and glycoprotein sorting in epithelial cells

In polarized epithelial cells measles virus (MV) is predominantly released at the apical cell surface, irrespective of the sorting of its two envelope glycoproteins F and H. It has been reported previously that the viral matrix (M) protein modulates the fusogenic capacity of the viral envelope glycoproteins. Here, extant MV mutants and chimeras were used to determine the role of M protein in the transport of viral glycoproteins and release of progeny virions in polarized epithelial CaCo2 cells. In the absence of M, envelope glycoproteins are sorted to the basolateral surface, suggesting that they possess intrinsic basolateral sorting signals. However, interactions of M with the glycoprotein cytoplasmic tails allow M-glycoprotein co-segregation to the apical surface, suggesting a vectorial function of M to retarget the glycoproteins for apical virion release. Whereas this may allow virus airway shedding, the intrinsic sorting of the glycoproteins to the basolateral surface may account for systemic host infection by allowing efficient cell-cell fusion.