Importance of Hydrophilic Hydration and Intramolecular Interactions in the Thermodynamics of Helix-Coil Transition and Helix-Helix Assembly in a Deca-Alanine Peptide

Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t. Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t. In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

Analysis of nonsynonymous SNPs in candidate genes that influence bovine temperament and evaluation of their effect in Brahman cattle

BACKGROUND

Temperament is an important production trait in cattle and multiple strategies had been developed to generate molecular markers to assist animal selection. As nonsynonymous single nucleotide polymorphisms are markers with the potential to affect gene functions, they could be useful to predict phenotypic effects. Genetic selection of less stress-responsive, temperamental animals is desirable from an economic and welfare point of view.

METHODS AND RESULTS

Two nonsynonymous single nucleotide polymorphisms identified in HTR1B and SLC18A2 candidate genes for temperament were analyzed in silico to determine their effects on protein structure. Those nsSNPs allowing changes in proteins were selected for a temperament association analysis in a Brahman population. Transversion effects on protein structure were evaluated in silico for each amino acid change model, revealing structural changes in the proteins of the HTR1B and SLC18A2 genes. The selected nsSNPs were genotyped in a Brahman population (n = 138), and their genotypic effects on three temperament traits were analyzed: exit velocity, pen score, and temperament score. Only the SNP rs209984404-HTR1B (C/A) showed a significant association (P = 0.0144) with pen score. The heterozygous genotype showed a pen score value 1.17 points lower than that of the homozygous CC genotype.

CONCLUSION

The results showed that in silico analysis could direct the selection of nsSNPs with the potential to change the protein. Non-synonymous single nucleotide polymorphisms causing structural changes and reduced protein stability were identified. Only rs209984404-HTR1B shows that the allele affecting protein stability was associated with the genotype linked to docility in cattle.

Toward a structural identification of metastable molecular conformations

Interpreting high-dimensional data from molecular dynamics simulations is a persistent challenge. In this paper, we show that for a small peptide, deca-alanine, metastable states can be identified through a neural net based on structural information alone. While processing molecular dynamics data, dimensionality reduction is a necessary step that projects high-dimensional data onto a low-dimensional representation that, ideally, captures the conformational changes in the underlying data. Conventional methods make use of the temporal information contained in trajectories generated through integrating the equations of motion, which forgoes more efficient sampling schemes. We demonstrate that EncoderMap, an autoencoder architecture with an additional distance metric, can find a suitable low-dimensional representation to identify long-lived molecular conformations using exclusively structural information. For deca-alanine, which exhibits several helix-forming pathways, we show that this approach allows us to combine simulations with different biasing forces and yields representations comparable in quality to other established methods. Our results contribute to computational strategies for the rapid automatic exploration of the configuration space of peptides and proteins.

Statistical Mechanical Design Principles for Coarse-Grained Interactions across Different Conformational Free Energy Surfaces

Systematic bottom-up coarse-graining (CG) of molecular systems provides a means to explore different coupled length and time scales while treating the molecular-scale physics at a reduced level. However, the configuration dependence of CG interactions often results in CG models with limited applicability for exploring the parametrized configurations. We propose a statistical mechanical theory to design CG interactions across different configurations and conditions. In order to span wide ranges of conformational space, distinct classical CG free energy surfaces for characteristic configurations are identified using molecular collective variables. The coupling interaction between different CG free energy surfaces can then be systematically determined by analogy to quantum mechanical approaches describing coupled states. The present theory can accurately capture the underlying many-body potentials of mean force in the CG variables for various order parameters applied to liquids, interfaces, and in principle proteins, uncovering the complex nature underlying the coupling interaction and imparting a new protocol for the design of predictive multiscale models.

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

A Structural Analysis of Proteinaceous Nanotube Cavities and Their Applications in Nanotechnology

Protein nanotubes offer unique properties to the materials science field that allow them to fulfill various functions in drug delivery, biosensors and energy storage. Protein nanotubes are chemically diverse, modular, biodegradable and nontoxic. Furthermore, although the initial design or repurposing of such nanotubes is highly complex, the field has matured to understand underlying chemical and physical properties to a point where applications are successfully being developed. An important feature of a nanotube is its ability to bind ligands via its internal cavities. As ligands of interest vary in size, shape and chemical properties, cavities have to be able to accommodate very specific features. As such, understanding cavities on a structural level is essential for their effective application. The objective of this review is to present the chemical and physical diversity of protein nanotube cavities and highlight their potential applications in materials science, specifically in biotechnology.

α-Helices propagating from stable nucleators exhibit unconventional thermal folding

Although the effect of thermal perturbations on protein structure has long been modeled in helical peptides, several details, such as the relation between the thermal stabilities of the propagating and nucleating segments of helices, remain elusive. We had earlier reported on the helix-nucleating propensities of covalent H-bond surrogate-constrained α-turns. Here, we analyze the thermal stabilities of helices that propagate along peptides appended to these α-helix nucleators using their NMR and far-UV CD spectra. Unconventional thermal folding of these helix models reveals that the helical fold in propagating backbones resists thermal perturbations as long as their nucleating template is intact. The threshold temperature of such resistance is also influenced by the extent of similarity between the natures of helical folds in the nucleating and propagating segments. Correlations between helicities and rigidities of helix-nucleating and helix-propagating segments reveal subtle interdependence, which explains cooperativity and residual helix formation during protein folding.

