Simulation of oriented NMR spectra: Combining molecular dynamics and chemical shift tensor calculations.

Solid state NMR is widely used to study the orientation and other structural features of proteins and peptides in lipid bilayers. Using data obtained by PISEMA (Polarization Inversion Spin Exchange at Magic Angle) experiments, periodic spectral patterns arise from well-aligned α-helical molecules. Significant problems in the interpretation of PISEMA spectra may arise for systems that do not form perfectly defined secondary structures, like α-helices, or the signal pattern is disturbed by molecular motion. Here, we present a new method that combines molecular dynamics simulation with tensorial orientational constraints (MDOC) and chemical shift tensor calculations for the simulation and interpretation of PISEMA-like spectra. The calculations include the spectra arising from non α-helical molecules and molecules with non-uniform intrinsic mobility. In a first step, dipolar or quadrupolar interaction tensors drive molecular rotations and reorientations to obtain the proper mean values as observed in corresponding NMR experiments. In a second step, the coordinate snapshots of the MDOC simulations are geometry optimized with the isotropic 15 N chemical shifts as constraints using Bond Polarization Theory (BPT) to provide reliable 15 N CS tensor data. The averaged dipolar 1 H-15 N couplings and the δzz tensor components can then be combined to simulate PISEMA patterns. We apply this method to the ß-helical peptide gramicidin A (gA) and demonstrate that this method enables the assignment of most PISEMA resonances. In addition, MDOC simulations provide local order parameters for the calculated sites. These local order parameters reveal large differences in backbone mobility between L- and D-amino acids of gA.

Residual dipolar couplings as a tool for structural analysis of ionic liquids.

An acrylonitrile/dimethylacrylamide cross-linked polymer could be swollen in different imidazolium ionic liquids. Mechanical compression of the obtained polymer gels inside an NMR tube allowed the measurement of residual dipolar couplings. Conformational analysis of the 1-methyl-3-butyl-imidazolium (BMIM) cation could be performed by including the measured RDCs as restraints in time-averaged molecular dynamics.

Structures Controlled by Entropy: The Flexibility of Strychnine as Example.

To study the flexibility of strychnine, we performed molecular dynamics simulations with orientational tensorial constraints (MDOC). Tensorial constraints are derived from nuclear magnetic resonance (NMR) interaction tensors, for instance, from residual dipolar couplings (RDCs). Used as orientational constraints, they rotate the whole molecule and molecular parts with low rotational barriers. Since the NMR parameters are measured at ambient temperatures, orientational constraints generate conformers that populate the whole landscape of Gibbs free energy. In MDOC, structures are populated that are not only controlled by energy but by the entropy term TΔS of the Gibbs free energy. In the case of strychnine, it is shown that ring conformers are populated, which has not been discussed in former investigations. These conformer populations are not only in accordance with RDCs but fulfill nuclear Overhauser effect (NOE)-derived distance constraints and 3JHH couplings as well.

Statistical evaluation of simulated NMR data of flexible molecules.

A new probability score-named χ-probability-is introduced for evaluating the fit of mixed NMR datasets to calculate molecular model ensembles, in order to answer challenging structural questions such as the determination of stereochemical configurations. Similar to the DP4 parameter, the χ-probability is based on Bayes theorem and expresses the probability that an experimental NMR dataset fits to a given individual within a finite set of candidate structures or configurations. Here, the χ-probability is applied to single out the correct configuration in four example cases, with increasing complexity and conformational mobility. The NMR data (which include RDCs, NOE distances and 3J couplings) are calculated from MDOC (Molecular Dynamics with Orientational Constraints) trajectories and are investigated against experimentally measured data. It is demonstrated that this approach singles out the correct stereochemical configuration with probabilities more than 98%, even for highly mobile molecules. In more demanding cases, a decisive χ-probability test requires that the datasets include high-quality NOE distances in addition to RDC values.

The simulation of NMR data of flexible molecules: sagittamide A as an example for MD simulations with orientational constraints.

