MicroRNAs inhibit the translation of target mRNAs on the endoplasmic reticulum in Arabidopsis.

Translation inhibition is a major but poorly understood mode of action of microRNAs (miRNAs) in plants and animals. In particular, the subcellular location where this process takes place is unknown. Here, we show that the translation inhibition, but not the mRNA cleavage activity, of Arabidopsis miRNAs requires ALTERED MERISTEM PROGRAM1 (AMP1). AMP1 encodes an integral membrane protein associated with endoplasmic reticulum (ER) and ARGONAUTE1, the miRNA effector and a peripheral ER membrane protein. Large differences in polysome association of miRNA target RNAs are found between wild-type and the amp1 mutant for membrane-bound, but not total, polysomes. This, together with AMP1-independent recruitment of miRNA target transcripts to membrane fractions, shows that miRNAs inhibit the translation of target RNAs on the ER. This study demonstrates that translation inhibition is an important activity of plant miRNAs, reveals the subcellular location of this activity, and uncovers a previously unknown function of the ER.

A century journey of organelles research in the plant endomembrane system.

We are entering an exciting century in the study of the plant organelles in the endomembrane system. Over the past century, especially within the past 50 years, tremendous advancements have been made in the complex plant cell to generate a much clearer and informative picture of plant organelles, including the molecular/morphological features, dynamic/spatial behavior, and physiological functions. Importantly, all these discoveries and achievements in the identification and characterization of organelles in the endomembrane system would not have been possible without: (1) the innovations and timely applications of various state-of-art cell biology tools and technologies for organelle biology research; (2) the continuous efforts in developing and characterizing new organelle markers by the plant biology community; and (3) the landmark studies on the identification and characterization of the elusive organelles. While molecular aspects and results for individual organelles have been extensively reviewed, the development of the techniques for organelle research in plant cell biology is less appreciated. As one of the ASPB Centennial Reviews on "organelle biology," here we aim to take a journey across a century of organelle biology research in plants by highlighting the important tools (or landmark technologies) and key scientists that contributed to visualize organelles. We then highlight the landmark studies leading to the identification and characterization of individual organelles in the plant endomembrane systems.

Autophagosome biogenesis and organelle homeostasis in plant cells.

Autophagy is one of the major highly inducible degradation processes in response to plant developmental and environmental signals. In response to different stimuli, cellular materials, including proteins and organelles, can be sequestered into a double membrane autophagosome structure either selectively or non-selectively. The formation of an autophagosome as well as its delivery into the vacuole involves complex and dynamic membrane processes. The identification and characterization of the conserved autophagy-related (ATG) proteins and their related regulators have greatly advanced our understanding of the molecular mechanism underlying autophagosome biogenesis and function in plant cells. Autophagosome biogenesis is tightly regulated by the coordination of multiple ATG and non-ATG proteins, and selective cargo recruitment. This review updates our current knowledge of autophagosome biogenesis, with special emphasis on the core molecular machinery that drives autophagosome formation, and autophagosome-organelle interactions under abiotic stress conditions.

Protein degrons and degradation: Exploring substrate recognition and pathway selection in plants.

Proteome composition is dynamic and influenced by many internal and external cues, including developmental signals, light availability, or environmental stresses. Protein degradation, in synergy with protein biosynthesis, allows cells to respond to various stimuli and adapt by reshaping the proteome. Protein degradation mediates the final and irreversible disassembly of proteins, which is important for protein quality control and to eliminate misfolded or damaged proteins, as well as entire organelles. Consequently, it contributes to cell resilience by buffering against protein or organellar damage caused by stresses. Moreover, protein degradation plays important roles in cell signaling, as well as transcriptional and translational events. The intricate task of recognizing specific proteins for degradation is achieved by specialized systems that are tailored to the substrate's physicochemical properties and subcellular localization. These systems recognize diverse substrate cues collectively referred to as "degrons", which can assume a range of structural configurations. They are molecular surfaces recognized by E3 ligases of the ubiquitin-proteasome system, but can also be considered as general features recognized by other degradation systems, including autophagy or even organellar proteases. Here we provide an overview of the newest developments in the field, delving into the intricate processes of protein recognition and elucidating the pathways through which they are recruited for degradation.

Subscaling of a cytosolic RNA binding protein governs cell size homeostasis in the multiple fission alga Chlamydomonas.

