Name: Building the System Coordinates for Molecular Dynamics Simulations
Uploaded: 2023-09-01T14:10:00
Description: Watch this Scientific Journal Video about Building the System Coordinates for Molecular Dynamics Simulations at JoVE.com

building the system coordinates for molecular dynamics simulations

आणविक गतिशीलता का उपयोग करके लिपिड बाइलेयर का अनुकरण और विश्लेषण

Lipids are major constituents of membranes, which provide boundaries for cells and enable intracellular compartmentalization1,2,3. Lipids are amphiphilic, with a polar head group and two hydrophobic fatty acid tails; these self-assemble into a bilayer to minimize contact of the hydrophobic chains with water3,4. Various combinations of hydrophilic head groups and hydrophobic tails result in different classes of lipids in biological membranes, such as glycerophospholipids, sphingolipids, and sterols (Figure 1)1,5,6. Glycerophospholipids are primary building blocks of eukaryotic cell membranes composed of glycerophosphate, long-chain fatty acids, and head groups of low-molecular weight7. Lipid nomenclature is based on differences in head groups; examples include phosphatidyl-choline (PC), phosphatidyl-ethanolamine (PE), phosphatidyl-serine (PS), phosphatidyl-glycerol (PG), phosphatidyl-inositol (PI), or the unmodified phosphatidic acid (PA)5,6. As for hydrophobic tails, the length and degree of saturation vary, along with the backbone structure. The possible combinations are numerous, resulting in thousands of lipid species in mammalian cells6. Changes in membrane lipid composition lead to different mechanical and structural membrane properties that impact the activity of both integral membrane proteins and peripheral proteins2,6.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/65712/65712fig01.jpg" /> 
Figure 1. Representative lipid structures. Fatty acid tails are shown in blue boxes, common lipid head groups in orange, and sample backbones in purple. <a href="https://www.jove.com/files/ftp_upload/65712/65712fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Lipids are active players in cellular processes, protein activation in signaling cascades, and healthy cell homeostasis8,9. Altered lipid dynamics are the result of infection or can be markers of pathogenesis of disease10,11,12,13,14,15. As barriers for the cell, the study of membrane lipids and their role in permeation of small molecules is of relevance for drug delivery systems and membrane disruption mechanisms16,17. Chemical diversity and different ratios of lipid species across organelles, tissues, and organisms give rise to complex membrane dynamics2. It is therefore important to retain these characteristics in modelling studies of lipid bilayers, especially when the goal of a study is to examine interactions of other biomolecules with the membrane. The lipid species to consider in a model depends on the organism and cellular compartment of interest. For example, PG lipids are important for electron transfer in photosynthetic bateria18, while phosphorylated inositol lipids (PIPs) are major players in plasma membrane (PM) dynamics and signaling cascades in mammalian cells19,20. Inside the cell, the PM, endoplasmic reticulum (ER), Golgi, and mitochondrial membranes contain unique lipid abundances that influence their function. For example, the ER is the hub for lipid biogenesis and transports cholesterol out to the PM and Golgi; it contains a high lipid diversity with an abundance of PC and PE, but low sterol content, which promotes membrane fluidity21,22,23,24. In contrast, the PM incorporates hundreds and even thousands of lipid species depending on the organism25, it contains high levels of sphingolipids and cholesterol that give it a characteristic rigidity compared to other membranes in the cell24. Leaflet asymmetry should be considered for membranes like the PM, which has an outer leaflet rich in sphingomyelin, PC, and cholesterol, and an inner leaflet rich in PE, PI and PS that are important for signaling cascades24. Finally, lipid diversity also prompts the formation of micro-domains that differ in packing and internal order, known as lipid rafts24,26; these exhibit lateral asymmetry, are hypothesized to play important roles in cellular signaling26, and are hard to study due to their transient nature.
Experimental techniques such as fluoroscopy, spectroscopy, and model membrane systems like giant unilamellar vesicles (GUVs) have been used to investigate interactions of biomolecules with membranes. However, the complex and dynamic nature of the components involved is difficult to capture with experimental methods alone. For example, there are limitations on the imaging of transmembrane domains of proteins, the complexity of membranes used in such studies, and the identification of intermediate or transient states during the process of interest27,28,29. Since the advent of molecular simulation of lipid monolayers and bilayers in the 1980s29, lipid-protein systems and their interactions can now be quantified at the molecular level. Molecular dynamics (MD) simulation is a common computational technique that predicts the movement of particles based on their intermolecular forces. An additive interaction potential describes the bonded and non-bonded interactions between particles of the system30. The set of parameters used to model these interactions is termed the simulation forcefield (FF). These parameters are obtained from ab initio computations, semi-empirical, and quantum mechanical calculations, and optimized to reproduced data from X-ray and electron diffraction experiments, NMR, infrared, Raman and neutron spectroscopy, amongst other methods31.
MD simulations can be used to study systems at various levels of resolution32,33,34. Systems that aim to characterize specific biomolecular interactions, hydrogen bonds, and other high-resolution details are studied with all-atom (AA) simulations. In contrast, coarse-grained (CG) simulations lump atoms into larger functional groups to reduce computational cost and examine larger scale dynamics33. Situated in between these two are united-atom (UA) simulations, where hydrogen atoms are combined with their respective heavy atoms to accelerate the computation33,35. MD simulations are a powerful tool for exploration of the dynamics of lipid membranes and their interactions with other molecules and can serve to provide molecular level mechanisms for processes of interest at the membrane interface. Additionally, MD simulations can serve to narrow down experimental targets and predict macromolecular properties of a given system based on microscopic interactions.
In brief, given a set of initial coordinates, velocities, and a set of conditions like constant temperature and pressure, positions and velocities of each particle are calculated through numerical integration of the interaction potential and Newton&#39;s Law of motion. This is repeated iteratively, thereby generating a simulation trajectory30. These calculations are performed with an MD engine; among several open-source packages, GROMACS36 is one of the most commonly used engines and the one we describe here. It also includes tools for analysis and constructing initial coordinates of systems to be simulated37. Other MD engines include NAMD38; CHARMM39, and AMBER40, which the user can select at their own discretion based on computational performance of a given system. It is critical to visualize the trajectories during the simulation as well as for analysis and interpretation of the results. A variety of tools are available; here we discuss visual molecular dynamics (VMD) that offers a wide range of features, including three dimensional (3-D) visualization with expansive drawing and coloring methods, volumetric data visualization, building, preparing, and analyzing trajectories of MD simulation systems, and trajectory-movie making with no limits on system size, if the memory is available41,42,43.
The accuracy of predicted dynamics between system components is directly influenced by the FF chosen for the propagation of the trajectory. Empirical FF parametrization efforts are pursued by few research groups. The most established and common FF for MD include CHARMM39, AMBER40, Martini44, OPLS45, and SIRAH46. The all-atom additive CHARMM36 (C36) force field47 is widely used for AA MD of membrane systems as it accurately reproduces experimental structural data. It was originally developed by the CHARMM community, and it is compatible with multiple MD engines like GROMACS and NAMD. Despite improvements across common FFs, there is a continuous effort to improve the parameter sets to allow predictions that closely reproduce experimental observables, driven by interests in particular systems of study48,49.
A challenge when simulating lipid membranes is determining the length of simulation trajectory. This is largely dependent on the metrics to be analyzed and the process that one aims to characterize. Typically, complex lipid mixtures require longer time to reach equilibrium, as more species must have enough time to diffuse on the membrane plane and reach a stable lateral organization. A simulation is said to be at equilibrium when the property of interest has reached a plateau and fluctuates about a constant value. It is common practice to obtain at least 100-200 ns of equilibrated trajectory to carry out appropriate statistical analysis on the properties and interactions of interest. It is common to run membrane-only simulations between 200-500 ns, depending on the complexity of the lipid mixture and research question. Protein-lipid interactions typically require longer simulation times, between 500-2000 ns. Some approaches to accelerate sampling and observable dynamics with membrane systems are: (i) the highly mobile membrane mimetic (HMMM) model, which replaces end carbons of lipids in the membrane with organic solvent to accelerate sampling50; and (ii) hydrogen mass repartitioning (HMR), which combines a fraction of the masses of heavy atoms within a system with those of hydrogen atoms to allow the use of a larger simulation timestep51.
The following protocol discusses a beginner-friendly approach to build, run, and analyze realistic membrane models using AA MD. Given the nature of MD simulations, multiple trajectories must be run to account for reproducibility and proper statistical analysis of the results. It is current practice to run at minimum three replicas per system of interest. Once the lipid species have been selected for the organism and process of interest, basic steps to build, run, and analyze a simulation trajectory of a membrane-only system are outlined and summarized in Figure 2.
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/65712/65712fig02.jpg" /> 
Figure 2. Schematic to run MD simulations. Orange boxes correspond to the three main steps described in the protocol. Underneath is the workflow of the simulation process. During system setup, the system containing the initial coordinates of a solvated membrane system is built with a system input generator like CHARMM-GUI Membrane Builder. After transferring the input files to a high-performance computing cluster, the simulation trajectory is propagated using an MD engine, such as GROMACS. Trajectory analysis can be done on the computer cluster or a local working station along with visualization. Analysis is then carried out, using either packages with built-in analysis code such as GROMACS and VMD, or by using Bash scripts or various Python libraries. <a href="https://www.jove.com/files/ftp_upload/65712/65712fig02large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>