The recently developed MDOC (Molecular Dynamics with Orientational Constraints) simulation is applied for the first time to a fully flexible molecule. MDOC simulations aim to single out the naturally existing configuration of molecules and to elucidate conformer populations. The performance of the method was first demonstrated on a well-studied test case, the five-membered ring lactone (α-methylene-γ-butyrolactone). In the case of sagittamide A, one-bond 1H-13C residual dipolar couplings (RDC) are used as orientational constraints that reorient the molecule or parts of it. In addition, NOE distances and 3J scalar couplings are used as constraints. Five possible configurations of sagittamide A (labelled a to e) are considered. One experimental RDC value per flexible unit was available and this was not sufficient to single out one valid configuration. The problem could be solved by including NOE distances as well as 3J couplings as complementary constraints into the MDOC simulations. In accordance with former investigations, we confirmed the configuration a for the natural product. A detailed analysis of conformers of the central chain of 6 chiral carbon atoms could be given by inspecting the MDOC trajectory. The relative abundance of these conformers is crucial in fulfilling all three sets of constraints.

Conformational Investigations in Flexible Molecules Using Orientational NMR Constraints in Combination with 3J-Couplings and NOE Distances.

The downscaling of NMR tensorial interactions, such as dipolar couplings, from tens of kilohertz to a few hertz in low-order media is the result of dynamics spanning several orders of magnitudes, including vibrational modes (~ns-fs), whole-molecule reorientation (~ns) and higher barrier internal conformational exchange (

Configuration determination by residual dipolar couplings: accessing the full conformational space by molecular dynamics with tensorial constraints.

Residual dipolar couplings (RDCs) and other residual anisotropic NMR parameters provide valuable structural information of high quality and quantity, bringing detailed structural models of flexible molecules in solution in reach. The corresponding data interpretation so far is directly or indirectly based on the concept of a molecular alignment tensor, which, however, is ill-defined for flexible molecules. The concept is typically applied to a single or a small set of lowest energy structures, ignoring the effect of vibrational averaging. Here, we introduce an entirely different approach based on time averaged molecular dynamics with dipolar couplings as tensorial orientational restraints that can be used to solve structural problems in molecules of any size without the need of introducing an explicit molecular alignment tensor into the computation. RDC restraints are represented by their full 3D interaction tensor in the laboratory frame, for which pseudo forces are calculated using a secular dipolar Hamiltonian as the target. The resulting rotational averaging of each individual tensorial restraint leads to structural ensembles that best fulfil the experimental data. Using one-bond RDCs, the approach has been implemented in the force field procedures of the program COSMOS and extensively tested. A concise theoretical introduction, including the special treatment of force fields for stable and fast MD simulations, as well as applications regarding configurational analyses of small to medium-sized organic molecules with different degrees of flexibility, is given. The observed results are discussed in detail.

Investigation of backbone dynamics and local geometry of bio-molecules using calculated NMR chemical shifts and anisotropies.

Prerequisite for chemical shift (CS) and CS tensor calculations are highly refined structures defining the molecular surroundings of the nuclei under study. Here, we present geometry optimizations with 13C and 15N CS constraints for large bio-molecules like peptides and proteins. The method discussed here provides both, refined structures and chemical shift tensors. Furthermore, since the experimental resonances of aligned systems are related to CS tensors, they strongly depend on the orientation and motion of molecules, their fragments, functional groups and moieties. For efficient CS calculations we apply a semi-empirical approach-the bond polarization theory (BPT). The BPT relies on linear bond polarization parameters and we present a new set of parameters based on ab initio second-order Møller-Plesset perturbation theory calculations. The new parametrization extends the applicability of the BPT approach to a wide range of organic molecules and bio-polymers. Here, the method has been applied to the protein ubiquitin and the membrane-active peptide gramicidin A (dimer) in oriented bilayers. The calculated 13C and 15N CS values of best-refined structures published until now gave a large scatter with respect to the experiment. It will be shown that BPT CS optimizations can reduce these errors to values near the experimental uncertainty. In combination with molecular dynamics with orientational constraints it is possible to study motional dynamics and BPT calculations can provide residual chemical shift anisotropies.

Molecular Dynamics with Orientational Tensorial Constraints: A New Approach to Probe the Torsional Angle Distributions of Small Rotationally Flexible Molecules.

The potential of residual dipolar couplings (RDCs) in conformational studies of small molecules is now widely recognized, but current theoretical approaches for their interpretation have several limitations and there is still the need for a general method to probe the torsional angle distributions applicable to any rotationally flexible molecule. Molecular dynamics simulations with RDC-based orientational tensorial constraints (MDOC), implemented in the software COSMOS, are presented here as a conceptually new strategy. For the cases of the fluorinated anti-inflammatory drug diflunisal and the disaccharide cellobiose, we demonstrate that MDOC simulations with one-bond RDCs as tensorial constraints unveil torsion distributions and allow the determination of relative configuration in the presence of rotational flexibility. The independence of the initial structure or any a priori assumption as well as the possibility to combine different experimental constraints represent features, which make the COSMOS software a promising tool for the investigation of torsional angle distributions of flexible molecules, regardless of their size and degree of freedom.

Insight into the chromophore of rhodopsin and its Meta-II photointermediate by 19F solid-state NMR and chemical shift tensor calculations.

19F nuclei are useful labels in solid-state NMR studies, since their chemical shift and tensor elements are very sensitive to the electrostatic and space-filling properties of their local environment. In this study we have exploited a fluorine substituent, strategically placed at the C-12-position of 11-cis retinal, the chromophore of visual rhodopsins. This label was used to explore the local environment of the chromophore in the ground state of bovine rhodopsin and its active photo-intermediate Meta II. In addition, the chemical shift and tensor elements of the chromophore in the free state in a membrane environment and the bound state in the protein were determined. Upon binding of the chromophore into rhodopsin and Meta II, the isotropic chemical shift changes in the opposite direction by +9.7 and -8.4 ppm, respectively. An unusually large isotropic shift difference of 35.9 ppm was observed between rhodopsin and Meta II. This partly originates in the light-triggered 11-cis to all-trans isomerization of the chromophore. The other part reflects the local conformational rearrangements in the chromophore and the binding pocket. These NMR data were correlated with the available X-ray structures of rhodopsin and Meta II using bond polarization theory. For this purpose hydrogen atoms have to be inserted and hereto a family of structures were derived that best correlated with the well-established 13C chemical shifts. Based upon these structures, a 12-F derivative was obtained that best corresponded with the experimentally determined 19F chemical shifts and tensor elements. The combined data indicate strong changes in the local environment of the C-12 position and a substantially different interaction pattern with the protein in Meta II as compared to rhodopsin.

1H line width dependence on MAS speed in solid state NMR - Comparison of experiment and simulation.

Recent developments in magic angle spinning (MAS) technology permit spinning frequencies of ≥100 kHz. We examine the effect of such fast MAS rates upon nuclear magnetic resonance proton line widths in the multi-spin system of ß-Asp-Ala crystal. We perform powder pattern simulations employing Fokker-Plank approach with periodic boundary conditions and 1H-chemical shift tensors calculated using the bond polarization theory. The theoretical predictions mirror well the experimental results. Both approaches demonstrate that homogeneous broadening has a linear-quadratic dependency on the inverse of the MAS spinning frequency and that, at the faster end of the spinning frequencies, the residual spectral line broadening becomes dominated by chemical shift distributions and susceptibility effects even for crystalline systems.

Molecular dynamics simulations on PGLa using NMR orientational constraints.

NMR data obtained by solid state NMR from anisotropic samples are used as orientational constraints in molecular dynamics simulations for determining the structure and dynamics of the PGLa peptide within a membrane environment. For the simulation the recently developed molecular dynamics with orientational constraints technique (MDOC) is used. This method introduces orientation dependent pseudo-forces into the COSMOS-NMR force field. Acting during a molecular dynamics simulation these forces drive molecular rotations, re-orientations and folding in such a way that the motional time-averages of the tensorial NMR properties are consistent with the experimentally measured NMR parameters. This MDOC strategy does not depend on the initial choice of atomic coordinates, and is in principle suitable for any flexible and mobile kind of molecule; and it is of course possible to account for flexible parts of peptides or their side-chains. MDOC has been applied to the antimicrobial peptide PGLa and a related dimer model. With these simulations it was possible to reproduce most NMR parameters within the experimental error bounds. The alignment, conformation and order parameters of the membrane-bound molecule and its dimer were directly derived with MDOC from the NMR data. Furthermore, this new approach yielded for the first time the distribution of segmental orientations with respect to the membrane and the order parameter tensors of the dimer systems. It was demonstrated the deuterium splittings measured at the peptide to lipid ratio of 1/50 are consistent with a membrane spanning orientation of the peptide.

A Toxoplasma gondii Gluconeogenic Enzyme Contributes to Robust Central Carbon Metabolism and Is Essential for Replication and Virulence.

The expression of gluconeogenic enzymes is typically repressed when glucose is available. The protozoan parasite Toxoplasma gondii utilizes host glucose to sustain high rates of intracellular replication. However, despite their preferential utilization of glucose, intracellular parasites constitutively express two isoforms of the gluconeogenic enzyme fructose 1,6-bisphosphatase (TgFBP1 and TgFBP2). The rationale for constitutive expression of FBPases in T. gondii remains unclear. We find that conditional knockdown of TgFBP2 results in complete loss of intracellular growth in vitro under glucose-replete conditions and loss of acute virulence in mice. TgFBP2 deficiency was rescued by expression of catalytically active FBPase and was associated with altered glycolytic and mitochondrial TCA cycle fluxes, as well as dysregulation of glycolipid, amylopectin, and fatty acid biosynthesis. Futile cycling between gluconeogenic and glycolytic enzymes may constitute a regulatory mechanism that allows T. gondii to rapidly adapt to changes in nutrient availability in different host cells.

Structural characterization of a peptoid with lysine-like side chains and biological activity using NMR and computational methods.

N-substituted glycine oligomers or peptoids with charged side chains are a novel class of cell penetrating peptide mimetics and have been shown to serve as drug delivery agents. Here, we investigated by NMR spectroscopy and quantum chemical calculations whether a Rhodamine B labelled peptoid [RhoB(Spiro)-Ahx]-[But](6A)NH(2) with lysine-like side chains adopts structural motifs similar to regular peptides. Due to a low chemical shift dispersion, high resolution structure determination with conventional NMR-derived distance restraints and J-couplings was not possible. Instead, a combined assignment and structure refinement strategy using the QM/MM force field COSMOS-NMR was developed to interpret the highly ambiguous chemical shift and distance constraints and obtain a medium resolution three-dimensional structural model. This allowed us to select for the all cis-amide conformation of the peptide with a pseudo-helical arrangement of extended side chains as a faithful representative structure of [RhoB(Spiro)-Ahx]-[But](6A)NH(2). We tested the biological activity of the peptoid by live-cell imaging, which showed that the cellular uptake of the peptoid was comparable to conventional cell-penetrating peptides.

Structure simulation with calculated NMR parameters - integrating COSMOS into the CCPN framework.

The Collaborative Computing Project for NMR (CCPN) has build a software framework consisting of the CCPN data model (with APIs) for NMR related data, the CcpNmr Analysis program and additional tools like CcpNmr FormatConverter. The open architecture allows for the integration of external software to extend the abilities of the CCPN framework with additional calculation methods. Recently, we have carried out the first steps for integrating our software Computer Simulation of Molecular Structures (COSMOS) into the CCPN framework. The COSMOS-NMR force field unites quantum chemical routines for the calculation of molecular properties with a molecular mechanics force field yielding the relative molecular energies. COSMOS-NMR allows introducing NMR parameters as constraints into molecular mechanics calculations. The resulting infrastructure will be made available for the NMR community. As a first application we have tested the evaluation of calculated protein structures using COSMOS-derived 13C Cα and Cß chemical shifts. In this paper we give an overview of the methodology and a roadmap for future developments and applications.

Rapid calculation of protein chemical shifts using bond polarization theory and its application to protein structure refinement.

Although difficult to analyze, NMR chemical shifts provide detailed information on protein structure. We have adapted the semi-empirical bond polarization theory (BPT) to protein chemical shift calculation and chemical shift driven protein structure refinement. A new parameterization for BPT amide nitrogen chemical shift calculation has been derived from MP2 ab initio calculations and successfully evaluated using crystalline tripeptides. We computed the chemical shifts of the small globular protein ubiquitin, demonstrating that BPT calculations can match the results obtained at the DFT level of theory at very low computational cost. In addition to the calculation of chemical shift tensors, BPT allows the calculation of chemical shift gradients and consequently chemical shift driven geometry optimizations. We applied chemical shift driven protein structure refinement to the conformational analysis of a set of Trypanosoma brucei (the causative agent of African sleeping sickness) tryparedoxin peroxidase Px III structures. We found that the interaction of Px III with its reaction partner Tpx seems to be governed by conformational selection rather than by induced fit.

Irregular structure of the HIV fusion peptide in membranes demonstrated by solid-state NMR and MD simulations.

To better understand peptide-induced membrane fusion at a molecular level, we set out to determine the structure of the fusogenic peptide FP23 from the HIV-1 protein gp41 when bound to a lipid bilayer. An established solid-state (19)F nuclear magnetic resonance (NMR) approach was used to collect local orientational constraints from a series of CF(3)-phenylglycine-labeled peptide analogues in macroscopically aligned membranes. Fusion assays showed that these (19)F-labels did not significantly affect peptide function. The NMR spectra were characteristic of well-behaved samples, without any signs of heterogeneity or peptide aggregation at 1:300 in 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC). We can conclude from these NMR data that FP23 has a well-defined (time-averaged) conformation and undergoes lateral diffusion in the bilayer plane, presumably as a monomer or small oligomer. Attempts to evaluate its conformation in terms of various secondary structures, however, showed that FP23 does not form any type of regular helix or ß-strand. Therefore, all-atom molecular dynamics (MD) simulations were carried out using the orientational NMR constraints as pseudo-forces to drive the peptide into a stable alignment and structure. The resulting picture suggests that FP23 can adopt multiple ß-turns and insert obliquely into the membrane. Such irregular conformation explains why the structure of the fusion peptide could not be reliably determined by any biophysical method so far.

