Selectivity and transport in aquaporins from molecular simulation studies.

Molecular selection, ion exclusion, and water permeation are well known regulatory mechanisms in aquaporin. Water permeability was found to be diverse in different subgroups of plasma membrane intrinsic proteins (PIPs), even though the residues surrounding the water holes remained the same across the subgroups. Upon homology modeling and structural comparison, a conserved Ala/Ile(Val) residue difference was identified in helix 2 that affected the conformation of the NPA region and consequently influenced the water permeability. The residue difference was found to be conservative within the two subgroups of PIPs in rice as well as in other plants. Functional tests further confirmed the prediction via site-directed mutagenesis where replacement of Ala(103) or Ala(102) in respective OsPIP1;1 or OsPIP1;3 with Val yielded 7.0- and 2.2-fold increases in water transportation, and substitution of Ile(98) or Val(95) in respective OsPIP2;3 or OsPIP2;7 with Ala resulted in 73 or 52% reduction of water transportation. Based on structural analyses and molecular dynamics simulations, we proposed that the difference in water permeability was attributed to the orientation variations of helix 2 that modified water-water and water-protein interactions.

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

The term molecular mechanics (MM) refers to the use of simple potential-energy functions (e.g., harmonic oscillator or Coulombic potentials) to model molecular systems. Molecular mechanics approaches are widely applied in molecular structure refinement, molecular dynamics (MD) simulations, Monte Carlo (MC) simulations, and ligand-docking simulations. Typically, molecular mechanics models consist of spherical atoms connected by springs which represent bonds. Internal forces experienced in the model structure are described using simple mathematical functions. For example, Hooke’s law is commonly used to describe bonded interactions, and the nonbonded atoms might be treated as inelastic hard spheres or may interact according to a Lennard-Jones potential. Using these simple models, a molecular dynamics simulation numerically solves Newton’s equations of motion, thus allowing structural fluctuations to be observed with respect to time. Dynamic simulation methods are widely used to obtain information on the time evolution of conformations of proteins and other biological macromolecules1–4 and also kinetic and thermodynamic information. Simulations can provide fine detail concerning the motions of individual particles as a function of time. They can be utilized to quantify the properties of a system at a precision and on a time scale that is otherwise inaccessible, and simulation is, therefore, a valuable tool in extending our understanding of model systems. Theoretical consideration of a system additionally allows one to investigate the specific contributions to a property through “computational alchemy”,5 that is, modifying the simulation in a way that is nonphysical but nonetheless allows a model’s characteristics to be probed. One particular example is the artificial conversion of the energy function from that representing one system to that of another during a simulation. This is an important technique in free-energy calculations.6 Thus, molecular dynamics simulations, along with a range of complementary computational approaches, have become valuable tools for investigating the basis of protein structure and function. This review offers an outline of the origin of molecular dynamics simulation for protein systems and how it has developed into a robust and trusted tool. This review then covers more recent advances in theory and an illustrative selection of practical studies in which it played a central role. The range of studies in which MD has played a considerable or pivotal role is immense, and this review cannot do justice to them; MD simulations of biomedical importance were recently reviewed.4 Particular emphasis will be placed on the study of dynamic aspects of protein recognition, an area where molecular dynamics has scope to provide broad and far-ranging insights. This review concludes with a brief discussion of the future potential offered to advancement of the biological and biochemical sciences and the remaining issues that must be overcome to allow the full extent of this potential to be realized. 1.1. Historical Background MD methods were originally conceived within the theoretical physics community during the 1950s. In 1957, Alder and Wainwright7 performed the earliest MD simulation using the so-called hard-sphere model, in which the atoms interacted only through perfect collisions. Rahman8 subsequently applied a smooth, continuous potential to mimic real atomic interactions. During the 1970s, as computers became more widespread, MD simulations were developed for more complex systems, culminating in 1976 with the first simulation of a protein9,10 using an empirical energy function constructed using physics-based first-principles assumptions. MD simulations are now widely and routinely applied and especially popular in the fields of materials science11,12 and biophysics. As will be discussed later in this review, a variety of experimental conditions may be simulated with modern theories and algorithms. The initial simulations only considered single molecules in vacuo. Over time, more realistic or at least biologically relevant simulations could be performed. This trend is continuing today. The initial protein MD simulation, of the small bovine pancreatic trypsin inhibitor (BPTI), covered only 9.2 ps of simulation time. Modern simulations routinely have so-called equilibration periods much longer than that, and production simulations of tens of nanoseconds are routine, with the first microsecond MD simulation being reported in 1998.13 In addition, the original BPTI simulation included only about 500 atoms rather than the 104-106 atoms that are common today. While much of this advancement results from an immense increase in availability of computing power, major theoretical and methodological developments also contribute significantly. The number of publications regarding MD theory and application of MD to biological systems is growing at an extraordinary pace. A single review cannot do justice to the recent applications of MD. Using data from ISI Web of Science, the authors estimate that during 2005 at least 800 articles will be published that discuss molecular dynamics and proteins. The historical counts are shown in Figure 1. Open in a separate window Figure 1 Articles matching ISI Web of Science query “TS=(protein) AND TS=(molecular dynamics)”.

