Parameterization of the proline analogue Aze (azetidine-2-carboxylic acid) for molecular dynamics simulations and evaluation of its effect on homo-pentapeptide conformations

Strategies to Control Human Health Risks Arising from Antibiotics in the Environment: Molecular Modification of QNs for Enhanced Plant-Microbial Synergistic Degradation

In the present work, a comprehensive screening and evaluation system was established to improve the plant-microbial synergistic degradation effects of QNs. The study included the construction of a 3D-QSAR model, the molecular modification, environmental friendliness and functional evaluation of drugs, degradation pathway simulation, and human health risk assessment. Molecular dynamics was applied to quantify the binding capacity of QNs toward the plant degradation enzyme (peroxidase) and microbial degradation enzymes (manganese peroxidase, lignin peroxidase, and laccase). The fuzzy comprehensive evaluation method was used in combination with the weighted average method for normalization and assigning equal weights to the plant and microbial degradation effect values of the QNs. Considering the synergistic degradation effect value as the dependent variable and the molecular information of the QNs as the independent variable, a 3D-QSAR model was constructed for the plant-microbial synergistic degradation effect of QNs. The constructed model was then employed to conduct the molecular modification, environmental friendliness and functional evaluation, degradation pathway simulation, and human health risk assessment of transformation products using pharmacokinetics and toxicokinetics. The results revealed that the synergistic degradation effect 3D-QSAR (CoMSIA) model exhibited good internal and external prediction ability, fitting ability, stability, and no overfitting phenomenon. Norfloxacin (NOR) was used as the target molecule in the molecular modification. A total of 35 NOR derivatives with enhanced plant-microbial synergistic degradation effect (1.32-21.51%) were designed by introducing small-volume, strongly electronegative, and hydrophobic hydrogen bond receptor groups into the active group of the norfloxacin structure. The environment-friendliness and the functionality of NOR were evaluated prior to and after the modification, which revealed seven environment-friendly FQs derivatives exhibiting moderate improvement in stability and bactericidal efficacy. The simulation of the NOR plant and microbial degradation pathways prior to and after the modification and the calculation of the reaction energy barrier revealed Pathway A (D-17 to D-17-2) and Pathway B (D-17 to D-17-4) as the most prone degradation pathways in plants and Pathway A (D-17 to D-17-1) and Pathway B (D-17 to D-17-4) as the most prone degradation pathways in microorganisms. This demonstrated that the degradation of the modified NOR derivatives was significantly enhanced, with the hydroxylation and piperazine ring substitution reaction playing an important role in the degradation process. Finally, the parameters, including hepatotoxicity, mutagenicity, and rodent carcinogenicity, among others, predicted using the pharmacokinetics and toxicokinetics analyses revealed a significant reduction in the human health risk associated with the modified NOR, along with a considerable reduction in the toxicity of its transformation products, implying that the human health risk associated with the transformation products was reduced remarkably. The present study provides a theoretical basis for novel ideas and evaluation programs for improving the plant-microbial synergistic degradation of the QNs antibiotics for source control and drug design, thereby reducing the residues of these antibiotics and the associated hazard in the complex plant-soil environment, ultimately decreasing the potential risks to human health.

Environmentally Friendly Fluoroquinolone Derivatives with Lower Plasma Protein Binding Rate Designed Using 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation

Comparative molecular similarity index analysis (CoMSIA) was used to establish a three-dimensional quantitative structure–activity relationship (3D-QSAR) model with structural parameters of quinolones as the independent variables and plasma protein binding rate (logf_b) as the dependent variable to predict the logf_b values of remaining quinolones in this study. In addition, the mono-substituted and bis-substituted reaction schemes that significantly influenced the plasma protein binding rate of quinolones were determined through an analysis of the 3D-QSAR contour maps. It was found that the replacement of small groups, hydrophobic groups, electronegative groups, or hydrogen bond acceptor groups at the substitution sites significantly reduce the logf_b values of quinolone derivatives. Furthermore, the mechanism of decrease in binding rate between trovafloxacin (TRO) derivatives and plasma protein was revealed qualitatively and quantitatively based on molecular docking and molecular dynamics simulation. After modification of the target molecule, 11 TRO derivatives with low plasma protein binding rates were screened (reduced by 0.50–24.18%). Compared with the target molecule, the molecular genotoxicity and photodegradability of the TRO derivatives was higher (genotoxicity increased by 4.89–21.36%, and photodegradability increased by 9.04–20.56%), and their bioconcentration was significantly lower (by 36.90–61.41%).

Neuropilin-1 peptide-like ligands with proline mimetics, tested using the improved chemiluminescence affinity detection method

