Search

Climate change is expected to significantly affect the interior climate of old, leaky buildings without HVAC systems. As a result, the items of cultural significance that are hosted indoors will experience new ambient conditions, which will affect their degradation. In the current research, the impact of climate change on the biological, mechanical, and chemical degradation of a cabinet and a storage trunk which are made of wood and have paintings on their outer surface is investigated. These two items are found in two different rooms of a historic timber building in Vestfold, Norway. Data from the REMO2015 driven by the global model MPI-ESM-LR are used in order to account for past, present, and future climate conditions. In addition, climate data from ERA5 reanalysis are used in order to assess the accuracy of the MPI-ES-LR_REMO2015 model results. Whole building hygrothermal simulations are employed to calculate the temperature and the relative humidity inside the rooms that host the items of interest. The transient hygrothermal condition and certain characteristics of the timber surfaces are used as inputs in models that describe their degradation. The biological degradation is examined by using i) the updated VTT mould model and ii) the Growing Degree Days (GDD) for temperature and humidity dependant insects. The mechanical deterioration is assessed by the method proposed by Mecklenburg et al. (1998). The concept of the Lifetime Multiplier (LM) is used in order to assess the chemical deterioration of the furniture. Results reveal a significant mechanical degradation risk and a very high chemical deterioration risk. The biodeterioration risk remains at acceptable levels. Moreover, it could be possible that the storage trunk would be damaged by certain insects in the future. It is then suggested that both items should be moved to a room with proper conditions in order to minimize their chemical and mechanical deterioration risk and extend their life span. Finally, the significance of implementing bias correction in the data from climate models is underlined.

International audience; Theoretical work suggests that collective spatiotemporal behaviour of integral membrane proteins (IMPs) should be modulated by boundary lipids sheathing their membrane anchors. Here, we show evidence for this prediction whilst investigating the mechanism for maintaining a steady amount of the active form of IMP Lck kinase (LckA) by Lck trans-autophosphorylation regulated by the phosphatase CD45. We used super-resolution microscopy, flow cytometry, and pharmacological and genetic perturbation to gain insight into the spatiotemporal context of this process. We found that LckA is generated exclusively at the plasma membrane, where CD45 maintains it in a ceaseless dynamic equilibrium with its unphosphorylated precursor. Steady LckA shows linear dependence, after an initial threshold, over a considerable range of Lck expression levels. This behaviour fits a phenomenological model of trans-autophosphorylation that becomes more efficient with increasing LckA. We then challenged steady LckA formation by genetically swapping the Lck membrane anchor with structurally divergent ones, such as that of Src or the transmembrane domains of LAT, CD4, palmitoylation-defective CD4 and CD45 that were expected to drastically modify Lck boundary lipids. We observed small but significant changes in LckA generation, except for the CD45 transmembrane domain that drastically reduced LckA due to its excessive lateral proximity to CD45. Comprehensively, LckA formation and maintenance can be best explained by lipid bilayer critical density fluctuations rather than liquid-ordered phase-separated nanodomains, as previously thought, with "like/unlike" boundary lipids driving dynamical proximity and remoteness of Lck with itself and with CD45.

The COVID-19 pandemic is caused by SARS-CoV-2. Currently, most of the research efforts towards the development of vaccines and antibodies against SARS-CoV-2 were mainly focused on the spike (S) protein, which mediates virus entry into the host cell by binding to ACE2. As the virus SARS-CoV-2 continues to spread globally, variants have emerged, characterized by multiple mutations of the S glycoprotein. Herein, we employed microsecond-long molecular dynamics simulations to study the impact of the mutations of the S glycoprotein in SARS-CoV-2 Variant of Concern 202012/01 (B.1.1.7), termed the “UK variant”, in comparison with the wild type, with the aim to decipher the structural basis of the reported increased infectivity and virulence. The simulations provided insights on the different dynamics of UK and wild-type S glycoprotein, regarding in particular the Receptor Binding Domain (RBD). In addition, we investigated the role of glycans in modulating the conformational transitions of the RBD. The overall results showed that the UK mutant experiences higher flexibility in the RBD with respect to wild type, this behavior might be correlated with the increased transmission reported for this variant. Our work also adds useful structural information on antigenic “hotspots” and epitopes targeted by neutralizing antibodies.

