Molecular dynamics of oligopeptides 6. A comparative study of poincare cross-section of monopeptide frames in media with different hydrophobicity

A comparative study of the molecular dynamics of natural amino acid residues and their closest homologues and isomers was carried out. Molecular dynamics protocols not interfering with the principle of equidistribution of energy with respect to degrees of freedom were used. Poincare cross-sections, auto- and cross-correlation of complex exponential curves as a function of dihedrons were considered. The classification of dynamic properties of conformational degrees of freedom in the series of amino acid residues was carried out.     

Molecular dynamics of human alpha-fetoprotein fragment LDSYQCT and its analogs at different dielectric constants

A comparative study of the conformation dynamics of the human alpha-fetoprotein fragment LDSYQCT and heptapeptides derived from it by point substitutions has revealed a significant influence of electrostatic interactions on the set of preferred conformations and dynamics of amino acid residues when the peptides with blocked termini are examined at ? = 1. Peptide flexibility rises when the termini are left free (charged). At ? = 10 or 80, the set of probable conformations for all residues expands to much the same extent, i.e., at higher permittivity of the medium the dynamic effects of amino acid changes are leveled off.     

A molecular dynamics simulation study of human alpha-fetoprotein-derived peptide and its analogues. Influence of dielectric constant

A comparative study of the conformation dynamics of a series of heptapeptides: the human alpha-fetoprotein fragment LDSYQCT and its seven analogues has been conducted. The effect of the dielectric constant of medium on the dynamics of heptapeptide conformation is considered. It is shown that electrostatic interactions have a marked effect on several accessible conformations and the dynamics of the behavior of amino acid residues.     

A comparative study of the interaction of Bacillus amyloliquefaciens ribonuclease (barnase) and Bacillus intermedius ribonuclease (binase) with barstar by brownian dynamics simulation

A comparative study of the interaction of two RNases (binase and barnase) with the polypeptide inhibitor barstar was performed by Brownian dynamics simulation. It was demonstrated that this method adequately reproduced the dependence of the association rate on the pH of solution as well as the effect of mutations at individual amino acid residues on the inhibition of barnase by barstar. Two types of energy-favorable binase-barstar encounter complexes were found. In type I complex, the amino acid residues of the binase active center are involved in formation of the complex; in type II complex, the active center remains free. It is suggested that temporary binding of free barstar into type II complex competes with the inhibition reaction. Presumably, this explains the decrease in the rate of binase inhibition by barstar as compared with the analogous reaction of barnase.     

Conformational dynamics of human α-fetoprotein-derived heptapeptide LDSYQCT analogs

Conformational dynamics of a biologically active fragment of alpha-fetoprotein, the heptapeptide LDSYQCT, and its analogs obtained by site-directed substitutions of amino acid residues were studied. The conformational dynamics of the peptide were conservative under the substitutions Y17F, Y17S, and D15E. Substitutions C19A and S16V resulted only in local changes in the dynamic behavior of the peptide. Chemical modification of cysteine (C19) or dimerization of the peptide by producing a disulfide bond between cysteine residues of two parallel peptide chains, as well as the substitutions C19G, C19S, Q18E, and D15N changed a set of possible conformations and dynamic behavior of all amino acid residues. The most significant changes were caused by substitution of uncharged amino acid residues by charged ones, and vice versa.     

Temperature adaptation of lactate dehydrogenase. Structural, functional and genetic aspects

Comparison of the primary structures of thermophilic, mesophilic and psychrophilic lactate dehydrogenase (LDH) reveals a multitude of temperature-related amino acid substitutions. In the substitutions amino acid residues occurring preferentially in thermophilic, mesophilic (psychrophilic) LDH were found. On this basis, amino acid residues could be classified in an order from typical thermophilic (thermostabilizing) to typical mesophilic (thermolabilizing, increasing dynamics of the enzyme molecule) residues. The temperature-dependent ratio between thermostabilizing and thermolabilizing amino acid residues forms the basis for the specific structural and functional properties of thermophilic or mesophilic LDH. It is interesting that there appears to be a relationship between this order from thermophilic to mesophilic amino acid residues and the type of bases coding for these individual residues in the translation step of protein biosynthesis. Temperature-related amino acid substitutions are based on temperature-related base substitutions. A possible mechanism of temperature adaptation of LDH through alternative selection of thermophilic and mesophilic amino acid residues at the level of tRNA (anticodon)-mRNA (codon) interactions is discussed. These temperature-adaptation processes are evolutionary events in which the evolution and structure of the genetic code are involved.     

