Research article Topical Sections

Investigation into the fungal diversity within different regions of the gastrointestinal tract of Panaque nigrolineatus, a wood-eating fish

  • Received: 27 June 2017 Accepted: 25 August 2017 Published: 04 September 2017
  • The Amazonian catfish, Panaque nigrolineatus have several physiological adaptions enabling the scraping and consumption of wood (xylivory), facilitating a detritivorous dietary strategy. Composed of lignocellulose, wood is a difficult substrate to degrade and as yet, it is unclear whether the fish obtains any direct nutritional benefits from wood ingestion and degradation. However, there are numerous systems that rely on microbial symbioses to provide energy and other nutritional benefits for host organisms via lignocellulose decomposition. While previous studies on the microbial community of P. nigrolineatus have focused upon the bacterial population, the role of fungi in lignocellulose degradation in the fish has not yet been examined. This study describes the detection of fungi within the fish gastrointestinal tract. Using next generation sequencing, the effects of diet on enteric fungal populations were examined in each gastrointestinal tract region. Fungal species were found to vary in different regions of the gastrointestinal tract as a function of diet. This study is the first to examine the fungal community in a xylivorous fish and results support the hypothesis that diet influences fungal distribution and diversity within the gastrointestinal tract of P. nigrolineatus.

    Citation: Caroline L. Marden, Ryan McDonald, Harold J. Schreier, Joy E.M. Watts. Investigation into the fungal diversity within different regions of the gastrointestinal tract of Panaque nigrolineatus, a wood-eating fish[J]. AIMS Microbiology, 2017, 3(4): 749-761. doi: 10.3934/microbiol.2017.4.749

    Related Papers:

    [1] Sunmi Lee, Chang Yong Han, Minseok Kim, Yun Kang . Optimal control of a discrete-time plant–herbivore/pest model with bistability in fluctuating environments. Mathematical Biosciences and Engineering, 2022, 19(5): 5075-5103. doi: 10.3934/mbe.2022237
    [2] Luis F. Gordillo . Optimal sterile insect release for area-wide integrated pest management in a density regulated pest population. Mathematical Biosciences and Engineering, 2014, 11(3): 511-521. doi: 10.3934/mbe.2014.11.511
    [3] Xinli Hu, Wenjie Qin, Marco Tosato . Complexity dynamics and simulations in a discrete switching ecosystem induced by an intermittent threshold control strategy. Mathematical Biosciences and Engineering, 2020, 17(3): 2164-2178. doi: 10.3934/mbe.2020115
    [4] Liping Wu, Zhongyi Xiang . A study of integrated pest management models with instantaneous and non-instantaneous impulse effects. Mathematical Biosciences and Engineering, 2024, 21(2): 3063-3094. doi: 10.3934/mbe.2024136
    [5] Teklebirhan Abraha, Fahad Al Basir, Legesse Lemecha Obsu, Delfim F. M. Torres . Farming awareness based optimum interventions for crop pest control. Mathematical Biosciences and Engineering, 2021, 18(5): 5364-5391. doi: 10.3934/mbe.2021272
    [6] Bing Liu, Gang Hu, Baolin Kang, Xin Huang . Analysis of a hybrid pest management model incorporating pest resistance and different control strategies. Mathematical Biosciences and Engineering, 2020, 17(5): 4364-4383. doi: 10.3934/mbe.2020241
    [7] Baolin Kang, Xiang Hou, Bing Liu . Threshold control strategy for a Filippov model with group defense of pests and a constant-rate release of natural enemies. Mathematical Biosciences and Engineering, 2023, 20(7): 12076-12092. doi: 10.3934/mbe.2023537
    [8] Yanyun Li, Juhua Liang . Models for determining the optimal switching time in chemical control of pest with pesticide resistance. Mathematical Biosciences and Engineering, 2021, 18(1): 471-494. doi: 10.3934/mbe.2021026
    [9] Zhenzhen Shi, Huidong Cheng, Yu Liu, Yanhui Wang . Optimization of an integrated feedback control for a pest management predator-prey model. Mathematical Biosciences and Engineering, 2019, 16(6): 7963-7981. doi: 10.3934/mbe.2019401
    [10] Yi Yang, Lirong Liu, Changcheng Xiang, Wenjie Qin . Switching dynamics analysis of forest-pest model describing effects of external periodic disturbance. Mathematical Biosciences and Engineering, 2020, 17(4): 4328-4347. doi: 10.3934/mbe.2020239
  • The Amazonian catfish, Panaque nigrolineatus have several physiological adaptions enabling the scraping and consumption of wood (xylivory), facilitating a detritivorous dietary strategy. Composed of lignocellulose, wood is a difficult substrate to degrade and as yet, it is unclear whether the fish obtains any direct nutritional benefits from wood ingestion and degradation. However, there are numerous systems that rely on microbial symbioses to provide energy and other nutritional benefits for host organisms via lignocellulose decomposition. While previous studies on the microbial community of P. nigrolineatus have focused upon the bacterial population, the role of fungi in lignocellulose degradation in the fish has not yet been examined. This study describes the detection of fungi within the fish gastrointestinal tract. Using next generation sequencing, the effects of diet on enteric fungal populations were examined in each gastrointestinal tract region. Fungal species were found to vary in different regions of the gastrointestinal tract as a function of diet. This study is the first to examine the fungal community in a xylivorous fish and results support the hypothesis that diet influences fungal distribution and diversity within the gastrointestinal tract of P. nigrolineatus.


    Abbreviations:

    AAMD: all-atom molecular dynamics; CGMD: coarse-grained molecular dynamics; MD: molecular dynamics; HIV: human immunodeficiency virus; PLA2: phospholipase A2; LDL: low-density lipoproteins; SR: scavenger receptor; AM: amphiphilic macromolecules; QSAR: quantitative structure-activity relationship; NP: nanoparticles; CPP: cell-penetrating peptide; DDS: drug delivery system; IgG: immunoglobulin G; SpA: Staphylococcus aureus protein A; RMSD: root mean square deviation; RMSF: root mean square fluctuation; GPCR: G protein-coupled receptor; POPE: 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine; POPC: 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine; DPPE: 1,2-Dipalmitoyl-sn-glycero-3-phosphoethanolamine; DPPC: 1,2-Dipalmitoyl-sn-glycero-3-phosphocholine

    1. Introduction

    Richard P. Feynman stated, “Everything we know is only some kind of approximation, because we know that we do not know all the laws as yet.” In contrast, simulations attempt to reveal nature’s secrets by applying approximations of nature’s laws. As the name suggests, coarse-grained molecular dynamics (CGMD) simulates the behaviors of atoms and molecules by audacious approximations. Nevertheless,when used effectively, CGMD simulations can explain the physicochemical nature or even predict the behavior of a biological system,which may be impossible experimentally. Because of this advantage, the amount of research that uses CGMD simulations has steadily increased in recent years (Figure 1).

    Figure 1. Increase in the number of CGMD-related papers. Pubmed was used to search for the numbers of articles published per year that included “coarse grained” and “molecular dynamics”.

    Although simulating a molecular system based on quantum chemistry is possible [1], molecular dynamics usually employs classical Newtonian physics to simulate molecular motions. The force between atoms is derived from bonded interactions, including 2-,3-,and 4-body interactions, and non-bonded interactions such as van der Waals and electrostatic interactions. All-atom molecular dynamics (AAMD) calculate the motions of every atom, including those of hydrogen, and thereby simulate the behavior of a system with considerable accuracy.

    However, high computational costs limit the spatial and temporal scales that can be used with AAMD simulations. In contrast, CGMD substantially reduces these costs by replacing multiple atoms with a larger, unified particle (coarse-grained atom); thus, the degrees of freedom of a system are limited. This makes the energy function smoother and allows the use of larger time steps (e.g., 30 fs) as compared with that of AAMD (e.g., 2 fs). By appropriately setting up a CGMD system, the number of simulations can be exponentially increased with a limited loss of accuracy.

    In this review, we discuss the variations, potentials,and limitations of CGMD simulations. Strategic designing of a molecular system of interest and applying appropriate analyses are the keys to success with CGMD simulations. The validity of CGMD simulations can be supported by comparing the simulation results with experimentally derived evidence. Here,we describe specific examples of CGMD applications and the methods used to correlate the results obtained from simulations with those obtained from experiments.

    2. Overview of CGMD

    2.1. Configuring a CGMD system

    In essence, coarse-graining of a molecular system is a process to describe the behaviors of molecules in a simpler way by discarding less essential effects. Historically,a simple model in which amino acid residues were placed in a two-dimensional square lattice,the Gō model,was introduced to investigate protein folding [2]. The Gō model takes into account only those interactions between residues that are present in the native protein structure. This strikingly simple theory was considerably successful in describing the folding mechanisms of real proteins. It inspired many researchers and, after various modifications, is still being applied.

    There are infinite ways to construct CGMD models. Variations of these models are derived from the differences in the natures of molecular systems. For example,there is no need to treat solvent molecules,such as water, explicitly when the target molecular event involves DNA self-assembly [3]. In contrast, incorporating solvent molecules, ions, and electrostatic interactions between charged molecules is necessary when the object of interest is the translocation of cationic nanoparticles through a lipid bilayer [4]. When proteins are the molecules of interest, the numbers and types of coarse-grained atoms assigned to represent each amino acid side chain are crucial factors that define a CGMD model [5].

    In more complex CGMD models, four to six coarse-grained atoms are allocated to one amino acid residue. Some CGMD models explicitly incorporate hydrogen bonding [6]. In contrast, only one pseudo atom is assigned for an entire subunit of a protein in a huge molecular system, such as a virus capsid.

    Some CGMD models are individually tailored for each study, whereas others are generalized for various molecular systems. Parameterization of atomic interactions (forcefield) requires expertise in physics and chemistry; thus, it is particularly difficult for researchers who are not expert in these fields to construct a universal CGMD model. However, there are some generalized CGMD models, including Martini [7, 8] and NAMD-CG builder [9, 10]. Among these,the Martini forcefield has recently gained popularity [4, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]. Because this CGMD model runs under AAMD software Gromacs [25] with a modified forcefield for coarse-grained atoms, the computational algorithm for CGMD simulation is exactly the same as that used for AAMD.

    The Martini forcefield incorporates preset parameters for general phospholipids and proteins. CGMD simulations can basically be run for lipid membranes and/or proteins with known three dimensional structures. However, there are barriers during actual simulations because some knowledge of several processes is required, including the construction of a molecular system and potential energy minimization. Applications of the Martini forcefield have increased after additional coarse-grained parameters were proposed for lipids, including cardiolipin [26] and glycolipids [22], carbohydrates such as monosaccharides, disaccharides, and oligosaccharides [27], DNA [23], and artificial carbon nanotubes [19]. Among CGMD models that have to be reparameterized each time and whose targets are restricted to specific types of molecules (e.g., peptides), Martini’s ready-to-use and highly applicable system is very attractive.

    In addition, the Martini forcefield can simulate the dielectric responses of solvent molecules. Coarse-grained water atoms in the Martini forcefield were originally apolar, similar to several CGMD models. A recently proposed polarized water model enabled simulations of realistic phenomena, such as electroporation of a lipid bilayer with voltage applied across a bilayer membrane [28].

    2.2. What can be done with CGMD simulations

    During a CGMD simulation, analysis methods and perturbations similar to those used with AAMD can be used. Basic quantities used in a CGMD analysis include indicators of structural alterations such as root mean square deviation (RMSD) and root mean square fluctuation (RMSF), tilt angles of lipids and proteins, area per lipid within a lipid membrane, and diffusion constants for molecules.

    Local pressure profiles [29], binding free energy, and potential of mean force are analyzed. In addition, voltage can be locally analyzed [4]. These energy profiles are important for large molecular systems, such as for protein-protein interactions.

    Regarding perturbations, the positions of every atom can be controlled during a CGMD simulation. Position restraints and force applications are used during steered molecular dynamics. These features are used for protein and ligand docking. Replacing amino acid residues (mutation) is possible for proteins, which is useful for comparing the results of simulations with those of experiments. In addition to temperature and pressure, the surface tension of a lipid bilayer is controllable during a CGMD simulation. Applying a stretch force to membrane proteins is achieved by incorporating surface tension.

    Membrane potential, the difference in electric potential between the inside and outside of a cell, is important for certain cellular functions. This voltage difference can be mimicked during a CGMD simulation by creating two compartments by modeling two parallel planes of lipid bilayers inside a simulation box, and then adding the appropriate types and numbers of ions into each compartment.

    In addition, pH effects can be mimicked during a CGMD simulation. The method used to alter the pH in a molecular system by modifying the protonation state of a titratable site, originally suggested for AAMD simulations [30], is also possible for CGMD simulations [26, 31]. pH alters the protonation states of aspartate, glutamate, and histidine in proteins [32]; thus, it changes a protein’s activities.

    Another feature that makes CGMD simulations very powerful is “reverse graining”, to reconstruct all atom geometry from coarse-grained molecular geometry. Reverse graining is achieved by a simulated annealing method, which efficiently determines the minimum potential energy [33, 34]. Reverse graining enables a strategy to first observe any large conformational changes in a molecular system by CGMD and then analyze the detailed molecular interactions by AAMD.

    2.3. Advantages and disadvantages

    If a molecular system is appropriately constructed, then high-speed CGMD simulations can be used with nearly the same accuracy as those achieved with AAMD simulations. An example of a CGMD simulation of a human immunodeficiency virus (HIV) membrane protein is shown (Figure 2). CGMD simulation predictions of binding free energies between proteins are more than 500 times faster than those with AAMD simulations, and with similar accuracy [13]. A CGMD simulation has the advantage of being fast; however, there are some disadvantages due to coarse graining.

    Figure 2. Example of a CGMD simulation. This molecular system consisted of one HIV-1 viral protein U (Vpu, PDB id: 2JPX), 30 POPEs, 30 POPCs, 30 DPPEs, 30 DPPCs, 30 galactosylceramides, 30 cholesterols, and 3000 CG water atoms. The total number of CG atoms was 5178. Lipids and water particles were initially randomly placed within a 12 × 12 × 12 nm box, after which the self-assembly of lipid molecules was observed. The position of the Vpu protein was constrained. The size of each picture corresponds to the size of a periodic boundary box. In the top views from 0.2 to 1.0 µs, water atoms were omitted for clarity. Simulation for 1 µs was performed using a personal laptop computer (2.1 GHz single core) with a computation time of about 9 h.

    A disadvantage with the Martini CGMD model is that alterations in the secondary structures of proteins are beyond the scope of simulations. If alpha helix parameters are assigned to the amino acid residues of a protein, then these parameters do not change throughout the simulations. This can potentially affect the accuracy of CGMD simulations.

    Stark et al. reported that the estimated protein-protein interaction between lysozyme and chymotrypsinogen with a CGMD simulation was greater than the experimentally derived value. They attributed this discrepancy to the intrinsic problems of CGMD parameters and proposed that down-scaling of van der Waals parameters could resolve this problem [11]. The key to a successful CGMD simulation is to set up the molecular system of interest based on sufficient knowledge regarding these issues.

    Although vast computational power enables more accurate AAMD to simulate larger molecular systems for a longer time [35, 36], CGMD simulations are still good at handling large molecular systems, such as a membrane tether consisting of 4 million particles [37], a fully solvated protein in a liposome [38], and for effectively sampling the configuration space of target molecules because of rapid calculation speed. An example of a conformational change in a mechanosensitive ion channel in response to membrane tension is shown (Figure 3), where a CGMD simulation can reproduce the dynamic actions of a protein and aid in explaining the observations derived from experiments [21]. Analysis of these simulation results leads to a better explanation of the experimentally observed phenomena based on the theories of physics and chemistry.

    Figure 3. Opening of a bacterial mechanosensitive channel of large conductance (MscL, PDB accession number: 2OAR) in response to lipid bilayer stretch obtained by CGMD simulations. An MscL channel is embedded in a lipid bilayer consisting of POPC lipids, which are shown in stick representation. The simulation box is filled with CG water atoms. The channel’s pore enlarged and the bilayer thickness decreased after applying stretch. Left: before applying tension, right: 4 µs after applying tension.

    3.1. Simulations of large molecular systems

    Protein folding is a problem that is scientifically interesting [39, 40] as well as medically important, given that protein misfolding is a cause of several diseases, such as Parkinson’s disease [41]. Borgia et al. demonstrated that similarities in amino acid sequences between adjacent domains induced protein misfolding in multidomain proteins using CGMD simulations [42]. They used a Gō-like model in which each amino acid residue was represented as a single unified atom, and either an attractive or a repulsive interaction was defined between any two residues [43]. Water molecules and ions were not incorporated, drastically decreasing the degrees of freedom in their simulation system. This is a pr actical choice for effectively sampling protein structures among a vast conformational space in protein folding studies.

    Protein aggregation, which is also a cause of certain diseases, including Alzheimer’s disease, prion diseases, and Type II diabetes, is a phenomenon related to protein folding in that it can be regarded as an alternate folding pathway [44]. Aggregation begins from the nucleation of proteins. This nucleation process has been elucidated using coarse-grained simulations, including CGMD. Nguyen and Hall reported that an amorphous aggregate of peptides was formed first, after which an ordered nucleus was formed, resulting in the lateral addition of β-sheets [45].

    In their CGMD model (designated PRIME), four pseudo atoms were mapped to each amino acid residue, hydrogen bonding was incorporated, a solvent was implicitly modeled, and discontinuous potential energy functions were employed [46]. Because the PRIME model allowed for alterations in the secondary structures of proteins, it was suitable for studying protein aggregation, as proteins undergo dynamic structural changes during this process. Discontinuous potential energy functions and implicit handling of solvents reduced the computational time and provided an advantage of sampling a vast conformational space.

    A capsid is the protein shell of a virus that can be simulated using CGMD simulations [47, 48, 49, 50]. Some types of viruses, including HIV, pose a threat to the life of their hosts. Thus, unraveling the processes involved in virus proliferation, including capsid formation, is important. In a huge molecular system that includes a virus capsid typically comprising hundreds of proteins, one CGMD atom unites a large number of amino acid residues (e.g., a protein subunit), and solvents are not incorporated. Although the definition of CGMD atoms (i.e., identifying CG sites) can be deduced from the structural requirements of virus capsids [51], this can also be achieved by AAMD simulations using X-ray crystallographic structures to parameterize constants, such as Lennard-Jones potential parameters [48]. The effects of the lengths of polyions on the formation of a huge virus capsid can be explained using highly coarse-grained models [50]. In this case, the effects of electrostatic interactions are naturally incorporated. In addition, because the effects of pH and salt concentration on capsid formation are thought to be non-negligible, in future applications these should be incorporated in CGMD models for more detailed simulations.

    CGMD simulations are also used for protein structure predictions. Although membrane proteins are important pharmacological targets, the three dimensional structures of most of these proteins remain unknown. Bucher et al. constructed a structural model for phospholipase A2 (PLA2) by combining homology modeling and CGMD simulations [18]. PLA2 is a membrane protein that releases fatty acids when it hydrolyzes phospholipids. These play important roles in intracellular signal transduction and inflammatory processes that are associated with Alzheimer’s disease, hypertensive heart failure, neurological diseases, and cancer.

    The appropriate localization of a protein within a lipid bilayer is important for a protein’s function. AAMD simulations require a large computational cost for the process of inserting a protein into a lipid bilayer, whereas CGMD simulations can handle this very efficiently. After the stabilization process of a protein within a lipid bilayer is simulated by CGMD, the coarse-grained geometry can be reverse grained to obtain all atom geometry. Then, extensive detailed analyses, including hydrogen bonding, can be achieved using AAMD simulations.

    Cholesterol is a crucial component of mammalian cells, as it determines the structural, thermodynamic, and mechanical properties of lipid membranes [52]. Flip-flop, namely the exchange of cholesterol molecules inside lipid bilayers, is important for efficient cholesterol tracking. This relatively time-consuming event (estimated half-time of <1 s) can be simulated using the Martini CGMD forcefield [53].

    3.2. Pharmacological/toxicological applications

    Applications of CGMD to in silico drug design are in progress. Lewis et al. used CGMD simulations to search for the optimal structure for an antiatherogenic agent [16]. The interactions between oxidized low-density lipoproteins (LDL) and scavenger receptors (SR) on the cell surfaces of macrophages are important phenomena during arterial stiffening. Amphiphilic macromolecules (AM) are known to competitively inhibit the interactions between oxidized LDL and SR [16].

    Because AMs can form highly variable structures, it is expensive to synthesize them chemically and to test their actions experimentally. Lewis et al. systematically generated different AM structures and estimated their actions using CGMD simulations. Using quantitative structure-activity relationships (QSAR) [54], a method for predicting the biological effects of a substance as per its chemical structure, they estimated the antiatherogenic effects of AMs based on their structures obtained from their CGMD simulations. AMs with high QSAR index values were experimentally confirmed to be effective because they exhibited strong interactions with SR and low oxidized LDL intake into macrophages.

    CGMD is suitable for searching the feature spaces of several agents. Even when using CGMD that can track the dynamics of a molecular system on a relatively large time scale, direct observations of a chemical reaction (AM-SR interactions in the example above) are usually impossible. Therefore, using ingenious methods like QSAR is a key component for applying CGMD to in silico drug design.

    If the target of a particular agent is located inside a cell, this agent should be able to pass through the lipid bilayer of a cell membrane. It has been shown that cationic nanoparticles (NPs) and cell-penetrating peptides (CPPs) are useful for delivering agents into cells. Determining the optimal size and the chemical composition of a carrier is crucial for developing a drug delivery system (DDS). Lin et al. suggested that the shapes of NPs and CPPs were important during their translocation through a lipid bilayer using CGMD simulations [4]. Furthermore, carbon nanotubes are used for DDS. Lee et al. used CGMD to analyze the processes involved when incorporating carbon nanotubes into a lipid bilayer [19].

    For gene therapy, DNA must be transferred into a cell. Lipofection is a method used to transfer DNA by coating DNA with lipids. Khalid et al. analyzed the interactions between DNA and lipids using CGMD [23].

    3.3. Protein-protein interactions

    A major application of CGMD is analyzing protein-protein interactions. The association and dissociation dynamics of an antibody (immunoglobulin G, IgG) and its affinity ligand Staphylococcus aureus protein A (SpA) were analyzed using CGMD [21]. The binding status of IgG with SpA was evaluated using an index of conformational change, root mean square deviation (RMSD), and the potential energy between these proteins. In addition, the detailed contributions of each amino acid residue to the potential energy could be obtained. Liu et al. reported that the dissociation dynamics in response to changes in pH consisted of four phases based on an analysis of the hydrophobic interactions and electrostatic interactions between charged residues [21]. To properly describe the protein-protein interactions involved, the effects of pH and Coulomb’s forces should be incorporated into a CGMD model, which in this case was the Martini forcefield.

    Oligomerization and self-aggregation of proteins are perfect targets for CGMD simulations because these molecular systems are relatively large and chemical actions play important roles within these systems. G protein-coupled receptors (GPCR) are important pharmacological targets, including photon receptor rhodopsins and cardiac adrenergic receptors. The functions of these proteins are significantly modified by the formation of oligomers within a lipid bilayer. Mondal et al. found that hydrophobic mismatches between a protein and lipids were crucial for the oligomerization of β-adrenergic receptors [15]. Another example is the protein Ras, which is associated with cancerous transformation of cells, as Ras forms clusters in a structure-dependent manner [17].

    4. Experimental Validation of CGMD Simulations

    We will now discuss methods used for comparing the results obtained from simulations with those experimentally obtained by roughly categorizing these into structure-oriented properties and dynamics-oriented properties. Length is a simple structural metric of a molecule. López et al. compared the thickness of a glycolipid membrane with that experimentally obtained [22]. Khalid et al. compared the persistence length (index of polymer hardness) of coarse-grained duplex DNA obtained from simulations with that experimentally obtained [23].

    The area per lipid is often used to analyze the properties of lipid membranes [22, 23, 26]. Lipid molecules can assume various phases, including liposomes, micelles, bilayer sheets, and hexagonal tubes, depending on environmental conditions like pH. Experimental observations of the hexagonal phase can be performed by freeze-fracture electron microscopy or X-ray diffraction. Dhalberg et al. compared the phase preferences and hexagonal spacing of lipid molecules obtained from simulations with those experimentally obtained [26]. Comparisons of the structural patterns of a molecular assembly discussed above are often used and applied to other cases, such as analyzing the pore diameters of ion channels and the distances between domains/molecules.

    In addition to the structural properties, dynamic properties are often used to compare the results from simulations and experiments. Arkhipov et al. compared viral capsid hardness by atomic force microscopy [47]. Liu et al. analyzed the binding affinities between two proteins under different salt concentrations and pH values and compared these with experimental results [21].

    With regard to the phase transitions of lipid molecules, the transition temperature is another property used for comparisons [26]. Khalid et al. confirmed that the concentration dependence of cationic lipids fit well with experimental results for DNA chain spacing [23]. Another sophisticated example for validating CGMD simulations is analyzing the solvent accessible surfaces of proteins [18]. This property can be measured by deuterium exchange mass spectrometry and compared with CGMD results.

    5. Conclusions

    There are certain biomolecular details that can only be revealed by CGMD simulations. CGMD is a powerful tool for simulating molecular systems with large spatial and temporal scales, where chemical interactions are crucial. Protein‑protein interactions, including self-assembly and oligomerization, protein‑lipid interactions, and even analyzing the effects of artificial substances, such as carbon nanotubes on biological systems are targets for CGMD simulations. CGMD simulation results can be validated by comparing the structural and dynamic properties of a molecular system with those experimentally obtained. Uniform CGMD methods to handle lipids, proteins, and nucleic acids are under development.

    As we have discussed in this article, applicable targets of CGMD simulations are rapidly expanding. However, CGMD alone is not sufficient to provide scientific and biological insights into the extensive molecular systems incorporated in a cell. A possible direction for the evolution of CGMD simulations is seamless cooperation with larger systems, such as the finite element method (FEM) based on continuum dynamics, and with smaller systems, such as AAMD. Approaches to associate simulations on different scales are currently in progress [55, 56, 57]. Combining CGMD and FEM is anticipated to be particularly useful for analyzing the mechanical properties of cellular components, such as the extracellular matrix, and their implications for cellular functions. Combinations of CGMD and AAMD will synergize each other’s strengths; namely, the accuracy of AAMD and the speed of CGMD. In conclusion, CGMD simulations will contribute to resolving problems in more varied situations than ever before.

    Acknowledgements

    The authors are very grateful to Mr. Kazuya Saruwatari at Okayama University Medical School for his CGMD simulations of the MscL protein.

    Conflicts of interest

    The authors declare that there were no conflicts of interest associated with this study.

    [1] Boddy L, Watkinson SC (1995) Wood decomposition, higher fungi, and their role in nutrient redistribution. Can J Bot 73: 1377–1383. doi: 10.1139/b95-400
    [2] Worrall JJ, Anagnost SE, Zabel RA, et al. (1997) Comparison of wood decay among diverse lignicolous fungi. Mycologia 89: 199–219. doi: 10.2307/3761073
    [3] Schwarze FWMR (2007) Wood decay under the microscope. Fungal Biol Rev 2: 133–170.
    [4] Mohebby B (2005) Attenuated total reflection infrared spectroscopy of white-rot decayed beech wood. Int Biodeter Biodegr 55: 247–251. doi: 10.1016/j.ibiod.2005.01.003
    [5] Boer W, Folman LB, Summerbell RC, et al. (2005) Living in a fungal world: Impact of fungi on soil bacterial niche development. FEMS Microbiol Rev 29: 795–811. doi: 10.1016/j.femsre.2004.11.005
    [6] Ingham RE, Trofymow JA, Ingham ER, et al. (1985) Interactions of bacteria, fungi, and their nematode grazers : Effects on nutrient cycling and plant growth. Ecol Monogr 55: 119–140. doi: 10.2307/1942528
    [7] Frey-Klett P, Burlinson P, Deveau A, et al. (2011) Bacterial-fungal interactions: Hyphens between agricultural, clinical, environmental, and food microbiologists. Microbiol Mol Biol R 75: 583–609. doi: 10.1128/MMBR.00020-11
    [8] Benner R, Newell SY, Maccubbin AE, et al. (1984) Relative contributions of bacteria and fungi to rates of degradation of lignocellulosic detritus in salt-marsh sediments. Appl Environ Microbiol 48: 36–40.
    [9] Shortle WC, Menge JA, Cowling EB (1978) Interaction of bacteria, decay fungi, and live sapwood in discoloration and decay of trees. Eur J Forest Pathol 8: 293–300. doi: 10.1111/j.1439-0329.1978.tb00642.x
    [10] Lang E, Kleeberg I, Zadrazil F (2000) Extractable organic carbon and counts of bacteria near the lignocellulose-soil interface during the interaction of soil microbiota and white rot fungi. Bioresource Technol 75: 57–65. doi: 10.1016/S0960-8524(00)00031-6
    [11] Romaní AM, Fischer H, Mille-lindblom C, et al. (2006) Interactions of bacteria and fungi on decomposing litter : Differential extracellular enzyme activities. Ecology 87: 2559–2569. doi: 10.1890/0012-9658(2006)87[2559:IOBAFO]2.0.CO;2
    [12] Lugtenberg B, Kamilova F (2009) Plant-growth-promoting Rhizobacteria. Annu Rev Microbiol 63: 541–556.
    [13] Clausen CA (1996) Bacterial associations with decaying wood: A review. Int Biodeter Biodegr 37: 101–107. doi: 10.1016/0964-8305(95)00109-3
    [14] Seidler RJ, Aho PE, Evans HJ, et al. (1972) Nitrogen fixation by bacterial isolates from decay in living white fir trees [Abies concolor (Gord. and Glend.) Lindl.]. J Gen Microbiol 73: 413–416. doi: 10.1099/00221287-73-2-413
    [15] Aho PE (1974) Distribution, enumeration, and identification of nitrogen-fixing bacteria associated with decay in living white fir trees. Phytopathology 64: 1413. doi: 10.1094/Phyto-64-1413
    [16] Bayer EA, Chanzy H, Lamed R, et al. (1998) Cellulose, cellulases and cellulosomes. Curr Opin Struc Biol 8: 548–557. doi: 10.1016/S0959-440X(98)80143-7
    [17] Watanabe H, Tokuda G (2010) Cellulolytic systems in insects. Annu Rev Entomol 55: 609–632. doi: 10.1146/annurev-ento-112408-085319
    [18] Breznak J, Burne A (1994) Role of microorganisms in the digestion of lignocellulose by termites. Annu Rev Entomol 39: 453–487.
    [19] Hall JB, Silver S (2009) Digestive system of the cow. Sciences-New York: 1–4.
    [20] Jami E, Mizrahi I (2012) Composition and similarity of bovine rumen microbiota across individual animals. PLoS One 7: 1–8.
    [21] Clements KD, Angert ER, Montgomery W, et al. (2014) Intestinal microbiota in fishes: What's known and what's not. Mol Ecol 23: 1891–1898. doi: 10.1111/mec.12699
    [22] Isbrücker IJH (1980) Classification and catalogue of the mailed Loricariidae (Pisces, Siluriformes). Verslagen en Technische Gegevens 22: 1–181.
    [23] Eschmeyer WN, Fricke R, Lann RVD, Catalog of fishes: Genera, species, references. Institute for Biodiversity Science & Sustainability, California Academy of Sciences, 2017. Available from: http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp.
    [24] German DP (2009) Inside the guts of wood-eating catfishes: Can they digest wood? J Comp Physiol B 179: 1011–1023. doi: 10.1007/s00360-009-0381-1
    [25] Lujan NK, German DP, Winemiller KO (2011) Do wood-grazing fishes partition their niche?: Morphological and isotopic evidence for trophic segregation in Neotropical Loricariidae. Funct Ecol 25: 1327–1338.
    [26] German DP, Bittong RA (2009) Digestive enzyme activities and gastrointestinal fermentation in wood-eating catfishes. J Comp Physiol B 179: 1025–1042. doi: 10.1007/s00360-009-0383-z
    [27] Araujo-Lima CA, Forsberg BR, Victoria R, et al. (1986) Energy sources for detritivorous fishes in the Amazon. Science 234: 1256–1258. doi: 10.1126/science.234.4781.1256
    [28] Nelson JA, Wubah DA, Whitmer ME, et al. (1999) Wood-eating catfishes of the genus Panaque: Gut microflora and cellulolytic enzyme activities. J Fish Biol 54: 1069–1082.
    [29] Di MN, Schwarzentruber P, Schenker M, et al. (2013) Microbial population dynamics in the faeces of wood-eating loricariid catfishes. Lett Appl Microbiol 56: 401–407. doi: 10.1111/lam.12061
    [30] McDonald R, Schreier HJ, Watts JEM (2012) Phylogenetic analysis of microbial communities in different regions of the gastrointestinal tract in Panaque nigrolineatus, a wood-eating fish. PLoS One 7: e48018. doi: 10.1371/journal.pone.0048018
    [31] Watts JEM, McDonald R, Daniel R, et al. (2013) Examination of a culturable microbial population from the gastrointestinal tract of the wood-eating loricariid catfish Panaque nigrolineatus. Diversity 5: 641–656. doi: 10.3390/d5030641
    [32] McDonald R, Zhang F, Watts JEM, et al. (2015) Nitrogenase diversity and activity in the gastrointestinal tract of the wood-eating catfish Panaque nigrolineatus. ISME J 9: 1–13. doi: 10.1038/ismej.2014.99
    [33] Yoshimizu M, Kimura T (1976) Study on the intestinal microflora of salmonids. Fish Pathol 10: 243–259. doi: 10.3147/jsfp.10.243
    [34] Ugajin M (1976) Studies on the taxonomy of major microflora on the intestinal contents of salmonoids. Bull Jpn Soc Sci Fish 45: 721–731.
    [35] Holben WE, Williams P, Saarinen M, et al. (2002) Phylogenetic analysis of intestinal microflora indicates a novel Mycoplasma phylotype in farmed and wild salmon. Microbial Ecol 44: 175–185. doi: 10.1007/s00248-002-1011-6
    [36] Desai AR, Links MG, Collins SA, et al. (2012) Effects of plant-based diets on the distal gut microbiome of rainbow trout (Oncorhynchus mykiss). Aquaculture 350: 134–142.
    [37] Kamei Y, Sakata T, Kakimoto D (1985) Microflora in the alimentary tract of tilapia: Characterization and distribution of anaerobic bacteria. J Gen Appl Microbiol 31: 115–124. doi: 10.2323/jgam.31.115
    [38] Wu S, Wang G, Angert ER, et al. (2012) Composition, diversity, and origin of the bacterial community in grass carp intestine. PLoS One 7: e30440. doi: 10.1371/journal.pone.0030440
    [39] Dİler Ö, Dİler A (1998) Quantitative and qualitative changes of the gastrointestinal microflora of Pike-perch (Stizostedion Lucioperca L. 1758) in Egirdir Lake. Turk J Vet Anim Sci 22: 325–328.
    [40] Roeselers G, Mittge EK, Stephens WZ, et al. (2011) Evidence for a core gut microbiota in the zebrafish. ISME J 5: 1595–1608. doi: 10.1038/ismej.2011.38
    [41] Wu S, Gao T, Zheng Y, et al. (2010) Microbial diversity of intestinal contents and mucus in yellow catfish (Pelteobagrus fulvidraco). Aquaculture 303: 1–7. doi: 10.1016/j.aquaculture.2009.12.025
    [42] Nieto TP, Toranzo AE, Barja JL (1984) Comparison between the bacterial flora associated with fingerling rainbow trout cultured in two different hatcheries in the North-West of Spain. Aquaculture 42: 193–206. doi: 10.1016/0044-8486(84)90100-5
    [43] Andlid T, Juárez RV, Gustafsson L (1995) Yeast colonizing the intestine of rainbow trout (Salmo gairdneri) and turbot (Scophtalmus maximus). Microbial Ecol 30: 321–334.
    [44] Andlid T, Vazquez-Juarez R, Gustafsson L (1998) Yeasts isolated from the intestine of rainbow trout adhere to and grow in intestinal mucus. Mol Mar Biol Biotechnol 7: 115–126.
    [45] Gatesoupe FJ (2007) Live yeasts in the gut: Natural occurrence, dietary introduction, and their effects on fish health and development. Aquaculture 267: 20–30. doi: 10.1016/j.aquaculture.2007.01.005
    [46] Laconi S, Pompei R, (2007) Study and characterization of intestinal yeasts of mullet (Mugil spp.) for potential probiotic use. J Food Agr Environ 5: 475–480.
    [47] Raggi P, Lopez P, Diaz A, et al. (2014) Debaryomyces hansenii and Rhodotorula mucilaginosa comprised the yeast core gut microbiota of wild and reared carnivorous salmonids, croaker and yellowtail. Environ Microbiol 16: 2791–2803. doi: 10.1111/1462-2920.12397
    [48] Gardes M, Bruns TD (1993) ITS primers with enhanced specificity for basidiomycetes-application to the identification of mycorrhizae and rusts. Mol Ecol 2: 113–118. doi: 10.1111/j.1365-294X.1993.tb00005.x
    [49] Ihrmark K, Bödeker ITM, Cruz-Martinez K, et al. (2012) New primers to amplify the fungal ITS2 region-evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol Ecol 82: 666–677. doi: 10.1111/j.1574-6941.2012.01437.x
    [50] Schloss PD, Westcott SL, Ryabin T, et al. (2009) Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75: 7537–7541. doi: 10.1128/AEM.01541-09
    [51] Edgar RC (2010) Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26: 2460–2461. doi: 10.1093/bioinformatics/btq461
    [52] Edgar RC, Haas BJ, Clemente JC, et al. (2011) UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27: 2194–2200. doi: 10.1093/bioinformatics/btr381
    [53] Fu L, Niu B, Zhu Z, et al. (2012) CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 28: 3150–3152. doi: 10.1093/bioinformatics/bts565
    [54] Koljalg U, Nilsson RH, Abarenkov K, et al. (2013) Towards a unified paradigm for sequence-based identification of fungi. Mol Ecol 22: 5271–5277. doi: 10.1111/mec.12481
    [55] Caporaso JG, Kuczynski J, Stombaugh J, et al. (2010) Qiime allows analysis of high-throughput community sequencing data. Nat Rev Microbiol 7: 335–336.
    [56] Chao A (1984) Nonparametric estimation of the number of classes in a population. Scand J Stat 11: 265–270.
    [57] Prosser JI (2010) Replicate or lie. Environ Microbiol 12: 1806–1810. doi: 10.1111/j.1462-2920.2010.02201.x
    [58] Hiscox J, Savoury M, Müller CT, et al. (2015) Priority effects during fungal community establishment in beech wood. ISME J 9: 1–15. doi: 10.1038/ismej.2014.99
    [59] Purahong W, Wubet T, Lentendu G, et al. (2016) Life in leaf litter: novel insights into community dynamics of bacteria and fungi during litter decomposition. Mol Ecol 25: 4059–4074. doi: 10.1111/mec.13739
    [60] Hervé V, Le RX, Uroz S, et al. (2014) Diversity and structure of bacterial communities associated with Phanerochaete chrysosporium during wood decay. Environ Microbiol 16: 2238–2252. doi: 10.1111/1462-2920.12347
    [61] Bengtsson G (2012) International association for ecology interactions between fungi, bacteria and beech leaves in a stream microcosm. Oecologia 89: 542–549.
    [62] Lee SS, Ha JK, Cheng KJ (2000) Relative contributions of bacteria, protozoa, and fungi to in vitro degradation of orchard grass cell walls and their interactions. Appl Environ Microbiol 66: 3807–3813. doi: 10.1128/AEM.66.9.3807-3813.2000
    [63] Mouzouras R, Jones EBG, Venkatasamy R, et al. (1986) Decay of wood by micro-organisms in marine environments. BWPA Annual Convention.
    [64] Björdal CG (2012) Evaluation of microbial degradation of shipwrecks in the Baltic Sea. Int Biodeter Biodegr 70: 126–140. doi: 10.1016/j.ibiod.2012.01.012
    [65] Alconada TM, Martinez MJ (1996) Purification and characterization of a beta-glucosidase from the phytopathogenic fungus Fusarium oxysporum f. sp. melonis.Lett Appl Microbiol 22: 106–110. doi: 10.1111/j.1472-765X.1996.tb01120.x
    [66] Meyer V, Andersen MR, Brakhage AA, et al. (2016) Current challenges of research on filamentous fungi in relation to human welfare and a sustainable bio-economy: a white paper. Fungal Biol Biotechnol 3: 6. doi: 10.1186/s40694-016-0024-8
    [67] Wang XC, Liu C, Huang L, et al. (2015) ITS1: A DNA barcode better than ITS2 in eukaryotes? Mol Ecol Resour 15: 573–586. doi: 10.1111/1755-0998.12325
    [68] Nilsson RH, Ryberg M, Kristiansson E, et al. (2006) Taxonomic reliability of DNA sequences in public sequences databases: A fungal perspective. PLoS One 1: e59. doi: 10.1371/journal.pone.0000059
    [69] Huber JA, Morrison HG, Huse SM, et al. (2009) Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure. Environ Microbiol 11: 1292–1302. doi: 10.1111/j.1462-2920.2008.01857.x
    [70] Fonseca VG, Nichols B, Lallias D, et al. (2012) Sample richness and genetic diversity as drivers of chimera formation in nSSU metagenetic analyses. Nucleic Acids Res 40: e66. doi: 10.1093/nar/gks002
    [71] Lindahl B, Nilsson RH, Tedersoo L, et al. (2013) Fungal community analysis by high-throughput sequencing of amplified markers-a user's guide. New Phytol 199: 288–299. doi: 10.1111/nph.12243
    [72] Lujan NK, Hidalgo M, Stewart DJ (2010) Revision of Panaque (Panaque), with descriptions of three new species from the Amazon Basin (Siluriformes, Loricariidae). Copeia 2010: 676–704. doi: 10.1643/CI-09-185
  • This article has been cited by:

    1. Petru Cardei, Constantina Chireceanu, A Differential Mathematical Model for Experiments to Determine the Efficacy of Treatments Against the Bean Weevil, 2022, 20, 2224-2902, 12, 10.37394/23208.2023.20.2
  • Reader Comments
  • © 2017 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(7303) PDF downloads(1070) Cited by(11)

Article outline

Figures and Tables

Figures(3)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog