Localized Soft Vibrational Modes and Coherent Structural Phase Transformations in Rutile TiO2 Nanoparticles under Negative Pressure

We study the effect of size on the vibrational modes and frequencies of nanoparticles, by applying a newly developed, robust, and efficient first-principles-based method that we present in outline. We focus on rutile TiO₂, a technologically important material whose bulk exhibits a softening of a transverse acoustic mode close to q=(12,12,14), which becomes unstable with the application of negative pressure. We demonstrate that, under these conditions, nanoparticles above a critical size exhibit unstable localized modes and we calculate their characteristic localization length and decomposition with respect to bulk phonons. We propose that such localized soft modes could initiate coherent structural phase transformations in small nanoparticles above a critical size.

Investigating the quasi-liquid layer on ice surfaces: a comparison of order parameters

Ice surfaces are characterized by pre-melted quasi-liquid layers (QLLs), which mediate both crystal growth processes and interactions with external agents. Understanding QLLs at the molecular level is necessary to unravel the mechanisms of ice crystal formation. Computational studies of the QLLs heavily rely on the accuracy of the methods employed for identifying the local molecular environment and arrangements, discriminating between solid-like and liquid-like water molecules. Here we compare the results obtained using different order parameters to characterize the QLLs on hexagonal ice (Ih) and cubic ice (Ic) model surfaces investigated with molecular dynamics (MD) simulations in a range of temperatures. For the classification task, in addition to the traditional Steinhardt order parameters in different flavours, we select an entropy fingerprint and a deep learning neural network approach (DeepIce), which are conceptually different methodologies. We find that all the analysis methods give qualitatively similar trends for the behaviours of the QLLs on ice surfaces with temperature, with some subtle differences in the classification sensitivity limited to the solid-liquid interface. The thickness of QLLs on the ice surface increases gradually as the temperature increases. The trends of the QLL size and of the values of the order parameters as a function of temperature for the different facets may be linked to surface growth rates which, in turn, affect crystal morphologies at lower vapour pressure. The choice of the order parameter can be therefore informed by computational convenience except in cases where a very accurate determination of the liquid-solid interface is important.

Combined Free-Energy Calculation and Machine Learning Methods for Understanding Ligand Unbinding Kinetics

The determination of drug residence times, which define the time an inhibitor is in complex with its target, is a fundamental part of the drug discovery process. Synthesis and experimental measurements of kinetic rate constants are, however, expensive and time consuming. In this work, we aimed to obtain drug residence times computationally. Furthermore, we propose a novel algorithm to identify molecular design objectives based on ligand unbinding kinetics. We designed an enhanced sampling technique to accurately predict the free-energy profiles of the ligand unbinding process, focusing on the free-energy barrier for unbinding. Our method first identifies unbinding paths determining a corresponding set of internal coordinates (ICs) that form contacts between the protein and the ligand; it then iteratively updates these interactions during a series of biased molecular dynamics (MD) simulations to reveal the ICs that are important for the whole of the unbinding process. Subsequently, we performed finite-temperature string simulations to obtain the free-energy barrier for unbinding using the set of ICs as a complex reaction coordinate. Importantly, we also aimed to enable the further design of drugs focusing on improved residence times. To this end, we developed a supervised machine learning (ML) approach with inputs from unbiased “downhill” trajectories initiated near the transition state (TS) ensemble of the string unbinding path. We demonstrate that our ML method can identify key ligand–protein interactions driving the system through the TS. Some of the most important drugs for cancer treatment are kinase inhibitors. One of these kinase targets is cyclin-dependent kinase 2 (CDK2), an appealing target for anticancer drug development. Here, we tested our method using two different CDK2 inhibitors for the potential further development of these compounds. We compared the free-energy barriers obtained from our calculations with those observed in available experimental data. We highlighted important interactions at the distal ends of the ligands that can be targeted for improved residence times. Our method provides a new tool to determine unbinding rates and to identify key structural features of the inhibitors that can be used as starting points for novel design strategies in drug discovery.