Novel Peptide-Polymer Conjugate with pH-Responsive Targeting/Disrupting Effects on Biomembranes

Conjugating polymers to peptides has become a new strategy of designing functional antitumor agents for their improved stability and enhanced activity. In this paper, a novel peptide-polymer conjugate PEPc-PMAA with pH responsiveness was designed and synthesized. The isoelectric point of PEPc was studied by dynamic light scattering for the targeting effect. Also, the transmittances of PMAA at different pHs were measured using an ultraviolet-visible spectrophotometer for determining the triggering pH of the disrupting effect. The results showed that PEPc-PMAA was hydrophilic under neutral conditions and changed to be amphiphilic composed of positively charged PEPc and hydrophobic PMAA under acidic conditions. The interactions between PEPc-PMAA and mimic cells were investigated by the measurements of membrane fluidity and cargo leakage from 1,2-dipalmitoyl-sn-glycerol-3-phosphocholine and 1,2-dipalmitoyl-sn-glycerol-3-phospho-(1-rac-glycerol) (DPPG) liposomes. It proved that PEPc-PMAA caused a distinct membrane disturbance of the DPPG liposome at pH 5.5, resulting in more serious cargo leakage. Because of its targeting and disrupting effects on negatively charged biomembranes under acidic conditions, PEPc-PMAA showed its good potential as an antitumor agent.

Thermodynamics of Hydration from the Perspective of the Molecular Quasichemical Theory of Solutions

The quasichemical organization of the potential distribution theorem, molecular quasichemical theory (QCT), enables practical calculations and also provides a conceptual framework for molecular hydration phenomena. QCT can be viewed from multiple perspectives: (a) as a way to regularize an ill-conditioned statistical thermodynamic problem; (b) as an introduction of and emphasis on the neighborship characteristics of a solute of interest; or (c) as a way to include accurate electronic structure descriptions of near-neighbor interactions in defensible statistical thermodynamics by clearly defining neighborship clusters. The theory has been applied to solutes of a wide range of chemical complexity, ranging from ions that interact with water with both long-ranged and chemically intricate short-ranged interactions, to solutes that interact with water solely through traditional van der Waals interations, and including water itself. The solutes range in variety from monatomic ions to chemically heterogeneous macromolecules. A notable feature of QCT is that, in applying the theory to this range of solutes, the theory itself provides guidance on the necessary approximations and simplifications that can facilitate the calculations. In this Perspective, we develop these ideas and document them with examples that reveal the insights that can be extracted using the QCT formulation.

Hydrophilic Interactions Dominate the Inverse Temperature Dependence of Polypeptide Hydration Free Energies Attributed to Hydrophobicity

We address the association of the hydrophobic driving forces in protein folding with the inverse temperature dependence of protein hydration, wherein stabilizing hydration effects strengthen with increasing temperature in a physiological range. All-atom calculations of the free energy of hydration of aqueous deca-alanine conformers, holistically including backbone and side-chain interactions together, show that attractive peptide-solvent interactions and the thermal expansion of the solvent dominate the inverse temperature signatures that have been interpreted traditionally as the hydrophobic stabilization of proteins in aqueous solution. Equivalent calculations on a methane solute are also presented as a benchmark for comparison. The present study calls for a reassessment of the forces that stabilize folded protein conformations in aqueous solutions and of the additivity of hydrophobic/hydrophilic contributions.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Orientational and Folding Thermodynamics via Electric Dipole Moment Restraining

The projection of molecular processes onto a small set of relevant descriptors, the so-called reaction coordinates or collective variables (CVs), is a technique nowadays routinely employed by the biomolecular simulation community. In this work, we implemented two CVs to manipulate the orientation (i.e., angle) (μ⃗_a) and magnitude (|μ⃗|) of the electric dipole moment. In doing so, we studied the thermodynamics of water orientation under the application of external voltages and the folding of two polypeptides at zero-field conditions. The projection of the free-energy [potential of mean force (PMF)] along water orientation defined an upper limit of around 0.3 V for irrelevant thermodynamic effects. On the other hand, sufficiently strong μ⃗_a restraints applied on 12-alanine (Ala₁₂) triggered structural effects because of the alignment of local dipoles; for lower restraints, a full-body rotation is achieved. The manipulation of |μ⃗| produced strong perturbations on the secondary structure of Ala₁₂, promoting an enhanced sampling to its configurational space. Rigorous free-energy calculations in the form of 2-D PMFs for deca-alanine showed the utility of |μ⃗| as a reaction coordinate to study folding in small α helices. As a whole, we propose that the manipulation of both components of the dipole moment, μ⃗_a and |μ⃗|, provides thermodynamics insights into the structural conformation and stability of biomolecules. These new CVs are implemented in the Colvars module, available for NAMD and LAMMPS.

