Main-chain Torsion Research Articles

Structures of Prostacyclin Synthase and Its Complexes with Substrate Analog and Inhibitor Reveal a Ligand-specific Heme Conformation Change

Prostacyclin synthase (PGIS) is a cytochrome P450 (P450) enzyme that catalyzes production of prostacyclin from prostaglandin H(2). PGIS is unusual in that it catalyzes an isomerization rather than a monooxygenation, which is typical of P450 enzymes. To understand the structural basis for prostacyclin biosynthesis in greater detail, we have determined the crystal structures of ligand-free, inhibitor (minoxidil)-bound and substrate analog U51605-bound PGIS. These structures demonstrate a stereo-specific substrate binding and suggest features of the enzyme that facilitate isomerization. Unlike most microsomal P450s, where large substrate-induced conformational changes take place at the distal side of the heme, conformational changes in PGIS are observed at the proximal side and in the heme itself. The conserved and extensive heme propionate-protein interactions seen in all other P450s, which are largely absent in the ligand-free PGIS, are recovered upon U51605 binding accompanied by water exclusion from the active site. In contrast, when minoxidil binds, the propionate-protein interactions are not recovered and water molecules are largely retained. These findings suggest that PGIS represents a divergent evolution of the P450 family, in which a heme barrier has evolved to ensure strict binding specificity for prostaglandin H(2), leading to a radical-mediated isomerization with high product fidelity. The U51605-bound structure also provides a view of the substrate entrance and product exit channels.

Crystal Structures and Site-directed Mutagenesis of a Mycothiol-dependent Enzyme Reveal a Novel Folding and Molecular Basis for Mycothiol-mediated Maleylpyruvate Isomerization

Mycothiol (MSH) is the major low molecular mass thiols in many Gram-positive bacteria such as Mycobacterium tuberculosis and Corynebacterium glutamicum. The physiological roles of MSH are believed to be equivalent to those of GSH in Gram-negative bacteria, but current knowledge of MSH is limited to detoxification of alkalating chemicals and protection from host cell defense/killing systems. Recently, an MSH-dependent maleylpyruvate isomerase (MDMPI) was discovered from C. glutamicum, and this isomerase represents one example of many putative MSH-dependent enzymes that take MSH as cofactor. In this report, fourteen mutants of MDMPI were generated. The wild type and mutant (H52A) MDMPIs were crystallized and their structures were solved at 1.75 and 2.05 A resolution, respectively. The crystal structures reveal that this enzyme contains a divalent metal-binding domain and a C-terminal domain possessing a novel folding pattern (alphabetaalphabetabetaalpha fold). The divalent metal-binding site is composed of residues His52, Glu144, and His148 and is located at the bottom of a surface pocket. Combining the structural and site-directed mutagenesis studies, it is proposed that this surface pocket including the metal ion and MSH moiety formed the putative catalytic center.

Is glycine a surrogate for a d‐amino acid in the collagen triple helix?

Collagen is the most abundant protein in animals. Every third residue in a collagen strand is a glycine with phi, psi = -70 degrees, 175 degrees. A recent computational study suggested that replacing these glycine residues with D-alanine or D-serine would stabilize the collagen triple helix. This hypothesis is of substantial importance, as the glycine residues in collagen constitute nearly 10% of the amino acid residues in humans. To test this hypothesis, we synthesized a series of collagen mimic peptides that contain one or more D-alanine or D-serine residues replacing the canonical glycine residues. Circular dichroism spectroscopy and thermal denaturation experiments indicated clearly that the substitution of glycine with D-alanine or D-serine greatly disfavors the formation of a triple helix. Host-guest studies also revealed that replacing a single glycine residue with D-alanine is more destabilizing than is its replacement with L-alanine, a substitution that results from a common mutation in patients with collagen-related diseases. These data indicate that the glycine residues in collagen are not a surrogate for a D-amino acid and support the notion that the main-chain torsion angles of a glycine residue in the native structure (especially, phi > 0 degrees ) are critical determinants for its beneficial substitution with a D-amino acid in a protein.