Extraction and high-throughput sequencing of oak heartwood DNA: Assessing the feasibility of genome-wide DNA methylation profiling

Tree ring features are affected by environmental factors and therefore are the basis for dendrochronological studies to reconstruct past environmental conditions. Oak wood often provides the data for these studies because of the durability of oak heartwood and hence the availability of samples spanning long time periods of the distant past. Wood formation is regulated in part by epigenetic mechanisms such as DNA methylation. Studies of the methylation state of DNA preserved in oak heartwood thus could identify epigenetic tree ring features informing on past environmental conditions. In this study, we aimed to establish protocols for the extraction of DNA, the high-throughput sequencing of whole-genome DNA libraries (WGS) and the profiling of DNA methylation by whole-genome bisulfite sequencing (WGBS) for oak (Quercus robur) heartwood drill cores taken from the trunks of living standing trees spanning the AD 1776-2014 time period. Heartwood contains little DNA, and large amounts of phenolic compounds known to hinder the preparation of high-throughput sequencing libraries. Whole-genome and DNA methylome library preparation and sequencing consistently failed for oak heartwood samples more than 100 and 50 years of age, respectively. DNA fragmentation increased with sample age and was exacerbated by the additional bisulfite treatment step during methylome library preparation. Relative coverage of the non-repetitive portion of the oak genome was sparse. These results suggest that quantitative methylome studies of oak hardwood will likely be limited to relatively recent samples and will require a high sequencing depth to achieve sufficient genome coverage.

Photopharmacology of Ion Channels through the Light of the Computational Microscope

The optical control and investigation of neuronal activity can be achieved and carried out with photoswitchable ligands. Such compounds are designed in a modular fashion, combining a known ligand of the target protein and a photochromic group, as well as an additional electrophilic group for tethered ligands. Such a design strategy can be optimized by including structural data. In addition to experimental structures, computational methods (such as homology modeling, molecular docking, molecular dynamics and enhanced sampling techniques) can provide structural insights to guide photoswitch design and to understand the observed light-regulated effects. This review discusses the application of such structure-based computational methods to photoswitchable ligands targeting voltage- and ligand-gated ion channels. Structural mapping may help identify residues near the ligand binding pocket amenable for mutagenesis and covalent attachment. Modeling of the target protein in a complex with the photoswitchable ligand can shed light on the different activities of the two photoswitch isomers and the effect of site-directed mutations on photoswitch binding, as well as ion channel subtype selectivity. The examples presented here show how the integration of computational modeling with experimental data can greatly facilitate photoswitchable ligand design and optimization. Recent advances in structural biology, both experimental and computational, are expected to further strengthen this rational photopharmacology approach.

Durability of Dolutegravir-Based Regimens: A 5-Year Prospective Observational Study

This study evaluates the frequency and causes of dolutegravir (DTG) discontinuation along 5 years of follow-up, in both antiretroviral treatment (ART)-naive and experienced people living with HIV (PLWH). This is a prospective multi-center cohort study enrolling PLWH on DTG from July 2014 until November 2020. DTG-durability was investigated using the Kaplan-Meier survival curve. The Cox proportional-hazards model was used for estimating the hazard ratio (HR) of DTG discontinuation for any cause, and for adverse events (AEs). Nine hundred sixty-three PLWH were included, 25.3% were women and 28.0% were ART-naive. Discontinuations for any causes were 10.1 [95% confidence interval (95% CI) 8.9-11.5] per 100 person-years, similar in most regimens, with the apparent exception of tenofovir alafenamide/emtricitabine+DTG (p < 0.0001). In the multivariable Cox regression model, non-Caucasian ethnicity, age ≥50 years, and lower estimated glomerular filtration rate (eGFR) were associated with a higher probability of DTG interruption. The incidence rate of virological failure was 0.4 (95% CI 0.2-0.7) per 100 person-years, while the estimated discontinuation rate for AEs was 4.0 (3.2-4.9) per 100 person-years. Thirty-four DTG interruptions were due to grade ≥3 events (10 central nervous system, 6 hypersensitivity, 3 renal, 3 myalgia/asthenia, 3 abdominal pain, 2 gastrointestinal, and 7 other events). People with lower body mass index, age ≥50 years, and lower eGFR were at higher risk of AEs, while dual combinations were protective (HR 0.41 compared with abacavir/lamivudine/DTG, 95% CI 0.22-0.77). In this prospective observational study, we found high DTG durability and a low rate of virological failures. Dual therapies seemed protective toward AEs and might be considered, when feasible, a suitable option to minimize drug interactions and improve tolerability.

