Major Conformation Research Articles

Characterization of a transitionally occupied state and thermal unfolding of domain 1.1 of factor of RNA polymerase from Bacillus subtilis.

factors are essential parts of bacterial RNA polymerase (RNAP) as they allow to recognize promotor sequences and initiate transcription. Domain 1.1 of vegetative factors occupies the primary channel of RNAP and also prevents binding of the factor to promoter DNA alone. Here, we show that domain 1.1 of Bacillus subtilis exists in more structurally distinct variants in dynamic equilibrium. The major conformation at room temperature is represented by a previously reported well-folded structure solved by nuclear magnetic resonance (NMR), but 4% of the protein molecules are present in a less thermodynamically favorable state. We show that this population increases with temperature and we predict its significant elevation at higher but still biologically relevant temperatures. We characterized the minor state of the domain 1.1 using specialized methods of NMR. We found that, in contrast to the major state, the detected minor state is partially unfolded. Its propensity to form secondary structure elements is especially decreased for the first and third helices, while the second helix and strand close to the C-terminus are more stable. We also analyzed thermal unfolding of the domain 1.1 and performed functional experiments with full length and its shortened version lacking domain 1.1 ( ). The results revealed that while full length increases transcription activity of RNAP with increasing temperature, transcription with remains constant. In summary, this study reveals conformational dynamics of domain 1.1 and provides a basis for studies of its interaction with RNAP and effects on transcription regulation.

Conformational switching and flexibility in cobalamin-dependent methionine synthase studied by small-angle X-ray scattering and cryoelectron microscopy

Cobalamin-dependent methionine synthase (MetH) catalyzes the synthesis of methionine from homocysteine and 5-methyltetrahydrofolate (CH3-H4folate) using the unique chemistry of its cofactor. In doing so, MetH links the cycling of S-adenosylmethionine with the folate cycle in one-carbon metabolism. Extensive biochemical and structural studies on Escherichia coliMetH have shown that this flexible, multidomain enzyme adopts two major conformations to prevent a futile cycle of methionine production and consumption. However, as MetH is highly dynamic as well as both a photosensitive and oxygen-sensitive metalloenzyme, it poses special challenges for structural studies, and existing structures have necessarily come from a "divide and conquer" approach. In this study, we investigate E. coli MetH and a thermophilic homolog from Thermus filiformis using small-angle X-ray scattering (SAXS), single-particle cryoelectron microscopy (cryo-EM), and extensive analysis of the AlphaFold2 database to present a structural description of the full-length MetH in its entirety. Using SAXS, we describe a common resting-state conformation shared by both active and inactive oxidation states of MetH and the roles of CH3-H4folate and flavodoxin in initiating turnover and reactivation. By combining SAXS with a 3.6-Å cryo-EM structure of the T. filiformis MetH, we show that the resting-state conformation consists of a stable arrangement of the catalytic domains that is linked to a highly mobile reactivation domain. Finally, by combining AlphaFold2-guided sequence analysis and our experimental findings, we propose a general model for functional switching in MetH.

Phosphorylation at Ser65 modulates ubiquitin conformational dynamics

Ubiquitin phosphorylation at Ser65 increases the population of a rare C-terminally retracted (CR) conformation. Transition between the Major and CR ubiquitin conformations is critical for promoting mitochondrial degradation. The mechanisms by which the Major and CR conformations of Ser65-phosphorylated (pSer65) ubiquitin interconvert, however, remain unresolved. Here, we perform all-atom molecular dynamics simulations using the string method with swarms of trajectories to calculate the lowest free-energy path between these two conformers. Our analysis reveals the existence of a Bent intermediate in which the C-terminal residues of the β5 strand shift to resemble the CR conformation, while pSer65 retains contacts resembling the Major conformation. This stable intermediate was reproduced in well-tempered metadynamics calculations but was less stable for a Gln2Ala mutant that disrupts contacts with pSer65. Lastly, dynamical network modeling reveals that the transition from the Major to CR conformations involves a decoupling of residues near pSer65 from the adjacent β1 strand.

Cyanobacterial compound Tolyporphine K as an inhibitor of Apo-PBP (penicillin-binding protein) in A. baumannii and its ADME assessment