Coordination of growth and division in eukaryotic cells is essential for populations of proliferating cells to maintain size homeostasis, but the underlying mechanisms that govern cell size have only been investigated in a few taxa. The green alga Chlamydomonas reinhardtii (Chlamydomonas) proliferates using a multiple fission cell cycle that involves a long G1 phase followed by a rapid series of successive S and M phases (S/M) that produces 2n daughter cells. Two control points show cell-size dependence: the Commitment control point in mid-G1 phase requires the attainment of a minimum size to enable at least one mitotic division during S/M, and the S/M control point where mother cell size governs cell division number (n), ensuring that daughter distributions are uniform. tny1 mutants pass Commitment at a smaller size than wild type and undergo extra divisions during S/M phase to produce small daughters, indicating that TNY1 functions to inhibit size-dependent cell cycle progression. TNY1 encodes a cytosolic hnRNP A-related RNA binding protein and is produced once per cell cycle during S/M phase where it is apportioned to daughter cells, and then remains at constant absolute abundance as cells grow, a property known as subscaling. Altering the dosage of TNY1 in heterozygous diploids or through mis-expression increased Commitment cell size and daughter cell size, indicating that TNY1 is a limiting factor for both size control points. Epistasis placed TNY1 function upstream of the retinoblastoma tumor suppressor complex (RBC) and one of its regulators, Cyclin-Dependent Kinase G1 (CDKG1). Moreover, CDKG1 protein and mRNA were found to over-accumulate in tny1 cells suggesting that CDKG1 may be a direct target of repression by TNY1. Our data expand the potential roles of subscaling proteins outside the nucleus and imply a control mechanism that ties TNY1 accumulation to pre-division mother cell size.

A plant-unique protein BLISTER coordinates with core retromer to modulate endosomal sorting of plasma membrane and vacuolar proteins.

The retromer is a heteromeric protein complex that localizes to endosomal membranes and drives the formation of endosomal tubules that recycle membrane protein cargoes. In plants, the retromer plays essential and canonical functions in regulating the transport of vacuolar storage proteins and the recycle of endocytosed plasma membrane proteins (PM); however, the mechanisms underlying the regulation of assembly, protein stability, and membrane recruitment of the plant retromer complex remain to be elucidated. In this study, we identify a plant-unique endosomal regulator termed BLISTER (BLI), which colocalizes and associates with the retromer complex by interacting with the retromer core subunits VPS35 and VPS29. Depletion of BLI perturbs the assembly and membrane recruitment of the retromer core VPS26-VPS35-VPS29 trimer. Consequently, depletion of BLI disrupts retromer-regulated endosomal trafficking function, including transport of soluble vacuolar proteins and recycling of endocytosed PIN-FORMED (PIN) proteins from the endosomes back to the PM. Moreover, genetic analysis in Arabidopsis thaliana mutants reveals BLI and core retromer interact genetically in the regulation of endosomal trafficking. Taken together, we identified BLI as a plant-specific endosomal regulator, which functions in retromer pathway to modulate the recycling of endocytosed PM proteins and the trafficking of soluble vacuolar cargoes.

Arabidopsis AUTOPHAGY-RELATED2 is essential for ATG18a and ATG9 trafficking during autophagosome closure.

As a fundamental metabolic pathway, autophagy plays important roles in plant growth and development, particularly under stress conditions. A set of autophagy-related (ATG) proteins is recruited for the formation of a double-membrane autophagosome. Among them, the essential roles of ATG2, ATG18, and ATG9 have been well established in plant autophagy via genetic analysis; however, the underlying molecular mechanism for ATG2 in plant autophagosome formation remains poorly understood. In this study, we focused on the specific role of ATG2 in the trafficking of ATG18a and ATG9 during autophagy in Arabidopsis (Arabidopsis thaliana). Under normal conditions, YFP-ATG18a proteins are partially localized on late endosomes and translocated to ATG8e-labeled autophagosomes upon autophagic induction. Real-time imaging analysis revealed sequential recruitment of ATG18a on the phagophore membrane, showing that ATG18a specifically decorated the closing edges and finally disassociated from the completed autophagosome. However, in the absence of ATG2, most of the YFP-ATG18a proteins are arrested on autophagosomal membranes. Ultrastructural and 3D tomography analysis showed that unclosed autophagosome structures are accumulated in the atg2 mutant, displaying direct connections with the endoplasmic reticulum membrane and vesicular structures. Dynamic analysis of ATG9 vesicles suggested that ATG2 depletion also affects the association between ATG9 vesicles and the autophagosomal membrane. Furthermore, using interaction and recruitment analysis, we mapped the interaction relationship between ATG2 and ATG18a, implying a possible role of ATG18a in recruiting ATG2 and ATG9 to the membrane. Our findings unveil a specific role of ATG2 in coordinating ATG18a and ATG9 trafficking to mediate autophagosome closure in Arabidopsis.