Cyclopropyl Substituents Transform the Viscosity-Sensitive BODIPY Molecular Rotor into a Temperature Sensor

A quantitative fluorescent probe that responds to changes in temperature is highly desirable for studies of biological environments, particularly in cellulo. Here, we report new cell-permeable fluorescence probes based on the BODIPY moiety that respond to environmental temperature. The new probes were developed on the basis of a well-established BODIPY-based viscosity probe by functionalization with cyclopropyl substituents at α and β positions of the BODIPY core. In contrast to the parent BODIPY fluorophore, α-cyclopropyl-substituted fluorophore displays temperature-dependent time-resolved fluorescence decays showing greatly diminished viscosity dependence, making it an attractive sensor to be used with fluorescence lifetime imaging microscopy (FLIM). We performed theoretical calculations that help rationalize the effect of the cyclopropyl substituents on the photophysical behavior of the new BODIPYs. In summary, we designed an attractive new quantitative FLIM-based temperature probe that can be used for temperature sensing in live cells.

A Single Mutation in the Outer Lipid-Facing Helix of a Pentameric Ligand-Gated Ion Channel Affects Channel Function Through a Radially-Propagating Mechanism

Pentameric ligand-gated ion channels (pLGICs) mediate fast synaptic transmission and are crucial drug targets. Their gating mechanism is triggered by ligand binding in the extracellular domain that culminates in the opening of a hydrophobic gate in the transmembrane domain. This domain is made of four α-helices (M1 to M4). Recently the outer lipid-facing helix (M4) has been shown to be key to receptor function, however its role in channel opening is still poorly understood. It could act through its neighboring helices (M1/M3), or via the M4 tip interacting with the pivotal Cys-loop in the extracellular domain. Mutation of a single M4 tyrosine (Y441) to alanine renders one pLGIC-the 5-HT3A receptor-unable to function despite robust ligand binding. Using Y441A as a proxy for M4 function, we here predict likely paths of Y441 action using molecular dynamics, and test these predictions with functional assays of mutant receptors in HEK cells and Xenopus oocytes using fluorescent membrane potential sensitive dye and two-electrode voltage clamp respectively. We show that Y441 does not act via the M4 tip or Cys-loop, but instead connects radially through M1 to a residue near the ion channel hydrophobic gate on the pore-lining helix M2. This demonstrates the active role of the M4 helix in channel opening.

Proline isomerization effects in the amyloidogenic protein β2-microglobulin

The protein β₂-microglobulin can aggregate in insoluble amyloid fibrils. By relying on extensive sampling simulations, we study the Pro32 isomerization as a possible triggering factor leading to structural modifications in β₂-m.

Time-Resolved Fluorescence Anisotropy of a Molecular Rotor Resolves Microscopic Viscosity Parameters in Complex Environments

Understanding viscosity in complex environments remains a largely unanswered question despite its importance in determining reaction rates in vivo. Here, time-resolved fluorescence anisotropy imaging (TR-FAIM) is combined with fluorescent molecular rotors (FMRs) to simultaneously determine two non-equivalent viscosity-related parameters in complex heterogeneous environments. The parameters, FMR rotational correlation time and lifetime, are extracted from fluorescence anisotropy decays, which in heterogeneous environments show dip-and-rise behavior due to multiple dye populations. Decays of this kind are found both in artificially constructed adiposomes and in live cell lipid droplet organelles. Molecular dynamics simulations are used to assign each population to nano-environments within the lipid systems. The less viscous population corresponds to the state showing an average 25° tilt to the lipid membrane normal, and the more viscous population to the state showing an average 55° tilt. This combined experimental and simulation approach enables a comprehensive description of the FMR probe behavior within viscous nano-environments in complex, biological systems.

