Search

An automated NMR chemical shift assignment algorithm was developed using multi-objective optimization techniques. The problem is modeled as a combinatorial optimization problem and its objective parameters are defined separately in different score functions. Some of the heuristic approaches of evolutionary optimization are employed in this problem model. Both, a conventional genetic algorithm and multi-objective methods, i.e., the non-dominated sorting genetic algorithms II and III (NSGA2 and NSGA3), are applied to the problem. The multi-objective approaches consider each objective parameter separately, whereas the genetic algorithm followed a conventional way, where all objectives are combined in one score function. Several improvement steps and repetitions on these algorithms are performed and their combinations are also created as a hyper-heuristic approach to the problem. Additionally, a hill-climbing algorithm is also applied after the evolutionary algorithm steps. The algorithms are tested on several different datasets with a set of 11 commonly used spectra. The test results showed that our algorithm could assign both sidechain and backbone atoms fully automatically without any manual interactions. Our approaches could provide around a 65% success rate and could assign some of the atoms that could not be assigned by other methods.

The automated assignment of NOESY cross peaks has become a fundamental technique for NMR protein structure analysis. A widely used algorithm for this purpose is implemented in the program CYANA. It has been used for a large number of structure determinations of proteins in solution but a systematic evaluation of its performance has not yet been reported. In this paper we systematically analyze the reliability of combined automated NOESY assignment and structure calculation with CYANA under a variety of conditions on the basis of the experimental NMR data sets of ten proteins. To evaluate the robustness of the algorithm, the original high-quality experimental data sets were modified in different ways to simulate the effect of data imperfections, i.e. incomplete or erroneous chemical shift assignments, missing NOESY cross peaks, inaccurate peak positions, inaccurate peak intensities, lower dimensionality NOESY spectra, and higher tolerances for the matching of chemical shifts and peak positions. The results show that the algorithm is remarkably robust with regard to imperfections of the NOESY peak lists and the chemical shift tolerances but susceptible to lacking or erroneous resonance assignments, in particular for nuclei that are involved in many NOESY cross peaks.

NMR structure calculation using NOE-derived distance restraints requires a considerable number of assignments of both backbone and sidechains resonances, often difficult or impossible to get for large or complex proteins. Pseudocontact shifts (PCSs) also play a well-established role in NMR protein structure calculation, usually to augment existing structural, mostly NOE-derived, information. Existing refinement protocols using PCSs usually either require a sizeable number of sidechain assignments or are complemented by other experimental restraints. Here, we present an automated iterative procedure to perform backbone protein structure refinements requiring only a limited amount of backbone amide PCSs. Already known structural features from a starting homology model, in this case modules of repeat proteins, are framed into a scaffold that is subsequently refined by experimental PCSs. The method produces reliable indicators that can be monitored to judge about the performance. We applied it to a system in which sidechain assignments are hardly possible, designed Armadillo repeat proteins (dArmRPs), and we calculated the solution NMR structure of YM4A, a dArmRP containing four sequence-identical internal modules, obtaining high convergence to a single structure. We suggest that this approach is particularly useful when approximate folds are known from other techniques, such as X-ray crystallography, while avoiding inherent artefacts due to, for instance, crystal packing.

Folding of G-protein coupled receptors (GPCRs) according to the two-stage model (Popot, J. L., and Engelman, D. M. (1990) Biochemistry 29, 4031–4037) is postulated to proceed in 2 steps: partitioning of the polypeptide into the membrane followed by diffusion until native contacts are formed. Herein we investigate conformational preferences of fragments of the yeast Ste2p receptor using NMR. Constructs comprising the first, the first two, and the first three transmembrane (TM) segments, as well as a construct comprising TM1–TM2 covalently linked to TM7 were examined. We observed that the isolated TM1 does not form a stable helix nor does it integrate well into the micelle. TM1 is significantly stabilized upon interaction with TM2, forming a helical hairpin reported previously (Neumoin, A., Cohen, L. S., Arshava, B., Tantry, S., Becker, J. M., Zerbe, O., and Naider, F. (2009) Biophys. J. 96, 3187–3196), and in this case the protein integrates into the hydrophobic interior of the micelle. TM123 displays a strong tendency to oligomerize, but hydrogen exchange data reveal that the center of TM3 is solvent exposed. In all GPCRs so-far structurally characterized TM7 forms many contacts with TM1 and TM2. In our study TM127 integrates well into the hydrophobic environment, but TM7 does not stably pack against the remaining helices. Topology mapping in microsomal membranes also indicates that TM1 does not integrate in a membrane-spanning fashion, but that TM12, TM123, and TM127 adopt predominantly native-like topologies. The data from our study would be consistent with the retention of individual helices of incompletely synthesized GPCRs in the vicinity of the translocon until the complete receptor is released into the membrane interior.

Lantibiotics are peptide-derived antibiotics that inhibit the growth of Gram-positive bacteria via interactions with lipid II and lipid II-dependent pore formation in the bacterial membrane. Due to their general mode of action the Gram-positive producer strains need to express immunity proteins (LanI proteins) for protection against their own lantibiotics. Little is known about the immunity mechanism protecting the producer strain against its own lantibiotic on the molecular level. So far, no structures have been reported for any LanI protein. We solved the structure of SpaI, a LanI protein from the subtilin producing strain Bacillus subtilis ATCC 6633. SpaI is a 16.8-kDa lipoprotein that is attached to the outside of the cytoplasmic membrane via a covalent diacylglycerol anchor. SpaI together with the ABC transporter SpaFEG protects the B. subtilis membrane from subtilin insertion. The solution-NMR structure of a 15-kDa biologically active C-terminal fragment reveals a novel fold. We also demonstrate that the first 20 N-terminal amino acids not present in this C-terminal fragment are unstructured in solution and are required for interactions with lipid membranes. Additionally, growth tests reveal that these 20 N-terminal residues are important for the immunity mediated by SpaI but most likely are not part of a possible subtilin binding site. Our findings are the first step on the way of understanding the immunity mechanism of B. subtilis in particular and of other lantibiotic producing strains in general.

To achieve efficient proton pumping in the light-driven proton pump bacteriorhodopsin, the protein must be tightly coupled to the retinal to rapidly convert retinal isomerization into protein structural rearrangements. Methyl group dynamics of bR embedded in lipid nanodiscs were determined in the dark-adapted state, and were found to be mostly well-ordered at the cytosolic side. Methyl groups in the M145A mutant of bR, which displays only 10% residual proton pumping activity, are less well ordered suggesting a link between side chain dynamics on the cytosolic side of the bR cavity and proton pumping activity. In addition, slow conformational exchange, attributed to low frequency motions of aromatic rings, was indirectly observed for residues on the extracellular side of the bR cavity. This may be related to reorganization of the water network. These observations provide a detailed picture of previously undescribed equilibrium dynamics on different time scales for ground-state bR.