A solution structure analysis reveals a bent collagen triple helix in the complement activation recognition molecule mannan-binding lectin.

Collagen triple helices are critical in the function of mannan-binding lectin (MBL), an oligomeric recognition molecule in complement activation. The MBL collagen regions form complexes with the serine proteases MASP-1 and MASP-2 in order to activate complement, and mutations lead to common immunodeficiencies. To evaluate their structure-function properties, we studied the solution structures of four MBL-like collagen peptides. The thermal stability of the MBL collagen region was much reduced by the presence of a GQG interruption in the typical (X-Y-Gly)n repeat compared to controls. Experimental solution structural data were collected using analytical ultracentrifugation and small angle X-ray and neutron scattering. As controls, we included two standard Pro-Hyp-Gly collagen peptides (POG)10-13, as well as three more peptides with diverse (X-Y-Gly)n sequences that represented other collagen features. These data were quantitatively compared with atomistic linear collagen models derived from crystal structures and 12,000 conformations obtained from molecular dynamics simulations. All four MBL peptides were bent to varying degrees up to 85o in the best-fit molecular dynamics models. The best-fit benchmark peptides (POG)n were more linear but exhibited a degree of conformational flexibility. The remaining three peptides showed mostly linear solution structures. In conclusion, the collagen helix is not strictly linear, the degree of flexibility in the triple helix depends on its sequence, and the triple helix with the GQG interruption showed a pronounced bend. The bend in MBL GQG peptides resembles the bend in the collagen of complement C1q and may be key for lectin pathway activation.

Sequence-independent recognition of the amyloid structural motif by GFP protein family.

Cnidarian fluorescent protein (FP) derivatives such as GFP, mCherry, and mEOS2 have been widely used to monitor gene expression and protein localization through biological imaging because they are considered functionally inert. We demonstrate that FPs specifically bind amyloid fibrils formed from many natural peptides and proteins. FPs do not bind other nonamyloid fibrillar structures such as microtubules or actin filaments and do not bind to amorphous aggregates. FPs can also bind small aggregates formed during the lag phase and early elongation phase of fibril formation and can inhibit amyloid fibril formation in a dose-dependent manner. These findings suggest caution should be taken in interpreting FP-fusion protein localization data when amyloid structures may be present. Given the pathological significance of amyloid-related species in some diseases, detection and inhibition of amyloid fibril formation using FPs can provide insights on developing diagnostic tools.

The volume changes of unfolding of dsDNA.

High hydrostatic pressure can have profound effects on the stability of biomacromolecules. The magnitude and direction (stabilizing or destabilizing) of this effect is defined by the volume changes in the system, ΔV. Positive volume changes will stabilize the starting native state, whereas negative volume changes will lead to the stabilization of the final unfolded state. For the DNA double helix, experimental data suggested that when the thermostability of dsDNA is below 50°C, increase in hydrostatic pressure will lead to destabilization; i.e., helix-to-coil transition has negative ΔV. In contrast, the dsDNA sequences with the thermostability above 50°C showed positive ΔV values and were stabilized by hydrostatic pressure. In order to get insight into this switch in the response of dsDNA to hydrostatic pressure as a function of temperature, first we further validated this trend using experimental measurements of ΔV for 10 different dsDNA sequences using pressure perturbation calorimetry. We also developed a computational protocol to calculate the expected volume changes of dsDNA unfolding, which was benchmarked against the experimental set of 50 ΔV values that included, in addition to our data, the values from the literature. Computation predicts well the experimental values of ΔV. Such agreement between computation and experiment lends credibility to the computation protocol and provides molecular level rational for the observed temperature dependence of ΔV that can be traced to the hydration. Difference in the ΔV value for A/T versus G/C basepairs is also discussed.

Effects of Hydrostatic Pressure on the Thermodynamics of CspB-Bs Interactions with the ssDNA Template.