Mutagenesis computer experiments in pentameric ligand-gated ion channels: the role of simulation tools with different resolution

Pentameric ligand-gated ion channels (pLGICs) are an important class of widely expressed membrane neuroreceptors, which play a crucial role in fast synaptic communications and are involved in several neurological conditions. They are activated by the binding of neurotransmitters, which trigger the transmission of an electrical signal via facilitated ion flux. They can also be activated, inhibited or modulated by a number of drugs. Mutagenesis electrophysiology experiments, with natural or unnatural amino acids, have provided a large body of functional data that, together with emerging structural information from X-ray spectroscopy and cryo-electron microscopy, are helping unravel the complex working mechanisms of these neuroreceptors. Computer simulations are complementing these mutagenesis experiments, with insights at various levels of accuracy and resolution. Here, we review how a selection of computational tools, including first principles methods, classical molecular dynamics and enhanced sampling techniques, are contributing to construct a picture of how pLGICs function and can be pharmacologically targeted to treat the disorders they are responsible for.

Trans- Cis Proline Switches in a Pentameric Ligand-Gated Ion Channel: How They Are Affected by and How They Affect the Biomolecular Environment

Pentameric ligand-gated ion channels (pLGICs) are important neuroreceptors, embedded in neuronal membranes, that mediate fast synaptic transmission. The molecular details of their working mechanisms have still to be fully unravelled due to their complexity and limited structural information available. Here we focus on a potential molecular switch in a prototypical pLGIC, the serotonin-activated 5-HT₃ receptor, consisting of the trans- cis isomerization of a proline at the interface between the extracellular and transmembrane domain. Mutagenesis electrophysiology experiments previously showed that if such isomerization could not take place, the channel would not open, but the hypothetical role of this mechanism as key to channel gating is still debated. We investigate this switch within the receptor with molecular dynamics and enhanced sampling simulations. We analyze how the isomerization free energy landscape is affected by the receptor environment in comparison to simplified models. Moreover, we reveal how the isomerization, in turn, affects the structural and electrostatic properties of the receptor at the extracellular-transmembrane domain interface, e.g., by tuning the ion selectivity filter.

Unravelling the Roles of Size, Ligands, and Pressure in the Piezochromic Properties of CdS Nanocrystals

Understanding the effects of pressure-induced deformations on the optoelectronic properties of nanomaterials is important not only from the fundamental point of view but also for potential applications such as stress sensors and electromechanical devices. Here, we describe the novel insights into these piezochromic effects gained from using a linear-scaling density functional theory framework and an electronic enthalpy scheme, which allow us to accurately characterize the electronic structure of CdS nanocrystals with a zincblende-like core of experimentally relevant size. In particular, we focus on unravelling the complex interplay of size and surface (phenyl) ligands with pressure. We show that pressure-induced deformations are not simple isotropic scaling of the original structures and that the change in HOMO-LUMO gap with pressure results from two competing factors: (i) a bulk-like linear increase due to compression, which is offset by (ii) distortions and disorder and, to a lesser extent, orbital hybridization induced by ligands affecting the frontier orbitals. Moreover, we observe that the main peak in the optical absorption spectra is systematically red-shifted or blue-shifted, as pressure is increased up to 5 GPa, depending on the presence or absence of phenyl ligands. These heavily hybridize the frontier orbitals, causing a reduction in overlap and oscillator strength, so that at zero pressure, the lowest energy transition involves deeper hole orbitals than in the case of hydrogen-capped nanocrystals; the application of pressure induces greater delocalization over the whole nanocrystals bringing the frontier hole orbitals into play and resulting in an unexpected red shift for the phenyl-capped nanocrystals, in part caused by distortions. In response to a growing interest in relatively small nanocrystals that can be difficult to accurately characterize with experimental techniques, this work exemplifies the detailed understanding of structure-property relationships under pressure that can be obtained for realistic nanocrystals with state-of-the-art first-principles methods and used for the characterization and design of devices based on these and similar nanomaterials.