The diffusion of five flavor organic molecules, including D-limonene, myrcene, ethyl hexanoate, 2-nonanone, and linalool in low density polyethylene (LDPE) film were investigated by combined experimental and molecular dynamics (MD) simulation studies. The diffusion coefficients derived from the prediction model, experimental determination, and MD simulation were compared, and the related microscopic diffusion mechanism was investigated. The effects of the presence of the flavor organic molecules on the diffusion of O2 and H2O in polyethylene (PE) were also studied by MD simulation. Results show that: The diffusion of five flavor molecules in LDPE is well followed to Fick’s second law by the immersion experiment; MD simulation indicates the dual mode diffusion mechanism of the flavor molecules is in LDPE; the diffusion coefficients from MD simulation are close to those from the experimental determination, but are slightly larger than those values; the presence of the flavor organic molecules hinders the diffusion of O2 and H2O, which can be well explained from the fraction of free volume (FFV) and interaction energy calculation results derived from MD simulation.

A combined study of experimental and molecular dynamics (MD) simulation methods is presented for hindered phenol AO‐80/nitrile‐butadiene rubber/poly(vinyl chloride) (AO‐80/NBR/PVC) composites with different AO‐80 contents to establish the microstructure‐damping property relations. MD simulation found that the AO‐80/NBR/PVC composite (abbreviated as AO‐80/NBVC) with an AO‐80 content of 99 phr had the largest hydrogen bonds (H‐bonds) and highest binding energy, indicating a good compatibility between NBR and AO‐80 and good damping performance of AO‐80/NBVC composites. Experimental results from SEM, DSC, and DMA were in good agreement with the MD simulation results. The tensile test results showed that the AO‐80/NBVC composite with an AO‐80 content of 99 phr had high tensile strength because of the strong H‐bonds of the composites and the disintegration and reintegration of the H‐bonds. The MD simulation technique proves to be a promising tool for the design and prediction of high damping properties of advanced composites in a microscopic view.

Wettability of organic-rich mudrocks has significant impacts on multiphase fluid flow and hydrocarbon recovery. Since kerogen can constitute a significant volumetric fraction of mudrocks, the wettability of kerogen can considerably affect the wettability of organic-rich mudrocks. In a recent project, we showed using a combination of experimental methods and molecular simulation studies that kerogen type, thermal maturity, and reservoir temperature conditions affect the contact angle of water droplet on kerogen surface. This indicates that the wettability of kerogen is variable with its geochemistry and reservoir conditions. However, these studies do not consider the impact of having a third phase hydrocarbon in the reservoir system. Moreover, wettability affects the distribution of fluids in the porous media and the relative fluid mobility of reservoir fluids. However, the impact of kerogen wettability on fluid mobility is yet to be understood at a molecular level. Therefore, in this research project, we plan to use molecular dynamics (MD) simulations to determine the impact of kerogen geochemistry and reservoir conditions on three phase contact angles of water and hydrocarbon on kerogen surface. We also use MD simulations to determine the impact of wettability on mobility of fluids in kerogen pores. We use realistic molecular models of kerogen. The chemical structure and composition of these molecular models are verified against experimental data for kerogen including Fourier transform infrared (FTIR) measurements, X-ray photoelectron spectroscopy (XPS) measurements and X-ray fluorescence (XRF) measurements. The results obtained from MD simulation will be verified against experientially measured three-phase contact angles on kerogen surface. Furthermore, wettability of kerogen surface impacts the interaction between the fluid and the kerogen pores and, therefore, affects the relative permeability/mobility of different fluids in the kerogen pore space. In this work, we simulate water and oil transport through the kerogen pores of varying hydrophilicities by using MD simulations to determine the impact of kerogen wettability on the fluid mobilities.

