The design of a peptide sequence to inhibit HIV replication: a search algorithm combining Monte Carlo and self-consistent mean field techniques

We developed a search algorithm combining Monte Carlo (MC) and self-consistent mean field techniques to evolve a peptide sequence that has good binding capability to the anticodon stem and loop (ASL) of human lysine tRNA species, tRNA^Lys3, with the ultimate purpose of breaking the replication cycle of human immunodeficiency virus-1. The starting point is the 15-amino-acid sequence, RVTHHAFLGAHRTVG, found experimentally by Agris and co-workers to bind selectively to hypermodified tRNA^Lys3. The peptide backbone conformation is determined via atomistic simulation of the peptide-ASL^Lys3 complex and then held fixed throughout the search. The proportion of amino acids of various types (hydrophobic, polar, charged, etc.) is varied to mimic different peptide hydration properties. Three different sets of hydration properties were examined in the search algorithm to see how this affects evolution to the best-binding peptide sequences. Certain amino acids are commonly found at fixed sites for all three hydration states, some necessary for binding affinity and some necessary for binding specificity. Analysis of the binding structure and the various contributions to the binding energy shows that: 1) two hydrophilic residues (asparagine at site 11 and the cysteine at site 12) “recognize” the ASL^Lys3 due to the VDW energy, and thereby contribute to its binding specificity and 2) the positively charged arginines at sites 4 and 13 preferentially attract the negatively charged sugar rings and the phosphate linkages, and thereby contribute to the binding affinity.     

Balancing exploration and exploitation in population‐based sampling improves fragment‐based de novo protein structure prediction

Conformational search space exploration remains a major bottleneck for protein structure prediction methods. Population‐based meta‐heuristics typically enable the possibility to control the search dynamics and to tune the balance between local energy minimization and search space exploration. EdaFold is a fragment‐based approach that can guide search by periodically updating the probability distribution over the fragment libraries used during model assembly. We implement the EdaFold algorithm as a Rosetta protocol and provide two different probability update policies: a cluster‐based variation (EdaRose_c) and an energy‐based one (EdaRose_en). We analyze the search dynamics of our new Rosetta protocols and show that EdaRose_c is able to provide predictions with lower C RMSD to the native structure than EdaRose_en and Rosetta AbInitio Relax protocol. Our software is freely available as a C++ patch for the Rosetta suite and can be downloaded from http://www.riken.jp/zhangiru/software/ . Our protocols can easily be extended in order to create alternative probability update policies and generate new search dynamics. Proteins 2017; 85:852–858. © 2016 Wiley Periodicals, Inc.     

Specificity of Phage Display Selected Peptides for Modified Anticodon Stem and Loop Domains of tRNA

Protein recognition of RNA has been studied using Peptide Phage Display Libraries, but in the absence of RNA modifications. Peptides from two libraries, selected for binding the modified anticodon stem and loop (ASL) of human tRNA^Lys3 having 2-thiouridine (s²U₃₄) and pseudouridine (Ψ₃₉), bound the modified human ASL^Lys3(s²U₃₄;Ψ₃₉) preferentially and had significant homology with RNA binding proteins. Selected peptides were narrowed to a manageable number using a less sensitive, but inexpensive assay before conducting intensive characterization. The affinity and specificity of the best binding peptide (with an N-terminal fluorescein) were characterized by fluorescence spectrophotometry. The peptide exhibited the highest binding affinity for ASL^Lys3(s²U₃₄;Ψ₃₉), followed by the hypermodified ASL^Lys3 (mcm⁵s²U₃₄;ms²t⁶A₃₇) and the unmodified ASL^Lys3, but bound poorly to singly modified ASL^Lys3 constructs (Ψ₃₉, ms²t⁶A_37, s²U₃₄), ASL^Lys1,2 (t⁶A₃₇) and Escherichia coli ASL^Glu (s²U₃₄). Thus, RNA modifications are potentially important recognition elements for proteins and can be targets for selective recognition by peptides.     