Amino acid sequence of Ondatra myoglobin (Ondatra zibethica). Characteristic features of myoglobins from semiaquatic animals

Based on the amino acid composition of globin, amino acid analysis and N-terminal sequencing of peptides as well as a comparative analysis of the primary structure of beaver, coypu rat and otter myoglobins with the use of the fingerprinting technique, the amino acid sequence of the major component of ondatra myoglobin including 153 amino acid residues was reconstructed. The results of a comparative analysis of the primary structure of myoglobin and the peculiarities of the functional morphology of myoglobins from semi-aquatic animals and sperm whale and the role of amino acid substitutions in the spatial structure of ondatra myoglobin are discussed.     

Evaluation of the protein interfaces that form an intermolecular four-helix bundle as studied by computer simulation

Rational design of protein surface is important for creating higher order protein structures, but it is still challenging. In this study, we designed in silico the several binding interfaces on protein surfaces that allow a de novo protein–protein interaction to be formed. We used a computer simulation technique to find appropriate amino acid arrangements for the binding interface. The protein–protein interaction can be made by forming an intermolecular four-helix bundle structure, which is often found in naturally occurring protein subunit interfaces. As a model protein, we used a helical protein, YciF. Molecular dynamics simulation showed that a new protein–protein interaction is formed depending on the number of hydrophobic and charged amino acid residues present in the binding surfaces. However, too many hydrophobic amino acid residues present in the interface negatively affected on the binding. Finally, we found an appropriate arrangement of hydrophobic and charged amino acid residues that induces a protein–protein interaction through an intermolecular four-helix bundle formation.     

Relationship between structure and amino acid sequence of strongly twisted and coiled β-hairpins in globular proteins

β-Hairpins are widespread in proteins, and it is possible to find them both within β-sheets and separately. In this work, a comparative analysis of amino acid sequences of β-strands within strongly twisted β-hairpins from different structural protein subclasses has been conducted. Strongly twisted and coiled β-hairpin generates in the space a right double helix out of β-strands that are connected by a loop region (connections). The frequencies of amino acid residues on the internal (concave) and external (convex) surfaces of strongly twisted β-hairpins have been determined (220 β-hairpins from nonhomologous proteins were studied). The concave surface of these β-hairpins is mainly generated by hydrophobic residues, while the convex surface by hydrophilic residues; accordingly, the alternation of hydrophobic internal and hydrophilic external residues is observed in their amino acid sequences. Amino acid residues of glycine and alanine (especially in places of the largest twisting of the strands) were anomalously frequently found in internal positions of strongly twisted and coiled β-hairpins. It was established that internal positions never contain the proline residues, while external positions in the twisting region contain them in a relatively large amount. It was demonstrated that at least one amino acid residue in α_L- or ε-conformation is required for generation of relatively short (up to 7 amino acid residues) connection. As a rule, these positions are occupied by glycines. Thus, not only the alternation of hydrophobic and hydrophilic amino acid residues, but also the presence of one or two glycine residues in the connection region and the excess of glycines and alanines in the places of the largest strand twisting on the concave surface, as well as the presence of prolines on the convex surface, are required to generate a strongly twisted and coiled β-hairpin.     

Molecular dynamics study of femtosecond events in photoactive yellow protein after photoexcitation of the chromophore

Molecular dynamics simulations were carried out to study what happens in a photoreceptor protein, photoactive yellow protein (PYP), immediately after the vertical transition of the chromophore from the ground to the excited state. A photon absorption simulation was performed to investigate the movement of amino acid residues upon photoexcitation. To calculate the excited state of the chromophore, SCF-CI calculation was carried out with INDO/S Hamiltonian. We observed that some amino acid residues have strong interactions with the chromophore. Most of these amino acid residues are conserved in PYPs from three different species of bacteria. This observation indicates the biological importance of these residues. Proteins 32:268–275, 1998. © 1998 Wiley-Liss, Inc.     

A model of the secondary structure and distribution of antigenic determinants in the VP1 protein of hepatitis A virus

