Crystal structure of Ser-22/Ile-25 form crambin confirms solvent, side chain substate correlations

Plant peptides - redefining an area of ribosomally synthesized and post-translationally modified peptides

Covering 1965 to February 2024Plants are prolific peptide chemists and are known to make thousands of different peptidic molecules. These peptides vary dramatically in their size, chemistry, and bioactivity. Despite their differences, all plant peptides to date are biosynthesized as ribosomally synthesized and post-translationally modified peptides (RiPPs). Decades of research in plant RiPP biosynthesis have extended the definition and scope of RiPPs from microbial sources, establishing paradigms and discovering new families of biosynthetic enzymes. The discovery and elucidation of plant peptide pathways is challenging due to repurposing and evolution of housekeeping genes as both precursor peptides and biosynthetic enzymes and due to the low rates of gene clustering in plants. In this review, we highlight the chemistry, biosynthesis, and function of the known RiPP classes from plants and recommend a nomenclature for the recent addition of BURP-domain-derived RiPPs termed burpitides. Burpitides are an emerging family of cyclic plant RiPPs characterized by macrocyclic crosslinks between tyrosine or tryptophan side chains and other amino acid side chains or their peptide backbone that are formed by copper-dependent BURP-domain-containing proteins termed burpitide cyclases. Finally, we review the discovery of plant RiPPs through bioactivity-guided, structure-guided, and gene-guided approaches.

Plant defensins: types, mechanism of action and prospects of genetic engineering for enhanced disease resistance in plants

Natural antimicrobial peptides have been shown as one of the important tools to combat certain pathogens and play important role as a part of innate immune system in plants and, also adaptive immunity in animals. Defensin is one of the antimicrobial peptides with a diverse nature of mechanism against different pathogens like viruses, bacteria and fungi. They have a broad function in humans, vertebrates, invertebrates, insects, and plants. Plant defensins primarily interact with membrane lipids for their biological activity. Several antimicrobial peptides (AMPs) have been overexpressed in plants for enhanced disease protection. The plants defensin peptides have been efficiently employed as an effective strategy for control of diseases in plants. They can be successfully integrated in plants genome along with some other peptide genes in order to produce transgenic crops for enhanced disease resistance. This review summarizes plant defensins, their expression in plants and enhanced disease resistance potential against phytopathogens.

Uncertainty, imprecision, and many-valued logics in protein bioinformatics

Understanding proteins, their structures, functions, mutual interactions, activity in cellular reactions, interactions with drugs, and expression in body cells is a key to efficient medical diagnosis, drug production, and treatment of patients. Machine learning and data exploration methods supported by many-valued logics allow to grasp the imprecision and uncertainties that naturally occur in proteins and other biomolecules. Many-valued logics, like Łukasiewicz logic or fuzzy logic, are non-classical logics that do not restrict the number of truth values to only two values of true or false, but they allow for a larger set of truth degrees. In this paper, we briefly review the use of many-valued logics, especially the fuzzy logic, in bioinformatics. Then, we focus on protein bioinformatics, and present selected applications of many-valued logics in the analysis of complex protein structures, including; (1) potential-based protein similarity searching, (2) matching proteins on the basis of secondary structures, (3) 3D protein structure alignment, (4) prediction of intrinsically disordered proteins, and (5) fuzzy querying in large collections of Big macromolecular Data. Results of presented studies show that the utilization of many-valued logics can enrich the investigations of protein molecules, in which uncertainty and imprecision are prevalent problems. The paper discusses all observed benefits brought by the application of many-valued logics in investigations related to selected protein analyzes carried out by the author.

Elucidating the fundamental forces in protein crystal formation: the case of crambin

This study demonstrates the feasibility of periodic all-electron hybrid density functional theory calculations in the description of protein crystals, using crambin as a test case.

