Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Does growth rate determine the rate of metabolism in shorebird chicks living in the Arctic?

We measured resting and peak metabolic rates (RMR and PMR, respectively) during development of chicks of seven species of shorebirds: least sandpiper (Calidris minutilla; adult mass 20-22 g), dunlin (Calidris alpina; 56-62 g), lesser yellowlegs (Tringa flavipes; 88-92 g), short-billed dowitcher (Limnodromus griseus; 85-112 g), lesser golden plover (Pluvialis dominicana; 150-156 g), Hudsonian godwit (Limosa haemastica; 205-274 g), and whimbrel (Numenius phaeopus; 380 g). We tested two opposing hypotheses: the growth rate-maturity hypothesis, which posits that growth rate in chicks is inversely related to functional maturity of tissues, and the fast growth rate-high metabolism hypothesis, which suggests that rapid growth is possible only with a concomitant increase in either RMR or PMR. We have found no evidence that chicks of shorebirds with fast growth rates have lower RMRs or lower PMRs, as would be predicted by the growth rate-maturity hypothesis, but our data suggested that faster-growing chest muscles resulted in increased thermogenic capacity, consistent with the fast growth-high metabolism hypothesis. The development of homeothermy in smaller species is a consequence primarily of greater metabolic intensities of heat-generating tissues. The maximum temperature gradient between a chick's body and environment that can be maintained in the absence of a net radiative load increased rapidly with body mass during development and was highest in least sandpipers and lowest among godwits. Chicks of smaller species could maintain a greater temperature gradient at a particular body mass because of their higher mass-specific maximum metabolic rates.     

Intrinsic Protein Fluorescence Interferes with Detection of Tear Glycoproteins in SDS-Polyacrylamide Gels Using Extrinsic Fluorescent Dyes

Intrinsic protein fluorescence may interfere with the visualization of proteins after SDS-polyacrylamide electrophoresis. In an attempt to analyze tear glycoproteins in gels, we ran tear samples and stained the proteins with a glycoprotein-specific fluorescent dye. The fluorescence detected was not limited to glycoproteins. There was strong intrinsic fluorescence of proteins normally found in tears after soaking the gels in 40% methanol plus 1-10% acetic acid and, to a lesser extent, in methanol or acetic acid alone. Nanograms of proteins gave visible native fluorescence and interfere with extrinsic fluorescent dye detection. Poly-L-lysine, which does not contain intrinsically fluorescent amino acids, did not fluoresce.     

Statins Do Not Reduce Atrial Fibrillation After Cardiac Valvular Surgery: A Single Centre Observational Study

Introduction

Statins may theoretically reduce postoperative atrial fibrillation (AF) in patients after cardiac valvular surgery due to preservation of endothelial function and anti-ischaemic, anti-inflammatory and anti-remodelling effects.

Methods

Two hundred seventy-two patients who underwent cardiac workup and subsequently cardiac valvular surgery without AF and concomitant coronary artery bypass grafting (CABG) at our hospital were selected. Preoperative drug use and postoperative AF were recorded. AF was defined as any episode of AF longer than 10 s. In addition, results from echocardiography and blood samples were retrieved.

Results

Baseline characteristics were as follows: mean age was 65 ± 11 years, 142 (52%) patients were male, 189 (70%) had undergone aortic valve surgery and the mean left ventricular ejection fraction was 57 ± 12%. Statins were used by 79 patients (29%). Statin users, more often, had a prior percutaneous coronary intervention (25% vs 9%, p < 0.001) or CABG (24% vs 4%, p < 0.001), diabetes mellitus (22% vs 5%, p < 0.001) and more often used β-blockers (51% vs 24%, p < 0.001). Patients in the non-statin group more often had surgery on more than one valve (10% vs 3%, p = 0.043) and had a higher cholesterol level (222 ± 48 vs 190 ± 43 mg/dl, p < 0.001). Postoperative AF occurred in 54% (43/79) of the patients with and in 55% (106/193) of the patients without statins (p = 0.941). There was also no difference in the timing of onset of AF or duration of hospital stay.

Conclusion

In this observational study, statin use was not associated with a reduced incidence of AF in patients after cardiac valvular surgery.     

Constitutive innate immunity is a component of the pace-of-life syndrome in tropical birds

We studied the relationship between one component of immune function and basal metabolic rate (BMR), an indicator of the 'pace-of-life syndrome', among 12 tropical bird species and among individuals of the tropical house wren (Troglodytes aedon), to gain insights into functional connections between life history and physiology. To assess constitutive innate immunity we introduced a new technique in the field of ecological and evolutionary immunology that quantifies the bactericidal activity of whole blood. This in vitro assay utilises a single blood sample to provide a functional, integrated measure of constitutive innate immunity. We found that the bactericidal activity of whole blood varied considerably among species and among individuals within a species. This variation was not correlated with body mass or whole-organism BMR. However, among species, bacteria killing activity was negatively correlated with mass-adjusted BMR, suggesting that species with a slower pace-of-life have evolved a more robust constitutive innate immune capability. Among individuals of a single species, the house wren, bacteria killing activity was positively correlated with mass-adjusted BMR, pointing to physiological differences in individual quality on which natural selection potentially could act.     