Quasi-Chemical Theory with Cluster Sampling from Ab Initio Molecular Dynamics: Fluoride (F-) Anion Hydration

Accurate predictions of the hydration free energy for anions typically has been more challenging than that for cations. Hydrogen bond donation to the anion in hydrated clusters such as F(H₂O) _n ^- can lead to delicate structures. Consequently, the energy landscape contains many local minima, even for small clusters, and these minima present a challenge for computational optimization. Utilization of cluster experimental results for the free energies of gas-phase clusters shows that even though anharmonic effects are interesting they need not be of troublesome magnitudes for careful applications of quasi-chemical theory to ion hydration. Energy-optimized cluster structures for anions can leave the central ion highly exposed, and application of implicit solvation models to these structures can incur more serious errors than those for metal cations. Utilizing cluster structures sampled from ab initio molecular dynamics simulations substantially fixes those issues.

Quasi-chemical theory of F-(aq): The "no split occupancies rule" revisited

We use ab initio molecular dynamics (AIMD) calculations and quasi-chemical theory (QCT) to study the inner-shell structure of F^-(aq) and to evaluate that single-ion free energy under standard conditions. Following the "no split occupancies" rule, QCT calculations yield a free energy value of -101 kcal/mol under these conditions, in encouraging agreement with tabulated values (-111 kcal/mol). The AIMD calculations served only to guide the definition of an effective inner-shell constraint. QCT naturally includes quantum mechanical effects that can be concerning in more primitive calculations, including electronic polarizability and induction, electron density transfer, electron correlation, molecular/atomic cooperative interactions generally, molecular flexibility, and zero-point motion. No direct assessment of the contribution of dispersion contributions to the internal energies has been attempted here, however. We anticipate that other aqueous halide ions might be treated successfully with QCT, provided that the structure of the underlying statistical mechanical theory is absorbed, i.e., that the "no split occupancies" rule is recognized.

Intramolecular Interactions Overcome Hydration to Drive the Collapse Transition of Gly15

Simulations and experiments show oligo-glycines, polypeptides lacking any side chains, can collapse in water. We assess the hydration thermodynamics of this collapse by calculating the hydration free energy at each of the end points of the reaction coordinate, here taken as the end-to-end distance (r) in the chain. To examine the role of the various conformations for a given r, we study the conditional distribution, P(R_g|r), of the radius of gyration for a given value of r. The free energy change versus R_g, -k_BT ln P(R_g|r), is found to vary more gently compared to the corresponding variation in the excess hydration free energy. Using this observation within a multistate generalization of the potential distribution theorem, we calculate a tight upper bound for the hydration free energy of the peptide for a given r. On this basis, we find that peptide hydration greatly favors the expanded state of the chain, despite primitive hydrophobic effects favoring chain collapse. The net free energy of collapse is seen to be a delicate balance between opposing intrapeptide and hydration effects, with intrapeptide contributions favoring collapse.

Solvation Thermodynamics of Oligoglycine with Respect to Chain Length and Flexibility

Oligoglycine is a backbone mimic for all proteins and is prevalent in the sequences of intrinsically disordered proteins. We have computed the absolute chemical potential of glycine oligomers at infinite dilution by simulation with the CHARMM36 and Amber ff12SB force fields. We performed a thermodynamic decomposition of the solvation free energy (ΔG(sol)) of Gly2-5 into enthalpic (ΔH(sol)) and entropic (ΔS(sol)) components as well as their van der Waals and electrostatic contributions. Gly2-5 was either constrained to a rigid/extended conformation or allowed to be completely flexible during simulations to assess the effects of flexibility on these thermodynamic quantities. For both rigid and flexible oligoglycine models, the decrease in ΔG(sol) with chain length is enthalpically driven with only weak entropic compensation. However, the apparent rates of decrease of ΔG(sol), ΔH(sol), ΔS(sol), and their elec and vdw components differ for the rigid and flexible models. Thus, we find solvation entropy does not drive aggregation for this system and may not explain the collapse of long oligoglycines. Additionally, both force fields yield very similar thermodynamic scaling relationships with respect to chain length despite both force fields generating different conformational ensembles of various oligoglycine chains.

Predicting unfolding thermodynamics and stable intermediates for alanine-rich helical peptides with the aid of coarse-grained molecular simulation

This report focuses on the molecular-level processes and thermodynamics of unfolding of a series of helical peptides using a coarse-grained (CG) molecular model. The CG model was refined to capture thermodynamics and structural changes as a function of temperature for a set of published peptide sequences. Circular dichroism spectroscopy (CD) was used to experimentally monitor the temperature-dependent conformational changes and stability of published peptides and new sequences introduced here. The model predictions were quantitatively or semi-quantitatively accurate in all cases. The simulations and CD results showed that, as expected, in most cases the unfolding of helical peptides is well described by a simply 2-state model, and conformational stability increased with increased length of the helices. A notable exception in a 19-residue helix was when two Ala residues were each replaced with Phe. This stabilized a partly unfolded intermediate state via hydrophobic contacts, and also promoted aggregates at higher peptide concentrations.