Potential energy landscape of a flexible water model: Equation of state, configurational entropy, and Adam-Gibbs relationship.

The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid-liquid critical point (LLCP) in the supercooled regime, at (Pc = 150 MPa, Tc = 190 K, and ρc = 1.04 g/cm3) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water and the corresponding location of the LLCP are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian, which allows us to calculate the configurational entropy of the system, Sconf. The Sconf of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam-Gibbs relation, which relates the diffusion coefficient D with Sconf, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high frequencies lead to unphysical (negative) entropy for q-TIP4P/F water when using classical statistical mechanics (yet, the PEL formalism can be applied successfully).

The Harmonic and Gaussian Approximations in the Potential Energy Landscape Formalism for Quantum Liquids.

The potential energy landscape (PEL) formalism has been used in the past to describe the behavior of classical low-temperature liquids and glasses. Here, we extend the PEL formalism to describe the behavior of liquids and glasses that obey quantum mechanics. In particular, we focus on the (i) harmonic and (ii) Gaussian approximations of the PEL, which have been commonly used to describe classical systems, and show how these approximations can be applied to quantum liquids/glasses. Contrary to the case of classical liquids/glasses, the PEL of quantum liquids is temperature-dependent, and hence, the main expressions resulting from approximations (i) and (ii) depend on the nature (classical vs quantum) of the system. The resulting theoretical expressions from the PEL formalism are compared with results from path-integral Monte Carlo (PIMC) simulations of a monatomic model liquid. In the PIMC simulations, every atom of the quantum liquid is represented by a ring-polymer. Our PIMC simulations show that at the local minima of the PEL (inherent structures, or IS), sampled over a wide range of temperatures and volumes, the ring-polymers are collapsed. This considerably facilitates the description of quantum liquids using the PEL formalism. Specifically, the normal modes of the ring-polymer system/quantum liquid at an IS can be calculated analytically if the normal modes of the classical liquid counterpart are known (as obtained, e.g., from classical MC or molecular dynamics simulations of the corresponding atomic liquid).

A continuum of amorphous ices between low-density and high-density amorphous ice.

Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.

The Importance of Nuclear Quantum Effects on the Thermodynamic and Structural Properties of Low-Density Amorphous Ice: A Comparison with Hexagonal Ice.

We study the nuclear quantum effects (NQE) on the thermodynamic properties of low-density amorphous ice (LDA) and hexagonal ice (Ih) at P = 0.1 MPa and T ≥ 25 K. Our results are based on path-integral molecular dynamics (PIMD) and classical MD simulations of H2O and D2O using the q-TIP4P/F water model. We show that the inclusion of NQE is necessary to reproduce the experimental properties of LDA and ice Ih. While MD simulations (no NQE) predict that the density ρ(T) of LDA and ice Ih increases monotonically upon cooling, PIMD simulations indicate the presence of a density maximum in LDA and ice Ih. MD and PIMD simulations also predict a qualitatively different T-dependence for the thermal expansion coefficient αP(T) and bulk modulus B(T) of both LDA and ice Ih. Remarkably, the ρ(T), αP(T), and B(T) of LDA are practically identical to those of ice Ih. The origin of the observed NQE is due to the delocalization of the H atoms, which is identical in LDA and ice Ih. H atoms delocalize considerably (over a distance ≈ 20-25% of the OH covalent-bond length) and anisotropically (preferentially perpendicular to the OH covalent bond), leading to less linear hydrogen bonds HB (larger HOO angles and longer OO separations) than observed in classical MD simulations.

Sc-HOPO: A Potential Construct for Use in Radioscandium-Based Radiopharmaceuticals.

