A better understanding of protein aggregation is bound to translate into critical advances in several areas, including the treatment of misfolded protein disorders and the development of self-assembling biomaterials for novel commercial applications. Because of its ubiquity and clinical potential, albumin is one of the best-characterized models in protein aggregation research; but its properties in different conditions are not completely understood. Here, we carried out all-atom molecular dynamics simulations of albumin to understand how electrostatics can affect the conformation of a single albumin molecule just prior to self-assembly. We then analyzed the tertiary structure and solvent accessible surface area of albumin after electrostatically triggered partial denaturation. The data obtained from these single protein simulations allowed us to investigate the effect of electrostatic interactions between two proteins. The results of these simulations suggested that hydrophobic attractions and counterion binding may be strong enough to effectively overcome the electrostatic repulsions between the highly charged monomers. This work contributes to our general understanding of protein aggregation mechanisms, the importance of explicit consideration of free ions in protein solutions, provides critical new insights about the equilibrium conformation of albumin in its partially denatured state at low pH, and may spur significant progress in our efforts to develop biocompatible protein hydrogels driven by electrostatic partial denaturation.

The alanine-based peptide Ac-XX(A) ₇ OO-NH ₂ , referred to as XAO (where X, A, and O denote diaminobutyric acid, alanine, and ornithine, respectively), has recently been proposed to possess a well defined polyproline II (P _II ) conformation at low temperatures. Based on the results of extensive NMR and CD investigations combined with theoretical calculations, reported here, we present evidence that, on the contrary, this peptide does not have any significant amount of organized P _II structure but exists in an ensemble of conformations with a distorted bend in the N- and C-terminal regions. The conformational ensemble was obtained by molecular dynamics/simulated annealing calculations using the amber suite of programs with time-averaged distance and dihedral-angle restraints obtained from rotating-frame nuclear Overhauser effect (ROE) volumes and vicinal coupling constants ³ J _HNΗα , respectively. The computed ensemble-averaged radius of gyration R _g (7.4 ± 1.0) Å is in excellent agreement with that measured by small-angle x-ray scattering (SAXS) whereas, if the XAO peptide were in the P _II conformation, R _g would be 11.6 Å. Depending on the pH, peptide concentration, and temperature, the CD spectra of XAO do or do not possess the maximum with positive ellipticity in the 217-nm region, which is characteristic of the P _II structure, reflecting a shifting conformational equilibrium rather than an all-or-none transition. The “P _II conformation” should, therefore, be considered as one of the accessible conformational states of individual amino acid residues in peptides and proteins rather than as a structure of most of the chain in the early stage of folding.

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH₂ in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

Recent improvements in the protein-structure prediction method developed in our laboratory, based on the thermodynamic hypothesis, are described. The conformational space is searched extensively at the united-residue level by using our physics-based UNRES energy function and the conformational space annealing method of global optimization. The lowest-energy coarse-grained structures are then converted to an all-atom representation and energy-minimized with the ECEPP/3 force field. The procedure was assessed in two recent blind tests of protein-structure prediction. During the first blind test, we predicted large fragments of α and α+β proteins [60–70 residues with C ^α rms deviation (rmsd) <6 Å]. However, for α+β proteins, significant topological errors occurred despite low rmsd values. In the second exercise, we predicted whole structures of five proteins (two α and three α+β, with sizes of 53–235 residues) with remarkably good accuracy. In particular, for the genomic target TM0487 (a 102-residue α+β protein from Thermotoga maritima ), we predicted the complete, topologically correct structure with 7.3-Å C ^α rmsd. So far this protein is the largest α+β protein predicted based solely on the amino acid sequence and a physics-based potential-energy function and search procedure. For target T0198, a phosphate transport system regulator PhoU from T. maritima (a 235-residue mainly α-helical protein), we predicted the topology of the whole six-helix bundle correctly within 8 Å rmsd, except the 32 C-terminal residues, most of which form a β-hairpin. These and other examples described in this work demonstrate significant progress in physics-based protein-structure prediction.