Solid-Phase Synthesis and NMR Structural Studies of the Marine Antibacterial Cyclic Tetrapeptide: Cyclo[GSPE]

Twelve-membered head-to-tail cyclic tetrapeptides (CTPs) are rigid molecules found in nature and possess a diverse range of biological activities. A possible reason may be due to their ability to adopt rigid conformations in solution mimicking reverse-turns. Reverse-turns are common structural motifs which serve as molecular recognition sites in many protein-receptor interactions. In this paper, we describe the solid-phase synthesis of the antibacterial cyclic tetrapeptide cyclo[Gly-Ser-Pro-Glu] (cyclo[GSPE]), first isolated from the Ruegeria strain of marine bacteria by Mitova et al. (J Nat Prod 67:1178–1181, 2004). Our NMR experiments in H₂O:D₂O:DMSO (18:1:1) revealed that it possessed three conformations in an approximate ratio of 4:2:1 based on NMR amide peak intensities. 2D NMR studies and computer calculations revealed that the major conformer adopted a reverse-turn conformation and have ω torsion angles twisted by up to 2°, with two transoid amide bonds between Gly-Ser, Pro-Glu and two cisoid amide bonds between Ser-Pro, Glu-Gly in a cis–trans-cis–trans (ctct) pattern. This supports previous reports that majority of CTPs adopt a ctct pattern when dissolved in hydrogen-bond disrupting solvents (Che and Marshall in J Med Chem 49:111–124, 2006 and references cited therein). An ensemble of ten lowest-energy-minimised 3D structures generated using XPLOR-NIH software revealed that cyclo[GSPE] possessed a rigid backbone ring scaffold. The remaining two minor conformers were present in quantities too low for NMR structural studies.     

Reductive methylation andpKa determination of the lysine side chains in calbindin D9k

The Lys residues in the 75-residue Ca²⁺-binding protein calbindin D_9k were reductively methylated with¹³C-enriched formaldehyde. The possible structural effects resulting from the chemical modification were critically investigated by comparing two-dimensional NMR spectra and the exchange rates of some of the amide protons of the native and the modified protein. Our results show that the protein retains its structure even though 10 Lys out of a total of 75 amino acid residues were modified. In the Ca²⁺- and apo-forms of the protein, the¹³C-methylated Lys residues can be detected with high sensitivity and resolution using two-dimensional (¹H,¹³C)-heteronuclear multiple quantum coherence (HMQC) NMR spectroscopy. ThepKa values of the individual Lys residues in Ca²⁺-calbindin D_9k and apo-calbindin D_9k were obtained by combiningpH titration experiments and (¹H,¹³C)-HMQC NMR spectroscopy. Each Lys residue in the Ca²⁺- and apo-forms of calbindin D_9k has a uniquepKa value. The LyspKa values in the calcium protein range from 9.3 to 10.9, while those in the apo-protein vary between 9.7 and 10.7. Although apo-calbindin D_9k has a very similar structure compared to Ca²⁺-calbindin D_9k, the removal of two Ca²⁺ ions from the protein leads to an increase of thepKa values of the Lys residues.     

Diffusion and sedimentation interaction parameters for measuring the second virial coefficient and their utility as predictors of protein aggregation

The concentration-dependence of the diffusion and sedimentation coefficients (k_D and k_s, respectively) of a protein can be used to determine the second virial coefficient (B₂), a parameter valuable in predicting protein-protein interactions. Accurate measurement of B₂ under physiologically and pharmaceutically relevant conditions, however, requires independent measurement of k_D and k_s via orthogonal techniques. We demonstrate this by utilizing sedimentation velocity (SV) and dynamic light scattering (DLS) to analyze solutions of hen-egg white lysozyme (HEWL) and a monoclonal antibody (mAb1) in different salt solutions. The accuracy of the SV-DLS method was established by comparing measured and literature B₂ values for HEWL. In contrast to the assumptions necessary for determining k_D and k_s via SV alone, k_D and ks were of comparable magnitudes, and solution conditions were noted for both HEWL and mAb1 under which 1), k_D and k_s assumed opposite signs; and 2), k_D ≥ k_s. Further, we demonstrate the utility of k_D and k_s as qualitative predictors of protein aggregation through agitation and accelerated stability studies. Aggregation of mAb1 correlated well with B₂, k_D, and k_s, thus establishing the potential for k_D to serve as a high-throughput predictor of protein aggregation.     