Transcriptional and Epigenetic Regulation of Autophagy in Plants.

Autophagy, a highly conserved quality control mechanism, is essential for maintaining cellular homeostasis and healthy growth of plants. Compared with extensive research in the cytoplasmic control of autophagy, studies regarding the nuclear events involved in the regulation of plant autophagy are just beginning to emerge. Accumulating evidence reveals a coordinated expression of plant autophagy genes in response to diverse developmental states and growth conditions. Here, we summarize recent progress in the identification of tightly controlled transcription factors and histone marks associated with the autophagic process in plants, and propose several modules, consisting of transcription regulators and epigenetic modifiers, as important nuclear players that could contribute to both short-term and long-term controls of plant autophagy at the transcriptional and post-transcriptional levels.

Anlotinib Combined With Toripalimab as First-Line Therapy for Unresectable Hepatocellular Carcinoma: A Prospective, Multicenter, Phase II Study.

BACKGROUND: For patients with unresectable hepatocellular carcinoma (HCC), the first-line therapeutic options are still relatively limited, and treatment outcomes remain poor. We aimed to assess the efficacy and safety of anlotinib combined with toripalimab as first-line therapy for unresectable HCC. METHODS: In this single-arm, multicenter, phase II study (ALTER-H-003), patients with advanced HCC without previous systemic anticancer therapy were recruited. Eligible patients were given anlotinib (12 mg on days 1-14) combined with toripalimab (240 mg on day 1) in a 3-week cycle. The primary endpoint was the objective response rate (ORR) by immune-related Response Evaluation Criteria in Solid Tumours (irRECIST)/RECIST v1.1 and modified RECIST (mRECIST). Secondary endpoints included disease control rate (DCR), duration of response (DoR), progression-free survival (PFS), overall survival (OS), and safety. RESULTS: Between January 2020 and Jul 2021, 31 eligible patients were treated and included in the full analysis set. At data cutoff (January 10, 2023), the ORR was 29.0% (95% CI: 12.1%-46.0%) by irRECIST/RECIST v1.1, and 32.3% (95% CI: 14.8%-49.7%) by mRECIST criteria, respectively. Confirmed DCR and median DoR by irRECIST/RECIST v1.1 and mRECIST criteria were 77.4 % (95% CI: 61.8%-93.0%) and not reached (range: 3.0-22.5+ months), respectively. Median PFS was 11.0 months (95% CI: 3.4-18.5 months) and median OS was 18.2 months (95% CI: 15.8-20.5 months). Of the 31 patients assessed for adverse events (AEs), the most common grade ≥ 3 treatment-related AEs were hand-foot syndrome (9.7%, 3/31), hypertension (9.7%, 3/31), arthralgia (9.7%, 3/31), abnormal liver function (6.5%, 2/31), and decreased neutrophil counts (6.5%, 2/31). CONCLUSIONS: Anlotinib combined with toripalimab showed promising efficacy and manageable safety in Chinese patients with unresectable HCC in the first-line setting. This combination therapy may offer a potential new therapeutic approach for patients with unresectable HCC.

Recent advances in plant endomembrane research and new microscopical techniques.

The endomembrane system consists of various membrane-bound organelles including the endoplasmic reticulum (ER), Golgi apparatus, trans-Golgi network (TGN), endosomes, and the lysosome/vacuole. Membrane trafficking between distinct compartments is mainly achieved by vesicular transport. As the endomembrane compartments and the machineries regulating the membrane trafficking are largely conserved across all eukaryotes, our current knowledge on organelle biogenesis and endomembrane trafficking in plants has mainly been shaped by corresponding studies in mammals and yeast. However, unique perspectives have emerged from plant cell biology research through the characterization of plant-specific regulators as well as the development and application of the state-of-the-art microscopical techniques. In this review, we summarize our current knowledge on the plant endomembrane system, with a focus on several distinct pathways: ER-to-Golgi transport, protein sorting at the TGN, endosomal sorting on multivesicular bodies, vacuolar trafficking/vacuole biogenesis, and the autophagy pathway. We also give an update on advanced imaging techniques for the plant cell biology research.

Indole Diterpene Derivatives from the Aspergillus flavus GZWMJZ-288, an Endophytic Fungus from Garcinia multiflora.