Molecular simulations of proteins have been usually accomplished through empirical or semi-empirical potentials, due to the large size and inherent complexity of these biological systems. On the other hand, a theoretical description of proteins based on quantum-mechanical methods would however provide an unbiased characterization of their electronic properties, possibly offering a link between these and the ultimate biological activity. Yet, such approaches have been historically hindered by the large amount of requested computational power. Here we demonstrate the feasibility of periodic all-electron density functional theory calculations in the description of the crystal of the protein crambin (46 aminoacids), which is determined with exceptional structural accuracy. We have employed the hybrid B3LYP functional, coupled to an empirical description of London interactions (D*) to simulate the crambin crystal with an increasing amount of lattice water molecules in the cell (up to 172H₂O per cell). The agreement with the experiment is good for both protein geometry and protein–water interactions. The energetics was computed to predict crystal formation energies, protein–water and protein–protein interaction energies. We studied the role of dispersion interactions which are crucial for holding the crambin crystal in place. B3LYP-D* electrostatic potential and dipole moment of crambin as well as the electronic charge flow from crambin to the solvating water molecules (0.0015e per H₂O) have also been predicted. These results proved that quantum-mechanical simulations of small proteins, both free and in their crystalline state, are now feasible in a reasonable amount of time, by programs capable of exploiting high performance computing architectures, allowing the study of protein properties not easily amenable through classical force fields.

Antimicrobial Peptides from Plants

Plant antimicrobial peptides (AMPs) have evolved differently from AMPs from other life forms. They are generally rich in cysteine residues which form multiple disulfides. In turn, the disulfides cross-braced plant AMPs as cystine-rich peptides to confer them with extraordinary high chemical, thermal and proteolytic stability. The cystine-rich or commonly known as cysteine-rich peptides (CRPs) of plant AMPs are classified into families based on their sequence similarity, cysteine motifs that determine their distinctive disulfide bond patterns and tertiary structure fold. Cystine-rich plant AMP families include thionins, defensins, hevein-like peptides, knottin-type peptides (linear and cyclic), lipid transfer proteins, α-hairpinin and snakins family. In addition, there are AMPs which are rich in other amino acids. The ability of plant AMPs to organize into specific families with conserved structural folds that enable sequence variation of non-Cys residues encased in the same scaffold within a particular family to play multiple functions. Furthermore, the ability of plant AMPs to tolerate hypervariable sequences using a conserved scaffold provides diversity to recognize different targets by varying the sequence of the non-cysteine residues. These properties bode well for developing plant AMPs as potential therapeutics and for protection of crops through transgenic methods. This review provides an overview of the major families of plant AMPs, including their structures, functions, and putative mechanisms.

Introducing DInaMo: A Package for Calculating Protein Circular Dichroism Using Classical Electromagnetic Theory

The dipole interaction model is a classical electromagnetic theory for calculating circular dichroism (CD) resulting from the π-π* transitions of amides. The theoretical model, pioneered by J. Applequist, is assembled into a package, DInaMo, written in Fortran allowing for treatment of proteins. DInaMo reads Protein Data Bank formatted files of structures generated by molecular mechanics or reconstructed secondary structures. Crystal structures cannot be used directly with DInaMo; they either need to be rebuilt with idealized bond angles and lengths, or they need to be energy minimized to adjust bond lengths and bond angles because it is common for crystal structure geometries to have slightly short bond lengths, and DInaMo is sensitive to this. DInaMo reduces all the amide chromophores to points with anisotropic polarizability and all nonchromophoric aliphatic atoms including hydrogens to points with isotropic polarizability; all other atoms are ignored. By determining the interactions among the chromophoric and nonchromophoric parts of the molecule using empirically derived polarizabilities, the rotational and dipole strengths are determined leading to the calculation of CD. Furthermore, ignoring hydrogens bound to methyl groups is initially explored and proves to be a good approximation. Theoretical calculations on 24 proteins agree with experiment showing bands with similar morphology and maxima.

Crystal structure of small protein crambin at 0.48 Å resolution

With the development of highly brilliant and extremely intense synchrotron X-ray sources, extreme high-resolution limits for biological samples are now becoming attainable. Here, a study is presented that sets the record in crystallographic resolution for a biological macromolecule. The structure of the small protein crambin was determined to 0.48 Å resolution on the PETRA II ring before its conversion to a dedicated synchrotron-radiation source. The results reveal a wealth of details in electron density and demonstrate the possibilities that are potentially offered by a high-energy source. The question now arises as to what the true limits are in terms of what can be seen at such high resolution. From what can be extrapolated from the results using crystals of crambin, this limit would be at approximately 0.40 Å, which approaches that for smaller compounds.

Performance of phased rotation, conformation and translation function: accurate protein model building with tripeptidic and tetrapeptidic fragments

