Structure and Dynamics of Coil-like Molecules by Residual Dipolar Couplings

Nuclear magnetic resonance (NMR) spectroscopy provides us with many means to study biological macromolecules in solution. Proteins in particular are the most intriguing targets for NMR studies. Protein functions are usually ascribed to specific three-dimensional structures but more recently tails, long loops and non-structural polypeptides have also been shown to be biologically active. Examples include prions, -synuclein, amylin and the NEF HIV-protein. However, conformational preferences in coil-like molecules are difficult to study by traditional methods. Residual dipolar couplings (RDCs) have opened up new opportunities; however their analysis is not trivial. Here we show how to interpret RDCs from these weakly structured molecules. The most notable residual dipolar couplings arise from steric obstruction effects. In dilute liquid crystalline media as well as in anisotropic gels polypeptides encounter nematogens. The shape of a polypeptide conformation limits the encounter with the nematogen. The most elongated conformations may come closest whereas the most compact remain furthest away. As a result there is slightly more room in the solution for the extended than for the compact conformations. This conformation-dependent concentration effect leads to a bias in the measured data. The measured values are not arithmetic averages but essentially weighted averages over conformations. The overall effect can be calculated for random flight chains and simulated for more realistic molecular models. Earlier there was an implicit thought that weakly structured or non-structural molecules would not yield to any observable residual dipolar couplings. However, in the pioneering study by Shortle and Ackerman RDCs were clearly observed. We repeated the study for urea-denatured protein at high temperature and also observed indisputably RDCs. This was very convincing to us but we could not possibly accept the proposed reason for the non-zero RDCs, namely that there would be some residual structure left in the protein that to our understanding was fully denatured. We proceeded to gain understanding via simulations and elementary experiments. In measurements we used simple homopolymers with only two labelled residues and we simulated the data to learn more about the origin of RDCs. We realized that RDCs depend on the position of the residue as well as on the length of the polypeptide. Investigations resulted in a theoretical model for RDCs from coil-like molecules. Later we extended the studies by molecular dynamics. Somewhat surprisingly the effects are small for non-structured molecules whereas the bias may be large for a small compact protein. All in all the work gave clear and unambiguous results on how to interpret RDCs as structural and dynamic parameters of weakly structured proteins.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Conformational study of valinomycin: a molecular dynamics approach

Valinomycin is a highly flexible cyclic dodecadepsipeptide that transports ions across membranes. Such a flexibility in the conformation is required for its biological function since it has to encounter a variety of environments and liganding state. Exploration of conformational space of this molecule is therefore important and is one of the objectives of the present study that has been carried out by means of high temperature Molecular Dynamics. Further, the stability of the known bracelet-like structure of the uncomplexed valinomycin and the inherent flexibility around this structure has been investigated. The uncomplexed form of valinomycin has been simulated at 75–100 K for 1 ns in order to elucidate the average conformational properties. An alanine-analog of valinomycin has been simulated under identical conditions in order to evaluate the effect of sidechain on the conformational properties, The studies confirm the effect of sidechain on conformational equilibrium.

NUPARM and NUCGEN: software for analysis and generation of sequence dependent nucleic acid structures

Software packages NUPARM and NUCGEN, are described, which can be used to understand sequence directed structural variations in nucleic acids, by analysis and generation of non-uniform structures. A set of local inter basepair parameters (viz. tilt, roll, twist, shift, slide and rise) have been defined, which use geometry and coordinates of two successive basepairs only and can be used to generate polymeric structures with varying geometries for each of the 16 possible dinucleotide steps. Intra basepair parameters, propeller, buckle, opening and the C6...C8 distance can also be varied, if required, while the sugar phosphate backbone atoms are fixed in some standard conformation ill each of the nucleotides. NUPARM can be used to analyse both DNA and RNA structures, with single as well as double stranded helices. The NUCGEN software generates double helical models with the backbone fixed in B-form DNA, but with appropriate modifications in the input data, it can also generate A-form DNA ar rd RNA duplex structures.

1-(2-Chloroacetyl)-3-methyl-2,6-diphenylpiperidin-4-one

