Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

Crystal structure and conformational analysis of angiotensinogen fragments

The tripeptide acetyl-L-prolyl-L-phenylalanyl-L-histidine crystallizes in the orthorhombic space group P2(1)2(1)2(1) with eight molecules in a unit cell of dimensions a = 9.028(2), b = 140.54(6) and c = 42.41(1)A. The structure has been solved by direct methods and refined to an R value of 0.056 for 2904 observed reflections. The molecule exists as a zwitterion with terminal (His)CO2- and (imidazole)H+ as charged groups. The two peptide molecules in the structure adopt a type I beta-turn with Pro and Phe as the corner residues. The main conformational difference between the two crystallographically independent molecules is seen to be in the histidine side-chain orientations. The molecules arrange themselves in sheets perpendicular to the c axis. All hydrophobic side chains lie on one side of the sheets thus generated, whereas the hydrophilic groups are located on the other side. An interesting feature of the crystal structure is the existence of a water layer between adjacent peptide sheets. The conformational study of the isolated Ac-His-Pro-Phe-His-MA using energy calculations gives a rather limited number of stable conformers. The most stable corresponds to a type I beta-turn stabilized through two hydrogen bonds, followed by a less stable type II beta-turn (delta E = 2.0 kcal) and a partly helical structure (delta E = 2.6 kcal).     

Crystal structure and conformational analysis of ampullosporin A.

Ampullosporin A is a 15-mer peptaibol type polypeptide that induces pigment formation by the fungus Phoma destructiva, forms voltage-dependent ion channels in membranes and exhibits hypothermic effects in mice. The structure of ampullosporin A has been determined by x-ray crystallography. This is the first three-dimensional (3D) structure of the peptaibol subfamily SF6. From the N-terminus to residue 13 the molecule adopts an approximate right-handed alpha-helical geometry, whereas a less regular structure pattern with beta-turn characteristics is found in the C-terminus. Even though ampullosporin A does not contain a single proline or hydroxyproline it is significantly bent. It belongs to both the shortest and the most strongly bent peptaibol 3D structures. The straight structure part encompasses residues Ac-Trp(1)-Aib(10) and is thus less extended than the alpha-helical subunit. The 3D structure of ampullosporin A is discussed in relation to other experimentally determined peptaibol structures and in the context of its channel-forming properties. As a part of this comparison a novel bending analysis based on a 3D curvilinear axis describing the global structural characteristics has been proposed and applied to all 3D peptaibol structures. A sampling of 2500 conformations using different molecular dynamics protocols yields, for the complete ampullosporin A structure, an alpha-helix as the preferred conformation in vacuo with almost no bend. This indicates that solvent or crystal effects may be important for the experimentally observed peptide backbone bending characteristics of ampullosporin A.     

Modeling ganglioside headgroups by conformational analysis and molecular dynamics

The conformations and dynamics of gangliosides GM1, GM2, 6-GM2 and GM4 have been studied by computational means, and the results compared to NMR data. Unconstrained conformational searches were run using the AMBER* force field augmented by MNDO derived parameters for the Neu5Ac anomeric torsion, the GB/SA water solvation model, and the MC/EM alogorithm; extended (10–12[emsp4 ]ns) dynamic simulations in GB/SA water were performed with the MC/SD protocol, and the stored structures were minimized. The overall mobility of the Neu5Ac2,3Gal linkage and the position of its minimum energy conformation have been shown to depend mainly on the presence or the absence of a GalNAc residue at the adjacent position. The best quantitative agreement with the available NOE data was achieved after minimization of the structures stored during the MC/SD dynamic runs. The latter protocol appears to reproduce satisfactorily the available experimental data, and can be used with confidence to build three-dimensional models of ganglioside headgroups.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.     

LeuT conformational sampling utilizing accelerated molecular dynamics and principal component analysis

Monoamine transporters (MATs) function by coupling ion gradients to the transport of dopamine, norepinephrine, or serotonin. Despite their importance in regulating neurotransmission, the exact conformational mechanism by which MATs function remains elusive. To this end, we have performed seven 250 ns accelerated molecular dynamics simulations of the leucine transporter, a model for neurotransmitter MATs. By varying the presence of binding-pocket leucine substrate and sodium ions, we have sampled plausible conformational states representative of the substrate transport cycle. The resulting trajectories were analyzed using principal component analysis of transmembrane helices 1b and 6a. This analysis revealed seven unique structures: two of the obtained conformations are similar to the currently published crystallographic structures, one conformation is similar to a proposed open inward structure, and four conformations represent novel structures of potential importance to the transport cycle. Further analysis reveals that the presence of binding-pocket sodium ions is necessary to stabilize the locked-occluded and open-inward conformations.     