The automatic building of protein structures with tripeptidic and tetrapeptidic fragments was investigated. The oligopeptidic conformers were positioned in the electron-density map by a phased rotation, conformation and translation function and refined by a real-space refinement. The number of successfully located fragments lay within the interval 75–95% depending on the resolution and phase quality. The overlaps of partially located fragments were analyzed. The correctly positioned fragments were connected into chains. Chains formed in this way were extended directly into the electron density and a sequence was assigned. In the initial stage of the model building the number of located fragments was between 60% and 95%, but this number could be increased by several cycles of reciprocal-space refinement and automatic model rebuilding. A nearly complete structure can be obtained on the condition that the resolution is reasonable. Computer graphics will only be needed for a final check and small corrections.

Analysis of protein hydration in ultrahigh-resolution structures of the SRP GTPase Ffh

Two new structures of the SRP GTPase Ffh have been determined at 1.1 A resolution and provide the basis for comparative examination of the extensive water structure of the apo conformation of these GTPases. A set of well defined water-binding positions have been identified in the active site of the two-domain ;NG' GTPase, as well as at two functionally important interfaces. The water hydrogen-bonding network accommodates alternate conformations of the protein side chains by undergoing local rearrangements and, in one case, illustrates binding of a solute molecule within the active site by displacement of water molecules without further disruption of the water-interaction network. A subset of the water positions are well defined in several lower resolution structures, including those of different nucleotide-binding states; these appear to function in maintaining the protein structure. Consistent arrangements of surface water between three different ultrahigh-resolution structures provide a framework for beginning to understand how local water structure contributes to protein-ligand and protein-protein binding in the SRP GTPases.

Plant gamma-thionins: novel insights on the mechanism of action of a multi-functional class of defense proteins

This review focuses on the first plant defense protein class described in literature, with growth inhibition activity toward pathogens. These peptides were named gamma-thionins or defensins, which are small proteins that can be classified into four main subtypes according to their specific functions. Gamma-thionins are small cationic peptides with different and special abilities. They are able to inhibit digestive enzymes or act against bacteria and/or fungi. Current research in this area focuses particularly these two last targets, being the natural crop plant defenses improved through the use of transgenic technology. Here, we will compare primary and tertiary structures of gamma-thionins and also will analyze their similarities to scorpion toxins and insect defensins. This last comparison offers some hypothesis for gamma-thionins mechanisms of action against certain pathogens. This specific area has benefited from the recent determination of many gamma-thionin structures. Furthermore, we also summarize molecular interactions between plant gamma-thionins and fungi receptors, which include membrane proteins and lipids, shedding some light over pathogen resistance. Researches on gamma-thionins targets could help on plant genetic improvement for production of increased resistance toward pathogens. Thus, positive results recently obtained for transgenic plants and future prospects in the area are also approached. Finally, gamma-thionins activity has also been studied for future drug development, capable of inhibit tumor cell growth in human beings.

Characterization and structural role of disulfide bonds in a highly knotted thionin from Pyrularia pubera

Disulfide bonds play a crucial role in the stabilization of the amphipathic folding of the diverse families of cysteine-rich antimicrobial peptides. The determination of cysteine pairings in these peptides has largely depended on sequence homology criteria, since the classical methods of disulfide bond characterization, which usually require proteolysis as a first step, encounter serious drawbacks derived from the tight folding and the presence of vicinal cysteines. We have chosen the Pyrularia pubera thionin, a 47-residue peptide with four internal disulfides and a remarkable resistance to most proteases, as a representative member of this type of cysteine-rich peptides and have shown that a combination of partial reduction and cyanylation readily allows the determination of its disulfide bonds. We have also studied by molecular dynamics and a combination of partial reduction and proteolysis the role of disulfide bonds in the stabilization of the tridimensional structure of this thionin and found a good agreement with our partial reduction data, suggesting that removal of only one disulfide bond is enough to significantly alter the folding of the peptide.

Proposal for molecular mechanism of thionins deduced from physico-chemical studies of plant toxins

