Molecular Structure Refinement Based on Residual Dipolar Couplings: A Comparison of the Molecular Rotational-Sampling Method with the Alignment-Tensor Approach.

In NMR experiments, residual dipolar couplings (RDCs) in a molecule can be measured by averaging the dipolar couplings (DCs) over the rotational motion of a molecule in an environment that induces a slight anisotropic orientation distribution of the molecule. Since the shape of the anisotropic distribution cannot be measured, it is standard practice to use a particular orientation distribution of the molecule with respect to the magnetic field, in the form of a so-called alignment tensor (AT), to calculate RDC-values for the molecule. Since the same alignment tensor is commonly used to calculate the different RDCs of a molecule, this approach rests on the assumption that the rotational motion of the molecule is decoupled from its internal motions and that the molecule is rigid. The validity of these two assumptions is investigated for a small, simple molecule, using a relatively rigid atomic interaction function or force field and a more flexible one. By simulating the molecule using an orientation-biasing force an anisotropic rotational distribution can be generated, for which RDCs can be obtained. Using these RDCs as target RDCs when applying one of the two approaches of structure refinement based on RDCs, it can be investigated how well the target RDCs are approximated in the RDC restraining and whether the corresponding nonuniform orientation distribution is reproduced. For the relatively rigid version of the molecule, the AT approach reproduces the target RDC-values, although the nonuniform orientation distribution of the angle θab,H between the vector r⃗ab connecting two atoms a and b in the molecule and the vector representing the direction of the magnetic field H⃗ as generated in the orientation-biasing simulation cannot be reproduced in the AT RDC-restraining simulation. For the relatively flexible version of the molecule, the AT approach fails to reproduce both the target RDC values and the nonuniform orientation distribution. For biomolecules with flexible parts, the application of the AT approach is thus not recommended. Instead, a method based on sampling of the rotational and internal degrees of freedom of the molecule should be applied in molecular structure determination or refinement based on measured RDCs.

A Method to Derive Structural Information on Molecules from Residual Dipolar Coupling NMR Data.

A method for structure refinement of molecules based on residual dipolar coupling (RDC) data is proposed. It calculates RDC values using rotational and molecule-internal configurational sampling instead of the common refinement procedure that is based on the approximation of the nonuniform rotational distribution of the molecule by a single alignment tensor representing the average nonuniformity of this distribution. Using rotational sampling, as is occurring in the experiment leading to observable RDCs, the method stays close to the experiment. It avoids the use of an alignment tensor and thus the assumption that the overall rotation of the molecule is decoupled from its internal motions and that the molecule be rigid. Two simple molecules, two-united-atomic ethane and a cyclooctane molecule with eight side chains, containing 24 united atoms, serve as the so-called "toy model" test systems. The method demonstrates the influence of molecular flexibility and force-field deficiencies on the outcome of structure refinement based on RDCs. For a molecule of a given size (number of atoms Nat), there must be a sufficiently large number NRDC of measured RDC values available to allow the restraining forces to bias the overall orientation distribution of the molecule. If the ratio NRDC/Nat gets too small, the RDC-restraining forces will either not be strong enough to change the overall rotational direction of the molecule such that the target RDC values are approximated well or will be so strong that they induce a local deformation of the molecule. In the latter case, the size or inertia of the molecule hinders a restraining-induced overall rotation and the internal structure of the molecule is not strong enough to avoid local deformation due to the restraining forces.

Molecular dynamics simulation or structure refinement of proteins: are solvent molecules required? A case study using hen lysozyme.

In protein simulation or structure refinement based on values of observable quantities measured in (aqueous) solution, solvent (water) molecules may be explicitly treated, omitted, or represented by a potential of mean-solvation-force term, depending on protein coordinates only, in the force field used. These three approaches are compared for hen egg white lysozyme (HEWL). This 129-residue non-spherical protein contains a variety of secondary-structure elements, and ample experimental data are available: 1630 atom-atom Nuclear Overhauser Enhancement (NOE) upper distance bounds, 213 3 J-couplings and 200 S2 order parameters. These data are used to compare the performance of the three approaches. It is found that a molecular dynamics (MD) simulation in explicit water approximates the experimental data much better than stochastic dynamics (SD) simulation in vacuo without or with a solvent-accessible-surface-area (SASA) implicit-solvation term added to the force field. This is due to the missing energetic and entropic contributions and hydrogen-bonding capacities of the water molecules and the missing dielectric screening effect of this high-permittivity solvent. Omission of explicit water molecules leads to compaction of the protein, an increased internal strain, distortion of exposed loop and turn regions and excessive intra-protein hydrogen bonding. As a consequence, the conformation and dynamics of groups on the surface of the protein, which may play a key role in protein-protein interactions or ligand or substrate binding, may be incorrectly modelled. It is thus recommended to include water molecules explicitly in structure refinement of proteins in aqueous solution based on nuclear magnetic resonance (NMR) or other experimentally measured data.