Background and AimsThe aim of this study was to apply machine learning models and a nomogram to differentiate critically ill from non-critically ill COVID-19 pneumonia patients.MethodsClinical symptoms and signs, laboratory parameters, cytokine profile, and immune cellular data of 63 COVID-19 pneumonia patients were retrospectively reviewed. Outcomes were followed up until Mar 12, 2020. A logistic regression function (LR model), Random Forest, and XGBoost models were developed. The performance of these models was measured by area under receiver operating characteristic curve (AUC) analysis.ResultsUnivariate analysis revealed that there was a difference between critically and non-critically ill patients with respect to levels of interleukin-6, interleukin-10, T cells, CD4+ T, and CD8+ T cells. Interleukin-10 with an AUC of 0.86 was most useful predictor of critically ill patients with COVID-19 pneumonia. Ten variables (respiratory rate, neutrophil counts, aspartate transaminase, albumin, serum procalcitonin, D-dimer and B-type natriuretic peptide, CD4+ T cells, interleukin-6 and interleukin-10) were used as candidate predictors for LR model, Random Forest (RF) and XGBoost model application. The coefficients from LR model were utilized to build a nomogram. RF and XGBoost methods suggested that Interleukin-10 and interleukin-6 were the most important variables for severity of illness prediction. The mean AUC for LR, RF, and XGBoost model were 0.91, 0.89, and 0.93 respectively (in two-fold cross-validation). Individualized prediction by XGBoost model was explained by local interpretable model-agnostic explanations (LIME) plot.ConclusionsXGBoost exhibited the highest discriminatory performance for prediction of critically ill patients with COVID-19 pneumonia. It is inferred that the nomogram and visualized interpretation with LIME plot could be useful in the clinical setting. Additionally, interleukin-10 could serve as a useful predictor of critically ill patients with COVID-19 pneumonia.

Acoustic emission is a well-known noninvasive methodology for the study of defects in materials but still rarely applied in the field of cultural heritage diagnostic. How alteration products and degradation processes affect the acoustic emission signal still is an open issue. The proposed study concerns the utilization of such techniques to investigate the relations existing between the moisture content and the typology and amount of crystallized salts with the acoustic emission features. This work focuses on spruce wood logs belonging to an historical warehouse in Trondheim, Norway, since no strictly scientific studies exist on the conservation conditions of these big wooden structures. The methodology, also involving the moisture content measurement of structures and the samplings of portions analyzed through vacuum microbalance, allowed identifying a clear relationship between the amount of water in logs as a function of their distance from the ground and variations in the amplitude of the acoustic emission signals.

The thyroid-stimulating hormone receptor (TSHR) has been indicated as a putative domestication gene in chicken. Comparison of WGS identified a variant in residue 558 of the transmembrane domain (TM) of TSHR, where the domestic chicken (GGD) presents an arginine, whereas the red jungle fowl (RJF) shares a conserved glycine with other vertebrates. This variant has been demonstrated to be associated with phenotypes that are important for domestication and related to thyroid regulation, such as less fearful behavior, reduced aggressive behavior and reduced dependence on seasonal reproduction in GGD as compared with RJF. By means of molecular dynamics simulations, we highlighted the structural and dynamic differences of variant Gly558Arg in the TSHR TM domain. Alterations in TM helix flexibility, structure and protein overall motion are described. The so-called 'arginine snorkeling' of residue 568 in GGD is observed and we hypothesize it as the originating force that produces the observed whole-protein perturbation in the helix bundle dynamics, capable of altering the TSHR signal transduction. The results are discussed in the context of their implications for a better understanding of biological mechanisms in chicken under control of the thyroid, such as body metabolism, as well as for their usefulness in biomedical research.

Given the enormous social and health impact of the pandemic triggered by severe acute respiratory syndrome 2 (SARS-CoV-2), the scientific community made a huge effort to provide an immediate response to the challenges posed by Coronavirus disease 2019 (COVID-19). One of the most important proteins of the virus is an enzyme, called 3CLpro or main protease, already identified as an important pharmacological target also in SARS and Middle East respiratory syndrome virus (MERS) viruses. This protein triggers the production of a whole series of enzymes necessary for the virus to carry out its replicating and infectious activities. Therefore, it is crucial to gain a deeper understanding of 3CLpro structure and function in order to effectively target this enzyme. All-atoms molecular dynamics (MD) simulations were performed to examine the different conformational behaviors of the monomeric and dimeric form of SARS-CoV-2 3CLpro apo structure, as revealed by microsecond time scale MD simulations. Our results also shed light on the conformational dynamics of the loop regions at the entry of the catalytic site. Studying, at atomic level, the characteristics of the active site and obtaining information on how the protein can interact with its substrates will allow the design of molecules able to block the enzymatic function crucial for the virus.

The conformational behavior of the wild-type endonucleases I-DmoI and two of its mutants has been studied in the presence and in the absence of DNA target sequences by means of extended molecular dynamics simulations. Our results show that in the absence of DNA, the three protein forms explore a similar essential conformational space, whereas when bound to the same DNA target sequence of 25 base pairs, they diversify and restrain the subspace explored. In addition, the differences in the essential subspaces explored by the residues near the catalytic site for both the bound and unbound forms are discussed in background of the experimental protein activity.