Many reports have suggested that NRP-1 acts as a co-receptor for VEGF-A165 and boosts tumour growth and metastasis. This NRP-1, due to its important role in tumour progression, triggered interest in the design of new molecules able to significantly inhibit NRP-1/VEGF-A165 interaction to suppress pathological angiogenesis. Our previous SAR studies of compounds, showing affinity for NRP-1, led us to develop branched peptides with general formula Lys(hArg)-AA2-AA3-Arg. Here, three series of analogues were synthesized, in which the middle fragment (AA2 and/or AA3) of initial sequences was substituted with unnatural Pro analogues with different rigidities and ring sizes. The synthesized compounds were screened for VEGF-A165 inhibitory activity on an improved assay (ELISA), which was selected based on our comparative inhibition study of the parent compounds, indicating that the method with chemiluminescence detection gives more accurate data. The results of affinity for NRP-1 and enzymatic stability of newly obtained compounds enabled the selection of new structures, showing a 2 and 4-fold lower IC50 value compared to parent peptides.

Amino acid substrates impose polyamine, eIF5A, or hypusine requirement for peptide synthesis

Whereas ribosomes efficiently catalyze peptide bond synthesis by most amino acids, the imino acid proline is a poor substrate for protein synthesis. Previous studies have shown that the translation factor eIF5A and its bacterial ortholog EF-P bind in the E site of the ribosome where they contact the peptidyl-tRNA in the P site and play a critical role in promoting the synthesis of polyproline peptides. Using misacylated Pro-tRNAPhe and Phe-tRNAPro, we show that the imino acid proline and not tRNAPro imposes the primary eIF5A requirement for polyproline synthesis. Though most proline analogs require eIF5A for efficient peptide synthesis, azetidine-2-caboxylic acid, a more flexible four-membered ring derivative of proline, shows relaxed eIF5A dependency, indicating that the structural rigidity of proline might contribute to the requirement for eIF5A. Finally, we examine the interplay between eIF5A and polyamines in promoting translation elongation. We show that eIF5A can obviate the polyamine requirement for general translation elongation, and that this activity is independent of the conserved hypusine modification on eIF5A. Thus, we propose that the body of eIF5A functionally substitutes for polyamines to promote general protein synthesis and that the hypusine modification on eIF5A is critically important for poor substrates like proline.

Docking and molecular dynamics simulations of the Fyn-SH3 domain with free and phospholipid bilayer-associated 18.5-kDa myelin basic protein (MBP)-Insights into a noncanonical and fuzzy interaction

The molecular details of the association between the human Fyn-SH3 domain, and the fragment of 18.5-kDa myelin basic protein (MBP) spanning residues S38-S107 (denoted as xα2-peptide, murine sequence numbering), were studied in silico via docking and molecular dynamics over 50-ns trajectories. The results show that interaction between the two proteins is energetically favorable and heavily dependent on the MBP proline-rich region (P93-P98) in both aqueous and membrane environments. In aqueous conditions, the xα2-peptide/Fyn-SH3 complex adopts a "sandwich""-like structure. In the membrane context, the xα2-peptide interacts with the Fyn-SH3 domain via the proline-rich region and the β-sheets of Fyn-SH3, with the latter wrapping around the proline-rich region in a form of a clip. Moreover, the simulations corroborate prior experimental evidence of the importance of upstream segments beyond the canonical SH3-ligand. This study thus provides a more-detailed glimpse into the context-dependent interaction dynamics and importance of the β-sheets in Fyn-SH3 and proline-rich region of MBP. Proteins 2017; 85:1336-1350. © 2017 Wiley Periodicals, Inc.

MyelStones: the executive roles of myelin basic protein in myelin assembly and destabilization in multiple sclerosis

The classic isoforms of myelin basic protein (MBP, 14–21.5 kDa) are essential to formation of the multilamellar myelin sheath of the mammalian central nervous system (CNS). The predominant 18.5-kDa isoform links together the cytosolic surfaces of oligodendrocytes, but additionally participates in cytoskeletal turnover and membrane extension, Fyn-mediated signalling pathways, sequestration of phosphoinositides and maintenance of calcium homoeostasis. All MBP isoforms are intrinsically disordered proteins (IDPs) that interact via molecular recognition fragments (MoRFs), which thereby undergo local disorder-to-order transitions. Their conformations and associations are modulated by environment and by a dynamic barcode of post-translational modifications, particularly phosphorylation by mitogen-activated and other protein kinases and deimination [a hallmark of demyelination in multiple sclerosis (MS)]. The MBPs are thus to myelin what basic histones are to chromatin. Originally thought to be merely structural proteins forming an inert spool, histones are now known to be dynamic entities involved in epigenetic regulation and diseases such as cancer. Analogously, the MBPs are not mere adhesives of compact myelin, but active participants in oligodendrocyte proliferation and in membrane process extension and stabilization during myelinogenesis. A central segment of these proteins is pivotal in membrane-anchoring and SH3 domain (Src homology 3) interaction. We discuss in the present review advances in our understanding of conformational conversions of this classic basic protein upon membrane association, including new thermodynamic analyses of transitions into different structural ensembles and how a shift in the pattern of its post-translational modifications is associated with the pathogenesis and potentially onset of demyelination in MS.

The effects of threonine phosphorylation on the stability and dynamics of the central molecular switch region of 18.5-kDa myelin basic protein