Investigating the linkage between disease‐causing amino acid variants and their effect on protein stability and binding

Single amino acid variations (SAV) occurring in human population result in natural differences between individuals or cause diseases. It is well understood that the molecular effect of SAV can be manifested as changes of the wild type characteristics of the corresponding protein, among which are the protein stability and protein interactions. Typically the effect of SAV on protein stability and interactions was assessed via the changes of the wild type folding and binding free energies. However, in terms of SAV affecting protein functionally and disease susceptibility, one wants to know to what extend the wild type function is perturbed by the SAV. Here it is demonstrated that relative, rather than the absolute, change of the folding and binding free energy serves as a good indicator for SAV association with disease. Using HumVar as a source for disease‐causing SAV and experimentally determined free energy changes from ProTherm and SKEMPI databases, correlation coefficients (CC) between the disease index and relative folding and binding probability indexes, respectively, was achieved. The obtained CCs demonstrated the applicability of the proposed approach and it served as good indicator for SAV association with disease. Proteins 2016; 84:232–239. © 2015 Wiley Periodicals, Inc.     

Archaea S‐layer nanotube from a “black smoker” in complex with cyclo‐octasulfur (S8) rings

Elemental sulfur exists primarily as an ring and serves as terminal electron acceptor for a variety of sulfur‐fermenting bacteria. Hyperthermophilic archaea from black smoker vents are an exciting research tool to advance our knowledge of sulfur respiration under extreme conditions. Here, we use a hybrid method approach to demonstrate that the proteinaceous cavities of the S‐layer nanotube of the hyperthermophilic archaeon Staphylothermus marinus act as a storage reservoir for cyclo‐octasulfur . Fully atomistic molecular dynamics (MD) simulations were performed and the method of multiconfigurational thermodynamic integration was employed to compute the absolute free energy for transferring a ring of elemental sulfur from an aqueous bath into the largest hydrophobic cavity of a fragment of archaeal tetrabrachion. Comparisons with earlier MD studies of the free energy of hydration as a function of water occupancy in the same cavity of archaeal tetrabrachion show that the sulfur ring is energetically favored over water.     

pKa predictions for proteins,RNAs, and DNAs with the Gaussian dielectric function using DelPhi pKa

We developed a Poisson‐Boltzmann based approach to calculate the values of protein ionizable residues (Glu, Asp, His, Lys and Arg), nucleotides of RNA and single stranded DNA. Two novel features were utilized: the dielectric properties of the macromolecules and water phase were modeled via the smooth Gaussian‐based dielectric function in DelPhi and the corresponding electrostatic energies were calculated without defining the molecular surface. We tested the algorithm by calculating values for more than 300 residues from 32 proteins from the PPD dataset and achieved an overall RMSD of 0.77. Particularly, the RMSD of 0.55 was achieved for surface residues, while the RMSD of 1.1 for buried residues. The approach was also found capable of capturing the large shifts of various single point mutations in staphylococcal nuclease (SNase) from ‐cooperative dataset, resulting in an overall RMSD of 1.6 for this set of pKa's. Investigations showed that predictions for most of buried mutant residues of SNase could be improved by using higher dielectric constant values. Furthermore, an option to generate different hydrogen positions also improves predictions for buried carboxyl residues. Finally, the calculations on two RNAs demonstrated the capability of this approach for other types of biomolecules. Proteins 2015; 83:2186–2197. © 2015 Wiley Periodicals, Inc.     

Exploration of conformational changes in lactose permease upon sugar binding and proton transfer through coarse‐grained simulations

