A single mutation from aspartate to glycine at position 614 has dominated all circulating variants of the severe acute respiratory syndrome coronavirus 2. D614G mutation induces structural changes in ...the spike (S) protein that strengthen the virus infectivity. Here, we use molecular dynamics simulations to dissect the effects of mutation and 630-loop rigidification on S-protein structure. The introduction of the mutation orders the 630-loop structure and thereby induces global structural changes toward the cryoelectron microscopy structure of the D614G S-protein. The ordered 630-loop weakens local interactions between the 614th residue and others in contrast to disordered structures in the wild-type protein. The mutation allosterically alters global interactions between receptor-binding domains, forming an asymmetric and mobile down conformation and facilitating transitions toward up conformation. The loss of salt bridge between D614 and K854 upon the mutation generally stabilizes S-protein protomer, including the fusion peptide proximal region that mediates membrane fusion. Understanding the molecular basis of D614G mutation is crucial as it dominates in all variants of concern, including Delta and Omicron.
The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of ...trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.
Conformational transitions in multidomain proteins are essential for biological functions. The Apo conformations are typically open and flexible, while the Holo states form more compact conformations ...stabilized by protein-ligand interactions. Unfortunately, the atomically detailed mechanisms for such open-closed conformational changes are difficult to be accessed experimentally as well as computationally. To simulate the transitions using atomistic molecular dynamics (MD) simulations, efficient conformational sampling algorithms are required. In this work, we propose a new approach based on generalized replica-exchange with solute tempering (gREST) for exploring the open-closed conformational changes in multidomain proteins. Wherein, selected surface charged residues in a target protein are defined as the solute region in gREST simulation and the solute temperatures are different in replicas and exchanged between them to enhance the domain motions. This approach is called gREST selected surface charged residues (gREST_SSCR) and is applied to the Apo and Holo states of ribose binding protein (RBP) in solution. The conformational spaces sampled with gREST_SSCR are much wider than those with the conventional MD, sampling open-closed conformational changes while maintaining RBP domains' stability. The free-energy landscapes of RBP in the Apo and Holo states are drawn along with twist and hinge angles of the two moving domains. The inter-domain salt-bridges that are not observed in the experimental structures are also important in the intermediate states during the conformational changes.
Proper balance between protein-protein and protein-water interactions is vital for atomistic molecular dynamics (MD) simulations of globular proteins as well as intrinsically disordered proteins ...(IDPs). The overestimation of protein-protein interactions tends to make IDPs more compact than those in experiments. Likewise, multiple proteins in crowded solutions are aggregated with each other too strongly. To optimize the balance, Lennard-Jones (LJ) interactions between protein and water are often increased about 10% (with a scaling parameter, λ = 1.1) from the existing force fields. Here, we explore the optimal scaling parameter of protein-water LJ interactions for CHARMM36m in conjunction with the modified TIP3P water model, by performing enhanced sampling MD simulations of several peptides in dilute solutions and conventional MD simulations of globular proteins in dilute and crowded solutions. In our simulations, 10% increase of protein-water LJ interaction for the CHARMM36m cannot maintain stability of a small helical peptide, (AAQAA)3 in a dilute solution and only a small modification of protein-water LJ interaction up to the 3% increase (λ = 1.03) is allowed. The modified protein-water interactions are applicable to other peptides and globular proteins in dilute solutions without changing thermodynamic properties from the original CHARMM36m. However, it has a great impact on the diffusive properties of proteins in crowded solutions, avoiding the formation of too sticky protein-protein interactions.
Conformational changes of proteins upon ligand binding are usually explained in terms of several mechanisms including the induced fit, conformational selection, or their mixtures. Due to the slow ...time scales, conventional molecular dynamics (cMD) simulations based on the atomistic models cannot easily simulate the open-to-closed conformational transition in proteins. In our previous study, we have developed an enhanced sampling scheme (generalized replica exchange with solute tempering selected surface charged residues: gREST_SSCR) for multidomain proteins and applied it to ligand-mediated conformational changes in the G134R mutant of ribose-binding protein (RBP
) in solution. The free-energy landscape (FEL) of RBP
in the presence of a ribose at the binding site included the open and closed states and two intermediates, open-like and closed-like forms. Only the open and open-like forms existed in the FEL without a ribose. In the current study, the coupling between the conformational changes and ligand binding is further investigated using coarse-grained MD, multiple atomistic cMD, and free-energy calculations. The ribose is easily dissociated from the binding site of wild-type RBP and RBP
in the cMD simulations starting from the open and open-like forms. In contrast, it is stable at the binding site in the simulations from the closed and closed-like forms. The free-energy calculations provide the binding affinities of different structures, supporting the results of cMD simulations. Importantly, cMD simulations from the closed-like structures reveal transitions toward the closed one in the presence of a bound ribose. On the basis of the computational results, we propose a molecular mechanism in which conformational selection and induced fit happen in the first and second halves of the open-to-closed transition in RBP, respectively.
Spike (S) protein is the primary antigenic target for neutralization and vaccine development for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It decorates the virus surface and ...undergoes large motions of its receptor binding domains (RBDs) to enter the host cell. Here, we observe Down, one-Up, one-Open, and two-Up-like structures in enhanced molecular dynamics simulations, and characterize the transition pathways via inter-domain interactions. Transient salt-bridges between RBD
and RBD
and the interaction with glycan at N343
support RBD
motions from Down to one-Up. Reduced interactions between RBD
and RBD
in one-Up induce RBD
motions toward two-Up. The simulations overall agree with cryo-electron microscopy structure distributions and FRET experiments and provide hidden functional structures, namely, intermediates along Down-to-one-Up transition with druggable cryptic pockets as well as one-Open with a maximum exposed RBD. The inherent flexibility of S-protein thus provides essential information for antiviral drug rational design or vaccine development.
We serendipitously found that chaperonin GroEL can hydrolyze
-nitrophenyl β-galactoside (ONPG), a well-known substrate of the enzyme β-galactosidase. The ONPG hydrolysis by GroEL follows typical ...enzyme kinetics. Our experiments and molecular docking studies suggest ONPG binding at the ATP binding site of GroEL.
The initial step to resolving the controversy regarding the repair mechanism of the pyrimidine (6-4) pyrimidone photoproduct DNA lesion is to determine the protonation state of the two conserved ...active site histidine residues (His365 and His369 in Drosophila melanogaster) in the (6-4) photolyase enzyme ((6-4) PHR). Both His residues were experimentally determined to be crucial for catalysis. Most previous theoretical studies assumed the presence of protonated His365 in the active site, which would transfer its proton to N3′ of the lesion as the first step in repair; however, other empirical/semiempirical pK a calculations suggested the presence of two neutral His residues in the active site. Here, we conduct molecular dynamics simulations (MD) of the (6-4) PHR/DNA complex on the basis of the X-ray crystal structure to investigate all combinations of the three possible protonation states of each His. The MD results show that HIP365 (both ND1 and NE2 protonated) and HID369 (only ND1 protonated) are the most probable protonation states in the active site. Furthermore, we employ quantum mechanics-cluster and quantum mechanical/molecular mechanical (QM/MM) approaches to investigate the protonation of N3′, starting with three plausible complexes resulting from MD (HIP365/HID369, HID365/HIE369, and HID365/HIP369). Surprisingly, protonation of the N3′ atom is found to be feasible starting from all three protonation states: protonated HIP365 transfers the proton via a barrierless step, and in the other cases of neutral HID365, the proton is transferred from O4′ of the lesion via His365 through an energy barrier of ∼80 kJ mol–1 in both complexes. These results might explain the conservation of repair activity over a wide range of pH values.
Ultraviolet (UV) light from the sun damages DNA by forming a cyclobutane pyrimidine dimer (CPD) and pyrimidine(6-4)pyrimidone photoproducts (6-4) PP. Photolyase (PHR) enzymes utilize near-UV/blue ...light for DNA repair, which is initiated by light-induced electron transfer from the fully reduced flavin adenine dinucleotide chromophore. Despite similar structures and repair mechanisms, the functions of PHR are highly selective; CPD PHR repairs CPD, but not (6-4) PP, and vice versa. In this study, we attempted functional conversion between CPD and (6-4) PHRs. We found that a triple mutant of (6-4) PHR is able to repair the CPD photoproduct, though the repair efficiency is 1 order of magnitude lower than that of wild-type CPD PHR. Difference Fourier transform infrared spectra for repair demonstrate the lack of secondary structural alteration in the mutant, suggesting that the triple mutant gains substrate binding ability while it does not gain the optimized conformational changes from light-induced electron transfer to the release of the repaired DNA. Interestingly, the (6-4) photoproduct is not repaired by the reverse mutation of CPD PHR, and eight additional mutations (total of 11 mutations) introduced into CPD PHR are not sufficient. The observed asymmetric functional conversion is interpreted in terms of a more complex repair mechanism for (6-4) repair, which was supported by quantum chemical/molecular mechanical calculation. These results suggest that CPD PHR may represent an evolutionary origin for photolyase family proteins.