Understanding the thermodynamic mechanisms of adaptation of biomacromolecules to high hydrostatic pressure can help shed light on how piezophilic organisms can survive at pressures reaching over 1000 atmospheres. Interaction of proteins with nucleic acids is one of the central processes that allow information flow encoded in the sequence of DNA. Here, we report the results of a study on the interaction of cold shock protein B from Bacillus subtilis (CspB-Bs) with heptadeoxythymine template (pDT7) as a function of temperature and hydrostatic pressure. Experimental data collected at different CspB-Bs:pDT7 ratios were analyzed using a thermodynamic linkage model that accounts for both protein unfolding and CspB-Bs:pDT7 binding. The global fit to the model provided estimates of the stability of CspB-Bs, ΔGProto, the volume change upon CspB-Bs unfolding, ΔVProt, the association constant for CspB-Bs:pDT7 complex, Kao, and the volume changes upon pDT7 single-stranded DNA (ssDNA) template binding, ΔVBind. The protein, CspB-Bs, unfolds with an increase in hydrostatic pressure (ΔVProt < 0). Surprisingly, our study showed that ΔVBind < 0, which means that the binding of CspB-Bs to ssDNA is stabilized by an increase in hydrostatic pressure. Thus, CspB-Bs binding to pDT7 represents a case of linked equilibrium in which folding and binding react differently upon an increase in hydrostatic pressure: protein folding/unfolding equilibrium favors the unfolded state, while protein-ligand binding equilibrium favors the bound state. These opposing effects set a "maximum attainable" pressure tolerance to the protein-ssDNA complex under given conditions.

Protein adaptation to high hydrostatic pressure: Computational analysis of the structural proteome.

Hydrostatic pressure has a vital role in the biological adaptation of the piezophiles, organisms that live under high hydrostatic pressure. However, the mechanisms by which piezophiles are able to adapt their proteins to high hydrostatic pressure is not well understood. One proposed hypothesis is that the volume changes of unfolding (ΔVTot ) for proteins from piezophiles is distinct from those of nonpiezophilic organisms. Since ΔVTot defines pressure dependence of stability, we performed a comprehensive computational analysis of this property for proteins from piezophilic and nonpiezophilic organisms. In addition, we experimentally measured the ΔVTot of acylphosphatases and thioredoxins belonging to piezophilic and nonpiezophilic organisms. Based on this analysis we concluded that there is no difference in ΔVTot for proteins from piezophilic and nonpiezophilic organisms. Finally, we put forward the hypothesis that increased concentrations of osmolytes can provide a systemic increase in pressure stability of proteins from piezophilic organisms and provide experimental thermodynamic evidence in support of this hypothesis.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes.

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Strategy for Stabilization of CutA1 Proteins Due to Ion-Ion Interactions at Temperatures of over 100 °C.

In order to elucidate the contribution of charged residues to protein stabilization at temperatures of over 100 °C, we constructed many mutants of the CutA1 protein ( EcCutA1) from Escherichia coli. The goal was to see if one can achieve the same stability as for a CutA1 from hyperthermophile Pyrococcus horikoshii that has the denaturation temperature near 150 °C. The hydrophobic mutant of EcCutA1 ( Ec0VV) with denaturation temperature ( Td) of 113.2 °C was used as a template for mutations. The highest Td of Ec0VV mutants substituted by a single charged residue was 118.4 °C. Multiple ion mutants were also constructed by combination of single mutants and found to have an increased thermostability. The highest stability of multiple mutants was a mutant substituted by nine charged residues that had a Td of 142.2 °C. To evaluate the energy of ion-ion interactions of mutant proteins, we used the structural ensemble obtained by a molecular dynamics simulation at 300 K. The Td of ionic mutants linearly increases with the increments of the computed energy of ion-ion interactions for ionic mutant proteins even up to the temperatures near 140 °C, suggesting that ion-ion interactions cumulatively contribute to the stabilization of a protein at high temperatures.

Predictive Model of Linear Antimicrobial Peptides Active against Gram-Negative Bacteria.

Antimicrobial peptides (AMPs) have been identified as a potential new class of anti-infectives for drug development. There are a lot of computational methods that try to predict AMPs. Most of them can only predict if a peptide will show any antimicrobial potency, but to the best of our knowledge, there are no tools which can predict antimicrobial potency against particular strains. Here we present a predictive model of linear AMPs being active against particular Gram-negative strains relying on a semi-supervised machine-learning approach with a density-based clustering algorithm. The algorithm can well distinguish peptides active against particular strains from others which may also be active but not against the considered strain. The available AMP prediction tools cannot carry out this task. The prediction tool based on the algorithm suggested herein is available on https://dbaasp.org.

Non-linearity of the collagen triple helix in solution and implications for collagen function.