The asymmetric unit of the title compound, C20H20ClNO2, contains two crystallographically independent molecules of similar geometry. The piperidine ring adopts a distorted boat conformation in both molecules, in which the N atom assumes an almost planar configuration.

Processing of SeMV polyproteins revisited

Processing of Sesbania mosaic virus (SeMV) polyprotein 2a and 2ab was reanalyzed in the view of the new genome organization of sobemoviruses. Polyprotein 2a when expressed in E coli, from the new cDNA clone, got cleaved at the earlier identified sites E325-T326, E402-T403 and E498-S499 to release protease, VPg, P10 and P8, respectively. Additionally, a novel cleavage was identified within the protease domain at position E132-S133, which was found to be essential for efficient polyprotein processing. Products, corresponding to cleavages identified in E. coli, were also detected in infected Sesbania leaves. Interestingly, though the sites are exactly the same in polyprotein 2ab, it got cleaved between Protease-VPg but not between VPg-RdRp. This indicates to a differential cleavage preference, governed probably by the conformation of 2ab. Also, the studies revealed that, in SeMV, processing is regulated by mode of cleavage and context of the cleavage site.

A Search for Energy Minimized Sequences of Proteins

In this paper, we present numerical evidence that supports the notion of minimization in the sequence space of proteins for a target conformation. We use the conformations of the real proteins in the Protein Data Bank (PDB) and present computationally efficient methods to identify the sequences with minimum energy. We use edge-weighted connectivity graph for ranking the residue sites with reduced amino acid alphabet and then use continuous optimization to obtain the energy-minimizing sequences. Our methods enable the computation of a lower bound as well as a tight upper bound for the energy of a given conformation. We validate our results by using three different inter-residue energy matrices for five proteins from protein data bank (PDB), and by comparing our energy-minimizing sequences with 80 million diverse sequences that are generated based on different considerations in each case. When we submitted some of our chosen energy-minimizing sequences to Basic Local Alignment Search Tool (BLAST), we obtained some sequences from non-redundant protein sequence database that are similar to ours with an E-value of the order of 10(-7). In summary, we conclude that proteins show a trend towards minimizing energy in the sequence space but do not seem to adopt the global energy-minimizing sequence. The reason for this could be either that the existing energy matrices are not able to accurately represent the inter-residue interactions in the context of the protein environment or that Nature does not push the optimization in the sequence space, once it is able to perform the function.

The crystal and molecular structure of benzyloxycarbonyl-a-amino-isobutyryl-L-prolyl-methylamide. The observation of an X2-Pro3 Type III b-turn

The crystal and molecular structure of N-benzyloxycarbonyl-a-aminoisobutyryl-L-prolyl methylamide, the amino terminal dipeptide fragment of alamethicin, has been determined using direct methods. The compound crystallizes in the orthorhombic system with the space group P212-21. Cell dimensions are a = 7.705 A, b = 11.365 A, and c = 21.904 A. The structure has been refined using conventional procedures to a final R factor of 0.054. The molecular structure possesses a 4 - 1 intramolecular N-H--0 hydrogen bond formed between the CO group of the urethane moiety and the NH group of the methylamide function. The peptide backbone adopts the type 111 P-turn conformation, with 42 = -51.0°, +* = -39.7",&j = -65.0', $3 = -25.4'. An unusual feature is the occurrence of the proline residue at position 3 of the P-turn. The observed structure supports the view that Aib residues initiate the formation of type 111 @-turn conformations. The pyrrolidine ring is puckered in Cy-exo fashion.

An incipient 310 helix in Piv-Pro-Pro-Ala-NHMe as a model for peptide folding

The molecular mechanism of helix nucleation in peptides and proteins is not yet understood and the question of whether sharp turns in the polypeptide backbone serve as nuclei for protein folding has evoked controversy1,2. A recent study of the conformation of a tetrapeptide containing the stereochemically constrained residue alpha-aminoisobutyric acid, both in solution and the solid state, yielded a structure consisting of two consecutive beta-turns, leading to an incipient 310 helical conformation3,4. This led us to speculate that specific tri- and tetra-peptide sequences may indeed provide a helical twist to the amino-terminal segment of helical regions in proteins and provide a nucleation site for further propagation. The transformation from a 310 helical structure to an alpha-helix should be facile and requires only small changes in the phi and psi conformational angles and a rearrangement of the hydrogen bonding pattern5. If such a mechanism is involved then it should be possible to isolate an incipient 310 helical conformation in a tripeptide amide or tetrapeptide sequence, based purely on the driving force derived from short-range interactions. We have synthesised and studied the model peptide pivaloyl-Pro-Pro-Ala-NHMe (compound I) and provide here spectroscopic evidence for a 310 helical conformation in compound I.