Introduction: Milk kefir, a fermenting milk made with kefir grains, has shown potential in promoting apoptosis, regulating the cell cycle, and reducing tumor growth in breast cancer cells. This study aimed to investigate the stability and potential of milk kefir metabolites as inhibitors of breast cancer growth by interacting with the estrogen receptor alpha (ER-α), a key protein involved in breast cancer cell proliferation. We used computational methods, specifically molecular docking simulations with AutoDock and molecular dynamics (MD) simulations with Gromacs, to analyze how these metabolites bind to ER-α. Methods: A combination of molecular docking and MD simulations was used to explore how metabolites derived from milk kefir interact with ER-α, a crucial target in breast cancer therapy. The methodology included multiple stages: preparation of target proteins, preparation and screening of the metabolites, geometry optimization, molecular docking, and MD simulations. Results: The molecular docking simulations of 43 metabolites revealed three promising candidates: 2-Methyl (S35), benzeneethanol (S42), and 2,6-dimethyl-4-heptanone (S54), with binding affinities (ΔG) of -5.08, -5.06, and -4.90 kcal/mol, respectively. MD simulations further showed that the selected metabolites stabilized the ER-α-metabolite complex, with the 2,6-dimethyl-4-heptanone (S54) metabolite demonstrating the most negative total MM-GBSA energy value (ΔG = -22.98 kcal/mol), indicating a strong and stable binding interaction. Conclusion: 2,6-Dimethyl-4-heptanone, a metabolite from milk kefir, showed promising potential as a candidate for further development as a breast cancer treatment, offering a novel alternative to conventional therapies.

Coarse-grained (CG) molecular dynamics (MD) simulations are performed to study the insertion of cyclic peptide nanotubes into cell membranes and to examine whether cyclic peptide nanotubes can function as an ion channel and thereby as an antibacterial agent. To do so, the two coarse-grained (CG) models for lipid molecules and for proteins developed by Marrink et al. (J. Phys. Chem. B 2004, 108, 750) and by Shih et al. (J. Phys. Chem. B 2006, 110, 3674), respectively, were extended and modified. These CG models were verified by performing CG MD and all-atom (AA) MD simulations for a cyclic peptide nanotube, 8 x cyclo[(-d-Ala-l-Glu-d-Ala-l-Gln-)2], in water and by comparing the results from the two simulations. Comparison between static and dynamic (water transport) properties obtained from both simulations shows good agreement. To study nanotube insertion, a CG cyclic peptide nanotube, 8 x cyclo[(-Trp-d-Leu-)4], was prepared above the surface of a CG DPPC lipid bilayer, restrained with constraints, and equilibrated, and then a series of CG MD simulations were carried out by lifting the constraints imposed on the nanotube. The CG MD simulations show that the cyclic peptide nanotube spontaneously inserts into and reorients inside the lipid bilayer. After insertion, the long axis of the cyclic peptide nanotube is aligned approximately perpendicular to the bilayer plane indicating that the nanotube can function as an ion channel and as an antibacterial agent. Tilt structures of the cyclic peptide nanotubes inside the lipid bilayer are found to be in agreement with experiment and earlier AA simulations. Lipid flip-flop, a migration of lipid molecules from one leaflet to the other leaflet of the lipid bilayer, is also observed from the CG MD simulations. Finally, the CG MD simulations reveal that a lipid headgroup can be inserted into the cyclic peptide nanotube. This process is confirmed by an AA MD simulation.

Gastric H(+),K(+)-ATPase is shown to transport 2 mol of H(+)/mol of ATP hydrolysis in isolated hog gastric vesicles. We studied whether the H(+) transport mechanism is due to charge transfer and/or transfer of hydronium ion (H(3)O(+)). From transport of [(18)O]H(2)O, 1.8 mol of water molecule/mol of ATP hydrolysis was found to be transported. We performed a molecular dynamics simulation of the three-dimensional structure model of the H(+),K(+)-ATPase alpha-subunit at E(1) conformation. It predicts the presence of a charge transfer pathway from hydronium ion in cytosolic medium to Glu-345 in cation binding site 2 (H(3)O(+)-Lys-164 -Gln-161-Glu-345). No charge transport pathway was formed in mutant Q161L, E345L, and E345D. Alternative pathways (H(3)O(+)-Gln-161-Glu-345) in mutant K164L and (H(3)O(+)-Arg-105-Gln-161-Gln-345) in mutant E345Q were formed. The H(+),K(+)-ATPase activity in these mutants reflected the presence and absence of charge transfer pathways. We also found charge transfer from sites 2 to 1 via a water wire and a charge transfer pathway (H(3)O(+)-Asn-794 -Glu-797). These results suggest that protons are charge-transferred from the cytosolic side to H(2)O in sites 2 and 1, the H(2)O comes from cytosolic medium, and H(3)O(+) in the sites are transported into lumen during the conformational transition from E(1)PtoE(2)P.

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.