A comparative analysis of amino acid sequence of the proteins VP1 of hepatitis A virus and poliovirus of the 1 type was carried out. A model is proposed of structural organization of VP1 of hepatitis A virus providing the presence of a bilayer core formed by 8 antiparallel beta-strands. Probable candidates for surface antigenic determinants are the amino acid sequences located in unordered fragments of the polypeptide chain (residues 101-106 and 115-125), and alpha-helical region (residues 127-135).     

Disorder predictors also predict backbone dynamics for a family of disordered proteins

Several algorithms have been developed that use amino acid sequences to predict whether or not a protein or a region of a protein is disordered. These algorithms make accurate predictions for disordered regions that are 30 amino acids or longer, but it is unclear whether the predictions can be directly related to the backbone dynamics of individual amino acid residues. The nuclear Overhauser effect between the amide nitrogen and hydrogen (NHNOE) provides an unambiguous measure of backbone dynamics at single residue resolution and is an excellent tool for characterizing the dynamic behavior of disordered proteins. In this report, we show that the NHNOE values for several members of a family of disordered proteins are highly correlated with the output from three popular algorithms used to predict disordered regions from amino acid sequence. This is the first test between an experimental measure of residue specific backbone dynamics and disorder predictions. The results suggest that some disorder predictors can accurately estimate the backbone dynamics of individual amino acids in a long disordered region.     

Comparative specificity of microbial acid proteinases for synthetic peptides. II. Effect of secondary interaction

To investigate the effect of “secondary interaction” on hydrolysis by various acid proteinases from molds and yeasts, synthetic peptides

amino acid residues) were used as substrates. Pepsin was used for the comparative study. These peptides were split at the peptide bonds indicated by the arrows, permitting examination of the effect of residue X distant by two or three amino acid residues from the hydrolytic site in the peptides. According to the system of Schechter and Berger (Biochem. Biophys. Res. Commun. 27; 157, 1967), the amino acid residues in peptide substrates were numbered P₁, P₂, etc. toward the N-terminal direction from the site of hydrolysis, and P₁′, P₂′, etc. toward the C-terminal direction. The results indicated that hydrolysis by these microbial enzymes is affected by at least six amino acid residues (P₁-P₃ and P₁′-P₃′) in peptide substrates, as is seen with pepsin. Elongation of the peptide chain with suitable amino acid residues from P₁ to P₂ or P₃ and from P₁′ to P₂′ or P₃′ in peptide substrates resulted in much or less increase of hydrolysis depending upon the species of the enzyme producers.     

A simple statistical method for estimating type-II (cluster-specific) functional divergence of protein sequences

Predicting functional amino acid residues in silico is important for comparative genomics. In this paper, we focus on the issue of how to statistically identify cluster-specific amino acid residues that are related to the functional divergence after gene duplication. We approach this problem using a framework based on site-specific shift of amino acid property (type-II functional divergence), as opposed to site-specific shift of evolutionary rate (type-I functional divergence). An efficient statistical procedure is implemented to facilitate the development of phylogenomic database for cluster-specific residues of large-scale protein families. Our method has the following features: 1) statistical testing of the type-II functional divergence and 2) the site-specific Bayesian profile to measure how amino acid residues contribute to type-II (cluster-specific) functional divergence. Consequently, one may obtain the posterior probability for "functional" cluster-specific residues. Case studies are presented and indicate that radical cluster-specific residues are responsible for most of inferred type-II functional divergence, whereas conserved cluster-specific residues appear less than even those imperfect radical cluster-specific residues to this type of functional divergence.     

Biomimetic design of affinity peptide ligands for human IgG based on protein A-IgG complex

Staphylococcal protein A (SpA) has been widely used as an affinity ligand for the purification of immunoglobulin G (IgG). Based on the affinity motif of SpA, we have herein developed a biomimetic design strategy for affinity peptide ligands of IgG. First, according to the distribution of the six hot spots of the SpA affinity motif determined previously, the number of residues that should be inserted into between the hot spots was determined. Cysteine was introduced as one of the middle inserted residues of the peptide for later immobilization. Then, amino acid location was performed to identify other amino acid residues for insertion, leading to the construction of a peptide library. The library was screened by using different molecular simulation protocols, resulting in the selection of 15 peptide candidates. Thereafter, molecular dynamics simulations were performed to validate the dynamics of the affinity interactions between the candidates and IgG, and 14 of them were found to keep high affinities. Finally, the affinity and specificity of the top one ligand FYWHCLDE were exemplified by protein chromatography and IgG purification. The results indicate that the design strategy was successful and the affinity peptide ligand for IgG is promising for application in antibody purifications.     