Escherichia coli lactose permease (LacY) actively transports lactose and other galactosides across cell membranes through lactose/H⁺ symport process. Lactose/H⁺ symport is a highly complex process that involves sugar translocation, H⁺ transfer, and large‐scale protein conformational changes. The complete picture of lactose/H⁺ symport is largely unclear due to the complexity and multiscale nature of the process. In this work, we develop the force field for sugar molecules compatible with PACE, a hybrid and coarse‐grained force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid. After validation, we implement the new force field to investigate the binding of a ‐d ‐galactopyranosyl‐1‐thio‐ ‐d ‐galactopyranoside (TDG) molecule to a wild‐type LacY. Results show that the local interactions between TDG and LacY at the binding pocket are consistent with the X‐ray experiment. Transitions from inward‐facing to outward‐facing conformations upon TDG binding and protonation of Glu269 have been achieved from ～5.5 µs simulations. Both the opening of the periplasmic side and closure of the cytoplasmic side of LacY are consistent with double electron–electron resonance and thiol cross‐linking experiments. Our analysis suggests that the conformational changes of LacY are a cumulative consequence of interdomain H‐bonds breaking at the periplasmic side, interdomain salt‐bridge formation at the cytoplasmic side, and the TDG orientational changes during the transition.     

Inferring the microscopic surface energy of protein–protein interfaces from mutation data

Mutations at protein–protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical values, yielding mean energies of ?48 cal mol^?1 Å^?2 for interactions between hydrophobic surfaces, ?51 to ?80 cal mol^?1 Å^?2 for surfaces of complementary charge, and 66–83 cal mol^?1 Å^?2 for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at ?24 cal mol^?1 Å^?2. Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid‐body interactions (r = 0.66). When used to rerank docking poses, it can place near‐native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50–95% of the cohesive energy. The model is available open‐source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/ . Proteins 2015; 83:640–650. © 2015 Wiley Periodicals, Inc.     

Ratio of ellipticities between 192 and 208 nm (R1): An effective electronic circular dichroism parameter for characterization of the helical components of proteins and peptides

Circular dichroism (CD) spectroscopy represents an important tool for characterization of the peptide and protein secondary structures that mainly arise from the conformational disposition of the peptide backbone in solution. In 1991 Manning and Woody proposed that, in addition to the signal intensity, the ratio between and ((R₂ ) ? [θ]₂₂₂/[θ]₂₀₈), along with and ((R₁ ) ? [θ]₁₉₂/[θ]₂₀₈), may be utilized towards identifying the peptide/protein conformation (especially 3₁₀‐ and α‐helices). However, till date the use of the ratiometric ellipticity component for helical structure analysis of peptides and proteins has not been reported. We studied a series of temperature dependent CD spectra of a thermally stable, model helical peptide and its related analogs in water as a function of added 2,2,2‐trifluoroethanol (TFE) in order to explore their landscape of helicity. For the first time, we have experimentally shown here that the R₁ parameter can characterize better the individual helices, while the other parameter R₂ and the signal intensity do not always converge. We emphasize the use of the R₁ ratio of ellipticities for helical characterization because of the common origin of these two bands (exciton splitting of the amide π→ π* transition in a helical polypeptide). This approach may become worthwhile and timely with the increasing accessibility of CD synchrotron sources.     

Nucleotide‐dependent structural fluctuations and regulation of microtubule‐binding affinity of KIF1A

Molecular motors such as kinesin regulate affinity to a rail protein during the ATP hydrolysis cycle. The regulation mechanism, however, is yet to be determined. To understand this mechanism, we investigated the structural fluctuations of the motor head of the single‐headed kinesin called KIF1A in different nucleotide states using molecular dynamics simulations of a Gō‐like model. We found that the helix at the microtubule (MT) binding site intermittently exhibits a large structural fluctuation when MT is absent. Frequency of this fluctuation changes systematically according to the nucleotide states and correlates strongly with the experimentally observed binding affinity to MT. We also showed that thermal fluctuation enhances the correlation and the interaction with the nucleotide suppresses the fluctuation of the helix . These results suggest that KIF1A regulates affinity to MT by changing the flexibility of the helix during the ATP hydrolysis process: the binding site becomes more flexible in the strong binding state than in the weak binding state. Proteins 2015; 83:809–819. © 2015 Wiley Periodicals, Inc.     

Rate limiting step of the allosteric activation of the bacterial adhesin FimH investigated by molecular dynamics simulations