Collagen adopts a characteristic supercoiled triple helical conformation which requires a repeating (Xaa-Yaa-Gly)n sequence. Despite the abundance of collagen, a combined experimental and atomistic modelling approach has not so far quantitated the degree of flexibility seen experimentally in the solution structures of collagen triple helices. To address this question, we report an experimental study on the flexibility of varying lengths of collagen triple helical peptides, composed of six, eight, ten and twelve repeats of the most stable Pro-Hyp-Gly (POG) units. In addition, one unblocked peptide, (POG)10unblocked, was compared with the blocked (POG)10 as a control for the significance of end effects. Complementary analytical ultracentrifugation and synchrotron small angle X-ray scattering data showed that the conformations of the longer triple helical peptides were not well explained by a linear structure derived from crystallography. To interpret these data, molecular dynamics simulations were used to generate 50 000 physically realistic collagen structures for each of the helices. These structures were fitted against their respective scattering data to reveal the best fitting structures from this large ensemble of possible helix structures. This curve fitting confirmed a small degree of non-linearity to exist in these best fit triple helices, with the degree of bending approximated as 4-17° from linearity. Our results open the way for further studies of other collagen triple helices with different sequences and stabilities in order to clarify the role of molecular rigidity and flexibility in collagen extracellular and immune function and disease.

Modulation of folding energy landscape by charge-charge interactions: linking experiments with computational modeling.

The kinetics of folding-unfolding of a structurally diverse set of four proteins optimized for thermodynamic stability by rational redesign of surface charge-charge interactions is characterized experimentally. The folding rates are faster for designed variants compared with their wild-type proteins, whereas the unfolding rates are largely unaffected. A simple structure-based computational model, which incorporates the Debye-Hückel formalism for the electrostatics, was used and found to qualitatively recapitulate the experimental results. Analysis of the energy landscapes of the designed versus wild-type proteins indicates the differences in refolding rates may be correlated with the degree of frustration of their respective energy landscapes. Our simulations indicate that naturally occurring wild-type proteins have frustrated folding landscapes due to the surface electrostatics. Optimization of the surface electrostatics seems to remove some of that frustration, leading to enhanced formation of native-like contacts in the transition-state ensembles (TSE) and providing a less frustrated energy landscape between the unfolded and TS ensembles. Macroscopically, this results in faster folding rates. Furthermore, analyses of pairwise distances and radii of gyration suggest that the less frustrated energy landscapes for optimized variants are a result of more compact unfolded and TS ensembles. These findings from our modeling demonstrates that this simple model may be used to: (i) gain a detailed understanding of charge-charge interactions and their effects on modulating the energy landscape of protein folding and (ii) qualitatively predict the kinetic behavior of protein surface electrostatic interactions.

Putting the Piezolyte Hypothesis under Pressure.

A group of small molecules that stabilize proteins against high hydrostatic pressure has been classified as piezolytes, a subset of stabilizing cosolutes. This distinction would imply that piezolytes counteract the effects of high hydrostatic pressure through effects on the volumetric properties of the protein. The purpose of this study was to determine if cosolutes proposed to be piezolytes have an effect on the volumetric properties of proteins through direct experimental measurements of volume changes upon unfolding of model proteins lysozyme and ribonuclease A, in solutions containing varying cosolute concentrations. Solutions containing the proposed piezolytes glutamate, sarcosine, and betaine were used, as well as solutions containing the denaturants guanidinium hydrochloride and urea. Changes in thermostability were monitored using differential scanning calorimetry whereas changes in volume were monitored using pressure perturbation calorimetry. Our findings indicate that increasing stabilizing cosolute concentration increases the stability and transition temperature of the protein, but does not change the temperature dependence of volume changes upon unfolding. The results suggest that the pressure stability of a protein in solution is not directly affected by the presence of these proposed piezolytes, and so they cannot be granted this distinction.

Gallic Acid Is an Antagonist of Semen Amyloid Fibrils That Enhance HIV-1 Infection.