Analysis of the specific interactions between the lectin domain of malectin and diglucosides.

The endoplasmic reticulum malectin is a highly conserved protein in the animal kingdom that has no counterpart so far in lower organisms. We recently determined the structure of its conserved domain and found a highly selective binding to Glc(2)Man(9)GlcNAc(2), an intermediate of N-glycosylation. In our quest for putative ligands during the initial characterization of the protein, we noticed that the malectin domain is highly specific for diglucosides but quite tolerant towards the linkage of the glucosidic bond. To understand the molecular requirements for the observed promiscuity of the malectin domain, here we analyze the binding to a range of diglucosides through comparison of the protein chemical shift perturbation patterns and the saturation transfer difference spectra of the ligands including two maltose-mimicking drugs. A comparison of the maltose-bound structure of the malectin domain with the complex of the native ligand nigerose reveals why malectin is able to tolerate such a diversity of ligands.

Two essential arginine residues in the T components of energy-coupling factor transporters.

Energy-coupling factor (ECF) transporters, a recently discovered class of importers of micronutrients, are composed of a substrate-specific transmembrane component (S component) and a conserved energy-coupling module consisting of a transmembrane protein (T component) and pairs of ABC ATPases (A proteins). Based on utilization of a dedicated (subclass I) or shared (subclass II) energy-coupling module, ECF systems fall into two subclasses. The T components are the least-characterized proteins of ECF importers, and their function is essentially unknown. Using RcBioN and LmEcfT, the T units of the subclass I biotin transporter (RcBioMNY) of a gram-negative bacterium and of the subclass II folate, pantothenate, and riboflavin transporters of a lactic acid bacterium, respectively, we analyzed the role of two strongly conserved short motifs, each containing an arginine residue. Individual replacement of the two Arg residues in RcBioN reduced ATPase activity, an indicator of the transporter function, by two-thirds without affecting the modular assembly of the RcBioMNY complex. A double Arg-to-Glu replacement destroyed the complex and abolished ATPase activity. The corresponding single mutation in motif II of LmEcfT, as well as a double mutation, led to loss of the T unit from the subclass II ECF transporters and inactivated these systems. A single Arg-to-Glu replacement in motif I, however, abolished vitamin uptake activity without affecting assembly of the modules. Our results indicate that the conserved motif I in T components is essential for intramolecular signaling and, in cooperation with motif II, for subunit assembly of modular ECF transporters.

Calculation of fluorine chemical shift tensors for the interpretation of oriented (19)F-NMR spectra of gramicidin A in membranes.

A semi-empirical method for the prediction of chemical shifts, based on bond polarization theory, has recently been introduced for (13)C. Here, we extended this approach to calculate the (19)F chemical shift tensors of fluorine bound to aromatic rings and in aliphatic CF(3) groups. For the necessary parametrization, ab initio chemical shift calculations were performed at the MP2 level for a set of fluorinated molecules including tryptophan. The bond polarization parameters obtained were used to calculate the (19)F chemical shift tensors for several crystalline molecules, and to reference the calculated values on a chemical shift scale relative to CFCl(3). As a first biophysical application, we examined the distribution of conformations of a (19)F-labeled tryptophan side chain in the membrane-bound ion channel peptide, gramicidin A. The fluorine chemical shift tensors were calculated from snapshots of a molecular dynamics simulation employing the (19)F-parametrized bond polarization theory. In this MD simulation, published (2)H quadrupolar and (15)N-(1)H dipolar couplings of the indole ring were used as orientational constraints to determine the conformational distribution of the 5F-Trp(13) side chain. These conformations were then used to interpret the spectra of (19)F-labeled gramicidin A in fluid and gel phase lipid bilayers.