The bacterial adhesin FimH is a model for the study of protein allostery because its structure has been resolved in multiple configurations, including the active and the inactive state. FimH consists of a pilin domain (PD) that anchors it to the rest of the fimbria and an allosterically regulated lectin domain (LD) that binds mannose on the surface of infected cells. Under normal conditions, the two domains are docked to each other and LD binds mannose weakly. However, in the presence of tensile force generated by shear the domains separate and conformational changes propagate across LD resulting in a stronger bond to mannose. Recently, the crystallographic structure of a variant of FimH has been resolved, called ${\text{FimH}}^{\text{FocH}}$, where PD contains 10 mutations near the inter-domain interface. Although the X-ray structures of FimH and ${\text{FimH}}^{\text{FocH}}$ are almost identical, experimental evidence shows that ${\text{FimH}}^{\text{FocH}}$ is activated even in the absence of shear. Here, molecular dynamics simulations combined with the Jarzynski equality were used to investigate the discrepancy between the crystallographic structures and the functional assays. The results indicate that the free energy barrier of the unbinding process between LD and PD is drastically reduced in ${\text{FimH}}^{\text{FocH}}$. Rupture of inter-domain hydrogen bonds involving R166 constitutes a rate limiting step of the domain separation process and occurs more readily in ${\text{FimH}}^{\text{FocH}}$ than FimH. In conclusion, the mutations in ${\text{FimH}}^{\text{FocH}}$ shift the equilibrium toward an equal occupancy of bound and unbound states for LD and PD by reducing a rate limiting step.     

Structure of iron saturated C‐lobe of bovine lactoferrin at pH 6.8 indicates a weakening of iron coordination

The bilobal lactoferrin is an approximately 76 kDa glycoprotein. It sequesters two Fe³⁺ ions together with two ions. The C‐terminal half (residues, Tyr342–Arg689, C‐lobe) of bovine lactoferrin (BLF) (residues Ala1–Arg689) was prepared by limited proteolysis using trypsin. Both C‐lobe and intact BLF were saturated to 100%. Both of them retained up to nearly 85% of iron at pH 6.5. At pH 5.0, C‐lobe retained 75% of iron whereas intact protein could retain only slightly more than 60%. At pH 4.0 both contained 25% iron and at pH 2.0 they were left with iron concentration of only 10%. The structure of iron saturated C‐lobe was determined at 2.79 Å resolution and refined to R_cryst and R_free factors of 0.205 and 0.273, respectively. The structure contains two crystallographically independent molecules, A and B. They were found to have identical structures with an r.m.s. shift of 0.5 Å for their C^α atoms. A high solvent content of 66% was observed in the crystals. The average value of an overall B‐factor was 68.0 Ų. The distance of 2.9 Å observed for the coordination bond between Fe³⁺ ion and N^e2 of His595 appeared to be considerably longer than the normally observed values of 1.9–2.2 Å. This indicated that the coordination bond involving His595 may be absent. Other coordination distances were observed in the range of 2.1–2.3 Å. Based on the present structure of iron saturated C‐lobe, it may be stated that His595 is the first residue to dissociate from ferric ion when the pH is lowered. Proteins 2016; 84:591–599. © 2016 Wiley Periodicals, Inc.     

Spectroscopic investigations of Er(3+) :CdO-Bi(2) O(3) -B(2) O(3) glasses

This article reports on the optical properties of Er³⁺ ions doped CdO–Bi₂O₃–B₂O₃ (CdBiB) glasses. The materials were characterized by optical absorption and emission spectra. By using Judd–Ofelt theory, the intensity parameters Ω_λ (λ = 2, 4, 6) and also oscillatory strengths were calculated from the absorption spectra. The results were used to compute the radiative properties of Er³⁺:CdBiB glasses. The concentration quenching and energy transfer from Yb³⁺–Er³⁺ were explained. The stimulated emission cross‐section, full width at half maximum (FWHM) and FWHM × values are also calculated for all the Er³⁺:CdBiB glasses. Copyright © 2011 John Wiley & Sons, Ltd.     

