Determination of three-dimensional solution structure of waglerin I,a toxin from Trimeresurus wagleri,using 2D-NMR and molecular dynamics simulation

The solution conformation of a synthetic snake venom toxin waglerin I, has been determined by using proton nuclear magnetic resonance spectroscopy. By y a combination of various two-dimensional NMR techniques, the ¹H-NMR spectrum of waglerin I was completely assigned. A set of 247 interproton distance restraints was derived from nuclear Overhauser enhancement (NOE) measurements. These NOE constraints, in addition to the 2 dihedral angle restraints (from coupling constant measurements) and 7 ω torsion angle restraints for prolines, formed the basis of three-dimensional structure determined by molecular dynamics techniques. The 19 structures that were obtained satisfy the experimental restraints, and display small deviation from idealized covalent geometry. Analysis of converged structures indicates that the toxin has no special secondary structure. In the solution structure of waglerin I, the central ring region is well defined but the N- and C-termini possesses more disorder.     

Use of photoaffinity crosslinking and molecular modeling to analyze the global architecture of ribonuclease P RNA.

Bacterial ribonuclease P (RNase P), an endonuclease involved in tRNA maturation, is a ribonucleoprotein containing a catalytic RNA. The secondary structure of this ribozyme is well established, but comparatively little is understood about its 3-D structure. In this analysis, orientation and distance constraints between elements within the Escherichia coli RNase P RNA-pre-tRNA complex were determined by intra- and intermolecular crosslinking experiments. A molecular mechanics-based RNA structure refinement protocol was used to incorporate the distance constraints indicated by crosslinking, along with the known secondary structure of RNase P RNA and the tertiary structure of tRNA, into molecular models. Seven different structures that satisfy the constraints equally well were generated and compared by superposition to estimate helix positions and orientations. Manual refinement within the range of conformations indicated by the molecular mechanics analysis was used to derive a model of RNase P RNA with bound substrate pre-tRNA that is consistent with the crosslinking results and the available phylogenetic comparisons.     

Improvement of RNA secondary structure prediction using RNase H cleavage and randomized oligonucleotides

RNA secondary structure prediction using free energy minimization is one method to gain an approximation of structure. Constraints generated by enzymatic mapping or chemical modification can improve the accuracy of secondary structure prediction. We report a facile method that identifies single-stranded regions in RNA using short, randomized DNA oligonucleotides and RNase H cleavage. These regions are then used as constraints in secondary structure prediction. This method was used to improve the secondary structure prediction of Escherichia coli 5S rRNA. The lowest free energy structure without constraints has only 27% of the base pairs present in the phylogenetic structure. The addition of constraints from RNase H cleavage improves the prediction to 100% of base pairs. The same method was used to generate secondary structure constraints for yeast tRNA^Phe, which is accurately predicted in the absence of constraints (95%). Although RNase H mapping does not improve secondary structure prediction, it does eliminate all other suboptimal structures predicted within 10% of the lowest free energy structure. The method is advantageous over other single-stranded nucleases since RNase H is functional in physiological conditions. Moreover, it can be used for any RNA to identify accessible binding sites for oligonucleotides or small molecules.     

Crystallographic evaluation of internal motion of human alpha-lactalbumin refined by full-matrix least-squares method