Antibiotic-resistant Acinetobacter baumannii, is a common pathogen found in hospital settings and has become nosocomial due to its high infection-causing tendency amongst ICU patients. The present study explores the cyanocompoundswhich were capable to inhibit the Penicillin Binding Protein of A. baumannii through molecular docking, ADMET, and molecular dynamicssimulation strategy. A database having structural and origin details was generated for 85 bioactive compounds in MS Excel. The 3-D structures weredownloaded from the PubChem database and minimized. The receptor protein was minimized and validated for structure correctness. The database was screened against the penicillin-binding protein of A. baumannii through PyRx software. The top 5 compounds including the control molecule werefurther redocked to the receptor molecule through Autodock Vina software. The molecule pose having the highest affinity was further subjected to 100ns MD- simulation and simultaneously the in-vitro activity of the methanol extract and hexane extract was checked through agar well diffusion assay.Docking studies indicate Tolyporphine K to be a lead molecule which was further assessed through Molecular dynamics and MM/PBSA. The in-silicoresults suggested that the protein-ligand complex was found to be stable over the 100 ns trajectory with a binding free energy of -8.56 Kcalmol-1. Theligand did not induce any major structural conformation in the protein moiety and was largely stabilized by hydrophobic interactions. The bioactivityscore and ADME properties of the compounds were also calculated. The in-vitro agar well diffusion assay showed a moderate zone of inhibition of12.33mm. The results indicate that the compound Tolyporphin- K could be a potential inhibitor of penicillin-binding protein in A. baumannii. Yet furtherwork needs to be done to have a more concrete basis for the pathway of inhibition. Communicated by Ramaswamy H. Sarma

Stereochemical Analysis of Trifluoroacetamide Derivatives Based on Through-Space 1H–19F Spin–Spin Couplings

In this study, the conformational properties of tertiary trifluoroacetamides in dibenzoazepine (1a and 1b) and benzodiazepine (2a and 2b) derivatives, which exist as equilibrated E- and Z-amide conformers in solution, were investigated by 1H and 19F NMR spectroscopy. In all cases, one of the methylene protons neighboring the nitrogen atom of the minor conformer showed a finely split pattern due to coupling with the trifluoromethyl fluorine atoms, as confirmed by 19F-decoupling experiments. One-dimensional (1D) and two-dimensional (2D) 1H-19F heteronuclear Overhauser spectroscopy (HOESY) experiments were performed to confirm whether these couplings are attributable to through-bond spin-spin couplings (TBCs) or through-space spin-spin couplings (TSCs). HOESY cross-peaks between CF3 (19F) and one of the CH2-N protons of the minor conformers indicate that the two nuclei are spatially close to each other, thus establishing the stereochemistry of the major (E-) and minor (Z-) conformers. The E-amide preferences of the trifluoroacetamides are consistent with the results of density functional theory calculations and X-ray crystallographic analyses. Furthermore, the otherwise incomprehensible 1H NMR spectra were accurately assigned using the HOESY-determined TSCs. The 1H NMR assignments of the E- and Z-methyl signals of N,N-dimethyl trifluoroacetamide, the simplest tertiary trifluoroacetamide, were revised for the first time in half a century.

Role of Metal Cofactor in Enhanced Thermal Stability of Azurin.

Metal cofactors are critical centers for different biochemical processes of metalloproteins, and often, this metal coordination renders additional structural stability. In this study, we explore the additional stability conferred by the copper ion on azurin by analyzing both the apo and holo forms using temperature replica exchange molecular dynamics (REMD) data. We find a 14 K decrease in denaturation temperature for apo (406 K) azurin relative to that of holo (420 K), indicating a copper ion-induced additional thermal stability for holo azurin. The unfolding of apo azurin begins with the melting of α-helix and β-sheet V, similar to that of holo form. β-Sheets IV, VII, and VIII are comparatively more stable than other β-strands and melt at higher temperatures. Similar to holo azurin, the strong hydrophobic interactions among the apolar residues in the protein core is the key factor that renders high stability to apo protein as well. We construct free energy surfaces at different temperatures to capture the major conformations along the unfolding basins of the protein. Using contact maps from different basins we show the changes in the interaction between different residues along the unfolding pathway. Furthermore, we compare the Cα root-mean-square fluctuations (Cα-RMSF) and B-factor of all residues of apo and holo forms to understand the flexibility of different regions. The concerted displacement of α-helix and β-sheets V and VI from the protein core is another distinction we observe for apo compared to the holo form, where β-sheet VI was relatively stable.