A new indole diterpene, 26-dihydroxyaflavininyl acetate (1), along with five known analogs (2-6) were isolated from the liquid fermentation of Aspergillus flavus GZWMJZ-288, an endophyte from Garcinia multiflora. The structures of these compounds were identified through NMR, MS, chemical reaction, and X-ray diffraction experiments. Enzyme inhibition activity screening found that compounds 1, 4, and 6 have a good binding affinity with NPC1L1, among which compound 6 exhibited a stronger binding ability than ezetimibe at a concentration of 10 µM. Moreover, compound 5 showed inhibitory activity against α-glucosidase with an IC50 value of 29.22 ± 0.83 µM, which is 13 times stronger than that of acarbose. The results suggest that these aflavinine analogs may serve as lead compounds for the development of drugs targeting NPC1L1 and α-glucosidase. The binding modes of the bioactive compounds with NPC1L1 and α-glucosidase were also performed through in silico docking studies.

MLKs kinases phosphorylate the ESCRT component FREE1 to suppress abscisic acid sensitivity of seedling establishment.

The FYVE domain protein required for endosomal sorting 1 (FREE1), which was previously identified as a plant-specific component of the endosomal sorting complex required for transport machinery, plays an essential role in endosomal trafficking. Moreover, FREE1 also functions as an important negative regulator in abscisic acid (ABA) signalling. Multiple phosphorylations and ubiquitination sites have been identified in FREE1, hence unveiling the factors involved in posttranslational regulation of FREE1 is critical for comprehensively understanding FREE1-related regulatory networks during plant growth. Here, we demonstrate that plant-specific casein kinase I members MUT9-like kinases 1-4 (MLKs 1-4)/Arabidopsis EL1-like 1-4 interact with and phosphorylate FREE1 at serine residue S582, thereby modulating the nuclear accumulation of FREE1. Consequently, mutation of S582 to non-phosphorylable residue results in reduced nuclear localization of FREE1 and enhanced ABA response. In addition, mlk123 and mlk134 triple mutants accumulate less FREE1 in the nucleus and display hypersensitive responses to ABA treatment, whereas overexpression of the nuclear-localized FREE1 can restore the ABA sensitivity of seedling establishment in mlks triple mutants. Collectively, our study demonstrates a previously unidentified function of MLKs in attenuating ABA signalling in the nucleus by regulating the phosphorylation and nuclear accumulation of FREE1.

AtNBR1 Is a Selective Autophagic Receptor for AtExo70E2 in Arabidopsis.

Selective autophagy is a subcellular process whereby cytoplasmic materials are selectively sequestered into autophagosomes for subsequent delivery to the vacuole for degradation and recycling. Arabidopsis (Arabidopsis thaliana) NBR1 (next to BRCA1 gene 1 protein; AtNBR1) has been proposed to function as a selective autophagy receptor in plants, whereby AtNBR1 anchors the ubiquitinated targets to autophagosomes for degradation. However, the specific cargos of AtNBR1 remain elusive. We previously showed that Arabidopsis exocyst subunit EXO70 family protein E2 (AtExo70E2), a marker for exocyst-positive organelle (EXPO), colocalized with the autophagosome marker Arabidopsis autophagy-related protein8 (AtATG8) and was delivered to the vacuole for degradation upon autophagic induction. Here, through multiple analyses, we demonstrate that AtNBR1 is a selective receptor for AtExo70E2 during autophagy in Arabidopsis. First, two novel loss-of-function nbr1 CRISPR mutants (nbr1-c1 and nbr1-c2) showed an early-senescence phenotype under short-day growth conditions. Second, during autophagic induction, the vacuolar delivery of AtExo70E2 or EXPO was significantly reduced in nbr1 mutants compared to wild-type plants. Third, biochemical and recruitment assays demonstrated that AtNBR1 specifically interacted and recruited AtExo70E2 or its EXPO to AtATG8-positive autophagosomes in a ubiquitin-associated (UBA)-independent manner during autophagy. Taken together, our data indicate that AtNBR1 functions as a selective receptor in mediating vacuolar delivery of AtExo70E2 or EXPO in a UBA-independent manner in plant autophagy.

Plant Mitophagy in Comparison to Mammals: What Is Still Missing?