PPARγ represents a key target for the treatment of type II diabetes and metabolic syndrome. Synthetic antidiabetic drugs activating PPARγ are accompanied by serious undesirable side effects related to their agonism. In the search for new PPARγ regulators, inhibitors of PPARγ phosphorylation mediated by CDK5 represent an opportunity for the development of an improved generation of anti-diabetic drugs acting through this nuclear receptor. We have employed a multi-disciplinary approach, including protein-protein docking, NMR, HDX, MD simulations and site-directed mutagenesis to investigate conformational changes in PPARγ that impair the ability of CDK5 to interact with PPARγ and hence inhibit PPARγ phosphorylation. Finally, we describe an alternate inhibition mechanism adopted by ligands bound far from the phosphorylation site.

We present and characterize a novel mutation in KCNJ10. Unlike previously reported EAST/SeSAME patients, our patient was heterozygous, and contrary to previous studies, mimicking the heterozygous state by coexpression resulted in loss of channel function. We report in the same patient co-occurrence of a KCNT1 mutation resulting in a more severe phenotype. This study provides new insights into the phenotypic spectrum and to the genotype-phenotype correlations associated with EAST/SeSAME and MMFSI.

The tau protein belongs to the category of intrinsically disordered proteins, which in their native state do not have an average stable structure and fluctuate between many conformations. In its physiological state, tau helps nucleating and stabilising the microtubules in the axons of the neurons. On the other hand, the same tau is involved in the development of Alzheimer disease, when it aggregates in paired helical filaments forming fibrils, which form insoluble tangles. The beginning of the pathological aggregation of tau has been attributed to a local transition of protein portions from random coil to a β-sheet. These structures would very likely be transient; therefore, we performed a molecular dynamics simulation of tau to gather information on the existence of segments of tau endowed with a secondary structure. We combined the results of our simulation with small-angle X-ray scattering experimental data to extract from the dynamics a set of most probable conformations of tau. The analysis of these ...

Episodic ataxia type 1 (EA1) is an autosomal dominant disorder characterized by continuous myokymia and episodic attacks of ataxia. Mutations in the gene KCNA1 that encodes the voltage-gated potassium channel Kv1.1 are responsible for EA1. In several brain areas, Kv1.1 coassembles with Kv1.4, which confers N-type inactivating properties to heteromeric channels. It is therefore likely that the rate of inactivation will be determined by the number of Kv1.4 inactivation particles, as set by the precise subunit stoichiometry. We propose that EA1 mutations affect the rate of N-type inactivation either by reduced subunit surface expression, giving rise to a reduced number of Kv1.1 subunits in heterotetramer Kv1.1-Kv1.4 channels, or by reduced affinity for the Kv1.4 inactivation domain. To test this hypothesis, quantified amounts of mRNA for Kv1.4 or Kv1.1 containing selected EA1 mutations either in the inner vestibule of Kv1.1 on S6 or in the transmembrane regions were injected into Xenopus laevis oocytes and the relative rates of inactivation and stoichiometry were determined. The S6 mutations, V404I and V408A, which had normal surface expression, reduced the rate of inactivation by a decreased affinity for the inactivation domain while the mutations I177N in S1 and E325D in S5, which had reduced subunit surface expression, increased the rate of N-type inactivation due to a stoichiometric increase in the number of Kv1.4 subunits.

The conformational behavior of the wild-type endonucleases I-DmoI and two of its mutants has been studied in the presence and in the absence of DNA target sequences by means of extended molecular dynamics simulations. Our results show that in the absence of DNA, the three protein forms explore a similar essential conformational space, whereas when bound to the same DNA target sequence of 25 base pairs, they diversify and restrain the subspace explored. In addition, the differences in the essential subspaces explored by the residues near the catalytic site for both the bound and unbound forms are discussed in background of the experimental protein activity.

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 x 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the 'quality' of the starting structure seems mainly to influence the simulation stability of the transmembrane beta-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal alpha-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a approximately 10 ns timescale.