Crystal structure and molecular dynamics simulation of ubiquitin-like domain of murine parkin

Parkin is the gene product identified as the major cause of autosomal recessive juvenile Parkinsonism (AR-JP). Parkin, a ubiquitin ligase E3, contains a unique ubiquitin-like domain in its N-terminus designated Uld which is assumed to be a interaction domain with the Rpn 10 subunit of 26S proteasome. To elucidate the structural and functional role of Uld in parkin at the atomic level, the X-ray crystal structure of murine Uld was determined and a molecular dynamics simulation of wild Uld and its five mutants (K27N, R33Q, R42P, K48A and V56E) identified from AR-JP patients was performed. Murine Uld consists of two alpha helices [Ile23-Arg33 (alpha1) and Val56-Gln57 (alpha2)] and five beta strands [Met1-Phe7 (beta1), Tyr11-Asp18 (beta2), Leu41-Phe45 (beta3), Lys48-Pro51 (beta4) and Ser65-Arg72 (beta5)] and its overall structure is essentially the same as that of human ubiquitin with a 1.22 A rmsd for the backbone atoms of residues 1-76; however, the sequential identity and similarity between both molecules are 32% and 63%, respectively. This close resemblance is due to the core structure built by same hydrogen bond formations between and within the backbone chains of alpha1 and beta1/2/5 secondary structure elements and by nearly the same hydrophobic interactions formed between the nonpolar amino acids of their secondary structures. The side chain NetaH of Lys27 on the alpha1 helix was crucial to the stabilization of the spatial orientations of beta3 and beta4 strands, possible binding region with Rpn 10 subunit, through three hydrogen bonds. The MD simulations showed the K27N and R33Q mutations increase the structural fluctuation of these beta strands including the alpha1 helix. Reversely, the V56E mutant restricted the spatial flexibility at the periphery of the short alpha2 helix by the interactions between the polar atoms of Glu56 and Ser19 residues. However, a large fluctuation of beta4 strand with respect to beta5 strand was induced in the R42P mutant, because of the impossibility of forming paired hydrogen bonds of Pro for Arg42 in wild Uld. The X-ray structure showed that the side chains of Asp39, Gln40 and Arg42 at the N-terminal periphery of beta3 strand protrude from the molecular surface of Uld and participate in hydrogen bonds with the polar residues of neighboring Ulds. Thus, the MD simulation suggests that the mutation substitution of Pro for Arg42 not only causes the large fluctuation of beta3 strand in the Uld but also leads to the loss of the ability of Uld to trap the Rpn 10 subunit. In contrast, the MD simulation of K48A mutant showed little influence on the beta3-beta4 loop structure, but a large fluctuation of Lys48 side chain, suggesting the importance of flexibility of this side chain for the interaction with the Rpn 10 subunit. The present results would be important in elucidating the impaired proteasomal binding mechanism of parkin in AR-JP.     

The dynamics,structure, and conformational free energy of proline-containing antifreeze glycoprotein

Recent NMR studies of the solution structure of the 14-amino acid antifreeze glycoprotein AFGP-8 have concluded that the molecule lacks long-range order. The implication that an apparently unstructured molecule can still have a very precise function as a freezing inhibitor seems startling at first consideration. To gain insight into the nature of conformations and motions in AFGP-8, we have undertaken molecular dynamics simulations augmented with free energy calculations using a continuum solvation model. Starting from 10 different NMR structures, 20 ns of dynamics of AFGP were explored. The dynamics show that AFGP structure is composed of four segments, joined by very flexible pivots positioned at alanine 5, 8, and 11. The dynamics also show that the presence of prolines in this small AFGP structure facilitates the adoption of the poly-proline II structure as its overall conformation, although AFGP does adopt other conformations during the course of dynamics as well. The free energies calculated using a continuum solvation model show that the lowest free energy conformations, while being energetically equal, are drastically different in conformations. In other words, this AFGP molecule has many structurally distinct and energetically equal minima in its energy landscape. In addition, conformational, energetic, and hydrogen bond analyses suggest that the intramolecular hydrogen bonds between the N-acetyl group and the protein backbone are an important integral part of the overall stability of the AFGP molecule. The relevance of these findings to the mechanism of freezing inhibition is discussed.     