Shift of the high-performance liquid chromatographic retention times of metabolites in relation to the original drug on an RP8 column with acidic mobile phase

The effect of the structural change in the metabolization of drugs on the HPLC retention time with an RP8 column with an acetonitrile–phosphate buffer (pH 2.3) as the mobile phase was investigated at model compound pairs of 29 functionalization reactions. A more or less typical region for T_M=log(k′_M/k′_D) was found for each of these reactions (with k′_M and k′_D being the capacity factors of the metabolite and the drug, respectively), which can be explained by an increase or a decrease of the hydrophilic properties caused by the structural change. This effect is superimposed by an essential influence of the unchanged part of the molecule and in some cases by special intramolecular interactions like the hydrogen bond. Despite the more complicated structure of real drugs the results obtained at the model compound pairs were confirmed for most of the 55 metabolite/drug pairs. The practical use of the T_M values as a support to distinguish between different metabolites in the HPLC-DAD analysis of intoxications is demonstrated with cases of poisoning with diphenhydramine, propafenone and methaqualone.     

Molecular dynamics simulations reveal initial structural and dynamic features for the A2AR as a result of ligand binding

G-protein-coupled receptors (GPCRs) are membrane proteins that have a wide variety of physiological roles. Adenosine receptors belong to the GPCR family. Adenosine receptors are implicated in many physiological disorders, such as Parkinson's disease, Huntington's disease, inflammatory and immune's disease and many others. Interestingly, crystal structures of the active and inactive conformations of the A₂-subtype adenosine receptor (A₂AR) have been solved. These two structures could be used to get insights about the conformational changes that occur during the process of activation/inactivation processes of this receptor. Therefore, two ligand-free simulations of the native active (PDB code: 3QAK) and inactive (PDB code: 3EML) conformations of the A₂AR and two halo-simulations were carried out to observe the initial conformational changes induced by coupling adenosine to the inactive conformation and caffeine to the active conformation. Furthermore, we constructed an A₂AR model that contained four thermostabilising mutations, L48A, T65A, Q89A and A54L, which had previously been determined to stabilise the bound conformation of the agonist, and we ran molecular dynamics simulations of this mutant to investigate how these point mutations might affect the inactive conformation of this receptor. This study provides insights about the initial structural and dynamic features that occur as a result of the binding of caffeine and adenosine in the active and inactive A₂AR structures, respectively, as well as the introduction of some mutations on the inactive structure of the A₂AR. Moreover, we provide useful and detailed information regarding structural features such as toggle switch and ionic lock during the activation/inactivation processes of this receptor.     

Light dependence of the decay of the proton gradient in broken chloroplasts

The initial rates and steady-state values of proton uptake by broken chloroplasts have been measured as functions of light intensity at various concentrations of chlorophyll, pyocyanine, supporting electrolyte, buffer, as well as pH and temperature. Kinetic analysis of the data shows that the rate of decay of proton gradient due to backward leakage depends on light intensity. Under steady illumination, the decay constant k_L is equal to k_D + mR₀, where R₀ is the initial rate of proton uptake which is a function of light intensity, k_D is the decay constant in the dark and m is a parameter which is independent of light intensity. Treatment of chloroplasts with lysolecithin, neutral detergent, 2,4-dinitrophenol, or valinomycin in the presence of K⁺ increases k_D without affecting m. Treatment with N,N′-dicyclohexylcarbodiimide or adenylyl imidodiphosphate under appropriate conditions decreases m without affecting k_D. Treatment with glutaraldehyde makes k_L independent of light intensity and hence m = 0. These results suggest that the light-dependent part (mR₀) of k_L is due to leakage of protons through the coupling factor (CF₁-CF₀) complex which can open or close depending on light intensity and that the light-independent part (k_D) of the decay constant k_L is due to proton leakage elsewhere.     

