An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints

In this work, we introduce an algorithm to compute the derivatives of physical observables along the constrained subspace when flexible constraints are imposed on the system (i.e., constraints in which the hard coordinates are fixed to configuration-dependent values). The presented scheme is exact, it does not contain any tunable parameter, and it only requires the calculation and inversion of a sub-block of the Hessian matrix of second derivatives of the function through which the constraints are defined. We also present a practical application to the case in which the sought observables are the Euclidean coordinates of complex molecular systems, and the function whose minimization defines the constraints is the potential energy. Finally, and in order to validate the method, which, as far as we are aware, is the first of its kind in the literature, we compare it to the natural and straightforward finite-differences approach in three molecules of biological relevance: methanol, N-methyl-acetamide and a tri-glycine peptideComment: 13 pages, 8 figures, published versio

Molecular dynamics simulations reveal that AEDANS is an inert fluorescent probe for the study of membrane proteins

Computer simulations were carried out of a number of AEDANS-labeled single cysteine mutants of a small reference membrane protein, M13 major coat protein, covering 60% of its primary sequence. M13 major coat protein is a single membrane-spanning, α-helical membrane protein with a relatively large water-exposed region in the N-terminus. In 10-ns molecular dynamics simulations, we analyze the behavior of the AEDANS label and the native tryptophan, which were used as acceptor and donor in previous FRET experiments. The results indicate that AEDANS is a relatively inert environmental probe that can move unhindered through the lipid membrane when attached to a membrane protein

A two-domain elevator mechanism for sodium/proton antiport

Sodium/proton (Na+/H+) antiporters, located at the plasma membrane in every cell, are vital for cell homeostasis1. In humans, their dysfunction has been linked to diseases, such as hypertension, heart failure and epilepsy, and they are well-established drug targets2. The best understood model system for Na+/H+ antiport is NhaA from Escherichia coli1, 3, for which both electron microscopy and crystal structures are available4, 5, 6. NhaA is made up of two distinct domains: a core domain and a dimerization domain. In the NhaA crystal structure a cavity is located between the two domains, providing access to the ion-binding site from the inward-facing surface of the protein1, 4. Like many Na+/H+ antiporters, the activity of NhaA is regulated by pH, only becoming active above pH 6.5, at which point a conformational change is thought to occur7. The only reported NhaA crystal structure so far is of the low pH inactivated form4. Here we describe the active-state structure of a Na+/H+ antiporter, NapA from Thermus thermophilus, at 3 Å resolution, solved from crystals grown at pH 7.8. In the NapA structure, the core and dimerization domains are in different positions to those seen in NhaA, and a negatively charged cavity has now opened to the outside. The extracellular cavity allows access to a strictly conserved aspartate residue thought to coordinate ion binding1, 8, 9 directly, a role supported here by molecular dynamics simulations. To alternate access to this ion-binding site, however, requires a surprisingly large rotation of the core domain, some 20° against the dimerization interface. We conclude that despite their fast transport rates of up to 1,500 ions per second3, Na+/H+ antiporters operate by a two-domain rocking bundle model, revealing themes relevant to secondary-active transporters in general

Is TEA an inhibitor for human Aquaporin-1?

Excessive water uptake through aquaporins can be life threatening, and disregulation of water permeability causes many diseases. Therefore, reversible aquaporin inhibitors are highly desired. In this paper, we identified the binding site for tetraethylammonium (TEA) of the membrane water channel aquaporin-1 by a combined molecular docking and molecular dynamics simulation approach. The binding site identified from docking studies was independently confirmed with an unbiased molecular dynamics simulation of an aquaporin tetramer embedded in a lipid membrane, surrounded by a 100-mM tetraethylammonium solution in water. A third independent assessment of the binding site was obtained by umbrella sampling simulations. These simulations, in addition, revealed a binding affinity of more than 17 kJ/mol, corresponding to an IC50 value of << 3 mM. Finally, we observed in our simulations a 50% reduction of the water flux in the presence of TEA, in agreement with water permeability measurements on aquaporin expressed in oocytes. These results confirm TEA as a putative lead for an aquaporin-1 inhibitor

The Dynamical Mechanism of Auto-Inhibition of AMP-Activated Protein Kinase

We use a novel normal mode analysis of an elastic network model drawn from configurations generated during microsecond all-atom molecular dynamics simulations to analyze the mechanism of auto-inhibition of AMP-activated protein kinase (AMPK). A recent X-ray and mutagenesis experiment (Chen, et al Nature 2009, 459, 1146) of the AMPK homolog S. Pombe sucrose non-fermenting 1 (SNF1) has proposed a new conformational switch model involving the movement of the kinase domain (KD) between an inactive unphosphorylated open state and an active or semi-active phosphorylated closed state, mediated by the autoinhibitory domain (AID), and a similar mutagenesis study showed that rat AMPK has the same auto-inhibition mechanism. However, there is no direct dynamical evidence to support this model and it is not clear whether other functionally important local structural components are equally inhibited. By using the same SNF1 KD-AID fragment as that used in experiment, we show that AID inhibits the catalytic function by restraining the KD into an unproductive open conformation, thereby limiting local structural rearrangements, while mutations that disrupt the interactions between the KD and AID allow for both the local structural rearrangement and global interlobe conformational transition. Our calculations further show that the AID also greatly impacts the structuring and mobility of the activation loop

A secondary RET mutation in the activation loop conferring resistance to vandetanib