Sequence evolution correlates with structural dynamics

Biochemical activity and core stability are essential properties of proteins, maintained usually by conserved amino acids. Structural dynamics emerged in recent years as another essential aspect of protein functionality. Structural dynamics enable the adaptation of the protein to binding substrates and to undergo allosteric transitions, while maintaining the native fold. Key residues that mediate structural dynamics would thus be expected to be conserved or exhibit coevolutionary patterns at least. Yet, the correlation between sequence evolution and structural dynamics is yet to be established. With recent advances in efficient characterization of structural dynamics, we are now in a position to perform a systematic analysis. In the present study, a set of 34 enzymes representing various folds and functional classes is analyzed using information theory and elastic network models. Our analysis shows that the structural regions distinguished by their coevolution propensity as well as high mobility are predisposed to serve as substrate recognition sites, whereas residues acting as global hinges during collective dynamics are often supported by conserved residues. We propose a mobility scale for different types of amino acids, which tends to vary inversely with amino acid conservation. Our findings suggest the balance between physical adaptability (enabled by structure-encoded motions) and chemical specificity (conferred by correlated amino acid substitutions) underlies the selection of a relatively small set of versatile folds by proteins.     

Structural and functional conservation profiles of novel cathepsin L-like proteins identified in the Drosophila melanogaster genome

Cathepsin L is a cysteine protease which degrades connective tissue proteins including collagen, elastin, and fibronectin. In this study, five well-characterized cathepsin L proteins from different arthropods were used as query sequences for the Drosophila genome database. The search yielded 10 cathepsin L-like sequences, of which eight putatively represent novel cathepsin L-like proteins. To understand the phylogenetic relationship among these cathepsin L-like proteins, a phylogenetic tree was constructed based on their sequences. In addition, models of the tertiary structures of cathepsin L were constructed using homology modeling methods and subjected to molecular dynamics simulations to obtain reasonable structure to understand its dynamical behavior. Our findings demonstrate that all of the potential Drosophila cathepsin L-like proteins contain at least one cathepsin propeptide inhibitor domain. Multiple sequence alignment and homology models clearly highlight the conservation of active site residues, disulfide bonds, and amino acid residues critical for inhibitor binding. Furthermore, comparative modeling indicates that the sequence/structure/function profiles and active site architectures are conserved.     

Cloning and sequencing of the gene for the cytoplasmic inorganic pyrophosphatase from the thermoacidophilic archaebacterium Thermoplasma acidophilum.

The gene (ppa) from the thermoacidophilic archaebacterium Thermoplasma acidophilum, encoding the cytoplasmic pyrophosphatase, has been cloned. Two degenerate oligonucleotide probes, synthesized according to the N-terminal amino acid sequence of the isolated protein, were used to screen subgenomic libraries. The DNA-derived amino acid sequence of the archaebacterial enzyme allows, for the first time, comparative studies of cytoplasmic pyrophosphatases to be extended to all three urkingdoms. The archaebacterial pyrophosphatase more closely resembles the eubacterial enzymes on the basis of sequence similarity and subunit size. The majority of amino acid residues considered to be essential for hydrolysis of pyrophosphate seem to have been conserved throughout evolution, as inferred from the results of an alignment of sequences from all three urkingdoms.     

Classification of Bacteria Based on the Biases of Terminal Amino Acid Residues

The frequencies of amino acid residues are known to be biased at both terminal regions of amino acid sequences deduced from bacterial genomic DNA. To investigate whether or not the features of biases of amino acid residues at the terminal regions are related to the bacterial phylogeny, we calculated the normalized amino acid compositions at both terminal regions, and used these compositions to classify 144 bacteria by hierarchical clustering analysis. Our results showed that most of these bacteria were classified into taxonomic classes by the hierarchical clustering analysis that was based on the normalized amino acid compositions at the N-terminal region. Therefore, we concluded that the features of biases of the N-terminal amino acid residues were related to the bacterial phylogeny.     

Molecular dynamics of zervamicin II and its analogs in water and methanol

A comparative study has been made of the molecular dynamics of zervamicin II (an antimicrobial peptide of the peptaibol group, which has channel-forming activity) in water and methanol. The influence of amino acid substitutions on the dynamics and stability of the peptide structure has been investigated. The amino acid sequence responsible for the absence of swivel motions in short peptaibols has been determined.