The Genetic Algorithm and the Conformational Search of Polypeptides and Proteins

Abstract

The genetic algorithm is a technique of function optimization derived from the principles of evolutionary theory. We have adapted it to perform conformational search on polypeptides and proteins. The algorithm was first tested on several small polypeptides and the 46 amino acid protein crambin under the AMBER potential energy function. The probable global minimum conformations of the polypeptides were located 90% of the time and a non-native conformation of crambin was located that was 150kcal/mol lower in potential energy than the minimized crystal structure conformation. Next, we used a knowledge-based potential function to predict the structures of melittin, pancreatic polypeptide, and crambin. A 2.31 Å ΔRMS conformation of melittin and a 5.33 Å ΔRMS conformation of pancreatic polypeptide were located by genetic algorithm-based conformational search under the knowledge-based potential function. Although the ΔRMS of pancreatic polypeptide was somewhat high, most of the secondary structure was correct. The secondary structure of crambin was predicted correctly, but the potential failed to promote packing interactions. Finally, we tested the packing aspects of our potential function by attempting to predict the tertiary structure of cytochrome b ₅₆₂ given correct secondary structure as a constraint. The final predicted conformation of cytochrome b ₅₆₂ was an almost completely extended continuous helix which indicated that the knowledge-based potential was useless for tertiary structure prediction. This work serves as a warning against testing potential functions designed for tertiary structure prediction on small proteins.     

Fast Protein Loop Sampling and Structure Prediction Using Distance-Guided Sequential Chain-Growth Monte Carlo Method

Loops in proteins are flexible regions connecting regular secondary structures. They are often involved in protein functions through interacting with other molecules. The irregularity and flexibility of loops make their structures difficult to determine experimentally and challenging to model computationally. Conformation sampling and energy evaluation are the two key components in loop modeling. We have developed a new method for loop conformation sampling and prediction based on a chain growth sequential Monte Carlo sampling strategy, called Distance-guided Sequential chain-Growth Monte Carlo (DiSGro). With an energy function designed specifically for loops, our method can efficiently generate high quality loop conformations with low energy that are enriched with near-native loop structures. The average minimum global backbone RMSD for 1,000 conformations of 12-residue loops is Å, with a lowest energy RMSD of Å, and an average ensemble RMSD of Å. A novel geometric criterion is applied to speed up calculations. The computational cost of generating 1,000 conformations for each of the x loops in a benchmark dataset is only about cpu minutes for 12-residue loops, compared to ca cpu minutes using the FALCm method. Test results on benchmark datasets show that DiSGro performs comparably or better than previous successful methods, while requiring far less computing time. DiSGro is especially effective in modeling longer loops (– residues).     

A competition smFRET assay to study ligand-induced conformational changes of the dengue virus protease

Ligand binding to proteins often is accompanied by conformational transitions. Here, we describe a competition assay based on single molecule Förster resonance energy transfer (smFRET) to investigate the ligand-induced conformational changes of the dengue virus (DENV) NS2B-NS3 protease, which can adopt at least two different conformations. First, a competitive ligand was used to stabilize the closed conformation of the protease. Subsequent addition of the allosteric inhibitor reduced the fraction of the closed conformation and simultaneously increased the fraction of the open conformation, demonstrating that the allosteric inhibitor stabilizes the open conformation. In addition, the proportions of open and closed conformations at different concentrations of the allosteric inhibitor were used to determine its binding affinity to the protease. The K_D value observed is in accordance with the IC₅₀ determined in the fluorometric assay. Our novel approach appears to be a valuable tool to study conformational transitions of other proteases and enzymes.     