Recent in vitro studies have demonstrated that amyloid fibrils found in semen from healthy and HIV-infected men, as well as semen itself, can markedly enhance HIV infection rates. Semen fibrils are made up of multiple naturally occurring peptide fragments derived from semen. The best characterized of these fibrils are SEVI (semen-derived enhancer of viral infection), made up of residues 248-286 of prostatic acidic phosphatase, and the SEM1 fibrils, made up of residues 86-107 of semenogelin 1. A small molecule screen for antagonists of semen fibrils identified four compounds that lowered semen-mediated enhancement of HIV-1 infectivity. One of the four, gallic acid, was previously reported to antagonize other amyloids and to exert anti-inflammatory effects. To better understand the mechanism by which gallic acid modifies the properties of semen amyloids, we performed biophysical measurements (atomic force microscopy, electron microscopy, confocal microscopy, thioflavin T and Congo Red fluorescence assays, zeta potential measurements) and quantitative assays on the effects of gallic acid on semen-mediated enhancement of HIV infection and inflammation. Our results demonstrate that gallic acid binds to both SEVI and SEM1 fibrils and modifies their surface electrostatics to render them less cationic. In addition, gallic acid decreased semen-mediated enhancement of HIV infection but did not decrease the inflammatory response induced by semen. Together, these observations identify gallic acid as a non-polyanionic compound that inhibits semen-mediated enhancement of HIV infection and suggest the potential utility of incorporating gallic acid into a multicomponent microbicide targeting both the HIV virus and host components that promote viral infection.

Applications of pressure perturbation calorimetry to study factors contributing to the volume changes upon protein unfolding.

BACKGROUND: Pressure perturbation calorimetry (PPC) is a biophysical method that allows direct determination of the volume changes upon conformational transitions in macromolecules. SCOPE OF THIS REVIEW: This review provides novel details of the use of PPC to analyze unfolding transitions in proteins. The emphasis is made on the data analysis as well as on the validation of different structural factors that define the volume changes upon unfolding. Four case studies are presented that show the application of these concepts to various protein systems. MAJOR CONCLUSIONS: The major conclusions are: 1. Knowledge of the thermodynamic parameters for heat induced unfolding facilitates the analysis of the PPC profiles. 2. The changes in the thermal expansion coefficient upon unfolding appear to be temperature dependent.3.Substitutions on the protein surface have negligible effects on the volume changes upon protein unfolding. 4. Structural plasticity of proteins defines the position dependent effect of amino acid substitutions of the residues buried in the native state. 5. Small proteins have positive volume changes upon unfolding which suggests difference in balance between the cavity/void volume in the native state and the hydration volume changes upon unfolding as compared to the large proteins that have negative volume changes. GENERAL SIGNIFICANCE: The information provided here gives a better understanding and deeper insight into the role played by various factors in defining the volume changes upon protein unfolding.

Linking computation and experiments to study the role of charge-charge interactions in protein folding and stability.

Over the past two decades there has been an increase in appreciation for the role of surface charge-charge interactions in protein folding and stability. The perception shifted from the belief that charge-charge interactions are not important for protein folding and stability to the near quantitative understanding of how these interactions shape the folding energy landscape. This led to the ability of computational approaches to rationally redesign surface charge-charge interactions to modulate thermodynamic properties of proteins. Here we summarize our progress in understanding the role of charge-charge interactions for protein stability using examples drawn from my own laboratory and touch upon unanswered questions.

Thermal expansivities of peptides, polypeptides and proteins as measured by pressure perturbation calorimetry.

The main goal of this work was to provide direct experimental evidence that the expansivity of peptides, polypeptides and proteins as measured by pressure perturbation calorimetry (PPC), can serve as a proxy to characterize relative compactness of proteins, especially the denatured state ensemble. This is very important as currently only small angle X-ray scattering (SAXS), intrinsic viscosity and, to a lesser degree, fluorescence resonance transfer (FRET) experiments are capable of reporting on the compactness of denatured state ensembles. We combined the expansivity measurements with other biophysical methods (far-UV circular dichroism spectroscopy, differential scanning calorimetry, and small angle X-ray scattering). Three case studies of the effects of conformational changes on the expansivity of polypeptides in solution are presented. We have shown that expansivity appears to be insensitive to the helix-coil transition, and appears to reflect the changes in hydration of the side-chains. We also observed that the expansivity is sensitive to the global conformation of the polypeptide chain and thus can be potentially used to probe hydration of different collapsed states of denatured or even intrinsically disordered proteins.

ProteinVolume: calculating molecular van der Waals and void volumes in proteins.