Evaluation of potential effects and genetic parameters in conformational limb defects in Pura Raza Española horses

Morphological limbs defects are the most frequent defects in horses. The aim of this study was to establish the prevalence, associated effects and genetic parameters of hock and knee defects in Pura Raza Española horses, using two different approaches. In Approach_1, the hock defects analysed were closed, open, convergent and divergent and the knee defects studied were buck, calf, bench and knock. Defects were classified into 3 levels: 0 no defect, 1 slight defect and 2 serious defects. Approach_2, which used a linear scale, in a pair of opposing defects, divided into 5 levels, from −2 to 2, where 0 corresponded to the absence of defects: Lateral view (closed/open) and rear view (convergent/divergent) hock defect, lateral view (buck/calf) and frontal view (bench/knock) knee defect. A total of 43,358 horses, with an average age of 5.07 years, were evaluated, with a prevalence of horses affected ranging from 3.31% (bench) to 74.12% (convergent). Genetic parameters were estimated using a Bayesian procedure with the BLUPF90 software. Heritability in Approach_1, ranged from 0.25 (bench) to 0.42 (divergent) and in Approach_2, from 0.18 (bench/knock) to 0.24 (convergent/divergent). The opposing defects may be related to different genes, and it is therefore better to study them as separate defects and not on the same linear scale. The highest positive genetic correlation was between calf vs knock (0.70). Our results imply that selection against limb defects is possible and would allow us to reduce the genetic risk of the horses’ offspring suffering from them. Highlights There high prevalence of hock and knee defects in PRE horses make the animals unsuitable for sport and cause a significant economic loss to the breeding herds. The evaluation of limb defects independently instead of the traditional linear scale, where the opposite defects are located at the extremes, has good results for its genetic evaluation. Using the hock and knee defects as a selection criteria in the PRE' s breeding program could contribute to the eradication of major conformation defects which make the animals unsuitable for sport.

Benchmarking the Drude Polarizable Force Field Using the r(GACC) Tetranucleotide.

Polarizable force fields, in particular, the Drude polarizable force field (Drude FF), may hold the key to more accurately modeling biomolecules with molecular dynamics simulations by explicitly accounting for atomic polarizability. Previous work has shown promising results in simulating duplex nucleic acids and protein structures with excellent agreement with experimental values. However, benchmarking the Drude polarizable force field with highly flexible, single-stranded structures has yet to be achieved. In this work, the r(GACC) tetranucleotide is simulated over a multimicrosecond time scale, starting with various different initial conformations. Despite the starting conformation, including starting from the expected dominant A-form major conformation, the experimental structural distribution is not matched. In fact, the major NMR conformation is never resampled. Instead, the r(GACC) tetranucleotide becomes stabilized in anomalous structures that are inconsistent with the NMR data and that favor base-pairing and electrostatic interactions over base stacking. These structures are maintained for lengthy time scales (>1 μs) themselves, suggesting a misbalance of forces in the Drude polarizable force field itself. This model system is suggestive of the fact that currently the Drude polarizable force field does not appear to produce the sensitive balance of forces required to accurately model other single-stranded or noncanonical RNA structures.

Distinct structure and gating mechanism in diverse NMDA receptors with GluN2C and GluN2D subunits.

N-methyl-D-aspartate (NMDA) receptors are heterotetramers comprising two GluN1 and two alternate GluN2 (N2A-N2D) subunits. Here we report full-length cryo-EM structures of the human N1-N2D di-heterotetramer (di-receptor), rat N1-N2C di-receptor and N1-N2A-N2C tri-heterotetramer (tri-receptor) at a best resolution of 3.0 Å. The bilobate N-terminal domain (NTD) in N2D intrinsically adopts a closed conformation, leading to a compact NTD tetramer in the N1-N2D receptor. Additionally, crosslinking the ligand-binding domain (LBD) of two N1 protomers significantly elevated the channel open probability (Po) in N1-N2D di-receptors. Surprisingly, the N1-N2C di-receptor adopted both symmetric (minor) and asymmetric (major) conformations, the latter further locked by an allosteric potentiator, PYD-106, binding to a pocket between the NTD and LBD in only one N2C protomer. Finally, the N2A and N2C subunits in the N1-N2A-N2C tri-receptor display a conformation close to one protomer in the N1-N2A and N1-N2C di-receptors, respectively. These findings provide a comprehensive structural understanding of diverse function in major NMDA receptor subtypes.