The low temperature form of human alpha-lactalbumin (HAL) was crystallized from a 2H2O solution and its structure was refined to the R value of 0.119 at 1.15 A resolution by the full-matrix least-squares method. Average estimated standard deviations of atomic parameters for non-hydrogen atoms were 0.038 A for coordinates and 0.044 A2 for anisotropic temperature factors (Uij). The magnitude of equivalent isotropic temperature factors (Ueqv) was highly correlated with the distance from the molecular centroid and fitted to a quadratic equation as a function of atomic coordinates. The atomic thermal motion was rather isotropic in the core region and the anisotropy increased towards the molecular surface. The statistical analysis revealed the out-of-plane motion of main-chain oxygen atoms, indicating that peptide groups are in rotational vibration around a Calpha.Calpha axis. The TLS model, which describes the rigid-body motion in terms of translation, libration, and screw motions, was adopted for the evaluation of the molecular motion and the TLS parameters were determined by the least-squares fit to Uij. The reproduced Ueqvcal from the TLS parameters was in fair agreement with observed Ueqv, but differences were found in regions of residues, 5-22, 44-48, 70-75, and 121-123, where Ueqv was larger than Ueqvcal because of large local motions. To evaluate the internal motion of HAL, the contribution of the rigid-body motion was determined to be 42.4 % of Ueqv in magnitude, which was the highest estimation to satisfy the condition that the Uijint tensors of the internal motion have positive eigen values. The internal motion represented with atomic thermal ellipsoids clearly showed local motions different from those observed in chicken-type lysozymes which have a backbone structure very similar to HAL. The result indicates that the internal motion is closely related to biological function of proteins.     

Refined crystal structure of deoxyhemoglobin S. I. Restrained least-squares refinement at 3.0-A resolution

The crystal structure of deoxyhemoglobin S has been refined at 3.0-A resolution using the Hendrickson-Konnert restrained least-squares method. Comparison with the structure of deoxyhemoglobin A reveals a hingelike movement of the beta-chain A helices, which are involved in molecular contacts, toward the EF corners of their respective subunits. This movement brings the amino termini of the beta-chains closer to the molecular dyad. The A helices remain alpha-helical throughout their entire lengths. No other major structural difference is found between deoxyhemoglobin A and deoxyhemoglobin S.     

Simulated annealing with restrained molecular dynamics using a flexible restraint potential: theory and evaluation with simulated NMR constraints.

A new functional representation of NMR-derived distance constraints, the flexible restraint potential, has been implemented in the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168) for molecular structure generation. In addition, flat-bottomed restraint potentials for representing dihedral angle and vicinal scalar coupling constraints have been introduced into CONGEN. An effective simulated annealing (SA) protocol that combines both weight annealing and temperature annealing is described. Calculations have been performed using ideal simulated NMR constraints, in order to evaluate the use of restrained molecular dynamics (MD) with these target functions as implemented in CONGEN. In this benchmark study, internuclear distance, dihedral angle, and vicinal coupling constant constraints were calculated from the energy-minimized X-ray crystal structure of the 46-amino acid polypeptide crambin (ICRN). Three-dimensional structures of crambin that satisfy these simulated NMR constraints were generated using restrained MD and SA. Polypeptide structures with extended backbone and side-chain conformations were used as starting conformations. Dynamical annealing calculations using extended starting conformations and assignments of initial velocities taken randomly from a Maxwellian distribution were found to adequately sample the conformational space consistent with the constraints. These calculations also show that loosened internuclear constraints can allow molecules to overcome local minima in the search for a global minimum with respect to both the NMR-derived constraints and conformational energy. This protocol and the modified version of the CONGEN program described here are shown to be reliable and robust, and are applicable generally for protein structure determination by dynamical simulated annealing using NMR data.     

Pseudo-structures for the 20 common amino acids for use in studies of protein conformations by measurements of intramolecular proton-proton distance constraints with nuclear magnetic resonance

"Pseudo-structures" of the 20 common amino acid residues are introduced for use in protein spatial structure determinations, which rely on the use of intramolecular proton-proton distance constraints determined by nuclear Overhauser effects as input for distance geometry calculations. The proposed structures satisfy requirements for the initial structural interpretation of the nuclear magnetic resonance data that arise from the absence of stereospecific assignments and/or limited spectral resolution for certain resonance lines. The pseudo-atoms used as reference points for the experimental distance constraints can be used in conjunction with the real amino acid structures representing the van der Waals' constraints on the spatial molecular structure, or with simplified models in order to reduce the computing time for the distance geometry calculations.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