Application of restrained minimization, simulated annealing and molecular dynamics simulations for the conformational analysis of oligosaccharides.

The purpose of the present study was to determine the confidence with which the small number of 1H NMR nuclear Overhauser effect (NOE) distance constraints measurable across glycosidic linkages in oligosaccharides could be used for solution conformational analysis. This was assessed by use of these constraints in restrained molecular mechanical minimization of the tetrasaccharide Gal beta 1----4(Fuc alpha 1----3)Glc-NAc beta 1----3Gal, a model compound of the Lewis-X antigenic determinant. This presents a particularly severe test case in view of extreme resonance overlap and a dearth of inter-residue distance constraints. It is concluded that these constraints, when used in conventional restrained minimization, result in the generation of 'virtual conformations' and local minima about glycosidic linkages. However, these restraints are nevertheless found to be useful in the initial stages of a conformational analysis strategy involving restrained minimization combined with dynamical simulated annealing to define more accurately the global minimum energy configuration, together with molecular dynamics simulation to explore conformational mobility about this minimum. Theoretical ROE values calculated over the time course of the MD simulation, using a formalism appropriate for the time scale of the internal motion, are compared with those obtained experimentally in the oligosaccharide.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Murine interleukin-3: structure, dynamics, and conformational heterogeneity in solution

Interleukin-3 (IL-3), a cytokine produced primarily by activated T-cells during immune responses, is a crucial regulator of allergic inflammation. The three-dimensional structure of murine IL-3 (mIL-3) has remained elusive owing to its poor solubility and strong tendency toward aggregation under solution conditions typically used for structural studies. Here we describe the solution properties and structure of mIL-3 determined by NMR using an engineered construct of mIL-3 (mIL-3(33-156)). mIL-3 adopts a four-helical bundle fold, typical of proteins belonging to the short-chain cytokine family, and features a core of highly conserved hydrophobic residues. While significant line broadening and peak disappearance were observed in NMR spectra at higher temperatures, there was no evidence for temperature-dependent changes of the oligomeric state of mIL-3(33-156). Further analysis of the temperature dependence of amide (1)H chemical shifts and backbone (15)N relaxation parameters, including (15)N relaxation dispersion, revealed the presence of significant conformational exchange and local conformational heterogeneity. Residues recently shown by mutagenesis to play key roles in β(IL-3) receptor recognition and activation, which are located within the α(A) and α(C) helices and aligned on one face of the mIL-3(33-156) structure, are relatively rigid. In contrast, pronounced conformational heterogeneity was observed for a cluster of residues located on the opposite side of mIL-3, which corresponds spatially to sites in the related cytokines human IL-3, IL-5, and GM-CSF that are known to mediate interactions with their respective α-receptor subunits. Such conformational heterogeneity may facilitate the interaction of mIL-3 with each of two naturally occurring mIL-3Rα isoforms, leading to structurally distinct high-affinity complexes.     

Vasopressin conformational fluctuations: a molecular dynamics study

Molecular dynamics simulations have been used to study the conformational fluctuations of the oligopeptide hormone vasopressin. Starting coordinates for these simulations were built upon the crystal structure of pressinoic acid, the cyclic ring moiety of vasopressin, recently determined by x-ray diffraction. Coordinates for the additional tripeptide “tail” of vasopressin were selected by arbitrary positioning of this segment using interactive computer graphics. Two such starting configurations were minimized to relax strains, and long dynamics simulations (20 and 40 ps) in vacuo were then conducted following extensive heating and equilibration sequences (36 ps). In these studies, vasopressin was found to undergo few substantial conformational changes at 300 K on the time scale simulated, in contrast to the results of a shorter previous simulation, but comparable structural transitions were observed during the equilibration periods. The pressinoic acid structure was found to be a reasonably stable possible conformation for vasopressin in vacuum on this time scale.     

Crystal structure of thrombin-ecotin reveals conformational changes and extended interactions