A Model for the Action of Vinblastine in Vivo

A model for the action of vinblastine (VLB) on cells multiplying exponentially in vivo with a generation time, T_G, has been derived. It is based on the assumption that cells attempting to pass through mitosis in the presence of VLB lose their proliferative capacity and that this lethal effect occurs only when the cells are exposed to a concentration of VLB which is above a critical value, C_k. The model leads to two predictions. First, that the percentage of cells surviving at any time after exposure to a dose, D, of VLB is 100% if D < D_k and decreases to 0% after a time, T_G, following a dose D ≥ D_k·2^{T G/T1/2}, where D_k represents the dose of VLB required to produce the concentration C_k, and T_1/2 is the half-life of the VLB in vivo. Second, that the time, T_G, at which the percentage of cells surviving an exposure to VLB, at doses greater than D_k·2^{U G/T1/2}, decreases to zero should be equal to the generation time of the cells. Both of these predictions were confirmed experimentally which indicates that the model adequately explains the action of VLB in vivo.     

Conformational calculations on the Ala14-Pro27 LCI segment of rabbit skeletal myosin

In order to define the conformational characteristics of a singular Ala¹⁴-Pro²⁷ segment in myosin LC1, conformational calculations were performed using the Simplex algorithm of Nelder and Mead (Computer J. 7 (1965) 308–313) in the ACME program proposed by Tournarie (J. Appl. Cryst. 6 (1973) 309–346). The (Ala-Pro) n = 1 unit was assigned a given conformation x; the conformation energy was then minimized for n = 1 to n = 7 by adjusting structural parameters (angle values). Similarly, 13 different possible conformations were optimized and compared, showing that a (β₂R)₇ conformation is favored by about 20 kcal per mol over the next most probable conformation (C₇R)₇. In the β₂R conformation, the (Ala-Pro)₇ segment is a wide helix, 15 Å in length and 8.65 Å in diameter, while the C₇R conformation results in a semi-extended structure of 25 Å long, with an approximate diameter of 6 Å. These characteristics are in agreement with available experimental data and putative functions of the LC1 N-terminus.     

Theoretical Conformational Analysis in the Determination of Productive Conformations of Substrates for Acetylcholinesterase and Butyrylcholinesterase

All the equilibrium conformations of 34 analogues of acetylcholine (ACh) with the general formula R-C(O)O-Alk-N⁺(CH₃)₃ are calculated by the method of molecular mechanics. In the series R-C(O)O-(CH₂)₂-N⁺(CH₃)₃, a reliable correlation is found between the molecular volume of the substrate and the rate of its hydrolysis by acetylcholinesterase (AChE); the absence of such a correlation is demonstrated for butyryl-cholinesterase (BChE). Theoretical conformational analysis confirms that the completely extended tt conformation of ACh is productive for the hydrolysis by AChE, which agrees with the results of X-ray analysis of AChE. AChE is shown to hydrolyze only those substrates that form equilibrium conformers compatible in the mutual arrangement of trimethylammonium group, carbonyl carbon, and carbonyl oxygen with the tt conformation of ACh; in this case, the rate of substrate hydrolysis depends on the total population of these conformers. A reliable correlation was found between the population of the semifolded (tg^?) conformation of the choline moiety of substrate molecules and rate of their BChE hydrolysis. In a series of CH₃-C(O)O-Alk-N⁺(CH₃)₃, the rate of BChE hydrolysis is demonstrated to depend on the total population of conformations compatible in the mutual arrangement of functionally important atoms with the tg^? conformation of ACh. The tg^? conformation of ACh is concluded to be productive for BChE hydrolysis. Similar orientations of the substrate molecules relative to the catalytic triads of both AChE and BChE are proven to coincide upon the substrate productive sorption in their active sites. It is hypothesized that the sorption stage is rate-limiting in cholinesterase hydrolysis and the enzyme hydrolyzes the ACh molecule in its energetically favorable conformation.     