On the use of intra-molecular distance and angle constraints to lengthen the time step in molecular and stochastic dynamics simulations of proteins.

Computer simulation of proteins in aqueous solution at the atomic level of resolution is still limited in time span and system size due to limited computing power available and thus employs a variety of time-saving techniques that trade some accuracy against computational effort. An example of such a time-saving technique is the application of constraints to particular degrees of freedom when integrating Newton's or Langevin's equations of motion in molecular dynamics (MD) or stochastic dynamics (SD) simulations, respectively. The application of bond-length constraints is standard practice in protein simulations and allows for a lengthening of the time step by a factor of three. Applying recently proposed algorithms to constrain bond angles or dihedral angles, it is investigated, using the protein trypsin inhibitor as test molecule, whether bond angles and dihedral angles involving hydrogen atoms or even stiff proper (torsional) dihedral angles as well as improper ones (maintaining particular tetrahedral or planar geometries) may be constrained without generating too many artificial side effects. Constraining the relative positions of the hydrogen atoms in the protein allows for a lengthening of the time step by a factor of two. Additionally constraining the improper dihedral angles and the stiff proper (torsional) dihedral angles in the protein does not allow for an increase of the MD or SD time step.

On the use of multiple-time-step algorithms to save computing effort in molecular dynamics simulations of proteins.

Computer simulation of proteins in aqueous solution at the atomic level of resolution is still limited in time span and system size due to limited computing power available and thus employs a variety of time-saving techniques that trade some accuracy against computational effort. Examples of such time-saving techniques are the application of constraints to particular degrees of freedom or the use of a multiple-time-step (MTS) algorithm distinguishing between particular forces when integrating Newton's equations of motion. The application of two types of MTS algorithms to bond-stretching forces versus the remaining forces in molecular dynamics (MD) simulations of a protein in aqueous solution or of liquid water is investigated and the results in terms of total energy conservation and the influence on various other properties are compared to those of MD simulations of the same systems using bond-length, and for water bond-angle, constraints. At comparable computational effort, the use of bond-length constraints in proteins leads to better energy conservation and less distorted properties than the two MTS algorithms investigated.

On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio-Molecular Systems: Overview and Perspective on Issues.

Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.

A method to apply bond-angle constraints in molecular dynamics simulations.

An algorithm to apply bond-angle constraints in molecular dynamics simulations of macromolecules or molecular liquids is presented. It uses Cartesian coordinates and determines the Lagrange multipliers required for maintaining the constraints iteratively. It constitutes an alternative to the use of only distance constraints (DCs) between particles to maintain a particular geometry. DCs are unsuitable to maintain particular, for example, linear or flat, geometries of molecules. The proposed algorithm can easily handle bond-length, bond-angle, and dihedral-angle constraints simultaneously, as when calculating a potential of mean force along a dihedral-angle degree of freedom.

Algorithms to apply dihedral-angle constraints in molecular or stochastic dynamics simulations.

Various algorithms to apply dihedral-angle constraints in molecular dynamics or stochastic dynamics simulations of molecular systems are presented, investigated, and tested. They use Cartesian coordinates and determine the Lagrangian multipliers necessary for maintaining the constraints iteratively. The most suitable algorithm to maintain a dihedral-angle constraint is numerically compared to the alternative to use distance constraints to this end. It can easily be used to obtain a potential of mean force along a dihedral-angle coordinate.

Saturation Mutagenesis by Efficient Free-Energy Calculation.

Single-point mutations in proteins can greatly influence protein stability, binding affinity, protein function or its expression per se. Here, we present accurate and efficient predictions of the free energy of mutation of amino acids. We divided the complete mutational free energy into an uncharging step, which we approximate by a third-power fitting (TPF) approach, and an annihilation step, which we approximate using the one-step perturbation (OSP) method. As a diverse set of test systems, we computed the solvation free energy of all amino acid side chain analogues and obtained an excellent agreement with thermodynamic integration (TI) data. Moreover, we calculated mutational free energies in model tripeptides and established an efficient protocol involving a single reference state. Again, the approximate methods agreed excellently with the TI references, with a root-mean-square error of only 3.6 kJ/mol over 17 mutations. Our combined TPF+OSP approach does show not only a very good agreement but also a 2-fold higher efficiency than full blown TI calculations.

Exploring genetic suppression interactions on a global scale.