लेखक स्पॉटलाइट: सेल झिल्ली बायोफिज़िक्स को आगे बढ़ाना - प्रयोगात्मक और कम्प्यूटेशनल दृष्टिकोण के माध्यम से इंटरैक्शन और चुनौतियों की खोज 

सिमुलेशन अध्ययन में जटिल लिपिड मिश्रण का उपयोग करते हुए यथार्थवादी झिल्ली मॉडलिंग

Cell membrane biophysics is rapidly growing. Common experimental approaches include fluorescence microscopy, neutron and X-ray scattering, mass spectroscopy, and atomic force microscopy. Computational approaches include phonological models and those based on statistical thermodynamics that help us understand atomistic or molecular level details of interactions occurring at the membrane interface and inside its hydrophobic core.When using molecular dynamic simulations, challenges include sampling events of interest, setting a simulations length, ensuring results converge and reproduce physical values, and accessing computational power. These challenges are particularly true for transient events at the membrane interface, such as the interaction of peripheral proteins with it or the aggregation of lipids and proteins that require large changes in confirmation. This protocol provides a beginner friendly step-by-step guide to start running molecular dynamic simulations of complex lipid membranes.There are many software alternatives to carry these simulations, and some packages have tutorials or manuals. We hope our protocol provides a concise foundation on realistic membrane modeling and tips on the considerations that influence the quality of results and these types of modeling studies. This protocol has allowed us to capture interactions between membrane lipids and other biomolecules that were not observed using pure or binary lipid mixtures.Many interactions at the membrane surface depend on the diversity of lipids on the membrane itself. Our models demonstrate the importance of incorporating appropriate lipid species to accurately explore biomolecular function within membranes.