Structure of the Sm Binding Site From Human U4 snRNA Derived From a 3 ns PME Molecular Dynamics Simulation

Abstract

A molecular dynamics simulation of the Sm binding site from human U4 snRNA was undertaken to determine the conformational flexibility of this region and to identify RNA conformations that were important for binding of the Sm proteins. The RNA was fully-solvated (>9,000 water molecules) and charge neutralized by inclusion of potassium ions. A three nanosecond MD simulation was conducted using AMBER with long-range electrostatic forces considered using the particle mesh Ewald summation method. The initial model of the Sm binding site region had the central and 3′ stem-loops that flanked the Sm site co-axial with one another, and with the single-stranded Sm binding site region ([I] conformation). During the course of the trajectory, the axes of the 3′ stem-loop, and later the central stem-loop, became roughly orthogonal from their original anti-parallel orientation. As these conformational changes occurred, the snRNA adopted first an [L] conformation, and finally a [U] conformation. The [U] conformation was more stable than either the [I] or [L] conformations, and persisted for the final 1 ns of the trajectory. Analysis of the structure resulting from the MD simulations revealed the bulged nucleotide, U₁₁₄, and the mismatched A₉₁-G₁₁₀ base pair provided distinctive structural features that may enhance Sm protein binding. Based on the results of the MD simulation and the available experimental data, we proposed a mechanism for the binding of the Sm protein sub-complexes to the snRNA. In this model, the D₁/D₂ and E/F/G Sm protein sub-complexes first bind the snRNA in the [U] conformation, followed by conformational re-arrangement to the [I] conformation and binding of the D₃/B Sm protein sub-complex.     

Machine Learning Methods for X-Ray Scattering Data Analysis from Biomacromolecular Solutions