A model of the plant-to-plant movement of aphids I. Description of the model

The plant-to-plant movement of insects in one of the factors determining the distribution of individuals in insect populations. In this report the movement of barley aphids was analyzed by a statistical model. The model is represented as the convolution of three probability functions:

1. The probability that s individuals are found on a plant at time t₀:Q(s);
2. The probability that i individuals leave the plant and remain on the ground from time t₀ to t₁:_sC_ipⁱq^s−i and p+q=1, where p and q are the proportions of individuals which do not leave a plant and which leave it once or more, respectively;
3. The probability that j individuals climb a plant between time t₀ to t₁ and stay there at time t₁:e^−λλ^j/λ^!, where λ is the mean of the individuals.

The probability that l individuals are located on a plant at time t₁ is represented by the following equation It was shown by simple experiments that the experimental populations were well fitted to the model.     

Role of active‐site residues Tyr55 and Tyr114 in catalysis and substrate specificity of Corynebacterium diphtheriae C‐S lyase

In recent years, there has been increased interest in bacterial methionine biosynthesis enzymes as antimicrobial targets because of their pivotal role in cell metabolism. C‐S lyase from Corynebacterium diphtheriae is a pyridoxal 5′‐phosphate‐dependent enzyme in the transsulfuration pathway that catalyzes the α,β‐elimination of sulfur‐containing amino acids, such as l ‐cystathionine, to generate ammonia, pyruvate, and homocysteine, the immediate precursor of L ‐methionine. In order to gain deeper insight into the functional and dynamic properties of the enzyme, mutants of two highly conserved active‐site residues, Y55F and Y114F, were characterized by UV‐visible absorbance, fluorescence, and CD spectroscopy in the absence and presence of substrates and substrate analogs, as well as by steady‐state kinetic studies. Substitution of Tyr55 with Phe apparently causes a 130‐fold decrease in at pH 8.5 providing evidence that Tyr55 plays a role in cofactor binding. Moreover, spectral data show that the mutant accumulates the external aldimine intermediate suggesting that the absence of interaction between the hydroxyl moiety and PLP‐binding residue Lys222 causes a decrease in the rate of substrate deprotonation. Mutation of Tyr114 with Phe slightly influences hydrolysis of l ‐cystathionine, and causes a change in substrate specificity towards l ‐serine and O‐acetyl‐l ‐serine compared to the wild type enzyme. These findings, together with computational data, provide useful insights in the substrate specificity of C‐S lyase, which seems to be regulated by active‐site architecture and by the specific conformation in which substrates are bound, and will aid in development of inhibitors. Proteins 2015; 83:78–90. © 2014 Wiley Periodicals, Inc.     

Thermodynamics of Aβ16–21 dissociation from a fibril: Enthalpy,entropy, and volumetric properties

Here, we provide insights into the thermodynamic properties of A dissociation from an amyloid fibril using all‐atom molecular dynamics simulations in explicit water. An umbrella sampling protocol is used to compute potentials of mean force (PMF) as a function of the distance ξ between centers‐of‐mass of the A peptide and the preformed fibril at nine temperatures. Changes in the enthalpy and the entropic energy are determined from the temperature dependence of these PMF(s) and the average volume of the simulation box is computed as a function of ξ. We find that the PMF at 310 K is dominated by enthalpy while the entropic energy does not change significantly during dissociation. The volume of the system decreases during dissociation. Moreover, the magnitude of this volume change also decreases with increasing temperature. By defining dock and lock states using the solvent accessible surface area (SASA), we find that the behavior of the electrostatic energy is different in these two states. It increases (unfavorable) and decreases (favorable) during dissociation in lock and dock states, respectively, while the energy due to Lennard‐Jones interactions increases continuously in these states. Our simulations also highlight the importance of hydrophobic interactions in accounting for the stability of A . Proteins 2015; 83:1963–1972. © 2015 Wiley Periodicals, Inc.