De Novo Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Chaenomeles speciosa (Sweet) Nakai Revealed the Existence of Two Structural Isomers.

As a valuable Chinese traditional medicinal species, Chaenomeles speciosa (Sweet) Nakai (C. speciosa) is a natural resource with significant economic and ornamental value. However, its genetic information is not well understood. In this study, the complete mitochondrial genome of C. speciosa was assembled and characterized to explore the repeat sequences, recombination events, rearrangements, and IGT, to predict RNA editing sites, and to clarify the phylogenetic and evolutionary relationship. The C. speciosa mitochondrial genome was found to have two circular chromosomes as its major conformation, with a total length of 436,464 bp and 45.2% GC content. The mitochondrial genome contained 54 genes, including 33 unique protein-coding genes, 18 tRNAs, and 3 rRNA genes. Seven pairs of repeat sequences involving recombination events were analyzed. Both the repeat pairs, R1 and R2, played significant roles in mediating the major and minor conformations. In total, 18 MTPTs were identified, 6 of which were complete tRNA genes. There were 454 RNA editing sites in the 33 protein-coding sequences predicted by the PREPACT3 program. A phylogenetic analysis based on 22 species of mitochondrial genomes was constructed and indicated highly conserved PCG sequences. Synteny analyses showed extensive genomic rearrangements in the mitochondrial genome of C. speciosa and closely related species. This work is the first to report the C. speciosa mitochondrial genome, which is of great significance for conducting additional genetic studies on this organism.

Conversion of Helix 1 into a Loop in Prion Protein Misfolding.

Cellular prion protein PrPC consists of three α-helices, one β-sheet, and an unstructured N-terminal domain. Misfolding of this protein into the scrapie form (PrPSc) increases dramatically the β-sheet content. H1 is the most stable helix on PrPC and contains an unusual number of hydrophilic amino acids. Its fate in PrPSc is not clear. We performed replica exchange molecular dynamics simulations on H1 alone, H1 together with an N-terminally flanking H1B1 loop and H1 in complex with other hydrophilic regions of the prion protein. In the presence of the H99SQWNKPSKPKTNMK113 sequence, H1 is almost completely converted to a loop structure stabilized by a network of salt bridges. On the other hand, H1 retains its helical structure alone or together with the other sequences considered in this study. We carried out an additional simulation by restraining the distance between the two ends of H1, mimicking a possible geometric restriction by the rest of the protein. Even though the loop was the major conformation, a significant amount of helical structure was also observed. This suggests that the interaction with H99SQWNKPSKPKTNMK113 is necessary for complete helix-to-loop conversion.

Light Absorption by Cinnamaldehyde Constituents of Biomass Burning Organic Aerosol Modeled Using Time-Dependent Density Functional Theory

Organic aerosol emitted from biomass burning absorbs visible radiation. However, the impact of this light absorption on the overall climate effects of atmospheric aerosol is not well known, partly due to variability in particle composition and absorptivity. Cinnamaldehydes, which consist of an aromatic ring with an unsaturated aldehyde substituent, are an important class of chromophores in light-absorbing organic aerosol, or brown carbon. Here, light absorption by a homologous series of three cinnamaldehydes─coumaraldehyde, coniferaldehyde, and sinapaldehyde─is modeled with time-dependent density functional theory (TD-DFT) calculations, in the gas and aqueous phases. Based on a survey of hydration and acid dissociation equilibria, the neutral aldehyde is expected to be the predominant form of each species in the atmospheric aqueous phase. These species have complicated conformational landscapes compared to many other brown carbon constituents, like rigid polycyclic aromatic hydrocarbons. For coumaraldehyde, coniferaldehyde, and sinapaldehyde, a total of 8, 26, and 18 conformers were located, respectively. For each species, most of the total population is accounted for by the four most-populated conformers. The relative contributions of the conformers to the total light absorption of the respective species are dictated more by differences in the relative free energies than by differences in the molar absorption coefficients. As the functionalization increases, the absorption is red-shifted. The peaks predicted in water agree well with experimental spectra of coniferaldehyde and sinapaldehyde. No conformers have vertical transitions in the visible spectral range, so absorption above 380 nm is due to the shoulders of transitions of major conformers at ultraviolet wavelengths. These results demonstrate the importance of exploring potential energy landscapes, determining conformer stability and absorptivity, to predict the light absorption of chromophores in brown carbon.