Three isotopes of scandium─43Sc, 44Sc, and 47Sc─have attracted increasing attention as potential candidates for use in imaging and therapy, respectively, as well as for possible theranostic use as an elementally matched pair. Here, we present the octadentate chelator 3,4,3-(LI-1,2-HOPO) (or HOPO), an effective chelator for hard cations, as a potential ligand for use in radioscandium constructs with simple radiolabeling under mild conditions. HOPO forms a 1:1 Sc-HOPO complex that was fully characterized, both experimentally and theoretically. [47Sc]Sc-HOPO exhibited good stability in chemical and biological challenges over 7 days. In healthy mice, [43,47Sc]Sc-HOPO cleared the body rapidly with no signs of demetalation. HOPO is a strong candidate for use in radioscandium-based radiopharmaceuticals.

Understanding Proton Transfer in Non-aqueous Biopolymers based on Helical Peptides: A Quantum Mechanical Study.

Histidine (an imidazole-based amino acid) is a promising building block for short aromatic peptides containing a proton donor/acceptor moiety. Previous studies have shown that polyalanine helical peptides substituted at regular intervals with histidine residues exhibit both structural stability as well as high proton affinity and high conductivity. Here, we present first-principle calculations of non-aqueous histidine-containing 310-, α- and π-helices and show that they are able to form hydrogen-bonded networks mimicking proton wires that have the ability to shuttle protons via the Grotthuss shuttling mechanism. The formation of these wires enhances the stability of the helices, and our structural characterizations confirm that the secondary structures are conserved despite distortions of the backbones. In all cases, the helices exhibit high proton affinity and proton transfer barriers on the order of 1~4 kcal/mol. Zero-point energy calculations suggest that for these systems, ground state vibrational energy can provide enough energy to cross the proton transport energy barrier. Additionally, ab initio molecular dynamics results suggests that the protons are transported unidirectionally through the wire at a rate of approximately 2 Å every 20 fs. These results demonstrate that efficient deprotonation-controlled proton wires can be formed using non-aqueous histidine-containing helical peptides.

Anomalous properties in the potential energy landscape of a monatomic liquid across the liquid-gas and liquid-liquid phase transitions.

As a liquid approaches the gas state, the properties of the potential energy landscape (PEL) sampled by the system become anomalous. Specifically, (i) the mechanically stable local minima of the PEL [inherent structures (IS)] can exhibit cavitation above the so-called Sastry volume, vS, before the liquid enters the gas phase. In addition, (ii) the pressure of the liquid at the sampled IS [i.e., the PEL equation of state, PIS(v)] develops a spinodal-like minimum at vS. We perform molecular dynamics simulations of a monatomic water-like liquid and verify that points (i) and (ii) hold at high temperatures. However, at low temperatures, cavitation in the liquid and the corresponding IS occurs simultaneously and a Sastry volume cannot be defined. Remarkably, at intermediate/high temperatures, the IS of the liquid can exhibit crystallization, i.e., the liquid regularly visits the regions of the PEL that belong to the crystal state. The model liquid considered also exhibits a liquid-liquid phase transition (LLPT) between a low-density and a high-density liquid (LDL and HDL). By studying the behavior of PIS(v) during the LLPT, we identify a Sastry volume for both LDL and HDL. The HDL Sastry volume marks the onset above which IS are heterogeneous (composed of LDL and HDL particles), analogous to points (i) and (ii) above. However, the relationship between the LDL Sastry volume and the onset of heterogeneous IS is less evident. We conclude by presenting a thermodynamic argument that can explain the behavior of the PEL equation of state PIS(v) across both the liquid-gas phase transition and LLPT.

Nuclear quantum effects on the dynamics and glass behavior of a monatomic liquid with two liquid states.