The Free Energy Landscape of GABA Binding to a Pentameric Ligand-Gated Ion Channel and Its Disruption by Mutations

Pentameric ligand-gated ion channels (pLGICs) of the Cys-loop superfamily are important neuroreceptors that mediate fast synaptic transmission. They are activated by the binding of a neurotransmitter, but the details of this process are still not fully understood. As a prototypical pLGIC, here we choose the insect resistance to dieldrin (RDL) receptor involved in resistance to insecticides and investigate the binding of the neurotransmitter GABA to its extracellular domain at the atomistic level. We achieve this by means of μ-sec funnel-metadynamics simulations, which efficiently enhance the sampling of bound and unbound states by using a funnel-shaped restraining potential to limit the exploration in the solvent. We reveal the sequence of events in the binding process from the capture of GABA from the solvent to its pinning between the charged residues Arg111 and Glu204 in the binding pocket. We characterize the associated free energy landscapes in the wild-type RDL receptor and in two mutant forms, where the key residues Arg111 and Glu204 are mutated to Ala. Experimentally these mutations produce nonfunctional channels, which is reflected in the reduced ligand binding affinities due to the loss of essential interactions. We also analyze the dynamical behavior of the crucial loop C, whose opening allows the access of GABA to the binding site and closure locks the ligand into the protein. The RDL receptor shares structural and functional features with other pLGICs; hence, our work outlines a valuable protocol to study the binding of ligands to pLGICs beyond conventional docking and molecular dynamics techniques.

Pressure-Induced Amorphization and a New High Density Amorphous Metallic Phase in Matrix-Free Ge Nanoparticles

Over the last two decades, it has been demonstrated that size effects have significant consequences for the atomic arrangements and phase behavior of matter under extreme pressure. Furthermore, it has been shown that an understanding of how size affects critical pressure-temperature conditions provides vital guidance in the search for materials with novel properties. Here, we report on the remarkable behavior of small (under ~5 nm) matrix-free Ge nanoparticles under hydrostatic compression that is drastically different from both larger nanoparticles and bulk Ge. We discover that the application of pressure drives surface-induced amorphization leading to Ge-Ge bond overcompression and eventually to a polyamorphic semiconductor-to-metal transformation. A combination of spectroscopic techniques together with ab initio simulations were employed to reveal the details of the transformation mechanism into a new high density phase-amorphous metallic Ge.

Insights into the binding of GABA to the insect RDL receptor from atomistic simulations: a comparison of models

The resistance to dieldrin (RDL) receptor is an insect pentameric ligand-gated ion channel (pLGIC). It is activated by the neurotransmitter γ-aminobutyric acid (GABA) binding to its extracellular domain; hence elucidating the atomistic details of this interaction is important for understanding how the RDL receptor functions. As no high resolution structures are currently available, we built homology models of the extracellular domain of the RDL receptor using different templates, including the widely used acetylcholine binding protein and two pLGICs, the Erwinia Chrysanthemi ligand-gated ion channel (ELIC) and the more recently resolved GluCl. We then docked GABA into the selected three dimensional structures, which we used as starting points for classical molecular dynamics simulations. This allowed us to analyze in detail the behavior of GABA in the binding sites, including the hydrogen bond and cation-π interaction networks it formed, the conformers it visited and the possible role of water molecules in mediating the interactions; we also estimated the binding free energies. The models were all stable and showed common features, including interactions consistent with experimental data and similar to other pLGICs; differences could be attributed to the quality of the models, which increases with increasing sequence identity, and the use of a pLGIC template. We supplemented the molecular dynamics information with metadynamics, a rare event method, by exploring the free energy landscape of GABA binding to the RDL receptor. Overall, we show that the GluCl template provided the best models. GABA forming direct salt-bridges with Arg211 and Glu204, and cation-π interactions with an aromatic cage including Tyr109, Phe206 and Tyr254, represents a favorable binding arrangement, and the interaction with Glu204 can also be mediated by a water molecule.