The classic isoforms of myelin basic protein (MBP) are essential for the formation and maintenance of myelin in the central nervous system of higher vertebrates. The protein is involved in all facets of the development, compaction, and stabilization of the multilamellar myelin sheath, and also interacts with cytoskeletal and signaling proteins. The predominant 18.5-kDa isoform of MBP is an intrinsically-disordered protein that is a candidate auto-antigen in the human demyelinating disease multiple sclerosis. A highly-conserved central segment within classic MBP consists of a proline-rich region (murine 18.5-kDa sequence -T92-P93-R94-T95-P96-P97-P98-S99-) containing a putative SH3-ligand, adjacent to a region that forms an amphipathic α-helix (P82-I90) upon interaction with membranes, or under membrane-mimetic conditions. The T92 and T95 residues within the proline-rich region can be post-translationally modified through phosphorylation by mitogen-activated protein (MAP) kinases. Here, we have investigated the structure of the α-helical and proline-rich regions in dilute aqueous buffer, and have evaluated the effects of phosphorylation at T92 and T95 on the stability and dynamics of the α-helical region, by utilizing four 36-residue peptides (S72-S107) with differing phosphorylation status. Nuclear magnetic resonance spectroscopy reveals that both the α-helical as well as the proline-rich regions are disordered in aqueous buffer, whereas they are both structured in a lipid environment (cf., Ahmed et al., Biochemistry 51, 7475-9487, 2012). Thermodynamic analysis of trifluoroethanol-titration curves monitored by circular dichroism spectroscopy reveals that phosphorylation, especially at residue T92, impedes formation of the amphipathic α-helix. This conclusion is supported by molecular dynamics simulations, which further illustrate that phosphorylation reduces the folding reversibility of the α-helix upon temperature perturbation and affect the global structure of the peptides through altered electrostatic interactions. The results support the hypothesis that the central conserved segment of MBP constitutes a molecular switch in which the conformation and/or intermolecular interactions are mediated by phosphorylation/dephosphorylation at T92 and T95.

Myelin management by the 18.5-kDa and 21.5-kDa classic myelin basic protein isoforms

The classic myelin basic protein (MBP) splice isoforms range in nominal molecular mass from 14 to 21.5 kDa, and arise from the gene in the oligodendrocyte lineage (Golli) in maturing oligodendrocytes. The 18.5-kDa isoform that predominates in adult myelin adheres the cytosolic surfaces of oligodendrocyte membranes together, and forms a two-dimensional molecular sieve restricting protein diffusion into compact myelin. However, this protein has additional roles including cytoskeletal assembly and membrane extension, binding to SH3-domains, participation in Fyn-mediated signaling pathways, sequestration of phosphoinositides, and maintenance of calcium homeostasis. Of the diverse post-translational modifications of this isoform, phosphorylation is the most dynamic, and modulates 18.5-kDa MBP's protein-membrane and protein-protein interactions, indicative of a rich repertoire of functions. In developing and mature myelin, phosphorylation can result in microdomain or even nuclear targeting of the protein, supporting the conclusion that 18.5-kDa MBP has significant roles beyond membrane adhesion. The full-length, early-developmental 21.5-kDa splice isoform is predominantly karyophilic due to a non-traditional P-Y nuclear localization signal, with effects such as promotion of oligodendrocyte proliferation. We discuss in vitro and recent in vivo evidence for multifunctionality of these classic basic proteins of myelin, and argue for a systematic evaluation of the temporal and spatial distributions of these protein isoforms, and their modified variants, during oligodendrocyte differentiation.

Proline: the distribution, frequency, positioning, and common functional roles of proline and polyproline sequences in the human proteome

Proline is an anomalous amino acid. Its nitrogen atom is covalently locked within a ring, thus it is the only proteinogenic amino acid with a constrained phi angle. Sequences of three consecutive prolines can fold into polyproline helices, structures that join alpha helices and beta pleats as architectural motifs in protein configuration. Triproline helices are participants in protein-protein signaling interactions. Longer spans of repeat prolines also occur, containing as many as 27 consecutive proline residues. Little is known about the frequency, positioning, and functional significance of these proline sequences. Therefore we have undertaken a systematic bioinformatics study of proline residues in proteins. We analyzed the distribution and frequency of 687,434 proline residues among 18,666 human proteins, identifying single residues, dimers, trimers, and longer repeats. Proline accounts for 6.3% of the 10,882,808 protein amino acids. Of all proline residues, 4.4% are in trimers or longer spans. We detected patterns that influence function based on proline location, spacing, and concentration. We propose a classification based on proline-rich, polyproline-rich, and proline-poor status. Whereas singlet proline residues are often found in proteins that display recurring architectural patterns, trimers or longer proline sequences tend be associated with the absence of repetitive structural motifs. Spans of 6 or more are associated with DNA/RNA processing, actin, and developmental processes. We also suggest a role for proline in Kruppel-type zinc finger protein control of DNA expression, and in the nucleation and translocation of actin by the formin complex.