Systematic mapping of the conformational landscape and dynamism of soluble fibrinogen

BackgroundFibrinogen is a soluble, multisubunit, and multidomain dimeric protein, which, upon its proteolytic cleavage by thrombin, is converted to insoluble fibrin, initiating polymerization that substantially contributes to clot growth. Fibrinogen contains numerous, transiently accessible “cryptic” epitopes for hemostatic and immunologic proteins, suggesting that fibrinogen exhibits conformational flexibility, which may play functional roles in its temporal and spatial interactions. Hitherto, there have been limited integrative approaches characterizing the solution structure and internal flexibility of fibrinogen. MethodsHere, utilizing a multipronged, biophysical approach involving 2 solution-based techniques, temperature-dependent hydrogen-deuterium exchange mass spectrometry and small angle X-ray scattering, corroborated by negative stain electron microscopy, we present a holistic, conformationally dynamic model of human fibrinogen in solution. ResultsOur data reveal 4 major and distinct conformations of fibrinogen accommodated by a high degree of internal protein flexibility along its central scaffold. We propose that the fibrinogen structure in the solution consists of a complex, conformational landscape with multiple local minima. This is further supported by the location of numerous point mutations that are linked to dysfibrinogenemia and posttranslational modifications, residing near the identified fibrinogen flexions. ConclusionThis work provides a molecular basis for the structural “dynamism” of fibrinogen that is expected to influence the broad swath of its functionally diverse macromolecular interactions and fine-tune the structural and mechanical properties of blood clots.

Dissecting movement of the transmembrane segments of non-gastric proton pump mutants with voltage-clamp fluorometry.

The Na+,K+-ATPase (NKA) extrudes 3 Na+ and imports 2 K+ across the plasmalemma of every animal cell. The non-gastric H+,K+- ATPase (ngHKA) exports one H+ and imports one K+ in the apical membrane of several epithelia. Both P-type pumps have ∼65% identity and nearly identical catalytic cycles driven by binding of different ions to two major conformations (E1, high affinity for Na+ or H+ and E2 higher affinity for K+). Recently, we reported the structures of the wild type (WT) ngHKA and its electrogenic mutant K794A and the NKA-like K794S/A797P/W940/R949C (SPWC) mutant, in the K+-occluded E2 state. We also reported SPWC's AMPCP-bound E1 Cryo-EM structure, which is nearly identical to that of the E1(3Na+) NKA. We obtained the WT (Na+ and AlF-ADP) and K794A (Na+ and AMPPCP) structures. While the WT structure was identical to E2P state, the AMPPCP-bound K794A structure with (presumably) 2 Na+ bound has a mixed conformation. The P and N domains as well as cytoplasmic portion of transmembrane segment TM4 shows E1-like conformation, but the A domain, TM1-TM3 and the luminal portion of TM4 took an E2-like conformation. To evaluate displacements of the moving TM segments in all three constructs, we introduced a single Cys residue in the loops between TM1-TM2, TM3-4, or TM5-6, labeled them with tetramethylrhodamine maleimide (TMRM) and evaluated currents and fluorescence changes under voltage clamp fluorometry in Xenopus oocytes. Consistent with the stable E2 state in WT, a TMRM introduced in TM1-TM2 showed minimal voltage dependent changes in fluorescence, while the voltage dependence of the TM1-TM2 fluorescent signals observed with K794A and SPWC mutants were progressively closer to signals observed in NKA. Studies with TMRM introduced at TM3-4 and TM5-6 positions are underway. Funded by NSF MCB-2003251.

Conformational diversity and environmental implications of trans-2-pentenal.