Genetic suppression occurs when the phenotypic defects caused by a mutation in a particular gene are rescued by a mutation in a second gene. To explore the principles of genetic suppression, we examined both literature-curated and unbiased experimental data, involving systematic genetic mapping and whole-genome sequencing, to generate a large-scale suppression network among yeast genes. Most suppression pairs identified novel relationships among functionally related genes, providing new insights into the functional wiring diagram of the cell. In addition to suppressor mutations, we identified frequent secondary mutations,in a subset of genes, that likely cause a delay in the onset of stationary phase, which appears to promote their enrichment within a propagating population. These findings allow us to formulate and quantify general mechanisms of genetic suppression.

Studying metal ion binding properties of a three-way junction RNA by heteronuclear NMR.

Self-splicing group II introns are highly structured RNA molecules, containing a characteristic secondary and catalytically active tertiary structure, which is formed only in the presence of Mg(II). Mg(II) initiates the first folding step governed by the κζ element within domain 1 (D1κζ). We recently solved the NMR structure of D1κζ derived from the mitochondrial group II intron ribozyme Sc.ai5γ and demonstrated that Mg(II) is essential for its stabilization. Here, we performed a detailed multinuclear NMR study of metal ion interactions with D1κζ, using Cd(II) and cobalt(III)hexammine to probe inner- and outer-sphere coordination of Mg(II) and thus to better characterize its binding sites. Accordingly, we mapped (1)H, (15)N, (13)C, and (31)P spectral changes upon addition of different amounts of the metal ions. Our NMR data reveal a Cd(II)-assisted macrochelate formation at the 5'-end triphosphate, a preferential Cd(II) binding to guanines in a helical context, an electrostatic interaction in the ζ tetraloop receptor and various metal ion interactions in the GAAA tetraloop and κ element. These results together with our recently published data on Mg(II) interaction provide a much better understanding of Mg(II) binding to D1κζ, and reveal how intricate and complex metal ion interactions can be.

Multiple Binding Poses in the Hydrophobic Cavity of Bee Odorant Binding Protein AmelOBP14.

In the first step of olfaction, odorants are bound and solubilized by small globular odorant binding proteins (OBPs) which shuttle them to the membrane of a sensory neuron. Low ligand affinity and selectivity at this step enable the recognition of a wide range of chemicals. Honey bee Apis mellifera's OBP14 (AmelOBP14) binds different plant odorants in a largely hydrophobic cavity. In long molecular dynamics simulations in the presence and absence of ligand eugenol, we observe a highly dynamic C-terminal region which forms one side of the ligand-binding cavity, and the ligand drifts away from its crystallized orientation. Hamiltonian replica exchange simulations, allowing exchanges of conformations sampled by the real ligand with those sampled by a noninteracting dummy molecule and several intermediates, suggest an alternative, quite different ligand pose which is adopted immediately and which is stable in long simulations. Thermodynamic integration yields binding free energies which are in reasonable agreement with experimental data.

Protonation-Dependent Base Flipping at Neutral pH in the Catalytic Triad of a Self-Splicing Bacterial Group†II Intron.

NMR spectroscopy has revealed pH-dependent structural changes in the highly conserved catalytic domain 5 of a bacterial group II intron. Two adenines with pK(a) values close to neutral pH were identified in the catalytic triad and the bulge. Protonation of the adenine opposite to the catalytic triad is stabilized within a G(syn)-AH(+) (anti) base pair. The pH-dependent anti-to-syn flipping of this G in the catalytic triad modulates the known interaction with the linker region between domains 2 and 3 (J23) and simultaneously the binding of the catalytic Mg(2+) ion to its backbone. Hence, this here identified shifted pK(a) value controls the conformational change between the two steps of splicing.

Processing of complex N-glycans in IgG Fc-region is affected by core fucosylation.

We investigated N-glycan processing of immunoglobulin G1 using the monoclonal antibody cetuximab (CxMab), which has a glycosite in the Fab domain in addition to the conserved Fc glycosylation, as a reporter. Three GlcNAc (Gn) terminating bi-antennary glycoforms of CxMab differing in core fucosylation (α1,3- and α1,6-linkage) were generated in a plant-based expression platform. These GnGn, GnGnF(3), and GnGnF(6) CxMab variants were subjected in vivo to further processing toward sialylation and GlcNAc diversification (bisected and branching structures). Mass spectrometry-based glycan analyses revealed efficient processing of Fab glycans toward envisaged structures. By contrast, Fc glycan processing largely depend on the presence of core fucose. A particularly strong support of glycan processing in the presence of plant-specific core α1,3-fucose was observed. Consistently, molecular modeling suggests changes in the interactions of the Fc carbohydrate chain depending on the presence of core fucose, possibly changing the accessibility. Here, we provide data that reveal molecular mechanisms of glycan processing of IgG antibodies, which may have implications for the generation of glycan-engineered therapeutic antibodies with improved efficacies.