Aggregation of apolar peptides in organic solvents. Concentration dependence of 1H-nmr parameters for peptide NH groups in 310 helical decapeptide fragment of suzukacillin

Peptide NH chemical shifts and their temperature dependences have been monitored as a function of concentration for the decapeptide, Boc-Aib-Pro-Val-Aib-Val-Ala-Aib-Ala-Aib-Aib-OMe in CDCl3 (0.001-0.06M) and (CD3)2SO (0.001-0.03M). The chemical shifts and temperature coefficients for all nine NH groups show no significant concentration dependence in (CD3)2SO. Seven NH groups yield low values of temperature coefficients over the entire range, while one yields an intermediate value. In CDCl3, the Aib(1) NH group shows a large concentration dependence of both chemical shift and temperature coefficient, in contrast to the other eight NH groups. The data suggest that in (CD3)2SO, the peptide adopts a 310 helical conformation and is monomeric over the entire concentration range. In CDCl3, the 310 helical peptide associates at a concentration of 0.01M, with the Aib(1) NH involved in an intermolecular hydrogen bond. Association does not disrupt the intramolecular hydrogen-bonding pattern in the decapeptide.

Ethyl 2-acetyl-3-(4-bromoanilino)butanoate

The title compound, C14H18BrNO3, adopts an extended conformation, with all of the main-chain torsion angles associated with the ester and amino groups close to trans. In the crystal, inversion dimers linked by pairs of N-H center dot center dot center dot O hydrogen bonds are observed.

Diethyl 4-(2,4-dichlorophenyl)-2,6-dimethyl 1,4-dihydropyridine-3,5-dicarboxylate

In the title compound, C19H21Cl2NO4, the dihydropyridine ring adopts a flattened boat conformation. The dichlorophenyl ring is oriented almost perpendicular to the planar part of the dihydropyridine ring [dihedral angle = 89.1 (1)degrees]. An intramolecular C-H center dot center dot center dot O hydrogen bond is observed. In the crystal structure, molecules are linked into chains along the b axis by N-H center dot center dot center dot O hydrogen bonds.

2-[2-(Hydroxymethyl)phenyl]-1-(1-naphthyl)ethanol

The molecular conformation of the title compound, C19H18O2, is stabilized by an intramolecular O-H-O hydrogen bond. In addition, intermolecular O-H-O interactions link the molecules into zigzag chains running along the c axis.

2-[(E)-(4-hydroxy-3-methoxybenzylidene)amino]-N-(2-methylphenyl)-4,5,6,7-tetrahydro-1-benzothiophene-3-carboxamide

The title compound, C24H24N2O3S, exhibits antifungal and antibacterial properties. The compound crystallizes with two molecules in the asymmetric unit, with one molecule exhibiting 'orientational disorder' in the crystal structure with respect to the cyclohexene ring. The o-toluidine groups in both molecules are noncoplanar with the respective cyclohexene-fused thiophene ring. In both molecules, there is an intramolecular N-H...N hydrogen bond forming a pseudo-six-membered ring which locks the molecular conformation and eliminates conformational flexibility. The crystal structure is stabilized by O-H...O hydrogen bonds; both molecules in the asymmetric unit form independent chains, each such chain consisting of alternating 'ordered' and 'disordered' molecules in the crystal lattice.

2,6-Bis(3-methoxyphenyl)-3-methylpiperidin-4-one

In the molecule of the title compound, C20H23NO3, the bulky methoxyphenyl substituents at the equatorial 2,6-positions crowd the vicinity of the equatorial amino H atom and prevent it from forming intermolecular hydrogen bonds. The piperidine ring adopts a distorted chair conformation.