Further refinement of the structure of yeast tRNAPhe.

We have refined the monoclinic crystal structure of yeast tRNA^Phe against a complete set of X-ray data at 2.5 Å resolution, using real-space refinement and a combination of energy minimisation and crystallographic least-squares. This refinement has allowed us to define the conformation of residue D16, and to make corrections to Y37 and A76. We have found an additional magnesium binding site (making a total of four), a number of water molecules, and a possible spermine molecule.A revised list of torsion angles is given.     

Three-dimensional structure of echistatin and dynamics of the active site

Summary The snake venom protein echistatin contains the cell recognition sequence Arg-Gly-Asp and is a potent inhibitor of platelet aggregation. The three-dimensional structure of echistatin and the dynamics of the active RGD site are presented. A set of structures was determined using the Distance Geometry method and subsequently refined by Molecular Dynamics and energy minimization. Disulfide pairings are suggested, based on violations of experimental constraints. The structures satisfy 230 interresidue distance constraints, derived from nuclear Overhauser effect measurements, five hydrogen-bonding constraints, and 21 torsional constraints from vicinal spin-spin coupling constants. The segment from Gly⁵ to Cys²⁰ and from Asp³⁰ to Asn⁴² has a well-defined conformation and the Arg-Gly-Asp sequence, which adopts a turn-like structure, is located at the apex of a nine-residue loop connecting the two strands of a distorted -sheet. The mobility of the Arg-Gly-Asp site has been quantitatively characterized by ¹⁵N relaxation measurements. The overall correlation time of echistatin was determined from fluorescence measurements, and was used in a model-free analysis to determine internal motional parameters. The active site has order parameters of 0.3–0.5, i.e., among the smallest values ever observed at the active site of a protein. Correlation of the flexible region of the protein as characterized by relaxation experiments and the NMR solution structures was made by calculating generalized order parameters from the ensemble of three-dimensional structures. The motion of the RGD site detected experimentally is more extensive than a simple RGD loop wagging motional model, suggested by an examination of superposed solution structures.     

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 ¹⁹F nuclear magnetic resonance (NMR) approach was used to collect local orientational constraints from a series of CF₃-phenylglycine-labeled peptide analogues in macroscopically aligned membranes. Fusion assays showed that these ¹⁹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.     

Transverse relaxation optimized 3D and 4D 15N/15N separated NOESY experiments of 15N labeled proteins

Distance constraints are an important complement to orientational constraints. While a high-resolution monomer structure of the ion channel forming polypeptide, gramicidin A, has been solved with 120 orientational constraints, the precise geometry of the dimer interface has not been characterized. Here, using both ¹³C and ¹⁵N labeled gramicidin A samples in hydrated phospholipid bilayers, both inter- and intramolecular distances have been measured with a recently developed simultaneous frequency and amplitude modulation (SFAM) solid-state NMR scheme. Using this approach ¹⁵N-¹³C₁ residual dipolar couplings across a hydrogen bond as small as 20 ± 2 Hz have been characterized. While such distances are on the order of 4.2 ± 0.2 Å, the spectroscopy is complicated by rapid global motion of the molecular structure about the bilayer normal and channel axis. Consequently, the nominal 40 Hz dipolar coupling is averaged depending on the orientation of the internuclear vector with respect to the motional axis. The intermolecular distance confirmed the previously described monomeric structure, while the intramolecular distance across the monomer–monomer interface defined this junction and confirmed the previous model of this interface.     

Methods for multiple-marker mapping of quantitative trait loci in half-sib populations