This study investigates the conformational intricacies of trans-2-pentenal (trans-2PA), a significant biogenic volatile organic compound. To unveil its potential implications in atmospheric chemistry and environmental pollution, we employ advanced infrared resonant vacuum ultraviolet mass-analysed threshold ionisation spectroscopy. Through this method, we identify the major conformers within trans-2PA, encompassing trans-s-trans (tt-) and trans-s-cis (tc-) structures with planar (cis) and non-planar (gauche) configurations introduced by a methyl group. In a pioneering spectroscopic examination, we analyze trans-2PA in both the neutral and cationic states. This approach allows us to gain a comprehensive understanding of its molecular behavior. Our conformer-specific vibrational spectra not only reveal the relative populations of the main conformers, notably tt-cis and tt-gauche conformers, but also shed light on atmospheric oxidation processes and lower tropospheric organic aerosol formation mechanisms. Our findings expand the understanding of the role of trans-2PA in environmental and biological contexts. Additionally, they contribute to a broader understanding of its influence on air quality, climate, and atmospheric dynamics. The collaboration between advanced experimental techniques and computational methods fortifies the scientific underpinning of this study, opening doors to further exploration in the realms of atmospheric chemistry and environmental science.

Assembly of the Complete Mitochondrial Genome of Gelsemium elegans Revealed the Existence of Homologous Conformations Generated by a Repeat Mediated Recombination.

Gelsemium elegans (G. elegans) is a Chinese medicinal plant with substantial economic and feeding values. There is a lack of detailed studies on the mitochondrial genome of G. elegans. In this study, the mitochondrial genome of G. elegans was sequenced and assembled, and its substructure was investigated. The mitochondrial genome of G. elegans is represented by two circular chromosomes of 406,009 bp in length with 33 annotated protein-coding genes, 15 tRNA genes, and three rRNA genes. We detected 145 pairs of repeats and found that four pairs of repeats could mediate the homologous recombination into one major conformation and five minor conformations, and the presence of conformations was verified by PCR amplification and Sanger sequencing. A total of 124 SSRs were identified in the G. elegans mitochondrial genome. The homologous segments between the chloroplast and mitochondrial genomes accounted for 5.85% of the mitochondrial genome. We also predicted 477 RNA potential editing sites and found that the nad4 gene was edited 38 times, which was the most frequent occurrence. Taken together, the mitochondrial genome of G. elegans was assembled and annotated. We gained a more comprehensive understanding on the genome of this medicinal plant, which is vital for its effective utilization and genetic improvement, especially for cytoplasmic male sterility breeding and evolution analysis in G. elegans.

Polymeric structure of the Cannabis sativa L. mitochondrial genome identified with an assembly graph model

Cannabis sativa L. belongs to the family Cannabaceae in Rosales. It has been widely used as medicines, building materials, and textiles. Elucidating its genome is critical for molecular breeding and synthetic biology study. Many studies have shown that the mitochondrial genomes (mitogenomes) and even chloroplast genomes (plastomes) had complex polymeric structures. Using the Nanopore sequencing platform, we sequenced, assembled, and analyzed its mitogenome and plastome. The resulting unitig graph suggested that the mitogenome had a complex polymeric structure. However, a gap-free, circular sequence was further assembled from the unitig graph. In contrast, a circular sequence representing the plastome was obtained. The mitogenome major conformation was 415,837 bp long, and the plastome was 153,927 bp long. To test if the repeat sequences promote recombination, which corresponds to the branch points in the structure, we tested the sequences around repeats by long-read mapping. Among 208 pairs of predicted repeats, the mapping results supported the presence of cross-over around 25 pairs of repeats. Subsequent PCR amplification confirmed the presence of cross-over around 15 of the 25 repeats. By comparing the mitogenome and plastome sequences, we identified 19 mitochondria plastid DNAs, including seven complete genes (trnW-CCA, trnP-UGG, psbJ, trnN-GUU, trnD-GUC, trnH-GUG, trnM-CAU) and nine gene fragments. Furthermore, the selective pressure analysis results showed that five genes (atp1, ccmB, ccmC, cox1, nad7) had 19 positively selected sites. Lastly, we predicted 28 RNA editing sites. A total of 8 RNA editing sites located in the coding regions were successfully validated by PCR amplification and Sanger sequencing, of which four were synonymous, and four were nonsynonymous. In particular, the RNA editing events appeared to be tissue-specific in C. sativa mitogenome. In summary, we have confirmed the major confirmation of C. sativa mitogenome and characterized its structural features in detail. These results provide critical information for future variety breeding and resource development for C. sativa.

De Novo Hybrid Assembly of the Salvia miltiorrhiza Mitochondrial Genome Provides the First Evidence of the Multi-Chromosomal Mitochondrial DNA Structure of Salvia Species.