Herbal plants are often used as alternative medicine because they contain active compounds for the treatment of diseases and disorders with minimal side effects, and are easily obtained from the surrounding environment. Some of them have antiviral activity. This study aimed to analyze the potential of phytochemical compounds in the leaf of “pangi” (Pangium edule Reinw) as antiSARS-CoV-2 using molecular docking study. The drug- and lead-likeness properties of the selected compounds were obtained from the Swiss ADME and admetSAR online server tools. Molecular dynamics (MD) simulation of the selected ligand was carried out to validate the stability of the interaction. The results suggested that pangi leaves contain three compounds with remarkable binding affinities with Mpro (main protease) and RBD (receptor binding domain) were (5.beta.) pregnane-3,20.beta.-diol, 14.alpha.,18.alpha.-[4-methyl-3-oxo-(1-oxa-4-azabutane-1,4-diyl)]-, diacetate (PD), ethyl cholate (EC), and bis(3,5,5-trimethylhexyl) phthalate. Because EC will be metabolized in the body into cholic acid (Cho), this compound was then docked and validated using MD simulation. The compound has the best free binding energy (ΔG) with SARS-CoV-2 (–7.1 kcal/ mol with Mpro and –6.0 kcal/mol with RBD). Moreover, the compound is bound strongly to the active cavity of Mpro on Thr24, Thr26, His41, and Cys145 residues. The MM-GBSA calculation showed that the interaction of Cho with Mpro was higher than with RBD. According to the RMSD (root mean square deviation), RMSF (root mean square fluctuation), the radius of gyration (Rg), and intermolecular hydrogen bond (H-bond) analysis obtained from 50 ns MD simulations, Cho formed stable interactions with Mpro and RBD. The finding of this study indicated that Cho showed good anti-SARS-CoV-2 activity. The potential of the compound to inhibit the virus can serve as a starting point in the process of developing COVID-19 therapeutic natural medicine.

Present knowledge obtained by molecular dynamics (MD) simulation studies regarding the dynamics of water, both in the vicinity of biological membranes and within the proteinaceous water channels, also known as aquaporins (AQPs), is reviewed. A brief general summary of the water models most extensively employed in MD simulations (SPC, SPC/E, TIP3P, TIP4P), indicating their most relevant pros and cons, is likewise provided. Structural considerations of water are also discussed, based on different order parameters, which can be extracted from MD simulations as well as from experiments. Secondly, the behaviour of water in the neighbourhood of membranes by means of molecular dynamics simulations is addressed. Consequently, the comparison with previous experimental evidence is pointed out. In living cells, water is transported across the plasma membrane through the lipid bilayer and the aforementioned AQPs, which motivates this review to focus mostly on MD simulation studies of water within AQPs. Relevant contributions explaining peculiar properties of these channels are discussed, such as selectivity and gating. Water models used in these studies are also summarised. Finally, based on the information presented here, further MD studies are encouraged.