The protease inhibitor ecotin fails to inhibit thrombin despite its broad specificity against serine proteases. A point mutation (M84R) in ecotin results in a 1.5 nM affinity for thrombin, 10(4) times stronger than that of wild-type ecotin. The crystal structure of bovine thrombin is determined in complex with ecotin M84R mutant at 2.5 A resolution. Surface loops surrounding the active site cleft of thrombin have undergone significant structural changes to permit inhibitor binding. Particularly, the insertion loops at residues 60 and 148 in thrombin, which likely mediate the interactions with macromolecules, are displaced when the complex forms. Thrombin and ecotin M84R interact in two distinct surfaces. The loop at residue 99 and the C-terminus of thrombin contact ecotin through mixed polar and nonpolar interactions. The active site of thrombin is filled with eight consecutive amino acids of ecotin and demonstrates thrombin's preference for specific features that are compatible with the thrombin cleavage site: negatively charged-Pro-Val-X-Pro-Arg-hydrophobic-positively charged (P1 Arg is in bold letters). The preference for a Val at P4 is clearly defined. The insertion at residue 60 may further affect substrate binding by moving its adjacent loops that are part of the substrate recognition sites.     

Degradation of (+)-catechin by Acinetobacter calcoaceticus MTC 127

Acinetobacter calcoaceticus MTC 127 was able to grow on catechin and protocatechuic acid (PCA) as sole carbon source. Cells induced with catechin oxidized catechin and PCA at rates higher than cells of uninduced cultures. Two aromatic compounds, PCA and phloroglucinol carboxylic acid (PGCA) were isolated from culture filtrate of cells grown in catechin and characterized by infrared spectrometry and high performance thin-layer chromatography. Moreover, A. calcoaceticus MTC 127 produced high levels of PCA compared to PGCA in the degradation of catechin. Based upon these results, a pathway for the degradation of (+)-catechin in A. calcoaceticus MTC 127 is proposed. Enzymes extracted from catechin-induced culture showed catechin oxygenase (cox) and protocatechuate 3,4-dioxygenase (pcd) activities. Catechin oxygenase was purified by column chromatography and SDS-PAGE analysis showed a single band with an apparent molecular weight of 47 kDa.     

Biotransformation of (+)-catechin into taxifolin by a two-step oxidation: primary stage of (+)-catechin metabolism by a novel (+)-catechin-degrading bacteria, Burkholderia sp. KTC-1, isolated from tropical peat

Peat contains various persistent compounds derived from plant materials. We isolated a novel (+)-catechin-degrading bacterium, Burkholderia sp. KTC-1 (KTC-1), as an example of a bacterium capable of degrading persistent aromatic compounds present in tropical peat. This bacterium was isolated by an enrichment technique and grew on (+)-catechin as the sole carbon source under acidic conditions. The reaction of a crude enzyme extract and a structural study of its products showed that (+)-catechin is biotransformed into taxifolin during the preliminary stages of its metabolism by KTC-1. HPLC analysis showed that this transformation occurs in two oxidation steps: 4-hydroxylation and dehydrogenation. Furthermore, both (+)-catechin 4-hydroxylanase and leucocyanidin 4-dehydrogenase were localized in the cytosol of KTC-1. This is the first report on biotransformation of (+)-catechin into taxifolin via leucocyanidin by an aerobic bacterium. We suggest that tropical peat could become a unique resource for microorganisms that degrade natural aromatic compounds.     

Synthesis and conformational analysis of N-glycopeptides. II. CD,molecular dynamics,and nmr spectroscopic studies on linear N-glycopeptides