A computational exploration of the interactions of the green tea polyphenol (-)-Epigallocatechin 3-Gallate with cardiac muscle troponin C

Thanks to its polyphenols and phytochemicals, green tea is believed to have a number of health benefits, including protecting from heart disease, but its mechanism of action at the molecular level is still not understood. Here we explore, by means of atomistic simulations, how the most abundant of the green tea polyphenols, (-)-Epigallocatechin 3-Gallate (EGCg), interacts with the structural C terminal domain of cardiac muscle troponin C (cCTnC), a calcium binding protein that plays an important role in heart contractions. We find that EGCg favourably binds to the hydrophobic cleft of cCTnC consistently with solution NMR experiments. It also binds to cCTnC in the presence of the anchoring region of troponin I (cTnI(34-71)) at the interface between the E and H helices. This appears to affect the strength of the interaction between cCTnC and cTnI(34-71) and also counter-acts the effects of the Gly159Asp mutation, related to dilated cardiomyopathy. Our simulations support the picture that EGCg interacting with the C terminal domain of troponin C may help in regulating the calcium signalling either through competitive binding with the anchoring domain of cTnI or by affecting the interaction between cCTnC and cTnI(34-71).

GABA binding to an insect GABA receptor: a molecular dynamics and mutagenesis study

RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-π interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.

Interacting Lewis-X carbohydrates in condensed phase: a first-principles molecular dynamics study

We performed first-principles molecular dynamics calculations at finite temperature, to study the interacting conformations of Lewis-X (LeX) trisaccharides in the crystalline phase. The calculated cell parameters and detailed atomic structure of the LeX molecule compare well to the experimental data obtained by X-ray diffraction. We identify and characterize the hydrogen-bond network, responsible for the mutual interaction of the LeX pairs, whereas we find the intramolecular conformation and stability to be mainly assured by dispersion forces. The relative contributions to the crystallization energy of the hydrogen bonds and of the dispersion forces are defined and quantified. From this study, candidate configurations for the fully hydrated, in vivo structures of homotypic LeX-LeX interactions at cell surfaces can be proposed. We discuss how these configurations could also be relevant for the adhesion and self-assembly of nanostructures.

Rapid desensitization of the rat α7 nAChR is facilitated by the presence of a proline residue in the outer β-sheet

The rat α7 nicotinic acetylcholine receptor (nAChR) has a proline residue near the middle of the β9 strand. The replacement of this proline residue at position 180 (P180) by either threonine (α7-P180T) or serine (α7-P180S) slowed the onset of desensitization dramatically, with half-times of ~930 and 700 ms, respectively, compared to 90 ms for the wild-type receptor. To investigate the importance of the hydroxyl group on the position 180 side-chains, the mutant receptors α7-P180Y and α7-P180F were studied and showed half-times of desensitization of 650 and 160 ms, respectively. While a position 180 side-chain OH group may contribute to the slow desensitization rates, α7-P180S and α7-P180V resulted in receptors with similar desensitization rates, suggesting that increased backbone to backbone H bonding expected in the absence of proline at position 180 would likely exert a great effect on desensitization. Single channel recordings indicated that for the α7-P180T receptor there was a significantly reduced closed time without any change in single channel conductance (as compared to wild-type). Kinetic simulations indicated that all changes observed for the mutant channel behaviour were reproduced by decreasing the rate of desensitization, and increasing the microscopic affinity to resting receptors. Molecular dynamics (MD) simulations on a homology model were used to provide insight into likely H bond interactions within the outer β-sheet that occur when the P180 residue is mutated. All mutations analysed increased about twofold the predicted number of H bonds between the residue at position 180 and the backbone of the β10 strand. Moreover, the α7-P180T and α7-P180S mutations also formed some intrastrand H bonds along the β9 strand, although H bonding of the OH groups of the threonine or serine side-chains was predicted to be infrequent. Our results indicate that rapid desensitization of the wild-type rat α7 nAChR is facilitated by the presence of the proline residue within the β9 strand.

