Proper treatment of nonbonded interactions
is essential for the
accuracy of molecular dynamics (MD) simulations, especially in studies
of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in
different MD simulation programs can result in disagreements with
published simulations performed with CHARMM due to differences in
the protocols used to treat the long-range and 1-4 nonbonded interactions.
In this study, we systematically test the use of the C36 lipid FF
in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of
Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested
to find the optimal simulation protocol to best match bilayer properties
of six lipids with varying acyl chain saturation and head groups.
MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine
(DPPC) bilayer were used to obtain the optimal protocol for each program.
MD simulations with all programs were found to reasonably match the
DPPC bilayer properties (surface area per lipid, chain order parameters,
and area compressibility modulus) obtained using the standard protocol
used in CHARMM as well as from experiments. The optimal simulation
protocol was then applied to the other five lipid simulations and
resulted in excellent agreement between results from most simulation
programs as well as with experimental data. AMBER compared least favorably
with the expected membrane properties, which appears to be due to
its use of the hard-truncation in the LJ potential versus a force-based
switching function used to smooth the LJ potential as it approaches
the cutoff distance. The optimal simulation protocol for each program
has been implemented in CHARMM-GUI. This protocol is expected to be
applicable to the remainder of the additive C36 FF including the proteins,
nucleic acids, carbohydrates, and small molecules.
CHARMM-GUI Membrane Builder,
http://www.charmm-gui.org/input/membrane, is a web-based user interface designed to interactively build all-atom protein/membrane or membrane-only systems for molecular dynamics simulation through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builderthat allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types such as phosphoinositides, cardiolipin, sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance-based algorithm for faster initial ion displacement, (5) CHARMM inputs for P21 image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: cardiolipin lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM-GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments.
Sequence-specific nucleases have been applied to engineer targeted modifications in polyploid genomes, but simultaneous modification of multiple homoeoalleles has not been reported. Here we use transcription activator-like effector nuclease (TALEN) and clustered, regularly interspaced, short palindromic repeats (CRISPR)-Cas9 (refs. 4,5) technologies in hexaploid bread wheat to introduce targeted mutations in the three homoeoalleles that encode MILDEW-RESISTANCE LOCUS (MLO) proteins. Genetic redundancy has prevented evaluation of whether mutation of all three MLO alleles in bread wheat might confer resistance to powdery mildew, a trait not found in natural populations. We show that TALEN-induced mutation of all three TaMLO homoeologs in the same plant confers heritable broad-spectrum resistance to powdery mildew. We further use CRISPR-Cas9 technology to generate transgenic wheat plants that carry mutations in the TaMLO-A1 allele. We also demonstrate the feasibility of engineering targeted DNA insertion in bread wheat through nonhomologous end joining of the double-strand breaks caused by TALENs. Our findings provide a methodological framework to improve polyploid crops.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.