Marine ladder-shaped polyether toxins are implicated in neurological symptoms of fish-borne food poisonings. The toxin gambierol, produced by the marine dinoflagellate Gambierdiscus toxicus, belongs to the group of ladder-shaped polyether toxins and inhibits Kv3.1 channels with nanomolar affinity through a mechanism of gating modification. Binding determinants for gambierol localize at the lipid-exposed interface of the pore forming S5 and S6 segments, suggesting that gambierol binds outside of the permeation pathway. To explore a possible involvement of the voltage-sensing domain (VSD), we made different chimeric channels between Kv3.1 and Kv2.1, exchanging distinct parts of the gating machinery. Our results showed that neither the electro-mechanical coupling nor the S1-S3a region of the VSD affect gambierol sensitivity. In contrast, the S3b-S4 part of the VSD (paddle motif) decreased gambierol sensitivity in Kv3.1 more than 100-fold. Structure determination by homology modeling indicated that the position of the S3b-S4 paddle and its primary structure defines the shape and∖or the accessibility of the binding site for gambierol, explaining the observed differences in gambierol affinity between the channel chimeras. Furthermore, these findings explain the observed difference in gambierol affinity for the closed and open channel configurations of Kv3.1, opening new possibilities for exploring the VSDs as selectivity determinants in drug design.

An autosomal dominant protein aggregate myopathy, characterized by high plasma creatine kinase and calsequestrin-1 (CASQ1) accumulation in skeletal muscle, has been recently associated with a missense mutation in CASQ1 gene. The mutation replaces an evolutionarily-conserved aspartic acid with glycine at position 244 (p.D244G) of CASQ1, the main sarcoplasmic reticulum (SR) Ca2+ binding and storage protein localized at the terminal cisternae of skeletal muscle cells. Here, immunocytochemical analysis of myotubes, differentiated from muscle-derived primary myoblasts, shows that sarcoplasmic vacuolar aggregations positive for CASQ1 are significantly larger in CASQ1-mutated cells than control cells. A strong co-immuno staining of both RyR1 and CASQ1 was also noted in the vacuoles of myotubes and muscle biopsies derived from patients. Electrophysiological recordings and sarcoplasmic Ca2+ measurements provide evidence for less Ca2+ release from the SR of mutated myotubes when compared to that of controls. These findings further clarify the pathogenic nature of the p.D244G variant and point out defects in sarcoplasmic Ca2+ homeostasis as a mechanism underlying this human disease, which could be distinctly classified as "CASQ1-couplonopathy"., peer-reviewed

Molecular dynamics simulations may be used to probe the interactions of membrane proteins with lipids and with detergents at atomic resolution. Examples of such simulations for ion channels and for bacterial outer membrane proteins are described. Comparison of simulations of KcsA (an α-helical bundle) and OmpA (a β-barrel) reveals the importance of two classes of side chains in stabilizing interactions with the head groups of lipid molecules: (i) tryptophan and tyrosine; and (ii) arginine and lysine. Arginine residues interacting with lipid phosphate groups play an important role in stabilizing the voltage-sensor domain of the KvAP channel within a bilayer. Simulations of the bacterial potassium channel KcsA reveal specific interactions of phosphatidylglycerol with an acidic lipid-binding site at the interface between adjacent protein monomers. A combination of molecular modelling and simulation reveals a potential phosphatidylinositol 4,5-bisphosphate-binding site on the surface of Kir6.2.

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