We propose a molecular model for phospholipid membrane lysis by the ubiquitous plant toxins called thionins. Membrane lysis constitutes the first major effect exerted by these toxins that initiates a cascade of cytoplasmic events leading to cell death. X-ray crystallography, solution nuclear magnetic resonance (NMR) studies, small angle X-ray scattering and fluorescence spectroscopy provide evidence for the mechanism of membrane lysis. In the crystal structures of two thionins in the family, alpha(1)- and beta-purothionins (MW: approximately 4.8 kDa), a phosphate ion and a glycerol molecule are modeled bound to the protein. (31)P NMR experiments on the desalted toxins confirm phosphate-ion binding in solution. Evidence also comes from phospholipid partition experiments with radiolabeled toxins and with fluorescent phospholipids. This data permit a model of the phospholipid-protein complex to be built. Further, NMR experiments, one-dimensional (1D)- and two-dimensional (2D)-total correlation spectroscopy (TOCSY), carried out on the model compounds glycerol-3-phosphate (G3P) and short chain phospholipids, supported the predicted mode of phospholipid binding. The toxins' high positive charge, which renders them extremely soluble (>300 mg/mL), and the phospholipid-binding specificity suggest the toxin-membrane interaction is mediated by binding to patches of negatively charged phospholipids [phosphatidic acid (PA) or phosphatidyl serine (PS)] and their subsequent withdrawal. The formation of proteolipid complexes causes solubilization of the membrane and its lysis. The model suggests that the oligomerization may play a role in toxin's activation process and provides insight into the structural principles of protein-membrane interactions.

Myoglobin cavities provide interior ligand pathway

The myoglobin protein binds oxygen and catalyzes NO oxidation. As a key model protein, its dynamics have been well studied by spectroscopy and by crystallography as well as by simulation. Nonetheless, visualization of the mechanism of movement of ligands within myoglobin has been difficult. Coordinates of the A1 and A3 taxonomic spectral states of myoglobin from the 1 A crystal structure (1a6g) are generated as consistent sets of correlated clusters of residues with A or B crystal alternates. Analysis of cavities in these A1 and A3 conformations clarifies the pathway of ligand motion from distal entry through interior movement to the proximal side of the heme. Cavities opened up by buried alternate conformations link the distal to the proximal side of the heme. Structural conservation highlights the relevance of this pathway to human neuroglobin. Cavity migration via myoglobin crystal alternates provides a specific link of protein structure to protein dynamics and protein function and demonstrates the relevance of substates (discrete disorder) to function for all proteins.

Solvent Interaction of a Hsp70 Chaperone Substrate-Binding Domain Investigated with Water−NOE NMR Experiments

The interaction of solvent of the substrate binding domain of the bacterial heat shock 70 chaperone protein DnaK was studied in its apo form and with bound hydrophobic substrate peptide, using refined nuclear magnetic resonance experiments. Distinct differences between the two states of the protein were observed. According to our data, the apo form interacts more extensively with solvent than the peptide-bound form. Significantly, the open hydrophobic substrate binding cleft of DnaK in the apo form is found to contain several molecules of water which are displaced by the binding of the hydrophobic substrate, the peptide NRLLLTG. The solvent in the hydrophobic cleft has a residence time longer than 400 ps. It is predicted that the displacement of this trapped water must contribute to the binding free energy of the natural hydrophobic substrates of this class of protein-folding chaperone proteins.

On the nature of a glassy state of matter in a hydrated protein: Relation to protein function

Diverse biochemical and biophysical experiments indicate that all proteins, regardless of size or origin, undergo a dynamic transition near 200 K. The cause of this shift in dynamic behavior, termed a "glass transition," and its relation to protein function are important open questions. One explanation postulated for the transition is solidification of correlated motions in proteins below the transition. We verified this conjecture by showing that crambin's radius of gyration (Rg) remains constant below approximately 180 K. We show that both atom position and dynamics of protein and solvent are physically coupled, leading to a novel cooperative state. This glassy state is identified by negative slopes of the Debye-Waller (B) factor vs. temperature. It is composed of multisubstate side chains and solvent. Based on generalization of Adam-Gibbs' notion of a cooperatively rearranging region and decrease of the total entropy with temperature, we calculate the slope of the Debye-Waller factor. The results are in accord with experiment.

Automatic assignment of NOESY cross peaks and determination of the protein structure of a new world scorpion neurotoxin using NOAH/DIAMOD