Exploring the binding site of the human muscarinic M3 receptor: Homology modeling and docking study

AbstractThe human muscarinic M3 receptor (hM3) and its interactions with selective agonists and antagonists were investigated by means of combined homology and docking approach. Also, two pharmacophoric models for the hM3 agonist and antagonist binding sites were proposed. The three‐dimensional (3D) structure of hM3 receptor was modeled based on the high‐resolution X‐ray structure of bovine rhodopsin from the Protein Data Bank (PDB). To validate the reliability of the model obtained, the main chain torsion angles phi (Ψ) and psi (Φ) were examined in a Ramachandran plot, and all omega angles were measured for peptidic bond planarity. The characteristics of the active site, the position, and the orientation of ligands in situ, as well as the binding modes of the representative agonists and antagonists, were analyzed by applying a molecular docking technique using the AutoDock 3.0.5 program. Specific interactions responsible for recognition of the hM3 receptor, like ionic bond formed between protonated amine of the ligands and the Asp3.6 side chain were identified. Structure–reactivity relationships have been explained by analyzing the 3D structure of the hM3 model and the ligand conformations resulted from molecular docking simulation. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007

The Ig Doublet Z1Z2: A Model System for the Hybrid Analysis of Conformational Dynamics in Ig Tandems from Titin

Titin is a gigantic elastic filament that determines sarcomere ultrastructure and stretch response in vertebrate muscle. It folds into numerous Ig and FnIII domains connected in tandem. Data on interdomain arrangements and dynamics are essential for understanding the function of this filament. Here, we report a mechanistic analysis of the conformational dynamics of two Ig domains from the N terminus of titin, Z1Z2, by using X-ray crystallography, SAXS, NMR relaxation data, and residual dipolar couplings in combination. Z1Z2 preferentially adopts semiextended conformations in solution, with close-hinge arrangements representing low-probability states. Although interdomain contacts are not observed, the linker appears to acquire moderate rigidity via small, local hydrophobic interactions. Thus, Z1Z2 constitutes an adaptable modular system with restricted dynamics. We speculate that its preexistent conformation contributes to the selective recruitment of the binding partner telethonin onto the repetitive surface of the filament. The structural interconversion of four Z1Z2 conformers is analyzed.

Crystal structure of phosphoribosylformylglycinamidine synthase II (smPurL) from Thermotoga maritima at 2.15 Å resolution

Crystal structure of phosphoribosylformylglycinamidine synthase II (smPurL) from <i>Thermotoga maritima</i> at 2.15 Å resolution

Characterization of β‐turn and Asx‐turns mimicry in a model peptide: Stabilization via CH · · · O interaction

The chemical synthesis and single-crystal X-ray diffraction analysis of a model peptide, Boc-Thr-Thr-NH2 (1) comprised of proteinogenic residues bearing an amphiphilic Cbeta -stereogenic center, has been described. Interestingly, the analysis of its molecular structure revealed the existence of a distinct conformation that mimics a typical beta-turn and Asx-turns, i.e., the two Thr residues occupy the left- and right-corner positions. The main-chain torsion angles of the N- and C-terminal residues i.e., semiextended: phi = -68.9 degrees , psi = 128.6 degrees ; semifolded: phi = -138.1 degrees , psi = 2.5 degrees conformations, respectively, in conjunction with a gauche- disposition of the obligatory C-terminus Thr CgammaH3 group, characterize the occurrence of the newly described beta-turn- and Asx-turns-like topology. The preferred molecular structure is suggested to be stabilized by an effective nonconventional main-chain to side-chain Ci=O . . . H--Cgamma(i+2)-type intraturn hydrogen bond. Noteworthy, the observed topology of the resulting 10-membered hydrogen-bonded ring is essentially similar to the one perceived for a classical beta-turn and the Asx-turns, stabilized by a conventional intraturn hydrogen bond. Considering the signs as well as magnitudes of the backbone torsion angles and the orientation of the central peptide bond, the overall mimicked topology resembles the type II beta-turn or type II Asx-turns. An analysis of Xaa-Thr sequences in high-resolution X-ray elucidated protein structures revealed the novel topology prevalence in functional proteins (unpublished). In view of indubitable structural as well as functional importance of nonconventional interactions in bioorganic and biomacromolecules, we intend to highlight the participation of Thr CgammaH in the creation of a short-range C=O . . . H--Cgamma -type interaction in peptides and proteins.