Modeling the Excited States of Biological Chromophores within Many-Body Green's Function Theory

First-principle many-body Green's function theory (MBGFT) has been successfully used to describe electronic excitations in many materials, from bulk crystals to nanoparticles. Here we assess its performance for the calculations of the excited states of biological chromophores. MBGFT is based on a set of Green's function equations, whose key ingredients are the electron's self-energy Σ, which is obtained by Hedin's GW approach, and the electron-hole interaction, which is described by the Bethe-Salpeter equation (BSE). The GW approach and the BSE predict orbital energies and excitation energies with high accuracy, respectively. We have calculated the low-lying excited states of a series of model biological chromophores, related to the photoactive yellow protein (PYP), rhodopsin, and the green fluorescent protein (GFP), obtaining a very good agreement with the available experimental and accurate theoretical data; the order of the excited states is also correctly predicted. MBGFT bridges the gap between time-dependent density functional theory and high-level quantum chemistry methods, combining the efficiency of the former with the accuracy of the latter: this makes MBGFT a promising method for studying excitations in complex biological systems.

Trans-cis switching mechanisms in proline analogues and their relevance for the gating of the 5-HT3 receptor

Trans-cis isomerization of a proline peptide bond is a potential mechanism to open the channel of the 5-HT(3) receptor. Here, we have used the metadynamics method to theoretically explore such a mechanism. We have determined the free energy surfaces in aqueous solution of a series of dipeptides of proline analogues and evaluated the free energy difference between the cis and trans isomers. These theoretical results were then compared with data from mutagenesis experiments, in which the response of the 5-HT(3) receptor was measured when the proline at the apex of the M2-M3 transmembrane domain loop was mutated. The strong correlation between the experimental and the theoretical data supports the existence of a trans-cis proline switch for opening the 5-HT(3) receptor ion channel.

Exploring the binding of serotonin to the 5-HT3 receptor by density functional theory

The 5-HT3 receptor is a typical ligand-gated ion channel of the Cys-loop superfamily, which is activated by binding of serotonin (5-HT). Models of the binding site of this protein reveal potential interactions between 5-HT and Tyr143, Tyr153, and Tyr234. Here we describe a series of ab initio calculations, based on density functional theory, to assess the effects of mutating these tyrosine residues on the binding of 5-HT. A series of mutations to these tyrosines, previously studied experimentally, were tested, and the binding energies compared with the available experimental data. Our results show that Tyr153 could form a hydrogen bond with the tertiary amine of 5-HT, and that mutation in this location revealed binding energies broadly in line with experimentally determined EC50s. Tyr143 could also form a hydrogen bond, but as EC50s do not relate to binding energies, it is unlikely that such a bond is formed here. Tyr234 is quite distinct in that it may interact with 5-HT via a mixed hydrogen bond/cation-pi interaction.

Ab initio studies of layering behavior of liquid sodium surfaces and interfaces

We have studied the liquid surface of sodium with extensive ab initio molecular dynamics simulations based on ensemble density-functional theory. We find clear evidence of layering in the direction perpendicular to the surface that persists to temperatures more than 100 K above the melting point. We also observe clear Friedel oscillations in the electronic density response to the presence of a surface, but their direct effect on atomic layering is ruled out. A careful finite-size effect analysis accompanies our results, showing that liquid slabs 20-25 A thick capture the essential details of the surface structure. We conclude that geometrical confinement is the common cause for layer formation, which is similar to what happens at a liquid-solid interface: at a free liquid surface, the rapid decay of the electronic density from the bulk liquid value to zero in the vapor forms a hard wall against which the atoms pack. Finally, we predict x-ray reflectivities from ab initio molecular dynamics data that include some of the large surface-normal wave vector-transfer regions that, for alkali metals, are not accessible to experiments.