Insights into structural features determining odorant affinities to honey bee odorant binding protein 14.

Molecular interactions between odorants and odorant binding proteins (OBPs) are of major importance for understanding the principles of selectivity of OBPs towards the wide range of semiochemicals. It is largely unknown on a structural basis, how an OBP binds and discriminates between odorant molecules. Here we examine this aspect in greater detail by comparing the C-minus OBP14 of the honey bee (Apis mellifera L.) to a mutant form of the protein that comprises the third disulfide bond lacking in C-minus OBPs. Affinities of structurally analogous odorants featuring an aromatic phenol group with different side chains were assessed based on changes of the thermal stability of the protein upon odorant binding monitored by circular dichroism spectroscopy. Our results indicate a tendency that odorants show higher affinity to the wild-type OBP suggesting that the introduced rigidity in the mutant protein has a negative effect on odorant binding. Furthermore, we show that OBP14 stability is very sensitive to the position and type of functional groups in the odorant.

Structure and conformational dynamics of the domain 5 RNA hairpin of a bacterial group II intron revealed by solution nuclear magnetic resonance and molecular dynamics simulations.

Nuclear magnetic resonance (NMR) nuclear Overhauser enhancement (NOE) data obtained for a 35-nucleotide RNA segment of a bacterial group II intron indicate a helical hairpin structure in which three parts, a terminal pentaloop, a bulge, and a G-A mismatch, display no Watson-Crick base pairing. The 668 NOE upper distance bounds for atom pairs are insufficient to uniquely determine the conformation of these segments. Therefore, molecular dynamics simulations including time-averaged distance restraints have been used to obtain a conformational ensemble compatible with the observed NMR data. The ensemble shows alternating hydrogen bonding patterns for the mentioned segments. In particular, in the pentaloop and in the bulge, the hydrogen bonding networks correspond to distinct conformational clusters that could not be captured by using conventional single-structure refinement techniques. This implies that, to obtain a realistic picture of the conformational ensemble of such flexible biomolecules, it is necessary to properly account for the conformational variability in the structure refinement of RNA fragments.

The structural stabilization of the κ three-way junction by Mg(II) represents the first step in the folding of a group II intron.

Folding of group II introns is characterized by a first slow compaction of domain 1 (D1) followed by the rapid docking of other domains to this scaffold. D1 compaction initiates in a small subregion encompassing the κ and ζ elements. These two tertiary elements are also the major interaction sites with domain 5 to form the catalytic core. Here, we provide the first characterization of the structure adopted at an early folding step and show that the folding control element can be narrowed down to the three-way junction with the κ motif. In our nuclear magnetic resonance studies of this substructure derived from the yeast mitochondrial group II intron Sc.ai5γ, we show that a high affinity Mg(II) ion stabilizes the κ element and enables coaxial stacking between helices d' and d'', favoring a rigid duplex across the three-way junction. The κ-element folds into a stable GAAA-tetraloop motif and engages in A-minor interactions with helix d'. The addition of cobalt(III)hexammine reveals three distinct binding sites. The Mg(II)-promoted structural rearrangement and rigidification of the D1 core can be identified as the first micro-step of D1 folding.

Characterization of metal ion-nucleic acid interactions in solution.

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

Can soaked-in scavengers protect metalloprotein active sites from reduction during data collection?

One of the first events taking place when a crystal of a metalloprotein is exposed to X-ray radiation is photoreduction of the metal centres. The oxidation state of a metal cannot always be determined from routine X-ray diffraction experiments alone, but it may have a crucial impact on the metal's environment and on the analysis of the structural data when considering the functional mechanism of a metalloenzyme. Here, UV-Vis microspectrophotometry is used to test the efficacy of selected scavengers in reducing the undesirable photoreduction of the iron and copper centres in myoglobin and azurin, respectively, and X-ray crystallography to assess their capacity of mitigating global and specific radiation damage effects. UV-Vis absorption spectra of native crystals, as well as those soaked in 18 different radioprotectants, show dramatic metal reduction occurring in the first 60 s of irradiation with an X-ray beam from a third-generation synchrotron source. Among the tested radioprotectants only potassium hexacyanoferrate(III) seems to be capable of partially mitigating the rate of metal photoreduction at the concentrations used, but not to a sufficient extent that would allow a complete data set to be recorded from a fully oxidized crystal. On the other hand, analysis of the X-ray crystallographic data confirms ascorbate as an efficient protecting agent against radiation damage, other than metal centre reduction, and suggests further testing of HEPES and 2,3-dichloro-1,4-naphtoquinone as potential scavengers.