Colorectal cancer (CRC) is the third leading cause of cancer-related deaths in both men and women. Toll-like receptor 5 (TLR5), an autoimmune signaling receptor that plays a role in cancer, can be exploited for the suppression of human colon cancer. Salmonella flagellin protein, a novel agonist of TLR5 activating downstream signaling, could be a basis for designing anticancer peptides. The three-dimensional crystal structure of TLR5 (PDB ID: 3J0A, Resolution = 26.0 Å) was optimized using the AMBER force field in the YASARA suit. In silico enzymatic digestion tool, PeptideCutter, was used to identify peptides from Salmonella flagellin, an agonist against human TLR5. The 3D structure of the peptides was generated using PEP-FOLD3. These peptides were screened against human TLR5 using shape complementarity principles based on the binding affinity and interactions with the active residue of TLR5 monomer, and the selected peptides were further validated by molecular dynamic (MD) simulation. In this study, we generated 42 peptides from Salmonella flagellin protein by in silico protein digestion. Then, based on a new hidden Markov model sub-optimal conformation sampling approach as well as the size of the fragments, we select 38 effective peptides from these 42 cleavages. These peptides were screened against the monomeric Xray structure of human TLR5 using shape complementarity principles. Based on the binding affinity and interactions with the active residue of TLR5 monomer (residues 294 and 366 of TLR5), nine top-scored peptides were selected for the initial molecular dynamic (MD) simulation. Among these peptides, Clv10, Clv17, and Clv28 showed high stability and less flexibility during MD simulation. A 1 μs MD simulation was performed on TLR5-Clv10, TLR-Clv17, and TLR5-Clv28 complexes to further analyze the stability, conformational changes, and binding mode (Clv10, Clv17, and Clv28). During this MD study, the peptides showed high salt bridges and ionic interactions with residue ASP294 and residue ASP366 throughout the simulation and remained in the concave of the human TLR5 monomer. The RMSD and Rg values showed that the peptide-protein complexes become stable after 200 ns of contraction and extraction. These findings can facilitate the rational design of selected peptides as an agonist of TLR5, which have antitumor activity, suppress colorectal cancer tumors, and can be used as promising candidates and novel agonists of TLR5.

HighlightsAccording to the design of the experiment (DoE), multivariate analysis models were used to optimize the critical process parameters combined with multi-objective optimization.Based on the optimized operating conditions, the MILT-HD method not only enhances the extraction efficiency from Amomi fructus but also reduces energy demands and CO2 emissions.Based on the density functional theoretical (DFT) and molecular dynamics (MD) simulations, the mechanisms for ionic liquids (ILs) to improve the extraction efficiency of essential oil was comprehensively revealed.In this paper, Amomi fructus (Latin) was used to explore the mechanism of ionic liquids (ILs) in improving the extraction efficiency of essential oils. Microwave assisted ionic liquid treatment followed by a hydro-distillation (MILT-HD) process for isolating Amomi fructus essential oil was optimized by multi-objective optimization. Under optimum operating conditions, the IL-assisted extraction method not only enhances extraction efficiency but also reduces energy demands and CO2 emissions. Since the hydrogen bond structure network of cellulose in the cell wall is an important reason for hindering diffusion of essential oils, the mechanism of ILs was explored by density functional theoretical (DFT) and molecular dynamics (MD) simulations. According to DFT calculations, ILs can facilitate the cleavage of cellulose chains and have strong non-covalent interactions with cellulose. Based on the MD simulations, the degree of destruction of the cellulose hydrogen bond structure was explored. According to the DFT and MD simulations, the ILs can significantly destroy cellulose structure, thereby promoting essential oil release from the plant. These results were confirmed by scanning electron microscopy (SEM) and Fourier transform infrared spectroscopy (FTIR). This work is conducive to better understand the MILT-HD process for isolating essential oil and comprehensively understand the mechanism of ILs in the extraction process.

We combine molecular dynamics (MD) and Monte Carlo (MC) simulations to study the properties of gas molecules confined between two hard walls of a microchannel or nanochannel. The coupling between MD and MC simulations is introduced by performing MD near the boundaries for accuracy and MC in the bulk because of the low computational cost. We characterize the influence of different densities and molecule sizes on the equilibrium properties of the gas in the microchannel. The effect of the particle size on the simulation results is very small in the case of a dilute gas and increases with the density. The hybrid MD-MC simulation method is validated by comparing the results for density and temperature profiles with those of pure MD and pure MC simulations. These results compare well for pure MD and pure MC, as well as hybrid MD-MC, both in the bulk and near the boundaries, when hard-sphere interactions are used. When Lennard-Jones potentials are used to accurately model the interactions between the gas and wall molecules instead, the results of pure MD simulations differ significantly from the pure MC simulations near the boundaries, but the results of the hybrid method compare well with the pure MD results near the wall, and with the pure MC and pure MD results in the middle of the channel. The hybrid method also very accurately simulates the interface between the MD and MC simulation domains. Comparisons between MD, MC, and hybrid MD-MC computational costs are outlined. The speedup when using 50% of the domain for MD simulations and 50% for MC simulations is very small compared to pure MD simulations times, but this speedup increases drastically for more realistic situations where the region near the wall is small compared to the bulk region.