Mitochondrial homeostasis refers to the balance of mitochondrial number and quality in a cell. It is maintained by mitochondrial biogenesis, mitochondrial fusion/fission, and the clearance of unwanted/damaged mitochondria. Mitophagy represents a selective form of autophagy by sequestration of the potentially harmful mitochondrial materials into a double-membrane autophagosome, thus preventing the release of death inducers, which can trigger programmed cell death (PCD). Recent advances have also unveiled a close interconnection between mitophagy and mitochondrial dynamics, as well as PCD in both mammalian and plant cells. In this review, we will summarize and discuss recent findings on the interplay between mitophagy and mitochondrial dynamics, with a focus on the molecular evidence for mitophagy crosstalk with mitochondrial dynamics and PCD.

TRAF Family Proteins Regulate Autophagy Dynamics by Modulating AUTOPHAGY PROTEIN6 Stability in Arabidopsis.

Eukaryotic cells use autophagy to recycle cellular components. During autophagy, autophagosomes deliver cytoplasmic contents to the vacuole or lysosome for breakdown. Mammalian cells regulate the dynamics of autophagy via ubiquitin-mediated proteolysis of autophagy proteins. Here, we show that the Arabidopsis thaliana Tumor necrosis factor Receptor-Associated Factor (TRAF) family proteins TRAF1a and TRAF1b (previously named MUSE14 and MUSE13, respectively) help regulate autophagy via ubiquitination. Upon starvation, cytoplasmic TRAF1a and TRAF1b translocated to autophagosomes. Knockout traf1a/b lines showed reduced tolerance to nutrient deficiency, increased salicylic acid and reactive oxygen species levels, and constitutive cell death in rosettes, resembling the phenotypes of autophagy-defective mutants. Starvation-activated autophagosome accumulation decreased in traf1a/b root cells, indicating that TRAF1a and TRAF1b function redundantly in regulating autophagosome formation. TRAF1a and TRAF1b interacted in planta with ATG6 and the RING finger E3 ligases SINAT1, SINAT2, and SINAT6 (with a truncated RING-finger domain). SINAT1 and SINAT2 require the presence of TRAF1a and TRAF1b to ubiquitinate and destabilize AUTOPHAGY PROTEIN6 (ATG6) in vivo. Conversely, starvation-induced SINAT6 reduced SINAT1- and SINAT2-mediated ubiquitination and degradation of ATG6. Consistently, SINAT1/SINAT2 and SINAT6 knockout mutants exhibited increased tolerance and sensitivity, respectively, to nutrient starvation. Therefore, TRAF1a and TRAF1b function as molecular adaptors that help regulate autophagy by modulating ATG6 stability in Arabidopsis.

SH3 Domain-Containing Protein 2 Plays a Crucial Role at the Step of Membrane Tubulation during Cell Plate Formation.

During cytokinesis in plants, trans-Golgi network-derived vesicles accumulate at the center of dividing cells and undergo various structural changes to give rise to the planar cell plate. However, how this conversion occurs at the molecular level remains elusive. In this study, we report that SH3 Domain-Containing Protein 2 (SH3P2) in Arabidopsis thaliana plays a crucial role in converting vesicles to the planar cell plate. SH3P2 RNAi plants showed cytokinesis-defective phenotypes and produced aggregations of vesicles at the leading edge of the cell plate. SH3P2 localized to the leading edge of the cell plate, particularly the constricted or curved regions of the cell plate. The BAR domain of SH3P2 induced tubulation of vesicles. SH3P2 formed a complex with dynamin-related protein 1A (DRP1A) and affected DRP1A accumulation to the cell plate. Based on these results, we propose that SH3P2 functions together with DRP1A to convert the fused vesicles to tubular structures during cytokinesis.

ATG9 regulates autophagosome progression from the endoplasmic reticulum in Arabidopsis.

Autophagy is a conserved pathway for bulk degradation of cytoplasmic material by a double-membrane structure named the autophagosome. The initiation of autophagosome formation requires the recruitment of autophagy-related protein 9 (ATG9) vesicles to the preautophagosomal structure. However, the functional relationship between ATG9 vesicles and the phagophore is controversial in different systems, and the molecular function of ATG9 remains unknown in plants. Here, we demonstrate that ATG9 is essential for endoplasmic reticulum (ER)-derived autophagosome formation in plants. Through a combination of genetic, in vivo imaging and electron tomography approaches, we show that Arabidopsis ATG9 deficiency leads to a drastic accumulation of autophagosome-related tubular structures in direct membrane continuity with the ER upon autophagic induction. Dynamic analyses demonstrate a transient membrane association between ATG9 vesicles and the autophagosomal membrane during autophagy. Furthermore, trafficking of ATG18a is compromised in atg9 mutants during autophagy by forming extended tubules in a phosphatidylinositol 3-phosphate-dependent manner. Taken together, this study provides evidence for a pivotal role of ATG9 in regulating autophagosome progression from the ER membrane in Arabidopsis.