The recently determined structure of a mammalian voltage-gated potassium channel has important implications for our understanding of voltage-sensing and gating mechanisms in channels. It is also the first crystal structure of an overexpressed eukaryotic membrane protein.

Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The β-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with “dimer-of-dimers” motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (

In this work, we investigated the structural and dynamic consequences of two substi- tutions, P58A and G36P, located in two different solvent-exposed loops of cytochrome c551. The re- sults show that both mutations affect regions that are distant from the site of mutation. Here, the two loops appear to be dynamically coupled to each other, because the substitution at one site affects the structure and the dynamics of the other site. However, the substitutions at Gly-36 and Pro-58 presented substantial differences, which were re- lated to the mechanical (rigidity and deformability) properties of the site surrounding the mutation. Although the P58A mutant conserved a significant dynamic similarity to the wild-type protein as the immediate surroundings of position 58 became more rigid, the G36P mutant, which had deformed its flexible surroundings, presented a dynamic behav- ior that was markedly different from that of the wild-type protein. These results suggest that pertur- bation of sterically isolated and flexible regions, such as solvent-exposed loops, can have strong dy- namic consequences on the protein as a whole, raising the possibility that these effects could in turn affect the stability or the function of the pro- tein. Proteins 2003;50:222-229. © 2002 Wiley-Liss, Inc.

In recent years, increased interest in the origin of protein thermal stability has gained attention both for its possible role in understanding the forces governing the folding of a protein and for the design of new highly stable engineered biocatalysts. To study the origin of thermostability, we have performed molecular dynamics simulations of two rubredoxins, from the mesophile Clostridium pasteurianum and from the hyperthermophile Pyrococcus furiosus. The simulations were carried out at two temperatures, 300 and 373 K, for each molecule. The length of the simulations was within the range of 6-7.2 ns. The rubredoxin from the hyperthermophilic organism was more flexible than its mesophilic counterpart at both temperatures; however, the overall flexibility of both molecules at their optimal growth temperature was the same, despite 59% sequence homology. The conformational space sampled by both molecules was larger at 300 K than at 373 K. The essential dynamics analysis showed that the principal overall motions of the two molecules are significantly different. On the contrary, each molecule showed similar directions of motion at both temperatures.

This chapter focuses on homology modeling and molecular dynamics (MD) simulations studies of ion channels for the range of single cell organisms from prokaryotes to eukaryotes. The glutamate receptor channels (GluRs) share some distant homology in their transmembrane (TM) domains with K channels but possess distinct extracellular ligand-binding domains for which several structures, of both bacterial and mammalian homologs, are known. It can be seen that molecular modeling and simulations can contribute to studies of ion channels in two respects. Modeling studies enable extrapolation from experimental structures of prokaryotic ion channels to molecular models of eukaryotic homologs, thus aiding design and interpretation of, for example, mutation experiments for dissecting structure-function relationships. Ion channel structures and ion channel models may also be used as the basis of multinanosecond MD simulations. Finally, it will become increasingly important to run multiple simulations on multiple channels to allow comparative analysis of simulation results, which in turn will enable the formulation of more general hypotheses concerning the relationship between the conformational dynamics of channel proteins and their physiological functions.

From a specific definition of the roto-translational (external) and intramolecular (internal) coordinates, a constrained dynamics algorithm is derived for removing the roto-translational motions during molecular dynamics simulations, within the leap-frog integration scheme. In the paper the theoretical basis of this new method and its statistical mechanical consistency are reported, together with two applications.

The salivary antimicrobial peptide histatin-5 is able to aggregate and fuse negatively charged small unilamellar vesicles, and this fusogenic activity is selectively induced by the presence of zinc ions. Circular dichroism spectroscopy shows that histatin-5, in the presence of negatively charged vesicles and zinc ions, undergoes a conformational change leading to the stabilization of an alpha-helical secondary structure. We attribute the specific action of the zinc ions to the presence of a consensus sequence, HEXXH, located in the C-terminal functional domain of histatin-5, a recognized zinc-binding motif in many proteins. Two-dimensional proton NMR spectroscopy of histatin-5 in a trifluoroethanol/water mixture (a membrane mimetic environment) has been performed and the results analyzed by means of distance geometry and restrained molecular dynamics simulations. Our results reveal that the peptide chain, including the Zn-binding consensus sequence corresponding to residues 15-19, is in a helicoidal conformation. Comparison of the chemical shifts of the individual amino acids in histatin-5 with those recently reported in other solvents indicates that trifluoroethanol/water has a structuring capability somewhere between water and dimethyl sulfoxide. The mechanism of action of this antimicrobial peptide is discussed on the basis of its structural characteristics with particular attention to the Zn-binding motif.

Gambierol is a ladder-shaped polyether toxin that acts as a gating modifier to inhibit Kv3.1 channels with nanomolar affinity. Binding determinants for gambierol have been identified at an interface between S5 and S6, located outside the permeation pathway. However, the high gambierol sensitivity of Kv3.1 channels could not be fully transplanted to the insensitive Kv2.1 channel by introducing the S5-S6 determinants. To explore whether also the voltage-sensing domain (VSD) is a determinant for gambierol sensitivity, we exchanged the complete VSD (S1-S4), parts of the VSD (the S1-S3a region and the Sb-S4 paddle), and the electromechanical coupling (L45+S6c) between Kv3.1 and Kv2.1. Our results show that the L45+S6c and the S1-S3a region did not alter the affinity of Kv3.1 channels for gambierol. In contrast, the distal part of the VSD, the S3b-S4 paddle, displayed a 100-fold decrease in affinity compared to WT Kv3.1. Since all VSD chimeras displayed similar biophysical properties and remained sensitive to well-known pore blockers, the loss in gambierol sensitivity in the S3b-S4 paddle chimera is most likely not the result of allosteric effects. Molecular-Dynamics simulations indicated that the S3b-S4 paddle motif resides in proximity of gambierol and that the structure and position of the VSD may regulate the space of the binding site between the pore domain (S5 and S6) and the gating machinery. Hence, our results suggest that the VSD, and especially the S3b-S4 paddle motif, contributes to the structure and/or the accessibility of the binding site for gambierol. (This research was supported by FWO grant G0433.12N to DJS and JT).

The inwardly-rectifying potassium channel Kir4.1 is a major player in the astrocyte-mediated regulation of [K(+)](o) in the brain, which is essential for normal neuronal activity and synaptic functioning. KCNJ10, encoding Kir4.1, has been recently linked to seizure susceptibility in humans and mice, and is a possible candidate gene for Autism Spectrum Disorders (ASD). In this study, we performed a mutational screening of KCNJ10 in 52 patients with epilepsy of "unknown cause" associated with impairment of either cognitive or communicative abilities, or both. Among them, 14 patients fitted the diagnostic criteria for ASD. We identified two heterozygous KCNJ10 mutations (p.R18Q and p.V84M) in three children (two unrelated families) with seizures, ASD, and intellectual disability. The mutations replaced amino acid residues that are highly conserved throughout evolution and were undetected in about 500 healthy chromosomes. The effects of mutations on channel activity were functionally assayed using a heterologous expression system. These studies indicated that the molecular mechanism contributing to the disorder relates to an increase in either surface-expression or conductance of the Kir4.1 channel. Unlike previous syndromic associations of genetic variants in KCNJ10, the pure neuropsychiatric phenotype in our patients suggests that the new mutations affect K(+) homeostasis mainly in the brain, by acting through gain-of-function defects. Dysfunction in astrocytic-dependent K(+) buffering may contribute to autism/epilepsy phenotype, by altering neuronal excitability and synaptic function, and may represent a new target for novel therapeutic approaches.

—Subdivision Surfaces provide a compact way to describe a smooth surface using a mesh model. They are widely used in 3D animation and nearly all modern modeling programs support them. In this work we describe a complete parallel pipeline for real-time interactive editing, processing and rendering of smooth surface primitives on the Cell BE. Our approach makes it possible to edit and render these high-order graphics primitives passing them directly to a parallel pipeline which tessellates them just before rendering. We describe a combination of algorithmic, architectural and back-end optimizations that enable us to render smooth subdivision surfaces in real time and to dynamically deform 3D models represented by subdivision surfaces.

Voltage-gated K+ channels (Kv) are composed of four subunits, each of which contains six trans-membrane domains (TMs), S1 through S6. The S1-S4 segments comprise the voltage-sensing domain (VS), which senses membrane potential and controls the gating of the pore domain (PD). Although still controversial, the voltage-sensing domain undergoes conformational changes within the membrane electric field, upon membrane depolarization, that is mechanically transferred via the S4-S5 linker to the intracellular gate of the channel. MD simulations of portions of PD and VS regions highlighted the importance of their flexibility for proper channel function. Nevertheless, a comprehensive description of the dynamics of both domains at atomic level has not been provided yet. Here we report the analysis of all-atom multiple molecular dynamics simulations (∼200 ns) of the entire Kv1.2/2.1 chimera, consisting of the α and β chain embedded in a 549 monomer POPC bilayer, and immersed in a box of 135K explicit SPC water molecules at 300 K. We used principal components analysis (PCA) of the Cα atomic fluctuations covariance matrix to analyze the essential subspace that characterizes the channel internal dynamics. Briefly, we observed an up to 4.5 A conformational drift of VS from its starting position. The average RMSF of the S3b-S4 domain was between 1.6 and 3.0 A. Relative to the pore region, i) the second principal component shows that all four VS domains fluctuate in a concerted manner and affect the flexibility of the intracellular gates; ii) the first principal component reveals that T1 domain moves approximately 3.5 A downwards, influencing the local structure and dynamics of the neighboring intracellular gate. Protein-lipids interactions are crucial for channel structure/function. Thus, the contributions of H-bonds and salt-bridges between channel atoms and lipid head-groups on global channel dynamics will be illustrated.