Watch this Scientific Journal Video about Building the System Coordinates for Molecular Dynamics Simulations at JoVE.com

Building the System Coordinates for Molecular Dynamics Simulations  

आणविक गतिशीलता सिमुलेशन के लिए सिस्टम निर्देशांक का निर्माण

शुरू करने के लिए, चार्म-जीआई पर जाने के लिए एक वेब ब्राउज़र का उपयोग करें। ऑर्ग वेबसाइट। फ़ाइलों का अपना पहला सेट बनाने से पहले अपना मुफ्त खाता बनाएं और सक्रिय करें।शीर्ष मेनू पर जाएं, इनपुट जनरेटर पर नेविगेट करें और स्क्रीन के बाईं ओर उपलब्ध विकल्पों में से झिल्ली बिल्डर का चयन करें। लिपिड बाइलेयर बनाने के लिए बाईलेयर बिल्डर का चयन करें, फिर केवल झिल्ली प्रणाली का चयन करें और भविष्य की पहुंच के लिए उत्पन्न जॉब आईडी को स्टोर करें। पृष्ठ के शीर्ष पर स्थित दृश्य संरचना पर क्लिक करके भवन प्रक्रिया के हर चरण के दौरान सिस्टम को विज़ुअलाइज़ करें।पैच आकार में अनुपलब्ध घटकों या त्रुटियों के लिए लगातार निरीक्षण करें। एकल घटक बाइलेयर बनाने की परवाह किए बिना विषम लिपिड विकल्प का चयन करें, फिर एक आयताकार बॉक्स प्रकार चुनें। हाइड्रेशन के लिए, प्रति लिपिड 45 पानी के अणुओं का चयन करें, जो पूरी तरह से हाइड्रेटेड बाइलेयर के लिए पर्याप्त है।अब लिपिड घटकों की संख्या के आधार पर एक्सवाई की लंबाई निर्धारित करें। पूर्व-नियोजित मॉडल के आधार पर प्रत्येक लिपिड प्रजातियों के लिए आवश्यक लिपिड की संख्या निर्धारित करें। एंडोप्लाज्मिक रेटिकुलम और यूकेरियोटिक कोशिकाओं का प्रतिनिधित्व करने के लिए, पीआई मॉडल में 336 डीओपीसी, 132 डीपीपीई, 60 कोलेस्ट्रॉल और 72 पीओपीआई लिपिड के संयोजन का उपयोग किया गया था, जबकि 330 डीओपीसी, 126 डीपीपीई, 54 कोलेस्ट्रॉल, 66 पीओपीआई और 24 डीओपीएस लिपिड का उपयोग पीआईपीएस मॉडल के लिए किया गया था।एक सममित झिल्ली संरचना बनाने के लिए लिपिड नाम के बगल में दो बक्से में ऊपरी और निचले पत्रक के लिए अणुओं की वांछित संख्या इनपुट करें। फिर लिपिड प्रजातियों की सूची के शीर्ष पर नेविगेट करें और सिस्टम जानकारी दिखाएं बटन पर क्लिक करके कार्रवाई शुरू करें। चार्मम-जीयूआई का उपयोग करके, सममित बाईलेयर में प्रत्येक पत्रक में लिपिड की कुल गणना को मान्य करें।सिमुलेशन के लिए तापमान और दबाव के लिए वांछित मान और पसंद के आणविक गतिशीलता इंजन के अनुसार सिमुलेशन फ़ाइलों के लिए सिंटैक्स सेट करने के लिए अगली स्क्रीन पर क्लिक करना जारी रखें। उदाहरण के लिए, GROMACS। परिणामी फ़ाइलें डाउनलोड करें और उन्हें कंप्यूटर क्लस्टर में स्थानांतरित करें।चयनित सॉफ़्टवेयर का उपयोग करें, उदाहरण के लिए, दृश्य आणविक गतिशीलता या PyMOL अंतिम सिस्टम की कल्पना करने के लिए।

Research

JoVE Journal

Biology

Building the System Coordinates for Molecular Dynamics Simulations

Realistic Membrane Modeling Using Complex Lipid Mixtures in Simulation Studies

Lipids are structural building blocks of cell membranes; lipid species vary across cell organelles and across organisms. This variety results in different ...