Strike a balance: Optimization of backbone torsion parameters of AMBER polarizable force field for simulations of proteins and peptides

Based on the AMBER polarizable model (ff02), we have re-optimized the parameters related to the main-chain (Phi, Psi) torsion angles by fitting to the Boltzmann-weighted average quantum mechanical (QM) energies of the important regions (i.e., beta, P(II), alpha(R), and alpha(L) regions). Following the naming convention of the AMBER force field series, this release will be called ff02pol.rl The force field has been assessed both by energetic comparison against the QM data and by the replica exchange molecular dynamics simulations of short alanine peptides in water. For Ace-Ala-Nme, the simulated populations in the beta, P(II) and alpha(R) regions were approximately 30, 43, and 26%, respectively. For Ace-(Ala)(7)-Nme, the populations in these three regions were approximately 24, 49, and 26%. Both were in qualitative agreement with the NMR and CD experimental conclusions. In comparison with the previous force field, ff02pol.rl demonstrated good balance among these three important regions. The optimized torsion parameters, together with those in ff02, allow us to carry out simulations on proteins and peptides with the consideration of polarization.

Crystal structure of alanyl‐tRNA synthetase editing‐domain homolog (PH0574) from a hyperthermophile, Pyrococcus horikoshii OT3 at 1.45 Å resolution

Crystal structure of alanyl‐tRNA synthetase editing‐domain homolog (PH0574) from a hyperthermophile, <i>Pyrococcus horikoshii</i> OT3 at 1.45 Å resolution

A Complete Set of NMR Chemical Shifts and Spin−Spin Coupling Constants forl-Alanyl-l-alanine Zwitterion and Analysis of Its Conformational Behavior

With the aid of labeling with stable isotopes ((15)N and (13)C) a complete set of chemical shifts and indirect spin-spin coupling constants was obtained for the zwitterionic form of L-alanyl-L-alanine in aqueous solution. Different sensitivities of the NMR parameters to the molecular geometry were discussed on the basis of comparison with ab initio (DFT) calculated values. An adiabatic two-dimensional vibrational wave function was constructed and used for determination of the main chain torsion angle dispersions and conformational averaging of the NMR shifts and coupling constants. The quantum description of the conformational dynamics based on the density functional theory and a polarizable continuum solvent model agrees reasonably with classical molecular dynamics simulations using explicit solvent. The results consistently evidence the presence of a single form in the aqueous solution with equilibrium main chain torsion angle values (psi = 147 degrees, varphi = -153 degrees), close to that one found previously in an X-ray study. Under normal temperature the torsion angles can vary by about 10 degrees around their equilibrium values, which leads, however, to minor corrections of the NMR parameters only. The main chain heavy atom chemical shifts and spin-spin coupling constants involving the alpha-carbon and hydrogen atoms appear to be most useful for the peptide structural predictions.

Mapping and Consensus Sequence Identification for Multiple Vinculin Binding Sites within the Talin Rod

The interaction between the cytoskeletal proteins talin and vinculin plays a key role in integrin-mediated cell adhesion and migration. Three vinculin binding sites (VBS1-3) have previously been identified in the talin rod using a yeast two-hybrid assay. To extend these studies, we spot-synthesized a series of peptides spanning all the alpha-helical regions predicted for the talin rod and identified eight additional VBSs, two of which overlap key functional regions of the rod, including the integrin binding site and C-terminal actin binding site. The talin VBS alpha-helices bind to a hydrophobic cleft in the N-terminal vinculin Vd1 domain. We have defined the specificity of this interaction by spot-synthesizing a series of 25-mer talin VBS1 peptides containing substitutions with all the commonly occurring amino acids. The consensus for recognition is LXXAAXXVAXX- VXXLIXXA with distinct classes of hydrophobic side chains at positions 1, 4, 5, 8, 9, 12, 15, and 16 required for vinculin binding. Positions 1, 8, 12, 15, and 16 require an aliphatic residue and will not tolerate alanine, whereas positions 4, 5, and 9 are less restrictive. These preferences are common to all 11 VBS sequences with a minor variation occurring in one case. A crystal structure of this variant VBS peptide in complex with the vinculin Vd1 domain reveals a subtly different mode of vinculin binding.