We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid-liquid phase transition and liquid-liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., Tg LDL(P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., Tg HDL(P) increases upon compression. (iii) NQE shift both Tg LDL(P) and Tg HDL(P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA-HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA-HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA-HDA transformation.

Evidence of a liquid-liquid phase transition in H[Formula: see text]O and D[Formula: see text]O from path-integral molecular dynamics simulations.

We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H[Formula: see text]O and D[Formula: see text]O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density [Formula: see text], isothermal compressibility [Formula: see text], and self-diffusion coefficients D(T) of H[Formula: see text]O and D[Formula: see text]O are in excellent agreement with available experimental data; the isobaric heat capacity [Formula: see text] obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H[Formula: see text]O and D[Formula: see text]O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H[Formula: see text]O and D[Formula: see text]O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H[Formula: see text]O, from PIMD simulations, is located at [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D[Formula: see text]O is estimated to be [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of [Formula: see text] for D[Formula: see text]O and, particularly, H[Formula: see text]O suggest that improved water models are needed for the study of supercooled water.

The role of high-density and low-density amorphous ice on biomolecules at cryogenic temperatures: a case study with polyalanine.

Experimental techniques, such as cryo-electron microscopy, require biological samples to be recovered at cryogenic temperatures (T ≈ 100 K) with water being in an amorphous ice state. However, (bulk) water can exist in two amorphous ices at P < 1 GPa, low-density amorphous (LDA) ice at low pressures and high-density amorphous ice (HDA) at high pressures; HDA is ≈20-25% denser than LDA. While fast/plunge cooling at 1 bar brings the sample into LDA, high-pressure cooling (HPC), at sufficiently high pressure, produces HDA. HDA can also be produced by isothermal compression of LDA at cryogenic temperatures. Here, we perform classical molecular dynamics simulations to study the effects of LDA, HDA, and the LDA-HDA transformation on the structure and hydration of a small peptide, polyalanine. We follow thermodynamic paths corresponding to (i) fast/plunge cooling at 1 bar, (ii) HPC at P = 400 MPa, and (iii) compression/decompression cycles at T = 80 K. While process (i) produced LDA in the system, path (iii) produces HDA. Interestingly, the amorphous ice produced in process (ii) is an intermediate amorphous ice (IA) with properties that fall in-between those of LDA and HDA. Remarkably, the structural changes in polyalanine are negligible at all conditions studied (0-2000 MPa, 80-300 K) even when water changes among the low and high-density liquid states as well as the amorphous solids LDA, IA, and HDA. The similarities and differences in the hydration of polyalanine vitrified in LDA, IA, and HDA are described. Since the studied thermodynamic paths are suitable for the cryopreservation of biomolecules, we also study the structure and hydration of polyalanine along isobaric and isochoric heating paths, which can be followed experimentally for the recovery of cryopreserved samples. Upon heating, the structure of polyalanine remains practically unchanged. We conclude with a brief discussion of the practical advantages of (a) using HDA and IA as a cryoprotectant environment (as opposed to LDA), and (b) the use of isochoric heating as a recovery process (as opposed to isobaric heating).

Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water.

We perform path-integral molecular dynamics (PIMD) simulations of H2O and D2O using the q-TIP4P/F model. Simulations are performed at P = 1 bar and over a wide range of temperatures that include the equilibrium (T≥ 273 K) and supercooled (210 ≤T < 273 K) liquid states of water. The densities of both H2O and D2O calculated from PIMD simulations are in excellent agreement with experiments in the equilibrium and supercooled regimes. We also evaluate important thermodynamic response functions, specifically, the thermal expansion coefficient αP(T), isothermal compressibility κT(T), isobaric heat capacity CP(T), and static dielectric constant ε(T). While these properties are in excellent [αP(T) and κT(T)] or semi-quantitative agreement [CP(T) and ε(T)] with experiments in the equilibrium regime, they are increasingly underestimated upon further cooling. It follows that the inclusion of nuclear quantum effects in PIMD simulations of (q-TIP4P/F) water is not sufficient to reproduce the anomalous large fluctuations in density, entropy, and electric dipole moment characteristic of supercooled water. It has been hypothesized that water may exhibit a liquid-liquid critical point (LLCP) in the supercooled regime at P > 1 bar and that such a LLCP generates a maximum in CP(T) and κT(T) at 1 bar. Consistent with this hypothesis and in particular, with experiments, we find a maximum in the κT(T) of q-TIP4P/F light and heavy water at T≈ 230-235 K. No maximum in CP(T) could be detected down to T≥ 210 K. We also calculate the diffusion coefficient D(T) of H2O and D2O using the ring-polymer molecular dynamics (RPMD) technique and find that computer simulations are in remarkable good agreement with experiments at all temperatures studied. The results from RPMD/PIMD simulations are also compared with the corresponding results obtained from classical MD simulations of q-TIP4P/F water where atoms are represented by single interacting sites. Surprisingly, we find minor differences in most of the properties studied, with CP(T), D(T), and structural properties being the only (expected) exceptions.

Serological study of brucellosis in Argentine Creole sheep / Estudio serológico de brucelosis en ovejas criollas de Argentina

Ovine cattle was introduced into America during the Spanish conquest with the second journey of Columbus to the Antilles and was disseminated throughout the region. In 1587, sheep were introduced into Argentina, later developing into the "Creole" breed. We selected 486 animals from different Argentine provinces with the aim of determining the serological status of brucellosis caused by Brucella melitensis and Brucella ovis. For the detection of antibodies against smooth Brucella spp., the Rose Bengal test (RBT) was performed as screening test while the serum agglutination test (SAT) and 2 mercapto-ethanol (2ME) were run as a confirmatory technique. Moreover, for the detection of antibodies against rough Brucella spp., we used the rapid slide agglutination test (RSAT) for screening and an indirect ELISA (IELISA) as confirmatory assay. This study showed that the total positive percentage of brucellosis due to B. ovis was 2.9%. Excluding the animals mixed with the Suffolk breed; seropositivity would be 0.6%. All animals tested negative for brucellosis caused by B. melitensis.

El ganado ovino fue introducido en América durante el segundo viaje de Colón a las Antillas y rápidamente se diseminó por la región. En 1587 entraron a la Argentina las primeras ovejas, estas se adaptaron a las condiciones naturales y se desarrolló así la denominada «oveja criolla¼. Se seleccionaron 486 animales de diferentes provincias para conocer el estado serológico de la enfermedad producida por Brucella melitensis y Brucella ovis. Para detectar anticuerpos anti-Brucella en fase lisa se utilizó la prueba de rosa de Bengala (RBT) como tamiz y como confirmatorias seroaglutinación lenta en tubo (SAT) y 2-mercaptoetanol (2ME). En fase rugosa se realizó como tamiz la prueba rápida de aglutinación (RSAT) y ELISA indirecto como prueba confirmatoria (IELISA). Este estudio mostró un porcentaje de seropositividad de brucelosis causada por B. ovis del 2,9%. Excluyendo los animales cruzados con la raza Suffolk, la seropositividad sería del 0,6%. Todos los animales fueron negativos a brucelosis causada por B. melitensis.

Nuclear quantum effects on the liquid-liquid phase transition of a water-like monatomic liquid.

Polyamorphic substances have the ability to exist in more than one liquid and/or glass states. Examples include water, silicon, and hydrogen. In many of these substances, nuclear quantum effects may become important in the proximity of the liquid-liquid and glass-glass transformation. Here, we study the nuclear quantum effects on a monatomic liquid that exhibits water-like anomalous properties and a liquid-liquid phase transition (LLPT) ending at a liquid-liquid critical point (LLCP). By performing path integral Monte Carlo simulations with different values of the Planck's constant h, we are able to explore how the location of the LLCP/LLPT in the P-T plane shifts as the system evolves from classical, h = 0, to quantum, h > 0. We find that, as the quantum nature of the liquid (as quantified by h) increases, and the atoms in the liquid become more delocalized, the LLCP pressure increases, the LLCP temperature decreases, and the LLCP volume remains constant. In addition, the crystallization temperature decreases with increasing h. For large values of h, the LLCP is not accessible due to rapid crystallization. The structure of the liquids studied at different values of h are also investigated.

Serological study of brucellosis in Argentine Creole sheep.

Ovine cattle was introduced into America during the Spanish conquest with the second journey of Columbus to the Antilles and was disseminated throughout the region. In 1587, sheep were introduced into Argentina, later developing into the "Creole" breed. We selected 486 animals from different Argentine provinces with the aim of determining the serological status of brucellosis caused by Brucella melitensis and Brucella ovis. For the detection of antibodies against smooth Brucella spp., the Rose Bengal test (RBT) was performed as screening test while the serum agglutination test (SAT) and 2 mercapto-ethanol (2ME) were run as a confirmatory technique. Moreover, for the detection of antibodies against rough Brucella spp., we used the rapid slide agglutination test (RSAT) for screening and an indirect ELISA (IELISA) as confirmatory assay. This study showed that the total positive percentage of brucellosis due to B. ovis was 2.9%. Excluding the animals mixed with the Suffolk breed; seropositivity would be 0.6%. All animals tested negative for brucellosis caused by B. melitensis.

Probing the Limits of Supramolecular G-Quadruplexes Using Atomistic Molecular Dynamics Simulations.

Guanosine and related derivatives self-assemble in the presence of cations like potassium into supramolecular G-quadruplexes (SGQs), where four guanine moieties form planar tetrads (T) that coaxially stack into columnar aggregates with broad size distributions. However, SGQs made from 8-aryl-2'-deoxyguanosine derivatives (8ArGs), form mostly octamers, or two-tetrad (2T)-SGQs, while some form dodecamers (3T-SGQs), or hexadecamers (4T-SGQs), and none reported to date form higher assemblies. A theoretical model that addresses the configurational space available for the multiple pathways available for 8ArGs to self-assemble into SGQs is used to frame a series of molecular dynamics simulations (MDS) with selected SGQs. Some key insights from this work include: (a) The predicted entropic costs are not significantly higher for SGQs with more subunits due to their hierarchical assembly pathways; (b) The multiple isomeric SGQs vary in the interfacial contacts between consecutive tetrads, due to their two distinct sides (head, h; tail, t), with the MDS supporting the predicted order of stability of hh > ht > tt for octamers. (c) Such order also applies to dodecamers and hexadecamers, but with context-dependent exceptions due to strong allosteric effects. (d) The main factor disfavoring the tt interface is the repulsive dipolar interactions between the O4' from ribose moieties on adjacent tetrads. (e) SGQs with 5 or more tetrads are disfavored because the attractive interactions are not large or strong enough to overcome the many repulsive forces resulting from the addition of further tetrads. We expect these findings provide some guidelines to enable the further development of SGQs into functional materials.

The gamma-butyrolactone receptors BulR1 and BulR2 of Streptomyces tsukubaensis: tacrolimus (FK506) and butyrolactone synthetases production control.

Streptomyces tsukubaensis is a well-established industrial tacrolimus producer strain, but its molecular genetics is very poorly known. This information shortage prevents the development of tailored mutants in the regulatory pathways. A region (named bul) contains several genes involved in the synthesis and control of the gamma-butyrolactone autoregulator molecules. This region contains ten genes (bulA, bulZ, bulY, bulR2, bulS2, bulR1, bulW, bluB, bulS1, bulC) including two γ-butyrolactone receptor homologues (bulR1, bulR2), two putative gamma-butyrolactone synthetase homologues (bulS1, bulS2) and two SARP regulatory genes (bulY, bulZ). Analysis of the autoregulatory element (ARE)-like sequences by electrophoretic mobility shift assays and footprinting using the purified BulR1 and BulR2 recombinant proteins revealed six ARE regulatory sequences distributed along the bul cluster. These sequences showed specific binding of both BulR1 (the gamma-butyrolactone receptor) and BulR2, a possible pseudo γ-butyrolactone receptor. The protected region in all cases covered a 28-nt sequence with a palindromic structure. Optimal docking area analysis of BulR1 proved that this protein can be presented as either monomer or dimer but not oligomers and that it binds to the conserved ARE sequence in both strands. The effect on tacrolimus production was analysed by deletion of the bulR1 gene, which resulted in a strong decrease of tacrolimus production. Meanwhile, the ΔbulR2 mutation did not affect the biosynthesis of this immunosuppressant.

Src activation by ß-adrenoreceptors is a key switch for tumour metastasis.

Noradrenaline can modulate multiple cellular functions important for cancer progression; however, how this single extracellular signal regulates such a broad array of cellular processes is unknown. Here we identify Src as a key regulator of phosphoproteomic signalling networks activated in response to beta-adrenergic signalling in cancer cells. These results also identify a new mechanism of Src phosphorylation that mediates beta-adrenergic/PKA regulation of downstream networks, thereby enhancing tumour cell migration, invasion and growth. In human ovarian cancer samples, high tumoural noradrenaline levels were correlated with high pSrc(Y419) levels. Moreover, among cancer patients, the use of beta blockers was significantly associated with reduced cancer-related mortality. Collectively, these data provide a pivotal molecular target for disrupting neural signalling in the tumour microenvironment.

Modeling nonaqueous proton wires built from helical peptides: biased proton transfer driven by helical dipoles.

We report gas-phase electronic structure calculations on helical peptides that act as scaffolds for imidazole-based hydrogen-bonding networks (proton wires). We have modeled various 21-residue polyalanine peptides substituted at regular intervals with histidines (imidazole-bearing amino acids), using a hybrid approach with a semiempirical method (AM1) for peptide scaffolds and density functional theory (B3LYP) for proton wires. We have computed energy landscapes including barriers for Grotthuss-shuttling-type proton motions though wires supported on 3(10)-, α- and π-helical structures, showing the 3(10)- and α-helices to be attractive targets in terms of high proton affinities, low Grotthuss shuttling barriers, and high stabilities. Moreover, bias forces provided by the helical dipole moments were found to promote unidirectional proton translocation.

Using a reduced dimensionality model to compute the thermodynamic properties of finite polypeptide aggregates.

By implementing a simple reduced dimensionality model to describe the interactions in finite systems composed of two seven-amino-acid peptides, the thermodynamic properties of ordered and disordered aggregates were computed. Within this model, the hydrophobicity of each amino acid was varied, and the stability of the systems computed. Accurate averages in the canonical ensemble were obtained using various replica exchange Monte Carlo algorithms. Low and high temperature regions were encountered where the ordered and disordered aggregates were stabilized. It was observed that as the degree of hydrophobicity increased, the stability of the aggregates increased, with a significant energetic stabilization obtained for the ordered aggregates. Upon decreasing the concentration of the solution, the stability of the amorphous aggregates increased when compared to the ordered systems.

Model for the peptide-free conformation of class II MHC proteins.

BACKGROUND: Major histocompatibility complex proteins are believed to undergo significant conformational changes concomitant with peptide binding, but structural characterization of these changes has remained elusive. METHODOLOGY/PRINCIPAL FINDINGS: Here we use molecular dynamics simulations and experimental probes of protein conformation to investigate the peptide-free state of class II MHC proteins. Upon computational removal of the bound peptide from HLA-DR1-peptide complex, the alpha50-59 region folded into the P1-P4 region of the peptide binding site, adopting the same conformation as a bound peptide. Strikingly, the structure of the hydrophobic P1 pocket is maintained by engagement of the side chain of Phe alpha54. In addition, conserved hydrogen bonds observed in crystal structures between the peptide backbone and numerous MHC side chains are maintained between the alpha51-55 region and the rest of the molecule. The model for the peptide-free conformation was evaluated using conformationally-sensitive antibody and superantigen probes predicted to show no change, moderate change, or dramatic changes in their interaction with peptide-free DR1 and peptide-loaded DR1. The binding observed for these probes is in agreement with the movements predicted by the model. CONCLUSION/SIGNIFICANCE: This work presents a molecular model for peptide-free class II MHC proteins that can help to interpret the conformational changes known to occur within the protein during peptide binding and release, and can provide insight into possible mechanisms for DM action.