The 3D NMR structures of the scorpion neurotoxin, CsE-v5, were determined from the same NOESY spectra with NOAH/DIAMOD, an automated assignment and 3D structure calculation software package, and with a conventional manual assignment combined with a distance geometry/simulated annealing (X-PLOR) refinement method. The NOESY assignments and the 3D structures obtained from the two independent methods were compared in detail. The NOAH/DIAMOD program suite uses feedback filtering and self-correcting distance geometry methods to automatically assign NOESY spectra and to calculate the 3D structure of a protein. NOESY cross peaks were automatically picked using a standard software package and combined with 74 manually assigned NOESY peaks to start the NOAH/DIAMOD calculations. After 63 NOAH/DIAMOD cycles, using REDAC procedures in the last 8 cycles, and final FANTOM constrained energy minimization, a bundle of 20 structures with the smallest target functions has a RMSD of 0.81 A for backbone atoms and 1.11 A for all heavy atoms to the mean structure. Despite some missing chemical shifts of side chain protons, 776 (including 74 manually assigned) of 1130 NOE peaks were unambiguously assigned, 150 peaks have more than one possible assignment compatible with the bundle structures, and only 30 peaks could not be assigned within the given chemical shift tolerance ranges in either the D1 or the D2 dimension. The remaining 174, mainly weak NOE peaks were not compatible with the final 20 best bundle structures at the last NOAH/DIAMOD cycle. The automatically determined structures agree well with the structures determined independently using the conventional method and the same NMR spectra, with the mean RMSD in well-defined regions of 0.84 A for bb and 1.48 A for all heavy atoms from residues 2-5, 18-26, 32-36, and 39-45. This study demonstrates the potential of the NOAH/DIAMOD program suite to automatically assign NMR data for proteins and determine their structure.

Accurate protein crystallography at ultra-high resolution: valence electron distribution in crambin

The charge density distribution of a protein has been refined experimentally. Diffraction data for a crambin crystal were measured to ultra-high resolution (0.54 A) at low temperature by using short-wavelength synchrotron radiation. The crystal structure was refined with a model for charged, nonspherical, multipolar atoms to accurately describe the molecular electron density distribution. The refined parameters agree within 25% with our transferable electron density library derived from accurate single crystal diffraction analyses of several amino acids and small peptides. The resulting electron density maps of redistributed valence electrons (deformation maps) compare quantitatively well with a high-level quantum mechanical calculation performed on a monopeptide. This study provides validation for experimentally derived parameters and a window into charge density analysis of biological macromolecules.

Water: foldase activity in catalyzing polypeptide conformational rearrangements

Polypeptide conformer interconversion in a low dielectric environment is shown to be highly dependent on water concentration. Water increases this rate by 10(3), apparently by catalyzing hydrogen bond exchange, and thereby presenting functional properties analogous to that of a foldase. This catalytic effect is demonstrated on the interconversion of a parallel gramicidin dimer into an antiparallel dimer. A Hill coefficient of 6.5 is observed, illustrating the highly cooperative nature of the process. Protein folding in nonpolar environs, such as the hydrophobic core of a protein or the hydrophobic domain of a lipid bilayer, may be contingent on and rate-limited by the scarcity of water.

Automated 2D NOESY assignment and structure calculation of Crambin(S22/I25) with the self-correcting distance geometry based NOAH/DIAMOD programs

The NOAH/DIAMOD program suite was used to automatically assign an experimental 2D NOESY spectrum of the 46 residue protein crambin(S22/I25), using feedback filtering and self-correcting distance geometry (SECODG). Automatically picked NOESY cross peaks were combined with 157 manually assigned peaks to start NOAH/DIAMOD calculations. At each cycle, DIAMOD was used to calculate an ensemble of 40 structures from these NOE distance constraints and random starting structures. The 10 structures with smallest target function values were analyzed by the structure-based filter, NOAH, and a new set of possible assignments was automatically generated based on chemical shifts and distance constraints violations. After 60 iterations and final energy minimization, the 10 structures with smallest target functions converged to 1.48 A for backbone atoms. Despite several missing chemical shifts, 426 of 613 NOE peaks were unambiguously assigned; 59 peaks were ambiguously assigned. The remaining 128 peaks picked automatically by FELIX are probably primarily noise peaks, with a few real peaks that were not assigned by NOAH due to the incomplete proton chemical shifts list.