Resistance to vandetanib, a type I RET kinase inhibitor, developed in a patient with metastatic lung adenocarcinoma harboring a CCDC6-RET fusion that initially exhibited a response to treatment. The resistant tumor acquired a secondary mutation resulting in a serine-to-phenylalanine substitution at codon 904 in the activation loop of the RET kinase domain. The S904F mutation confers resistance to vandetanib by increasing the ATP affinity and autophosphorylation activity of RET kinase. A reduced interaction with the drug is also observed in vitro for the S904F mutant by thermal shift assay. A crystal structure of the S904F mutant reveals a small hydrophobic core around F904 likely to enhance basal kinase activity by stabilizing an active conformer. Our findings indicate that missense mutations in the activation loop of the kinase domain are able to increase kinase activity and confer drug resistance through allosteric effects

Walks4work: Rationale and study design to investigate walking at lunchtime in the workplace setting

Background: Following recruitment of a private sector company, an 8week lunchtime walking intervention was implemented to examine the effect of the intervention on modifiable cardiovascular disease risk factors, and further to see if walking environment had any further effect on the cardiovascular disease risk factors. Methods. For phase 1 of the study participants were divided into three groups, two lunchtime walking intervention groups to walk around either an urban or natural environment twice a week during their lunch break over an 8week period. The third group was a waiting-list control who would be invited to join the walking groups after phase 1. In phase 2 all participants were encouraged to walk during their lunch break on self-selecting routes. Health checks were completed at baseline, end of phase 1 and end of phase 2 in order to measure the impact of the intervention on cardiovascular disease risk. The primary outcome variables of heart rate and heart rate variability were measured to assess autonomic function associated with cardiovascular disease. Secondary outcome variables (Body mass index, blood pressure, fitness, autonomic response to a stressor) related to cardiovascular disease were also measured. The efficacy of the intervention in increasing physical activity was objectively monitored throughout the 8-weeks using an accelerometer device. Discussion. The results of this study will help in developing interventions with low researcher input with high participant output that may be implemented in the workplace. If effective, this study will highlight the contribution that natural environments can make in the reduction of modifiable cardiovascular disease risk factors within the workplace. © 2012 Brown et al.; licensee BioMed Central Ltd

Unrelenting plasmons

Following a brief historic introduction to plasmons, their useful properties and early applications, we highlight some of the key advances in the field over the past decade. We then discuss new directions for the future, such as the use of 2D materials and strong coupling phenomena, which are likely to shape the field over the next ten years. For centuries, metals were employed in optical applications only as mirrors and gratings. New vistas opened up in the late 1970s and early 1980s with the discovery of surface-enhanced Raman scattering and the use of surface plasmon (SP) resonances for sensing. However, it was not until the 1990s, with the appearance of accurate and reliable nanofabrication techniques, that plasmonics blossomed1. Initially, the attention focused on the exploitation of SPs (collective electronic oscillations at the surface of metals) for sensing, subwavelength waveguiding and extraordinary optical transmission2. Since then, the scientific and technological interest in SPs has expanded. Correspondingly, as illustrated in Fig. 1, the number of publications in the field has increased in a steady exponential fashion for more than two decades, and the momentum driving plasmonics research looks set to continue (...

Folding of the apolipoprotein A1 driven by the salt concentration as a possible mechanism to improve cholesterol trapping

The folding of the cholesterol trapping apolipoprotein A1 in aqueous solution at increasing ionic strength is studied using atomically detailed molecular dynamics simulations. We calculate various structural properties to characterize the conformation of the protein, such as the radius of gyration, the radial distribution function and the end to end distance. Additionally we report information using tools specifically tailored for the characterization of proteins, such as the mean smallest distance matrix and the Ramachandran plot. We find that two qualitatively different configurations of this protein are preferred, one where the protein is extended, and one where it forms loops or closed structures. It is argued that the latter promote the association of the protein with cholesterol and other fatty acids.Comment: 14 pages, 6 figures. To appear in "Selected Topics of Computational and Experimental Fluid Mechanics", Springer, J. Klapp, G. Ru\'iz, A. Medina, A. L\'opez & L. Di G. Sigalotti (eds.), 201

Dimer Formation Enhances Structural Differences between Amyloid β-Protein (1–40) and (1–42): An Explicit-Solvent Molecular Dynamics Study

Amyloid -protein (A) is central to the pathology of Alzheimer's disease. A 5% difference in the primary structure of the two predominant alloforms, A and A, results in distinct assembly pathways and toxicity properties. Discrete molecular dynamics (DMD) studies of A and A assembly resulted in alloform-specific oligomer size distributions consistent with experimental findings. Here, a large ensemble of DMD–derived A and A monomers and dimers was subjected to fully atomistic molecular dynamics (MD) simulations using the OPLS-AA force field combined with two water models, SPCE and TIP3P. The resulting all-atom conformations were slightly larger, less compact, had similar turn and lower -strand propensities than those predicted by DMD. Fully atomistic A and A monomers populated qualitatively similar free energy landscapes. In contrast, the free energy landscape of A dimers indicated a larger conformational variability in comparison to that of A dimers. A dimers were characterized by an increased flexibility in the N-terminal region D1-R5 and a larger solvent exposure of charged amino acids relative to A dimers. Of the three positively charged amino acids, R5 was the most and K16 the least involved in salt bridge formation. This result was independent of the water model, alloform, and assembly state. Overall, salt bridge propensities increased upon dimer formation. An exception was the salt bridge propensity of K28, which decreased upon formation of A dimers and was significantly lower than in A dimers. The potential relevance of the three positively charged amino acids in mediating the A oligomer toxicity is discussed in the light of available experimental data