BACKGROUND: Voids and cavities in the native protein structure determine the pressure unfolding of proteins. In addition, the volume changes due to the interaction of newly exposed atoms with solvent upon protein unfolding also contribute to the pressure unfolding of proteins. Quantitative understanding of these effects is important for predicting and designing proteins with predefined response to changes in hydrostatic pressure using computational approaches. The molecular surface volume is a useful metric that describes contribution of geometrical volume, which includes van der Waals volume and volume of the voids, to the total volume of a protein in solution, thus isolating the effects of hydration for separate calculations. RESULTS: We developed ProteinVolume, a highly robust and easy-to-use tool to compute geometric volumes of proteins. ProteinVolume generates the molecular surface of a protein and uses an innovative flood-fill algorithm to calculate the individual components of the molecular surface volume, van der Waals and intramolecular void volumes. ProteinVolume is user friendly and is available as a web-server or a platform-independent command-line version. CONCLUSIONS: ProteinVolume is a highly accurate and fast application to interrogate geometric volumes of proteins. ProteinVolume is a free web server available on http://gmlab.bio.rpi.edu . Free-standing platform-independent Java-based ProteinVolume executable is also freely available at this web site.

Crystal and NMR structures of a Trp-cage mini-protein benchmark for computational fold prediction.

To provide high-resolution X-ray crystallographic structures of a peptide with the Trp-cage fold, we prepared a cyclized version of this motif. Cyclized Trp-cage is remarkably stable and afforded two crystal forms suitable for X-ray diffraction. The resulting higher resolution crystal structures validate the prior NMR models and provide explanations for experimental observations that could not be rationalized by NMR structural data, including the structural basis for the increase in fold stability associated with motif cyclization and the manner in which a polar serine side chain is accommodated in the hydrophobic interior. A hexameric oligomer of the cyclic peptide is found in both crystal forms and indicates that under appropriate conditions, this minimized system may also serve as a model for protein-protein interactions.

Structural characterization of semen coagulum-derived SEM1(86-107) amyloid fibrils that enhance HIV-1 infection.

SEM1(86-107) is a 22-residue peptide corresponding to residues 86-107 in the semenogelin I protein. SEM1(86-107) is an abundant component of freshly liquefied semen and forms amyloid fibrils capable of enhancing HIV infection. To probe the factors affecting fibril formation and gain a better understanding of how differences in pH between semen and vaginal fluid affect fibril stability, this study determined the effect of pH on SEM1(86-107) fibril formation and dissociation. The SEM1(86-107) fibril structure (i.e., residues that comprise the fibrillar core) was also probed using hydrogen-deuterium exchange mass spectrometry (HDXMS) and hydroxyl radical-mediated protein modification. The average percent exposure to hydroxyl radical-mediated modification in the SEM1(86-107) fibrils was determined without requiring tandem mass spectrometry spectral acquisition or complete separation of modified peptides. It was found that the residue exposures calculated from HDXMS and hydroxyl radical-mediated modification were similar. These techniques demonstrated that three regions of SEM1(86-107) comprise the amyloid fibril core and that positively charged residues are exposed, suggesting that electrostatic interactions between SEM1(86-107) and HIV or the cell surface may be responsible for mediating HIV infection enhancement by the SEM1(86-107) fibrils.

The role of cross-chain ionic interactions for the stability of collagen model peptides.

The contribution of ionic interactions to the stability of the collagen triple helix was studied using molecular dynamics (MD) simulations and biophysical methods. To this end, we examined the stability of a host-guest collagen model peptide, Ac-GPOGPOGPYGXOGPOGPO-NH2, substituting KGE, KGD, EGK, and DGK for the YGX sequence. All-atom, implicit solvent MD simulations show that the fraction of cross-chain ionic interactions formed is different, with the most pronounced in the KGE and KGD sequences, and the least in the DGK sequence. To test whether the fraction of cross-chain ionic interactions correlates with the stability, experimental measurements of thermostability were done using differential scanning calorimetry and circular dichroism spectroscopy. It was found that the melting temperature is very similar for KGE and KGD peptides, whereas the EGK peptide has lower thermostability and the DGK peptide is the least thermostable. A novel, to our knowledge, computational protocol termed temperature-scan MD was applied to estimate the relative stabilities of the peptides from MD simulations. We found an excellent correlation between transition temperatures obtained from temperature-scan MD and those measured experimentally. These results suggest the importance of cross-chain ionic interactions for the stability of collagen triple helix and confirm the utility of MD simulations in predicting interactions and stability in this system.