2H-NMR study and molecular dynamics simulation of the location, alignment, and mobility of pyrene in POPC bilayers

The alignment of pyrene in a 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer was investigated using two different approaches, namely solid-state (2)H-NMR spectroscopy and molecular dynamics (MD) simulations. Quadrupolar splittings from (2)H-NMR spectra of deuterated pyrene-d(10) in an oriented lipid bilayer give information about the orientation of C-D bonds with respect to the membrane normal. From MD simulations, geometric information is accessible via trajectories. By defining molecular and bond order parameters, the data from MD trajectories and NMR spectra can be compared straightforwardly. To ensure that the results from both methods are comparable, parameters of the experimental and the simulation setup were chosen to be as similar as possible. From simulations, we saw that pyrene prefers a position inside the lipid membrane near the headgroups and has no tendency to diffuse from one monolayer of the membrane to the other. The results from simulation and NMR show that the normal of the molecular plane is aligned nearly perpendicular to the bilayer normal. The long axis of pyrene lies preferentially parallel to the bilayer normal within a range of +/-30 degrees . The results from the two different methods are remarkably consistent. The good agreement can be explained by the fact that the different kind of motions of a pyrene molecule are already averaged within a few nanoseconds, which is the timescale covered by the MD simulation.     

Structural investigation of syringomycin-E using molecular dynamics simulation and NMR

Syringomycin-E (SR-E) is a cyclic lipodepsinonapeptide produced by certain strains of the bacterium Pseudomonas syringae pv. syringae. It shows inhibitory effects against many fungal species, including human pathogens. Its primary biological target is the plasma membrane, where it forms channels comprised of at least six SR-E molecules. The high-resolution structure of SR-E and the structure of the channels are currently not known. In this paper, we investigate in atomic detail the molecular features of SR-E in water by NMR and in water and octane by molecular dynamics simulation (MD). We built a model of the peptide and examined its structure in water and octane in 200 ns MD simulations both with and without distance restraints derived from NMR NOE data. The resulting trajectories show good agreement with the measured NOEs and circular dichroism data from the literature and provide atomistic models of SR-E that are an important step toward a better understanding of the antifungal and antibacterial activity of this peptide.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Molecular dynamics simulations of the Apo-, Holo-, and acyl-forms of Escherichia coli acyl carrier protein

Acyl carrier protein (ACP) is an essential co-factor protein in fatty acid biosynthesis that shuttles covalently bound fatty acyl intermediates in its hydrophobic pocket to various enzyme partners. To characterize acyl chain-ACP interactions and their influence on enzyme interactions, we performed 19 molecular dynamics (MD) simulations of Escherichia coli apo-, holo-, and acyl-ACPs. The simulations were started with the acyl chain in either a solvent-exposed or a buried conformation. All four short-chain (< or = C10) and one long-chain (C16) unbiased acyl-ACP MD simulation show the transition of the solvent-exposed acyl chain into the hydrophobic pocket of ACP, revealing its pathway of acyl chain binding. Although the acyl chain resides inside the pocket, Thr-39 and Glu-60 at the entrance stabilize the phosphopantetheine linker through hydrogen bonding. Comparisons of the different ACP forms indicate that the loop region between helices II and III and the prosthetic linker may aid in substrate recognition by enzymes of fatty acid synthase systems. The MD simulations consistently show that the hydrophobic binding pocket of ACP is best suited to accommodate an octanoyl group and is capable of adjusting in size to accommodate chain lengths as long as decanoic acid. The simulations also reveal a second, novel binding mode of the acyl chains inside the hydrophobic binding pocket directed toward helix I. This study provides a detailed dynamic picture of acyl-ACPs that is in excellent agreement with available experimental data and, thereby, provides a new understanding of enzyme-ACP interactions.     

Seasonal redistribution of immune function in a migrant shorebird: annual-cycle effects override adjustments to thermal regime

Throughout the annual cycle, demands on competing physiological systems change, and animals must allocate resources to maximize fitness. Immune function is one such system and is important for survival. Yet detailed empirical data tracking immune function over the entire annual cycle are lacking for most wild animals. We measured constitutive immune indices once a month for a year on captive red knots (Calidris canutus). We also examined temperature as an environmental contributor to immune variation by manipulating ambient temperature to vary energy expenditure. To identify relationships among immune indices, we performed principal-component analysis. We found significant repeatability in immune indices over the annual cycle and covariation of immune indices within and among individuals. This covariation suggests immune strategies as individual traits among individuals and the use of different immune strategies during different annual-cycle stages within individuals. Over the annual cycle, both higher-cost phagocyte-based immunity and lower-cost lymphocyte-based immunity were high during mass change, but there was a clear shift toward lower-cost lymphocyte-based immunity during peak molt. Experimental manipulation of temperature had little effect on annual variation in immune function. This suggests that other environmental factors, such as food availability and disease, should also be examined in the future.