Prefusion spike protein conformational changes are slower in SARS-CoV-2 than in SARS-CoV-1.

Within the last 2 decades, severe acute respiratory syndrome coronaviruses 1 and 2 (SARS-CoV-1 and SARS-CoV-2) have caused two major outbreaks; yet, for reasons not fully understood, the coronavirus disease 2019 pandemic caused by SARS-CoV-2 has been significantly more widespread than the 2003 SARS epidemic caused by SARS-CoV-1, despite striking similarities between these two viruses. The SARS-CoV-1 and SARS-CoV-2 spike proteins, both of which bind to host cell angiotensin-converting enzyme 2, have been implied to be a potential source of their differential transmissibility. However, the mechanistic details of prefusion spike protein binding to angiotensin-converting enzyme 2 remain elusive at the molecular level. Here, we performed an extensive set of equilibrium and nonequilibrium microsecond-level all-atom molecular dynamics simulations of SARS-CoV-1 and SARS-CoV-2 prefusion spike proteins to determine their differential dynamic behavior. Our results indicate that the active form of the SARS-CoV-2 spike protein is more stable than that of SARS-CoV-1 and the energy barrier associated with the activation is higher in SARS-CoV-2. These results suggest that not only the receptor-binding domain but also other domains such as the N-terminal domain could play a crucial role in the differential binding behavior of SARS-CoV-1 and SARS-CoV-2 spike proteins.

cpSRP43 Is Both Highly Flexible and Stable: Structural Insights Using a Combined Experimental and Computational Approach.

The novel multidomain protein, cpSRP43, is a unique subunit of the post-translational chloroplast signal recognition particle (cpSRP) targeting pathway in higher plants. The cpSRP pathway is responsible for targeting and insertion of light-harvesting chlorophyll a/b binding proteins (LHCPs) to the thylakoid membrane. Upon emergence into the stroma, LHCPs form a soluble transit complex with the cpSRP heterodimer, which is composed of cpSRP43 and cpSRP54. cpSRP43 is irreplaceable as a chaperone to LHCPs in their translocation to the thylakoid membrane and remarkable in its ability to dissolve aggregates of LHCPs without the need for external energy input. In previous studies, cpSRP43 has demonstrated significant flexibility and interdomain dynamics. In this study, we explore the structural stability and flexibility of cpSRP43 using a combination of computational and experimental techniques and find that this protein is concurrently highly stable and flexible. In addition to microsecond-level unbiased molecular dynamics (MD), biased MD simulations based on system-specific collective variables are used along with biophysical experimentation to explain the basis of the flexibility and stability of cpSRP43, showing that the free and cpSRP54-bound cpSRP43 has substantially different conformations and conformational dynamics.

Transient local secondary structure in the intrinsically disordered C-term of the Albino3 insertase.

Albino3 (Alb3) is an integral membrane protein fundamental to the targeting and insertion of light-harvesting complex (LHC) proteins into the thylakoid membrane. Alb3 contains a stroma-exposed C-terminus (Alb3-Cterm) that is responsible for binding the LHC-loaded transit complex before LHC membrane insertion. Alb3-Cterm has been reported to be intrinsically disordered, but precise mechanistic details underlying how it recognizes and binds to the transit complex are lacking, and the functional roles of its four different motifs have been debated. Using a novel combination of experimental and computational techniques such as single-molecule fluorescence resonance energy transfer, circular dichroism with deconvolution analysis, site-directed mutagenesis, trypsin digestion assays, and all-atom molecular dynamics simulations in conjunction with enhanced sampling techniques, we show that Alb3-Cterm contains transient secondary structure in motifs I and II. The excellent agreement between the experimental and computational data provides a quantitatively consistent picture and allows us to identify a heterogeneous structural ensemble that highlights the local and transient nature of the secondary structure. This structural ensemble was used to predict both the inter-residue distance distributions of single molecules and the apparent unfolding free energy of the transient secondary structure, which were both in excellent agreement with those determined experimentally. We hypothesize that this transient local secondary structure may play an important role in the recognition of Alb3-Cterm for the LHC-loaded transit complex, and these results should provide a framework to better understand protein targeting by the Alb3-Oxa1-YidC family of insertases.

Binding affinity estimation from restrained umbrella sampling simulations.

