Computing Individual Area per Head Group Reveals Lipid Bilayer Dynamics.

Lipid bilayers express a range of phases from solid-like to gel-like to liquid-like as a function of temperature and lipid surface concentration. The area occupied per lipid head group serves as one useful indicator of the bilayer phase, in conjunction with the two-dimensional radial distribution function (i.e., structure factor) within the bilayer. Typically, the area per head group is determined by dividing the bilayer area equally among all head groups. Such an approach is less satisfactory for a multicomponent set of diverse lipids. In this work, area determination is performed on a lipid-by-lipid basis by attributing to a lipid the volume that surrounds each atom. Voronoi tessellation provides this division of the interfacial region on a per-atom basis. The method is applied to a multicomponent system of water, NaCl, and 19 phospholipid types that was devised recently [Langmuir2022, 38, 9481-9499] as a computational representation of the Gram-positive Staphylococcus aureus phospholipid bilayer. Results demonstrate that lipids and water molecules occupy similar extents of area within the interfacial region; ascribing area only to head groups implicitly incorporates assumptions about head group hydration. Results further show that lipid tails provide non-negligible contributions to area on the membrane side of the bilayer-water interface. Results for minimum and maximum area of individual lipids reveal that spontaneous fluctuations displace head groups more than 10 Šfrom the interfacial region during an NPT simulation at 310 K, leading to a zero contribution to total area at some times. Total area fluctuations and fluctuations per individual lipid relax with a correlation time of ∼10 ns. The method complements density profile as an approach to quantify the structure and dynamics of computational lipid bilayers.

Generation and Computational Characterization of a Complex Staphylococcus aureus Lipid Bilayer.

Studies indicate a crucial cell membrane role in the antibiotic resistance of Staphylococcus aureus. To simulate its membrane structure and dynamics, a complex molecular-scale computational representation of the S. aureus lipid bilayer was developed. Phospholipid types and their amounts were optimized by reverse Monte Carlo to represent characterization data from the literature, leading to 19 different phospholipid types that combine three headgroups [phosphatidylglycerol, lysyl-phosphatidylglycerol (LPG), and cardiolipin] and 10 tails, including iso- and anteiso-branched saturated chains. The averaged lipid bilayer thickness was 36.7 Å, and area per headgroup was 67.8 Å2. Phosphorus and nitrogen density profiles showed that LPG headgroups tended to be bent and oriented more parallel to the bilayer plane. The water density profile showed that small amounts reached the membrane center. Carbon density profiles indicated hydrophobic interactions for all lipids in the middle of the bilayer. Bond vector order parameters along each tail demonstrated different C-H ordering even within distinct lipids of the same type; however, all tails followed similar trends in average order parameter. These complex simulations further revealed bilayer insights beyond those attainable with monodisperse, unbranched lipids. Longer tails often extended into the opposite leaflet. Carbon at and beyond a branch showed significantly decreased ordering compared to carbon in unbranched tails; this feature arose in every branched lipid. Diverse tail lengths distributed these disordered methyl groups throughout the middle third of the bilayer. Distributions in mobility and ordering reveal diverse properties that cannot be obtained with monodisperse lipids.

Computational Study of Helical and Helix-Hinge-Helix Conformations of an Anti-Microbial Peptide in Solution by Molecular Dynamics and Vibrational Analysis.

Many classical antimicrobial peptides adopt an amphipathic helical structure at a water-membrane interface. Prior studies led to the hypothesis that a hinge near the middle of a helical peptide plays an important role in facilitating peptide-membrane interactions. Here, dynamics and vibrations of a designed hybrid antimicrobial peptide LM7-2 in solution were simulated to investigate its hinge formation. Molecular dynamics simulation results on the basis of the CHARMM36 force field showed that the α-helix LM7-2 bent around two or three residues near the middle of the peptide, stayed in a helix-hinge-helix conformation for a short period of time, and then returned to a helical conformation. High-resolution computational vibrational techniques were applied on the LM7-2 system when it has α-helical and helix-hinge-helix conformations to understand how this structural change affects its inherent vibrations. These studies concentrated on the calculation of frequencies that correspond to backbone amide bands I, II, and III: vibrational modes that are sensitive to changes in the secondary structure of peptides and proteins. To that end, Fourier transforms were applied to thermal fluctuations in C-N-H angles, C-N bond lengths, and C═O bond lengths of each amide group. In addition, instantaneous all-atom normal mode analysis was applied to monitor and detect the characteristic amide bands of each amide group within LM7-2 during the MD simulation. Computational vibrational results indicate that shapes and frequencies of amide bands II and especially III were altered only for amide groups near the hinge. These methods provide high-resolution vibrational information that can complement spectroscopic vibrational studies. They assist in interpreting spectra of similar systems and suggest a marker for the presence of the helix-hinge-helix motif. Moreover, radial distribution functions indicated an increase in the probability of hydrogen bonding between water and a hydrogen atom connected to nitrogen (HN) in such a hinge. The probability of intramolecular hydrogen bond formation between HN and an amide group oxygen atom within LM7-2 was lower around the hinge. No correlation has been found between the presence of a hinge and hydrogen bonds between amide group oxygen atoms and the hydrogen atoms of water molecules. This result suggests a mechanism for hinge formation wherein hydrogen bonds to oxygen atoms of water replace intramolecular hydrogen bonds as the peptide backbone folds.

Planarity and out-of-plane vibrational modes of tryptophan and tyrosine in biomolecular modeling.