The energetic profile of an ion translated along the axis of an ion channel should reveal whether the structure corresponds to a functionally open or closed state of the channel. In this study, we explore the combined use of Poisson-Boltzmann electrostatic calculations and evaluation of van der Waals interactions between ion and pore to provide an initial appraisal of the gating state of a channel. This approach is exemplified by its application to the bacterial inward rectifier potassium channel KirBac3.1, where it reveals the closed gate to be formed by a ring of leucine (L124) side chains. We have extended this analysis to a comparative survey of gating profiles, including model hydrophobic nanopores, the nicotinic acetylcholine receptor, and a number of potassium channel structures and models. This enables us to identify three gating regimes, and to show the limitation of this computationally inexpensive method. For a (closed) gate radius of 0.4 nm < R < 0.8 nm, a hydrophobic gate may be present. For a gate radius of 0.2 nm < R < 0.4 nm, both electrostatic and van der Waals interactions will contribute to the barrier height. Below R = 0.2 nm, repulsive van der Waals interactions are likely to dominate, resulting in a sterically occluded gate. In general, the method is more useful when the channel is wider; for narrower channels, the flexibility of the protein may allow otherwise-unsurmountable energetic barriers to be overcome.

Ion channels are proteins that form pores of nanoscopic dimensions in cell membranes. As a consequence of advance in protein crystallography we now know the three-dimensional structures of a number of ion channels. However, X-ray diffraction techniques yield an essentially static (time- and space-averaged) structure of an ion channel, in an environment often somewhat distantly related to that which the protein experiences when in a cell membrane. Thus, additional techniques are required to fully understand the relationship between channel structure and function. Potassium (K) channels (Yellen, 2002) provide an opportunity to explore the relationship between membrane protein structure, dynamics, and function. Furthermore, K channels are of considerable physiological and biomedical interest. They regulate K + ion flux across cell membranes. K channel regulation is accomplished by a conformational change that allows the protein to switch between two alternative (closed vs. open) conformations, a process known as gating. Gating is an inherently dynamic process that cannot be fully characterized by static structures alone. The elucidation of the structures of several K + channels (Mackinnon, 2003

Molecular dynamics simulations may be used to probe the interactions of membrane proteins with lipids and with detergents at atomic resolution. Examples of such simulations for ion channels and for bacterial outer membrane proteins are described. Comparison of simulations of KcsA (an alpha-helical bundle) and OmpA (a beta-barrel) reveals the importance of two classes of side chains in stabilizing interactions with the head groups of lipid molecules: (i) tryptophan and tyrosine; and (ii) arginine and lysine. Arginine residues interacting with lipid phosphate groups play an important role in stabilizing the voltage-sensor domain of the KvAP channel within a bilayer. Simulations of the bacterial potassium channel KcsA reveal specific interactions of phosphatidylglycerol with an acidic lipid-binding site at the interface between adjacent protein monomers. A combination of molecular modelling and simulation reveals a potential phosphatidylinositol 4,5-bisphosphate-binding site on the surface of Kir6.2.

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

Further experimental and computational studies are required before we reach a complete structural understanding of the mechanisms of voltage-sensing and voltage-gating of Kv channels. Progress is likely to be made using a range of indirect techniques to complement the structural data that have come from X-ray crystallography. In particular, spectroscopic [7xMolecular architecture of the KvAP voltage-dependent K+ channel in a lipid bilayer. Cuello, L.G., Cortes, D.M., and Perozo, E. Science. 2004; 306: 491–495Crossref | PubMed | Scopus (173)See all References[7] and chemical modification [8xMolecular mechanism of voltage sensor movements in a potassium channel. Elliott, D.J.S., Neale, E.J., Aziz, Q., Dunham, J.P., Munsey, T.S., Hunter, M., and Sivaprasadarao, A. EMBO J. 2005; in pressSee all References[8] studies of Kv channels in situ in lipid bilayers will help to resolve the structure of the resting (closed state) channel. These methods will then have to be deployed in combination with a transmembrane voltage difference to reveal the change in sensor conformation and orientation during the transition from the closed to the open state.It is likely that computational methods will be used to integrate data from these diverse sources in a molecular mechanism. Given the proposed changes in conformation of the voltage sensor upon activation of Kv channels, it is important to characterize the intrinsic flexibility of this domain. Molecular dynamics simulations offer one possibility for exploring the conformational dynamics of the sensor domain (our unpublished data).A number of other techniques may also yield information on the location of the voltage sensor relative to the membrane. For example, toxins that interact with the sensor, for example Vstx1 [9xA membrane-access mechanism of ion channel inhibition by voltage sensor toxins from spider venom. Lee, S.Y. and MacKinnon, R. Nature. 2004; 430: 232–235Crossref | PubMed | Scopus (169)See all References[9], provide valuable probes. Their location – and hence by extension the location of the voltage-sensor – relative to a lipid bilayer may be established via molecular simulations. Similarly, electron microscopy [10xElectron microscopic analysis of KvAP voltage-dependent K+ channels in an open conformation. Jiang, Q.X., Wang, D.N., and MacKinnon, R. Nature. 2004; 430: 806–810Crossref | PubMed | Scopus (82)See all References[10] may reveal the overall shape of Kv molecules trapped in different conformational states. Thus, by combining information from these disparate sources, a complete mechanism of Kv voltage-sensing and gating may emerge.

In a previous article (Zbilut et al., Biophys J 2003;85:3544-3557), we demonstrated how an aggregation versus folding choice could be approached considering hydrophobicity distribution and charge. In this work, our aim is highlighting the mutual interaction of charge and hydrophobicity distribution in the aggregation process. Use was made of two different peptides, both derived from a transmembrane protein (amyloid precursor protein; APP), namely, Abeta(1-28) and Abeta(1-40). Abeta(1-28) has a much lower aggregation propensity than Abeta(1-40). The results obtained by means of molecular dynamics simulations show that, when submitted to the most "aggregation-prone" environment, corresponding to the isoelectric point and consequently to zero net charge, both peptides acquire their maximum flexibility, but Abeta(1-40) has a definitely higher conformational mobility than Abeta(1-28). The absence of a hydrophobic "tail," which is the most mobile part of the molecule in Abeta(1-40), is the element lacking in Abeta(1-28) for obtaining a "fully aggregating" phenotype. Our results suggest that conformational flexibility, determined by both hydrophobicity and charge effect, is the main mechanistic determinant of aggregation propensity.

Although the use of surgical staples is a well established practice in intestinal tract surgery on adults, their role in biliodigestive anastomoses in adults and children has been more limited. The Authors describe a 12-year-old girl with a type-IV choledochal cyst, who was successfully treated with cyst excision and Roux-en-Y hepaticojejunostomy created with a surgical stapler.

Background. Bladder autoaugmentation uses partial detrusorectomy to create a diverticular bulge in the bladder mucosa. This technique has eliminated certain serious complications of cystoplasty with gastrointestinal tissues (e.g., fluid/electrolyte/acid-base imbalances, mucous hypersecretion), but the exposed mucosa is subject to fibrosis and, sometimes, to perforation, which can annul the benefits of surgery. Methods. We have developed an original technique based on traditional autoaugmentation with protection of the herniated mucosa by split-thickness pedunculated rectus abdominis muscle flaps that are sutured to the incised margins of the detrusor. Preliminary testing was done on 30 adult Wistar rats. A control group of 15 rats underwent laparotomy alone. Bladder capacity was measured via suprapubic cystography before and after (4 weeks, 8 weeks, 1 year) surgery, just before sacrifice. Sections of the reconstructed bladder were examined histologically. Results. Twenty-three bladder-augmented rats and 13 controls survived. In the experimental group, bladder capacity increased by 38% (mean). None of the rats experienced urinary retention, although one developed bladder stones. Histology revealed no pathologic changes (other than chronic inflammatory infiltrates at suture sites) in the mucosa, detrusor, or muscle flaps, which were all viable and well integrated by the fourth postoperative week. There were no signs of mucosal or muscle fibrosis. Conclusions. Preliminary results in a rat model suggest that this new technique can produce an enlarged bladder that is fully functional and less vulnerable to fibrotic retraction and rupture. Residual contractility in the muscle flaps might theoretically be exploited to facilitate paraphysiologic micturition.

The effects of core-packing on the structure, function and mechanics of the RNA-binding 4-helix-bundle Rop have been studied by molecular dynamics simulations. The structural, dynamical and geometrical properties of the Rop homodimer, (formed by the antiparallel juxtaposition of two helix-turn-helix motifs), have been compared with those of three protein variants described by Munson et al. (Protein Sci, 5:1584-1593, 1996), where the core of the native protein has been systematically repacked using a two-amino acid alphabet: Ala(2)Leu(2)-8, Ala(2)Leu(2)-8-rev, and Leu(2)Ala(2)-8. The results showed that it was possible to readily distinguish the inactive protein Leu(2)Ala(2)-8 from the other functionally active systems based on tertiary and quaternary structure criteria. Structural properties such as native secondary structure content did not correlate with biological activity. Biological activity was related in part to the relative arrangement of the residues within the binding site. But, more global aspects, related to the overall topology of the helical bundle, accounted for the small functional differences between Ala(2)Leu(2)-8 and Ala(2)Leu(2)-8-rev. Mechanically, the 4-helix-bundle absorbed core mutations by altering the local structure at the sequence termini and in the turns that join the two helices of each monomer, and by changing the overall orientation and separation of the extremely rigid helices. Proteins 1999;36:436-446.

The 28-residue peptide thymosin alpha1 was studied by circular dichroism and two-dimensional NMR. Circular dichroism indicates that thymosin alpha1 in water solution does not assume a preferred conformation, while in the presence of small unilamellar vesicles of dimiristoylphosphatidylcholine and dimiristoylphosphatidic acid (10:1) and in sodium dodecyl sulphate, it assumes a partly structured conformation. Presence of zinc ions produces similar effects. In a more hydrophobic environment like a solution of a mixed solvent water-2,2,2 trifluoroethanol, it adopts a structured conformation. NMR spectra indicated that in this mixture as solvent, thymosin alpha1 has a structure characterized by two regions. A beta-turn is present between residue 5 and residue 8, while the region between residues 17 and 24 shows an alpha helix conformation. These changes of conformation in different environments may be considered structural requirements in the steps of its interaction with the lymphocyte membrane. In fact, these conformational changes may correspond to the first event of the mechanism of lymphocyte activation in the immune response modulation by thymosin alpha1.