The protein-ligand binding affinity quantifies the binding strength between a protein and its ligand. Computer modeling and simulations can be used to estimate the binding affinity or binding free energy using data- or physics-driven methods or a combination thereof. Here we discuss a purely physics-based sampling approach based on biased molecular dynamics simulations. Our proposed method generalizes and simplifies previously suggested stratification strategies that use umbrella sampling or other enhanced sampling simulations with additional collective-variable-based restraints. The approach presented here uses a flexible scheme that can be easily tailored for any system of interest. We estimate the binding affinity of human fibroblast growth factor 1 to heparin hexasaccharide based on the available crystal structure of the complex as the initial model and four different variations of the proposed method to compare against the experimentally determined binding affinity obtained from isothermal titration calorimetry experiments.

Mechanistic Picture for Monomeric Human Fibroblast Growth Factor 1 Stabilization by Heparin Binding.

Human fibroblast growth factor (FGF) 1 or hFGF1 is a member of the FGF family that is involved in various vital processes such as cell proliferation, cell differentiation, angiogenesis, and wound healing. hFGF1, which is associated with low stability in vivo, is known to be stabilized by binding heparin sulfate, a glycosaminoglycan that aids the protein in the activation of its cell surface receptor. The poor thermal and proteolytic stability of hFGF1 and the stabilizing role of heparin have long been observed experimentally; however, the mechanistic details of these phenomena are not well understood. Here, we have used microsecond-level equilibrium molecular dynamics (MD) simulations to quantitatively characterize the structural dynamics of monomeric hFGF1 in the presence and absence of heparin hexasaccharide. We have observed a conformational change in the heparin-binding pocket of hFGF1 that occurs only in the absence of heparin. Several intramolecular interactions were also identified within the heparin-binding pocket that form only when hFGF1 interacts with heparin. The loss of both intermolecular and intramolecular interactions in the absence of heparin plausibly leads to the observed conformational change. This conformational transition results in increased flexibility of the heparin-binding pocket and provides an explanation for the susceptibility of apo hFGF1 to proteolytic degradation and thermal instability. This study provides a glimpse into mechanistic details of the heparin-mediated stabilization of hFGF1 and encourages the use of microsecond-level MD in studying the effect of binding on protein structure and dynamics. In addition, the observed differential behavior of hFGF1 in the absence and presence of heparin provides an example, where microsecond-level all-atom MD simulations are necessary to see functionally relevant biomolecular phenomena that otherwise will not be observed on sub-microsecond time scales.

Differential Dynamic Behavior of Prefusion Spike Proteins of SARS Coronaviruses 1 and 2.

The coronavirus spike protein, which binds to the same human receptor in both SARS-CoV-1 and 2, has been implied to be a potential source of their differential transmissibility. However, the mechanistic details of spike protein binding to its human receptor remain elusive at the molecular level. Here, we have used an extensive set of unbiased and biased microsecond-level all-atom molecular dynamics (MD) simulations of SARS-CoV-1 and 2 spike proteins to determine the differential dynamic behavior of prefusion spike protein structure in the two viruses. Our results indicate that the active form of the SARS-CoV-2 spike protein is more stable than that of SARS-CoV-1 and the energy barrier associated with the activation is higher in SARS-CoV-2. Our results also suggest that not only the receptor binding domain (RBD) but also other domains such as the N-terminal domain (NTD) could play a role in the differential binding behavior of SARS-CoV-1 and 2 spike proteins.

Characterization of the structural forces governing the reversibility of the thermal unfolding of the human acidic fibroblast growth factor.

Human acidic fibroblast growth factor (hFGF1) is an all beta-sheet protein that is involved in the regulation of key cellular processes including cell proliferation and wound healing. hFGF1 is known to aggregate when subjected to thermal unfolding. In this study, we investigate the equilibrium unfolding of hFGF1 using a wide array of biophysical and biochemical techniques. Systematic analyses of the thermal and chemical denaturation data on hFGF1 variants (Q54P, K126N, R136E, K126N/R136E, Q54P/K126N, Q54P/R136E, and Q54P/K126N/R136E) indicate that nullification of charges in the heparin-binding pocket can significantly increase the stability of wtFGF1. Triple variant (Q54P/K126N/R136E) was found to be the most stable of all the hFGF1 variants studied. With the exception of triple variant, thermal unfolding of wtFGF1 and the other variants is irreversible. Thermally unfolded triple variant refolds completely to its biologically native conformation. Microsecond-level molecular dynamic simulations reveal that a network of hydrogen bonds and salt bridges linked to Q54P, K126N, and R136E mutations, are responsible for the high stability and reversibility of thermal unfolding of the triple variant. In our opinion, the findings of the study provide valuable clues for the rational design of a stable hFGF1 variant that exhibits potent wound healing properties.