In this paper we consider the detection of individual loci controlling quantitative traits of interest (quantitative trait loci or QTLs) in the large half-sib family structure found in some species. Two simple approaches using multiple markers are proposed, one using least squares and the other maximum likelihood. These methods are intended to provide a relatively fast screening of the entire genome to pinpoint regions of interest for further investigation. They are compared with a more traditional single-marker least-squares approach. The use of multiple markers is shown to increase power and has the advantage of providing an estimate for the location of the QTL. The maximum-likelihood and the least-squares approaches using multiple markers give similar power and estimates for the QTL location, although the likelihood approach also provides estimates of the QTL effect and sire heterozygote frequency. A number of assumptions have been made in order to make the likelihood calculations feasible, however, and computationally it is still more demanding than the least-squares approach. The least-squares approach using multiple markers provides a fast method that can easily be extended to include additional effects.     

3D structure of bovine pancreatic ribonuclease A in aqueous solution: An approach to tertiary structure determination from a small basis of1H NMR NOE correlations

Summary A method is proposed to generate initial structures in cases where the distance geometry method may fail, such as when the set of¹H NMR NOE-based distance constraints is small in relation to the size of the protein. The method introduces an initial correlation between the and backbone angles (based on empirical observations) which is relaxed in later stages of the calculation. The obtained initial structures are refined by well-established methods of energy minimization and restrained molecular dynamics. The method is applied to determine the solution structure of Ribonuclease A (124 residues) from a NOE basis consisting of 467 NOE cross-correlations (97 intra-residue, 206 sequential, 23 medium-range and 141 long-range) obtained at 360 MHz. The global shape and backbone overall fold of the eight final refined structures are close to those shown by the crystal structure. A meaningful difference in the positioning of the catalytically important His¹¹⁹ side chain in the solution and crystal structures has been detected.     

The structure and antiviral activity of ribavirin analogs. III. Molecular and crystal structure of 1-(2-hydroxyethoxymethyl)-1,2,4-triazole-carboxamide

The crystal and molecular structure of a ribavirin acyclic analogue, 1-(2-hydroxyethoxymethyl)-1,2,4-triazole-3-carboxamide, has been determined by X-ray diffraction method. The space group is P1, unit cell parameters: a = 5,237, b = 6,960, c = 11,483 A, alpha = 93,89, beta = 97,43, gamma = 94,26 degrees; Z = 2. The structure was solved by the direct method and refined by least-squares procedure to R = 3.7%. Two molecular conformers statistically coexist in the unit cell, differing in the hydroxyethoxymethyl group conformation. Trans-conformation about O4'-C4' bond and gauche about C4'-C5' bond are observed in both molecules. C1'-O4' bond is approximately perpendicular to the aglicon.     

Crystal and molecular structure of the inhibitor eglin from leeches in complex with subtilisin Carlsberg

The crystal structure of the molecular complex of eglin, a serine proteinase inhibitor from leeches, with subtilisin Carlsberg has been determined at 2.0 A resolution by the molecular replacement method. The complex has been refined by restrained-parameter least-squares. The present crystallographic R factor (Formula: see text) is 0.183. Eglin is a member of the potato inhibitor 1 family, a group of serine proteinase inhibitors lacking disulfide bonds. Eglin shows strong structural homology to CI-2, a related inhibitor from barley seeds. The structure of subtilisin Carlsberg in this complex is very similar to the known structure from barley seeds. The structure of subtilisin Carlsberg in this complex is very similar to the known structure of subtilisin novo, despite changes of 84 out of 274 amino acids.     

Crystallographic refinement of yeast aspartic acid transfer RNA

