252 research outputs found
Recommended from our members
On the Role of Dewetting Transitions in Host-Guest Binding Free Energy Calculations.
We use thermodynamic integration (TI) and explicit solvent molecular dynamics (MD) simulation to estimate the absolute free energy of host-guest binding. In the unbound state, water molecules visit all of the internally accessible volume of the host, which is fully hydrated on all sides. Upon binding of an apolar guest, the toroidal host cavity is fully dehydrated; thus, during the intermediate λ stages along the integration, the hydration of the host fluctuates between hydrated and dehydrated states. Estimating free energies by TI can be especially challenging when there is a considerable difference in hydration between the two states of interest. We investigate these aspects using the popular TIP3P and TIP4P water models. TI free energy estimates through MD largely depend on water-related interactions, and water dynamics significantly affect the convergence of binding free energy calculations. Our results indicate that wetting/dewetting transitions play a major role in slowing the convergence of free energy estimation. We employ two alternative approaches-one analytical and the other empirically based on actual MD sampling-to correct for the standard state free energy. This correction is sizable (up to 4 kcal/mol), and the two approaches provide corrections that differ by about 1 kcal/mol. For the system considered here, the TIP4P water model combined with an analytical correction for the standard state free energy provides higher overall accuracy. This observation might be transferable to other systems in which water-related contributions dominate the binding process
The radical cation of bacteriochlorophyll b. A liquid-phase endor and triple resonance study
The previous termradical cationnext term of bacterioehlorophyll b (BChl b) is investigated by ENDOR and TRIPLE resonance in liquid solution. The experimental hyperfine coupling constants, ten proton and three nitrogen couplings, are compared with the predictions from advanced molecular-orbital calculations (RHF INDO/SP). The detailed picture obtained of the spin density distribution is a prerequisite for the investigation of the primary electron donor previous termradical cationnext term in BChl b containing photosynthetic bacteria
PHOTOCHEMICAL RING-OPENING IN meso-CHLORINATED CHLOROPHYLLS
Irradiation of 20-chloro-chlorophylls of the a-type with visible light produces long-wavelength shifted photoproducts, which transform in the dark to linear tetrapyrroles (bile pigments). The possible significance for chlorophyll degradation is discussed
Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics
Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ∼1.5 μs simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans–gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2–3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems
Simulation vs. Reality: A Comparison of In Silico Distance Predictions with DEER and FRET Measurements
Site specific incorporation of molecular probes such as fluorescent- and nitroxide spin-labels into biomolecules, and subsequent analysis by Förster resonance energy transfer (FRET) and double electron-electron resonance (DEER) can elucidate the distance and distance-changes between the probes. However, the probes have an intrinsic conformational flexibility due to the linker by which they are conjugated to the biomolecule. This property minimizes the influence of the label side chain on the structure of the target molecule, but complicates the direct correlation of the experimental inter-label distances with the macromolecular structure or changes thereof. Simulation methods that account for the conformational flexibility and orientation of the probe(s) can be helpful in overcoming this problem. We performed distance measurements using FRET and DEER and explored different simulation techniques to predict inter-label distances using the Rpo4/7 stalk module of the M. jannaschii RNA polymerase. This is a suitable model system because it is rigid and a high-resolution X-ray structure is available. The conformations of the fluorescent labels and nitroxide spin labels on Rpo4/7 were modeled using in vacuo molecular dynamics simulations (MD) and a stochastic Monte Carlo sampling approach. For the nitroxide probes we also performed MD simulations with explicit water and carried out a rotamer library analysis. Our results show that the Monte Carlo simulations are in better agreement with experiments than the MD simulations and the rotamer library approach results in plausible distance predictions. Because the latter is the least computationally demanding of the methods we have explored, and is readily available to many researchers, it prevails as the method of choice for the interpretation of DEER distance distributions
Weak temperature dependence of P (+) H A (-) recombination in mutant Rhodobacter sphaeroides reaction centers
International audienceIn contrast with findings on the wild-type Rhodobacter sphaeroides reaction center, biexponential P (+) H A (-) → PH A charge recombination is shown to be weakly dependent on temperature between 78 and 298 K in three variants with single amino acids exchanged in the vicinity of primary electron acceptors. These mutated reaction centers have diverse overall kinetics of charge recombination, spanning an average lifetime from ~2 to ~20 ns. Despite these differences a protein relaxation model applied previously to wild-type reaction centers was successfully used to relate the observed kinetics to the temporal evolution of the free energy level of the state P (+) H A (-) relative to P (+) B A (-) . We conclude that the observed variety in the kinetics of charge recombination, together with their weak temperature dependence, is caused by a combination of factors that are each affected to a different extent by the point mutations in a particular mutant complex. These are as follows: (1) the initial free energy gap between the states P (+) B A (-) and P (+) H A (-) , (2) the intrinsic rate of P (+) B A (-) → PB A charge recombination, and (3) the rate of protein relaxation in response to the appearance of the charge separated states. In the case of a mutant which displays rapid P (+) H A (-) recombination (ELL), most of this recombination occurs in an unrelaxed protein in which P (+) B A (-) and P (+) H A (-) are almost isoenergetic. In contrast, in a mutant in which P (+) H A (-) recombination is relatively slow (GML), most of the recombination occurs in a relaxed protein in which P (+) H A (-) is much lower in energy than P (+) H A (-) . The weak temperature dependence in the ELL reaction center and a YLH mutant was modeled in two ways: (1) by assuming that the initial P (+) B A (-) and P (+) H A (-) states in an unrelaxed protein are isoenergetic, whereas the final free energy gap between these states following the protein relaxation is large (~250 meV or more), independent of temperature and (2) by assuming that the initial and final free energy gaps between P (+) B A (-) and P (+) H A (-) are moderate and temperature dependent. In the case of the GML mutant, it was concluded that the free energy gap between P (+) B A (-) and P (+) H A (-) is large at all times
Using Selectively Applied Accelerated Molecular Dynamics to Enhance Free Energy Calculations
Accelerated molecular dynamics (aMD) has been shown to enhance conformational space sampling relative to classical molecular dynamics; however, the exponential reweighting of aMD trajectories, which is necessary for the calculation of free energies relating to the classical system, is oftentimes problematic, especially for systems larger than small poly peptides. Here, we propose a method of accelerating only the degrees of freedom most pertinent to sampling, thereby reducing the total acceleration added to the system and improving the convergence of calculated ensemble averages, which we term selective aMD. Its application is highlighted in two biomolecular cases. First, the model system alanine dipeptide is simulated with classical MD, all-dihedral aMD, and selective aMD, and these results are compared to the infinite sampling limit as calculated with metadynamics. We show that both forms of aMD enhance the convergence of the underlying free energy landscape by 5-fold relative to classical MD; however, selective aMD can produce improved statistics over all-dihedral aMD due to the improved reweighting. Then we focus on the pharmaceutically relevant case of computing the free energy of the decoupling of oseltamivir in the active site of neuraminidase. Results show that selective aMD greatly reduces the cost of this alchemical free energy transformation, whereas all-dihedral aMD produces unreliable free energy estimates
Conformational Basis for Asymmetric Seeding Barrier in Filaments of Three- and Four-Repeat Tau
*S Supporting Information ABSTRACT: Tau pathology in Alzheimer’s disease is intimately linked to the deposition of proteinacious filaments, which akin to infectious prions, have been proposed to spread via seeded conversion. Here we use double electron−electron resonance (DEER) spectroscopy in combination with extensive computational analysis to show that filaments of three- (3R) and four-repeat (4R) tau are conformationally distinct. Distance measurements between spin labels in the third repeat, reveal tau amyloid filaments as ensembles of known β-strand−turn−β-strand U-turn motifs. Whereas filaments seeded with 3R tau are structurally homogeneous, filaments seeded with 4R tau are heterogeneous, composed of at least three distinct conformers. These findings establish a molecular basis for the seeding barrier between different tau isoforms and offer a new powerful approach for investigating the composition and dynamics of amyloid fibril ensembles
Assessment of the Food Habits of the Moroccan Dorcas Gazelle in M’Sabih Talaa, West Central Morocco, Using the trnL Approach
Food habits of the Moroccan dorcas gazelle, Gazella dorcas massaesyla, previously investigated in the 1980s using microhistological fecal analysis, in the M’Sabih Talaa Reserve, west central Morocco, were re-evaluated over three seasons (spring, summer and autumn 2009) using the trnL approach to determine the diet composition and its seasonal variation from fecal samples. Taxonomic identification was carried out using the identification originating from the database built from EMBL and the list of plant species within the reserve. The total taxonomic richness in the reserve was 130 instead of 171 species in the 1980s. The diet composition revealed to be much more diversified (71 plant taxa belonging to 57 genus and 29 families) than it was 22 years ago (29 identified taxa). Thirty-four taxa were newly identified in the diet while 13 reported in 1986–87 were not found. Moroccan dorcas gazelle showed a high preference to Acacia gummifera, Anagallis arvensis, Glebionis coronaria, Cladanthus arabicus, Diplotaxis tenuisiliqua, Erodium salzmannii, Limonium thouini, Lotus arenarius and Zizyphus lotus. Seasonal variations occurred in both number (40–41 taxa in spring-summer and 49 taxa in autumn vs. respectively 23–22 and 26 in 1986–1987) and taxonomic type of eaten plant taxa. This dietary diversification could be attributed either to the difference in methods of analysis, trnL approach having a higher taxonomic resolution, or a potential change in nutritional quality of plants over time
- …