The Role of a Crystallographically Unresolved Cytoplasmic Loop in Stabilizing the Bacterial Membrane Insertase YidC2.

YidC, a bacterial member of the YidC/Alb3/Oxa1 insertase family, mediates membrane protein assembly and insertion. Cytoplasmic loops are known to have functional significance in membrane proteins such as YidC. Employing microsecond-level molecular dynamics (MD) simulations, we show that the crystallographically unresolved C2 loop plays a crucial role in the structural dynamics of Bacillus halodurans YidC2. We have modeled the C2 loop and used all- atom MD simulations to investigate the structural dynamics of YidC2 in its apo form, both with and without the C2 loop. The C2 loop was found to stabilize the entire protein and particularly the C1 region. C2 was also found to stabilize the alpha-helical character of the C-terminal region. Interestingly, the highly polar or charged lipid head groups of the simulated membranes were found to interact with and stabilize the C2 loop. These findings demonstrate that the crystallographically unresolved loops of membrane proteins could be important for the stabilization of the protein despite the apparent lack of structure, which could be due to the absence of the relevant lipids to stabilize them in crystallographic conditions.

An Effective Electric Dipole Model for Voltage-induced Gating Mechanism of Lysenin.

Lysenin is a pore-forming toxin, which self-inserts open channels into sphingomyelin containing membranes and is known to be voltage regulated. The mechanistic details of its voltage gating mechanism, however, remains elusive despite much recent efforts. Here, we have employed a novel combination of experimental and computational techniques to examine a model for voltage gating, that is based on the existence of an "effective electric dipole" inspired by recent reported structures of lysenin. We support this mechanism by the observations that (i) the charge-reversal and neutralization substitutions in lysenin result in changing its electrical gating properties by modifying the strength of the dipole, and (ii) an increase in the viscosity of the solvent increases the drag force and slows down the gating. In addition, our molecular dynamics (MD) simulations of membrane-embedded lysenin provide a mechanistic picture for lysenin conformational changes, which reveals, for the first time, the existence of a lipid-dependent bulge region in the pore-forming module of lysenin, which may explain the gating mechanism of lysenin at a molecular level.

ALS-causing mutations in profilin-1 alter its conformational dynamics: A computational approach to explain propensity for aggregation.

Profilin-1 (PFN1) is a 140-amino-acid protein with two distinct binding sites-one for actin and one for poly-L-proline (PLP). The best-described function of PFN1 is to catalyze actin elongation and polymerization. Thus far, eight DNA mutations in the PFN1 gene encoding the PFN1 protein are associated with human amyotrophic lateral sclerosis (ALS). We and others recently showed that two of these mutations (Gly118Val or G118V and Cys71Gly or C71G) cause ALS in rodents. In vitro studies suggested that Met114Thr and Thr109Met cause the protein to behave abnormally and cause neurotoxicity. The mechanism by which a single amino acid change in human PFN1 causes the degeneration of motor neurons is not known. In this study, we investigated the structural perturbations of PFN1 caused by each ALS-associated mutation. We used molecular dynamics simulations to assess how these mutations alter the secondary and tertiary structures of human PFN1. Herein, we present our in silico data and analysis on the effect of G118V and T109M mutations on PFN1 and its interactions with actin and PLP. The substitution of valine for glycine reduces the conformational flexibility of the loop region between the α-helix and ß-strand and enhances the hydrophobicity of the region. Our in silico analysis of T109M indicates that this mutation alters the shape of the PLP-binding site and reduces the flexibility of this site. Simulation studies of PFN1 in its wild type (WT) and mutant forms (both G118V and T109M mutants) revealed differential fluctuation patterns and the formation of salt bridges and hydrogen bonds between critical residues that may shed light on differences between WT and mutant PFN1. In particular, we hypothesize that the flexibility of the actin- and PLP-binding sites in WT PFN1 may allow the protein to adopt slightly different conformations in its free and bound forms. These findings provide new insights into how each of these mutations in PFN1 might increase its propensity for misfolding and aggregation, leading to its dysfunction.