AMBER-based hybrid force field for conformational sampling of polypeptides

AMBER-based hybrid force field was developed, with mixing the mainchain torsion energies of AMBER parm94 and parm96 force fields, with introducing a weighting factor to control the mixing ratio. First, dependence of an adiabatic Φ– Ψ potential–energy map of alanine dipeptide on the weighting factor was investigated, and the map was compared with an ab initio Φ– Ψ map. Second, we did enhanced conformational sampling of two six-residue peptides known to form a helical or turn conformation in solution. The hybrid force field reproduced the ab initio energy better than either parm94 or parm96, and exhibited the secondary-structure preference depending on the amino acid sequence.

Monte Carlo Backbone Sampling for Nucleic Acids Using Concerted Rotations Including Variable Bond Angles

An efficient concerted rotation algorithm for use in Monte Carlo statistical mechanics simulations of nucleic acids is reported. The corresponding algorithm “concerted rotations with flexible bond angles” (CRA) for sampling polypeptides was found to be superior to local moves that included only flexible dihedral angles by allowing exploration of a larger conformational space as well as facilitating backbone transitions. The performance of the present CRA algorithm for polynucleotides is compared to two alternatives, a simple update of main-chain torsion angles and a previously reported, concerted rotation algorithm with fixed bond angles and a mix of flexible and rigid main-chain dihedral angles. The test system is a 12 base-pair duplex B-form DNA helix, and the performance comparisons are made for the system both in a vacuum and with continuum GB/SA solvation. The results demonstrate the superior efficiency of the CRA method over the alternatives.

The crystal structure of placental growth factor in complex with domain 2 of vascular endothelial growth factor receptor-1.

Placental growth factor (PlGF) is a member of the vascular endothelial growth factor (VEGF) family and plays an important role in pathological angiogenic events. PlGF exerts its biological activities through binding to VEGFR1, a receptor tyrosine kinase that consists of seven immunoglobulin-like domains in its extracellular portion. Here we report the crystal structure of PlGF bound to the second immunoglobulin-like domain of VEGFR1 at 2.5 A resolution and compare the complex to the closely related structure of VEGF bound to the same receptor domain. The two growth factors, PlGF and VEGF, share a sequence identity of approximately 50%. Despite this moderate sequence conservation, they bind to the same binding interface of VEGFR1 in a very similar fashion, suggesting that both growth factors could induce very similar if not identical signaling events.

A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations.

Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third-generation point-charge all-atom force field for proteins. Following the earlier approach of Cornell et al., the charge set was obtained by fitting to the electrostatic potentials of dipeptides calculated using B3LYP/cc-pVTZ//HF/6-31G** quantum mechanical methods. The main-chain torsion parameters were obtained by fitting to the energy profiles of Ace-Ala-Nme and Ace-Gly-Nme di-peptides calculated using MP2/cc-pVTZ//HF/6-31G** quantum mechanical methods. All other parameters were taken from the existing AMBER data base. The major departure from previous force fields is that all quantum mechanical calculations were done in the condensed phase with continuum solvent models and an effective dielectric constant of epsilon = 4. We anticipate that this force field parameter set will address certain critical short comings of previous force fields in condensed-phase simulations of proteins. Initial tests on peptides demonstrated a high-degree of similarity between the calculated and the statistically measured Ramanchandran maps for both Ace-Gly-Nme and Ace-Ala-Nme di-peptides. Some highlights of our results include (1) well-preserved balance between the extended and helical region distributions, and (2) favorable type-II poly-proline helical region in agreement with recent experiments. Backward compatibility between the new and Cornell et al. charge sets, as judged by overall agreement between dipole moments, allows a smooth transition to the new force field in the area of ligand-binding calculations. Test simulations on a large set of proteins are also discussed.