Dual roles of an Arabidopsis ESCRT component FREE1 in regulating vacuolar protein transport and autophagic degradation.

Protein turnover can be achieved via the lysosome/vacuole and the autophagic degradation pathways. Evidence has accumulated revealing that efficient autophagic degradation requires functional endosomal sorting complex required for transport (ESCRT) machinery. However, the interplay between the ESCRT machinery and the autophagy regulator remains unclear. Here, we show that FYVE domain protein required for endosomal sorting 1 (FREE1), a recently identified plant-specific ESCRT component essential for multivesicular body (MVB) biogenesis and plant growth, plays roles both in vacuolar protein transport and autophagic degradation. FREE1 also regulates vacuole biogenesis in both seeds and vegetative cells of Arabidopsis. Additionally, FREE1 interacts directly with a unique plant autophagy regulator SH3 domain-containing protein2 and associates with the PI3K complex, to regulate the autophagic degradation in plants. Thus, FREE1 plays multiple functional roles in vacuolar protein trafficking and organelle biogenesis as well as in autophagic degradation via a previously unidentified regulatory mechanism of cross-talk between the ESCRT machinery and autophagy process.

Modeling structural transitions from the periplasmic-open state of lactose permease and interpretations of spin label experiments.

Lactose permease of E. coli (LacY) is a secondary active transporter (SAT) that belongs to the major facilitator superfamily (MFS). Experimental structures of the cytoplasmic-open and more recently occluded-like structure have been determined, however, the crystal structure of LacY in the periplasmic-open state is still not available. The periplasmic-open LacY structure is important for understanding complete proton/sugar transport process of LacY as well as other similar SAT proteins. Previously, a structural model of periplasmic-open LacY has been obtained through a two-step hybrid implicit-explicit (IM-EX) simulation method (JMB404: 506). Molecular dynamics simulations are performed to further test the IM-EX model for the periplasmic-open LacY with ßß-(Galp)2 in a lipid membrane. The comparison of the calculated pore radii to the data of the crystal structure indicates that the IM-EX model of LacY remains periplasmic-open in E269-protonated states. The neighbor residue distance change based on Cα are very similar in simulation results, but they are significantly different in double electron-electron resonance (DEER) experimental data, which motivates us to perform the molecular dynamics dummy spin-label (MDDS) simulations to test the effect of spin labels (size and internal flexibility) on DEER spin label distance measurements. The MDDS simulation results show that the orientation and movement of the spin labels significantly affect the residue pair distance measurement. DEER data alone may not provide an accurate guide for predicting protein structures. MDDS simulations can be applied to analyze the distance distribution due to spin labels and also aid in proper interpretation of DEER experimental data.

An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations.

Previous MD simulations of six phosphocholine (PC) lipid bilayers demonstrated the accuracy of the CHARMM36 force field (C36FF) for PC bilayer simulation at varied temperatures (BBA-Biomembranes, 1838 (2014): 2520-2529). In this work, we further examine the accuracy of C36FF over a wide temperature range for a broader range of lipid types such as various head groups (phosphatidic acid (PA), PC, phosphoethanolamine (PE), phosphoglycerol (PG), and phosphoserine (PS)), and tails (saturated, mono-, mixed- and poly-unsaturated acyl chains with varied chain lengths). The structural properties (surface area per lipid (SA/lip), overall bilayer thickness, hydrophobic thickness, headgroup-to-headgroup thickness, deuterium order parameter (SCD), and spin-lattice relaxation time (T1)) obtained from simulations agree well with nearly all available experimental data. Our analyses indicate that PS lipids have the most inter-lipid hydrogen bonds, while PG lipids have the most intra-lipid hydrogen bonds, which play the main role in their low SA/lip in PS lipids and low thicknesses in PG lipids, respectively. PS, PE, and PA lipids have the largest contact clusters with on average 5-8 lipids per cluster, while PC and PG have clusters of 4 lipids based on a cutoff distance of 6.5Å. PS lipids have much slower lipid wobble (i.e., higher correlation time) than other head groups at a given temperature as the hydrogen bonded network significantly reduces a lipid's mobility, and the rate of lipid wobble increases dramatically as temperature increases. These in-depth analyses facilitate further understanding of lipid bilayers at the atomic level.