Tryptophan and tyrosine are amino acids that play significant roles in the folding processes of proteins at water-membrane interfaces because of their amphipathic heteroaromatic rings. Employing appropriate heteroaromatic molecular structures is essential for obtaining accurate dynamics and predictive capabilities in molecular simulations of these amino acids. In this study, molecular dynamics simulations that applied the most recent version of the CHARMM36 force field were conducted on aqueous solutions of tryptophan and of tyrosine. Geometric analysis and dynamics quantified how aromatic rings deviated from planar structures and exhibited out-of-plane fluctuations. Radial distribution functions showed possible biological significance because the extent of ring planarity slightly affected local water concentrations near aromatic rings. Instantaneous all-atom normal mode analysis (NMA) and Fourier transformation of time autocorrelation functions of out-of-plane displacements were applied to study out-of-plane vibrations of atoms in these rings. The NMA started with minimum energy configurations and then averaged over fluctuations in aqueous solution. The frequencies and frequency patterns that were obtained for tryptophan and tyrosine with CHARMM36 differed from literature reports of Raman spectra, infrared spectra, and frequencies calculated using quantum mechanics, with some out-of-plane modes found at higher frequencies. Effects of imposing improper torsion potentials and changing torsion angle force constants were investigated for all atoms in the rings of tryptophan and tyrosine. Results show that these coarse force field variations only affect planarity and out-of-plane vibrations of atoms within the rings, and not other vibrations. Although increasing improper torsion force constants reduced deviations from aromatic ring planarity significantly, it increased out-of-plane mode frequencies. Reducing torsion angle force constants (with and without improper torsions) shifted modes to lower frequencies. A combination of decreasing most torsion angle force constants for ring atoms in both amino acids and including improper torsion forces attained frequencies and frequency patterns for out-of-plane normal modes that were more similar to the literature spectra. These force field variations decreased the extents of out-of-plane vibrations within the heteroaromatic rings of tryptophan, especially around the nitrogen atom in the ring, but not within the heteroaromatic ring of tyrosine. Conclusions were unaffected by the peptide endgroup, water, or simulation ensemble.

The unique N-terminal insert in the ribosomal protein, phosphoprotein P0, of Tetrahymena thermophila: Bioinformatic evidence for an interaction with 26S rRNA.

Phosphoprotein P0 (P0) is part of the stalk complex of the eukaryotic large ribosomal subunit necessary for recruiting elongation factors. While the P0 sequence is highly conserved, our group noted a 15-16 residue insert exclusive to the P0s of ciliated protists, including Tetrahymena thermophila. We hypothesized that this insert may have a function unique in ciliated protists, such as stalk regulation via phosphorylation of the insert. Almost no mention of this insert exists in the literature, and although the T. thermophila ribosome has been crystallized, there is limited structural data for Tetrahymena's P0 (TtP0) and its insert. To investigate the structure and function of the TtP0 insert, we performed in silico analyses. The TtP0 sequence was scanned with phosphorylation site prediction tools to detect the likelihood of phosphorylation in the insert. TtP0's sequence was also used to produce a homology model of the N-terminal domain of TtP0, including the insert. When the insert was modeled in the context of the 26S rRNA, it associated with a region identified as expansion segment 7B (ES7B), suggesting a potential functional interaction between ES7B and the insert in T. thermophila. We were not able to obtain sufficient data to determine whether a similar relationship exists in other ciliated protists. This study lays the groundwork for future experimental studies to verify the presence of TtP0 insert/ES7 interactions in Tetrahymena, and to explore their functional significance during protein synthesis.

Collagen gel formation in the presence of a carbon nanobrush.

Type I, bovine skin collagen was allowed to gel in the presence of various concentrations of a carbon nanotube material covered with a polystyrene/polyaniline copolymer, called a carbon nanobrush (CNB). The rate of collagen gelation was enhanced by the presence of the CNB in a dose dependent manner. The extent of collagen gelation was due to the concentration of collagen and not the amount of CNB. Collagen D-periodicity, and average fibril diameter were unchanged by the CNB material as seen in transmission electron micrographs. Gel tensile strength was reduced by the presence of the CNB in a dose related manner. The collagen-CNB mixture may have a role in the repair and reconstruction of wounds or degenerated connective tissue.

Homeoviscous response of Clostridium pasteurianum to butanol toxicity during glycerol fermentation.

Clostridium pasteurianum ATCC 6013 achieves high n-butanol production when glycerol is used as the sole carbon source. In this study, the homeoviscous membrane response of C. pasteurianum ATCC 6013 has been examined through n-butanol challenge experiments. Homeoviscous response is a critical aspect of n-butanol tolerance and has not been examined in detail for C. pasteurianum. Lipid membrane compositions were examined for glycerol fermentations with n-butanol production, and during cell growth in the absence of n-butanol production, using gas chromatography-mass spectrometry (GC-MS) and proton nuclear magnetic resonance ((1)H-NMR). Membrane stabilization due to homeoviscous response was further examined by surface pressure-area (π-A) analysis of membrane extract monolayers. C. pasteurianum was found to exert a homeoviscous response that was comprised of an increase lipid tail length and a decrease in the percentage of unsaturated fatty acids with increasing n-butanol challenge. This led to a more rigid or stable membrane that counteracted n-butanol fluidization. This is the first report on the changes in the membrane lipid composition during n-butanol production by C. pasteurianum ATCC 6013, which is a versatile microorganism that has the potential to be engineered as an industrial n-butanol producer using crude glycerol.