The comprehensive structural analysis reported herein of eight N-glycopeptides, in three different solvents, is based on quantitative CD experiments, homonuclear nuclear Overhauser effect measurements, and molecular dynamics (MD) calculations. Although several orientations of the two amide planes attached to the carbohydrate pyranose ring are possible, according to NOE, CD data, and MD simulations, of all of the glycopeptide models, regardless of the type of the carrier peptide, only one dominant conformer population was found. This conformer is characterized by a nearly trans orientation of the CH and NH hydrogens of both acetamido groups. This finding is in perfect agreement with x-ray crystallographic data on the solid state conformation of the 1-N-acetyl- and 1-N-(β-aspartoyl)-2-acetamido-2-deoxy-β-D -glucopyranosylamine. The precise identification of this dominant conformer of N-glycopeptides in solution was the major question addressed herein by the structural analyses. A “CD additivity” experiment was carried out using an equimolar solution of Boc-Pro-Asp-NHCH₃ and l-N-acetyl-3,4,6-tri-O-acetyl-2-acetamido-2-deoxy-β-D -gluco-pyranosylamine at ambient temperature in acetonitrile. The CD spectrum obtained from the equimolar solution of the above two molecules (the “spectroscopic sum”) was identical with the CD curve obtained from the algebraic summation of the individually recorded CD spectra of the peptide and the carbohydrate moiety (“mathematical sum”). The global picture of the CD spectral analyses of the eight parent peptides with the eight N-glycopeptides revealed that in trifluoroethanol and acetonitrile, the side-chain modification of the Asn models (natural N-glycopeptide analogues) by N-glycosylation has a significant effect on the conformation of the carrier peptide, resulting in a decrease in the original type I β-turn content. Simultaneously, the type II β-turn conformational percentage increased to ≈ 20%. Such a conformational ratio change seems to be larger than the expected errors arising from the CD analyses, and agrees with the results of MD calculations. N-glycosylation of Asn residues causes perturbations, not only through the covalent bond, but also through specific hydrogen bonds between the backbone and side chain atoms. CD spectroscopy, augmented by efficient CD curve deconvolution techniques, has proved to be a useful tool for studying multicomponent conformer mixtures of small linear peptides in solution and changes of conformational equilibria caused by N-glycosylation. © 1993 John Wiley & Sons, Inc.     

Crystal and molecular structure of cyclohepta-amylose dodecahydrate

β-Cyclodextrin (cyclohepta-amylose, β-CD) is a torus-shaped, cyclic heptasaccharide consisting of (1→4)-linked α-d-glucopyranosyl residues. It is able to form inclusion complexes with small molecules in aqueous solution because of its annular aperture (width, 6.2 Å). β-Cyclodextrin dodecahydrate, the “empty” β-CD, crystallises from water in space group P2₁, with cell constants a = 21.29(2), b = 10.33(1), c = 15.10(2) Å, and β = 112.3(5)°. A total of 5189 X-ray counter-data were collected on a four-circle diffractometer. The crystal structure was solved on the basis of the highly isomorphous β-CD · 2HI · 8H₂O adduct, and the atomic parameters were refined by the full matrix, least-squares method to R = 7.3% for all data. The crystal structure belongs to the cage type. The β-CD macrocycle exists in an open, circular conformation stabilised by intramolecular hydrogen-bonds between HO-2 and HO-3 of adjacent glucosyl residues; four of the seven HO-6 groups are in the favoured (?)gauche orientation with respect to O-5, two are in the (+)gauche orientation, and one is disordered over these two orientations. The 6.5 water molecules within the cavity are distributed over 8 sites and display extensive thermal motion which is probably correlated with statistical disorder.     

Crystal and molecular structure of V-anhydrous amylose

The conformation and crystalline packing of V-anhydrous amylose has been investigated by a combination of linked atom model building and X-ray diffraction analysis. The unit cell, the P2₁2₁2₁ space group, the left-handed sixfold helical conformation with all O(6) in gt rotational positions, and the intrahelical O(2)---O(3) and O(2)---O(6) hydrogen bonds are substantially in agreement with previous studies. A new model for packing of the chains in the unit cell and the presence of crystallographic water is proposed. Packing appears to be stabilized by corner-to-center chain O(2)---O(2) hydrogen bonds. The nature of the transition from the amylose–DMSO complex to V_a-amylose was considered and it is shown that the transition involves translation of the amylose chains parallel to the a and b unit cell axes with only slight changes in the orientation of the helix. No significant conformational changes result from the transition.     

Crystal and molecular structure of Destruxin B

The cyclic hexadepsipeptide mycotoxin Destruxin B, produced by Metarrhizium anisopliae, crystallizes in the orthorhombic space group P212121, with a = 11.010(2)A, b = 14.679(5)A, c = 21.273(7)A and Z = 4. The structure was solved by direct methods and refined by least-squares technique to a final unweighted R value of 0.051, for 3361 reflections with I greater than 3 sigma (I). The backbone of the peptide is asymmetric and is made of 5 trans peptide and ester units and 1 cis peptide unit. The backbone conformation of this cyclic depsipeptide is very similar to that of Roseotoxin B, an analogous mycotoxin produced by Trichothecium roseum. The conformation in the crystalline state also correlates well with the solution conformation, as reported from proton n.m.r. studies. The crystal packing is directed by van der Waals contacts.