Ramachandran plot on the web

A graphics package has been developed to display the main chain torsion angles phi, psi (phi, Psi); (Ramachandran angles) in a protein of known structure. In addition, the package calculates the Ramachandran angles at the central residue in the stretch of three amino acids having specified the flanking residue types. The package displays the Ramachandran angles along with a detailed analysis output. This software is incorporated with all the protein structures available in the Protein Databank.

Assessing equilibration and convergence in biomolecular simulations

If molecular dynamics simulations are used to characterize the folding of peptides or proteins, a wide range of conformational states needs to be sampled. This study reports an analysis of peptide simulations to identify the best methods for assessing equilibration and sampling in these systems where there is significant conformational disorder. Four trajectories of a beta peptide in methanol and four trajectories of an alpha peptide in water, each of 5 ns in length, have been studied. Comparisons have also been made with two 50-ns trajectories of the beta peptide in methanol. The convergence rates of quantities that probe both the extent of conformational sampling and the local dynamical properties have been characterized. These include the numbers of hydrogen bonds populated, clusters identified, and main chain torsion angle transitions in the trajectories. The relative equilibrium rates of different quantities are found to vary significantly between the two systems studied reflecting both the differences in peptide primary structure and the different solvents used. A cluster analysis of the simulation trajectories is identified as a very effective method for judging the convergence of the simulations. This is particularly the case if the analysis includes a comparison of multiple trajectories calculated for the same system from different starting structures.

Prediction of nonregular secondary structures of proteins based on wavelet analysis

The secondary structures of proteins fall into two classes: regular structure and nonregular structure. Helices and sheets are termed “regular” structures because their residues have repeating main-chain torsion angles, and their backbone N-H and C-O groups are arranged in a periodic pattern of hydrogen bonding. In contrast, the remaining structures with nonrepeating backbone torsion angles are called nonregular secondary structures. In this note, we performed an extensive sequence analysis of nonregular secondary structures and showed that these nonregular parts could be effectively predicted by continuous wavelet transform.

An Unusual Chain Conformation: s(6*4/3) Helices in Form II of Isotactic Poly(vinylcyclopentane)

The crystalline form II of isotactic poly(vinylcyclopentane) was studied on the basis of fiber X-ray diffraction and of potential energy calculations of both isolated and packed chains. The hypothesis of a s(6*4/3) conformation of the chain has been investigated in detail. According to our interpretation, s(6*4/3) helices, which represent an anomalous case among isotactic poly(α-olefins), result from the superposition of a basically s(2*4/1) conformation of the backbone chain and of side groups following one another with a regular sequence of three conformations, mainly different for their orientation. Slight deformations of the main chain torsion angles from s(2*4/1) symmetry are found to be relevant to improve the data fitting. The resulting structure is indicative of a quite singular polymorphism, mainly ascribable to the different conformations accessible to the side groups.

Phosphorylation disrupts the central helix in Op18/stathmin and suppresses binding to tubulin.

Protein phosphorylation represents a ubiquitous control mechanism in living cells. The structural prerequisites and consequences of this important post-translational modification, however, are poorly understood. Oncoprotein 18/stathmin (Op18) is a globally disordered phosphoprotein that is involved in the regulation of the microtubule (MT) filament system. Here we document that phosphorylation of Ser63, which is located within a helix initiation site in Op18, disrupts the transiently formed amphipathic helix. The phosphoryl group reduces tubulin binding 10-fold and suppresses the MT polymerization inhibition activity of Op18's C-terminal domain. Op18 represents an example where phosphorylation occurs within a regular secondary structural element. Together, our findings have implications for the prediction of phosphorylation sites and give insights into the molecular behavior of a globally disordered protein.