Small-angle x-ray scattering (SAXS) of biological macromolecules in solutions is a widely employed method in structural biology. SAXS patterns include information about the overall shape and low-resolution structure of dissolved particles. Here, we describe how to transform experimental SAXS patterns to feature vectors and how a simple k-nearest neighbor approach is able to retrieve information on overall particle shape and maximal diameter (D_max) as well as molecular mass directly from experimental scattering data. Based on this transformation, we develop a rapid multiclass shape-classification ranging from compact, extended, and flat categories to hollow and random-chain-like objects. This classification may be employed, e.g., as a decision block in automated data analysis pipelines. Further, we map protein structures from the Protein Data Bank into the classification space and, in a second step, use this mapping as a data source to obtain accurate estimates for the structural parameters (D_max, molecular mass) of the macromolecule under study based on the experimental scattering pattern alone, without inverse Fourier transform for D_max. All methods presented are implemented in a Fortran binary DATCLASS, part of the ATSAS data analysis suite, available on Linux, Mac, and Windows and free for academic use.     

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D₃ dopamine receptor (D₃R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D₃R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.     

Paramagnetic relaxation enhancement‐assisted structural characterization of a partially disordered conformation of ubiquitin

Nuclear magnetic resonance (NMR) is a powerful tool to study three‐dimensional structures as well as protein conformational fluctuations in solution, but it is compromised by increases in peak widths and missing signals. We previously reported that ubiquitin has two folded conformations, N₁ and N₂ and plus another folded conformation, I, in which some amide group signals of residues 33–41 almost disappeared above 3 kbar at pH 4.5 and 273 K. Thus, well‐converged structural models could not be obtained for this region owing to the absence of distance restraints. Here, we reexamine the problem using the ubiquitin Q41N variant as a model for this locally disordered conformation, I. We demonstrate that the variant shows pressure‐induced loss of backbone amide group signals at residues 28, 33, 36, and 39–41 like the wild‐type, with a similar but smaller effect on CαH and CβH signals. In order to characterize this I structure, we measured paramagnetic relaxation enhancement (PRE) under high pressure to obtain distance restraints, and calculated the structure assisted by Bayesian inference. We conclude that the more disordered I conformation observed at pH 4.0, 278 K, and 2.5 kbar largely retained the N₂ conformation, although the amide groups at residues 33–41 have more heterogeneous conformations and more contact with water, which differ from the N₁ and N₂ states. The PRE‐assisted strategy has the potential to improve structural characterization of proteins that lack NMR signals, especially for relatively more open and hydrated protein conformations.     

X-Pro peptides. A theoretical study of the hydrogen bonded conformations of (α-aminoisobutyryl-l-prolyl)n sequences

Intramolecularly hydrogen bonded conformations of (Aib-Pro)_n sequences have been analysed theoretically. Both 4→1 (C₁₀ and 3→1 (C₇ hydrogen bonded regular structures are shown to be stereochemically feasible. Conformational energies for the helical structures have been estimated using classical potential energy methods. Both C₁₀ and C₇ conformations have very similar energies. Pyrrolidine ring puckering has a pronounced effect on the energies, and only Cγ-endo puckered Pro residues can be accommodated. The theoretical calculations using spectroscopic data suggest that the recently proposed novel 3₁₀ helical conformation for benzyloxycarbonyl(Aib-Pro)₄-methyl ester is in solution, is indeed energetically and stereochemically favourable.     

Cloning, purification, and characterization of galactomannan-degrading enzymes from Myceliophthora thermophila

Genes of β-mannosidase 97 kDa, GH family 2 (bMann9), β-mannanase 48 kDa, GH family 5 (bMan2), and α-galactosidase 60 kDa, GH family 27 (aGal1) encoding galactomannan-degrading glycoside hydrolases of Myceliophthora thermophila C1 were successfully cloned, and the recombinant enzymes were purified to homogeneity and characterized. bMann9 displays only exo-mannosidase activity, the K _m and k _cat values are 0.4 mM and 15 sec^?1 for p-nitrophenyl-β-D-mannopyranoside, and the optimal pH and temperature are 5.3 and 40°C, respectively. bMann2 is active towards galac-tomannans (GM) of various structures. The K _m and k _cat values are 1.3 mg/ml and 67 sec^?1 for GM carob, and the optimal pH and temperature are 5.2 and 69°C, respectively. aGal1 is active towards p-nitrophenyl-α-D-galactopyranoside (PNPG) as well as GM of various structures. The K _m and k _cat values are 0.08 mM and 35 sec^?1 for PNPG, and the optimal pH and temperature are 5.0 and 60°C, respectively.     

Conformational studies of per-O-trimethylsilyl derivatives of D-fructoses and oligosaccharides containing β-D-fructofuranose residues by 220- and 300-MHz P.M.R. spectroscopy

The interpretation of 220- and 300-MHz P.M.R. spectra and the accurate chemical shifts and coupling constants of a number of per-O-trimethylsilyl-(TMS-) D-fructose derivatives and TMS-oligosaccharides containing β-D-fructofuranose residues are presented. On the basis of calculations with an adapted Karplus equation it is concluded that TMS-α- and -β-D-fructopyranose occur in the ²C₅(D) chair conformation whereas the D-glucopyranose rings in the oligosaccharides adopt the usual ⁴C₁(D) chair conformation. The structure of the latter units is very similar to that of TMS-α-_D-glucopyranose. The ⁴E(D) envelope and ⁴T₅(D) twist are the principal conformations of the D-fructofuranose rings. The conformation of the furanose ring depends on the number and kind of monosaccharide units attached thereto. The calculated, preferred conformation of the C-5-CH₂OTMS group of the D-fructofuranose moieties correlates with the time-averaged displacement of C-4 above the plane of C-2, C-3, and O-5.