Structural relaxations in electronically excited poly(para-phenylene)

Structural relaxations in electronically excited poly(para-phenylene) are studied using many-body perturbation theory and density-functional-theory methods. A sophisticated description of the electron-hole interaction is required to describe the excitonic energies, but the associated structural relaxations can be obtained quite accurately within a constrained density-functional-theory approach. We find that the structural relaxations in the low-energy excitonic states extend over about eight monomers, leading to an energy reduction of 0.22 eV and a Stokes shift of 0.40 eV.

Molecular epidemiology study of exogenous reinfection in an area with a low incidence of tuberculosis

In geographical areas with a low incidence of tuberculosis, recurrent tuberculosis is generally due to reactivation of the disease. However, the relative contribution of tuberculosis reinfection increases in parallel with the incidence of disease and is likely to depend on the epidemiological context: factors such as the spread of multidrug resistance, human immunodeficiency virus (HIV) infection, and immigration from developing countries could modify disease transmission in areas at low risk for tuberculosis. A molecular epidemiology study was performed in Lombardy, Northern Italy, where the incidence of tuberculosis is 17.5 cases per 100,000 persons. A total of 2,452 cases of culture-confirmed tuberculosis in 2,127 patients were studied. A group of 32 patients (1.5%), each of whom had two episodes of tuberculosis with cure as the outcome of the first episode and with more than 6 months between the two episodes, were studied by means of restriction fragment length polymorphism DNA fingerprinting analysis. For 5 of the 32 patients (16%), the DNA fingerprinting patterns of Mycobacterium tuberculosis strains responsible for the second episode did not match those of the corresponding isolates of the first episode, indicating exogenous reinfection. Two of these patients developed multidrug-resistant tuberculosis during the second episode, and in three cases the isolates belonged to clusters of M. tuberculosis strains spreading in the community. A fourfold-increased risk for reinfection was observed in immigrant patients compared to Italian subjects. In contrast, a higher risk of relapse rather than reinfection was evidenced in HIV-positive subjects and in patients infected with multidrug-resistant tuberculosis. Episodes of tuberculosis reinfection in areas with a low incidence of tuberculosis are rare compared to those in high-incidence geographical regions. In populations that have immigrated from high-risk areas, reinfection may represent a considerable contributor to the rate of recurrent tuberculosis. This finding emphasizes the importance of containing the spread of epidemic strains in close communities, in order to prevent changes in global tuberculosis trends for developed countries.

PCR-hybridization assay for Mycobacterium avium complex: optimization of detection in peripheral blood from humans

We evaluated the sensitivity of a DNA amplification test for the detection of Mycobacterium avium in blood samples using different blood components and different DNA extraction methods. M. avium-inoculated blood samples were processed to obtain separate blood components: peripheral blood mononuclear cells (PBMCs), polymorphonuclear cells (PMNCs), and whole-blood sodium dodecyl sulfate (SDS)-lysate pellets. The sensitivity for the detection of the lowest mycobacterial load (1 CFU/ml) was significantly greater (P < 0.01) with DNA extracted from SDS-lysate pellets than with DNA extracted from PBMCs or PMNCs. Subsequently, DNA extraction methods based on guanidine NaOH, and proteinase were compared. The sensitivity of the guanidine-based method was significantly greater (P < 0.01) than those of the others.

Ab initio molecular dynamics with a classical pressure reservoir: simulation of pressure-induced amorphization in a Si35H36 cluster

We present a new constant-pressure ab initio molecular dynamics method suitable for studying, e.g., pressure-induced structural transformations in finite nonperiodic systems such as clusters. We immerse an ab initio treated cluster into a model classical liquid, described by a soft-sphere potential, which acts as a pressure reservoir. The pressure is varied by tuning the parameter of the liquid potential. We apply the method to a Si35H36 cluster, which undergoes a pressure-induced amorphization at approximately 35 GPa, and remains in a disordered state even upon pressure release.