The structure of yeast transfer RNA aspartic acid has been refined in one crystal form to 3 Å resolution using the restrained least-squares method of Hendrickson and Konnert and real-space fitting using the FRODO program of Jones. The final Crystallographic discrepancy index R is 23.5% for 4585 reflections with magnitudes twice their standard deviations between 10 and 3 Å. With lower occupancies for some residues of the D-loop, the phosphate U1, and the base U33, the R-factor is 22.3%. The adaptation of the restrained least-squares program for nucleic acids and the progress of the refinement are described. The conformations are analysed with respect to stereochemistry and folding of the backbone. The contacts and hydrogen bonds of the secondary structure are compared with those of yeast tRNA^Phe. The presence of only four bases in the variable loop, instead of five as in yeast tRNA^Phe, leads to a rotation of residue 48 and a lateral movement of residue 46. These two rearrangements induce different environments for [U8 … A14] … A21 as well as for A9 and G45. Otherwise, all tertiary contacts observed in yeast tRNA^Phe are present in yeast tRNA^Asp, except for the absence of hydrogen-bonding between G18 of the D-loop and C56 of the T-loop. The presence of anticodon triplet pairing leads to a distribution of temperature factors different from that observed in yeast tRNA^Phe with a stabilization of the AC stem-and-loop and a destabilization of the T and D-loops. We are inclined to suggest that the labilization of the interactions between the T and D-loops is a consequence of the interaction of the anticodon triplets of symmetry-related molecules through hydrogen bonding, which mimics the interaction between the anticodon and its cognate codon on the messenger RNA.     

The crystal and molecular structure of a 3:2 mixture of laminarabiose and O-α-d-glucopyranosyl-(1→3)-β-d-glucopyranose

The crystal and molecular structure of a 3:2 mixture of laminarabiose and 3-O-α-d-glucopyranosyl-β-d-glucopyranose has been determined by X-ray diffraction. The crystal belongs to the monoclinic system, space group P2, a 14.778(1), b 4.794(1), c 10.516(1) Å and β 98.10(1)^c, Dm 1.54 g. cm^-3, Z 2. The structure was solved by the direct method and refined by the block-diagonal, least-squares procedure to R 0.057 for 1034 observed reflections. Difference synthesis showed all hydrogen atoms and indicated a partial (～39%), random substitution of the β anomer molecules by the α anomer molecules, which are accompanied by water molecules on the crystallographic two-fold axis (～19%). The molecule shows a conformation, different from the fully-extended one, which is stabilized by an intramolecular hydrogen-bond between O-1-H and O-5 [2.786(7) Å]. The ring-to-ring conformation can be described as (Φ, Ψ)=(27.9·–37.5·), according to the definition of Sathyanarayana and Rao, and it is located in the comparatively low-energy region of the energy-contour diagram of laminarabiose. Four intermolecular hydrogen-bonds hold molecules together to form infinite sheets, which are approximately parallel to the ab-plane and linked by additional hydrogen-bonds in the c-direction.     

The horizontal magnetic dance of the honeybee is compatible with a single-domain ferromagnetic magnetoreceptor

Although honeybees are able to sense the geomagnetic field, very little is known about the method in which they are able to detect it. The recent discovery of biochemically precipitated magnetite (Fe₃O₄) in bees, however, suggests the possibility that they might use a simple compass organelle for magnetoreception. If so, their orientation accuracy ought to be related to the accuracy of the compass, e.g. it should be poor in weak background fields and enhanced in strong fields. When dancing to the magnetic directions on a horizontal honeycomb, bees clearly show this type of alignment behavior. A least-squares fit between the expected alignment of a compass and this horizontal dance data is consistent with this hypothesis, and implies that the receptors have magnetic moments of 5 × 10^?13 emu, or magnetite volumes near 10^?15 cm³. Additional considerations suggests that these crystals are slightly sub-spherical and single-domain in size, held symmetrically in their receptors, and have a magnetic orientation energy of approximately to 6 kT in the geomagneticfield. A model of a magnetite-based magnetoreceptor consistent with these constraints is discussed.     

Estimation of protein-ligand binding parameters from contiuous ultrafiltration results.

A method of analyzing the result of continuous ultrafiltration experiments to obtain protein-ligand binding parameters is presented. This method employs a nonlinear least-squares regression algorithm coupled with a model of protein-ligand binding which alloww the computation of free ligand concentrations, and a second-order Runge-Kutta method to integrate free concentrations with respect to collected ultrafiltrate. The approach is general and effectively removes the constraints on maximum fraction size imposed by other methods.