Salvia miltiorrhiza has been an economically important medicinal plant. Previously, an S. miltiorrhiza mitochondrial genome (mitogenome) assembled from Illumina short reads, appearing to be a single circular molecule, has been published. Based on the recent reports on the plant mitogenome structure, we suspected that this conformation does not accurately represent the complexity of the S. miltiorrhiza mitogenome. In the current study, we assembled the mitogenome of S. miltiorrhiza using the PacBio and Illumina sequencing technologies. The primary structure of the mitogenome contained two mitochondrial chromosomes (MC1 and MC2), which corresponded to two major conformations, namely, Mac1 and Mac2, respectively. Using two approaches, including (1) long reads mapping and (2) polymerase chain reaction amplification followed by Sanger sequencing, we observed nine repeats that can mediate recombination. We predicted 55 genes, including 33 mitochondrial protein-coding genes (PCGs), 3 rRNA genes, and 19 tRNA genes. Repeat analysis identified 112 microsatellite repeats and 3 long-tandem repeats. Phylogenetic analysis using the 26 shared PCGs resulted in a tree that was congruent with the phylogeny of Lamiales species in the APG IV system. The analysis of mitochondrial plastid DNA (MTPT) identified 16 MTPTs in the mitogenome. Moreover, the analysis of nucleotide substitution rates in Lamiales showed that the genes atp4, ccmB, ccmFc, and mttB might have been positively selected. The results lay the foundation for future studies on the evolution of the Salvia mitogenome and the molecular breeding of S. miltiorrhiza.

Structural dynamics of the N-terminal SH2 domain of PI3K in its free and CD28-bound states.

Protein conformational changes with fluctuations are fundamental aspects of protein-protein interactions (PPIs); understanding these motions is required for the rational design of PPI-regulating compounds. Src homology 2 (SH2) domains are commonly found in adapter proteins involved in signal transduction and specifically bind to consensus motifs of proteins containing phosphorylated tyrosine (pY). Here, we analysed the interaction between the N-terminal SH2 domain (nSH2) of the regulatory subunit in phosphoinositide 3-kinase (PI3K) and the cytoplasmic region of the T-cell co-receptor, CD28, using NMR and molecular dynamics (MD) simulations. First, we assigned the backbone signals of nSH2 on 1 H-15 N heteronuclear single quantum coherence spectra in the absence or presence of the CD28 phosphopeptide, SDpYMNMTPRRPG. Chemical shift perturbation experiments revealed allosteric changes at the BC loop and the C-terminal region of nSH2 upon CD28 binding. NMR relaxation experiments showed a conformational exchange associated with CD28 binding in these regions. The conformational stabilisation of the C-terminal region correlated with the regulation of PI3K catalytic function. Further, using 19 F- and 31 P-labelled CD28 phosphopeptide, we analysed the structural dynamics of CD28 and demonstrated that the aromatic ring of the pY residue fluctuated between multiple conformations upon nSH2 binding. Our MD simulations largely explained the NMR results and the structural dynamics of nSH2 and CD28 in both bound and unbound states. Notably, in addition to its major conformation, we detected a minor conformation of nSH2 in the CD28 bound state that may explain the allosteric conformational change in the BC loop.

Mechanism of voltage gating in the voltage-sensing phosphatase Ci-VSP

Conformational changes in voltage-sensing domains (VSDs) are driven by the transmembrane electric field acting on the protein charges. Yet, the overall energetics and detailed mechanism of this process are not fully understood. Here, we determined free energy and displacement charge landscapes as well as the major conformations visited during a complete functional gating cycle in the isolated VSD of the phosphatase Ci-VSP (Ci-VSD) comprising four transmembrane helices (segments S1 to S4). Molecular dynamics simulations highlight the extent of S4 movements. In addition to the crystallographically determined activated "Up" and resting "Down" states, the simulations predict two Ci-VSD conformations: a deeper resting state ("down-minus") and an extended activated ("up-plus") state. These additional conformations were experimentally probed via systematic cysteine mutagenesis with metal-ion bridges and the engineering of proton conducting mutants at hyperpolarizing voltages. The present results show that these four states are visited sequentially in a stepwise manner during voltage activation, each step translocating one arginine or the equivalent of ∼1 e0 across the membrane electric field, yielding a transfer of ∼3 e0 charges in total for the complete process.