Simulation of the ligand-leaving process of the human heat shock protein.

The human heat shock protein plays a critical role in various diseases and is an important target for pharmacological modulation. Simulation of conformational changes and free energy profiles of the human heat shock protein derived by the ligand-leaving process is a challenging issue. In this work, steered molecular dynamics simulation was adopted to simulate the ligand-leaving process. Two composite systems of heat shock protein NHSP90 and small molecules 6FJ and 6G7 are selected as research objects. The free energy during the leaving of ligand small molecules is calculated using conventional molecular dynamics simulation, steered molecular dynamics simulation (SMD), and the umbrella sampling method. We found that the a slower pulling velocity (0.001 nm ns-1) will result in 2.19 kcal mol-1, and the umbrella sampling method gives a value of 3.26 kcal mol-1 for the free energy difference for the two systems, which reasonably agrees with experimental results. A faster-pulling velocity (0.01 nm ns-1) leads to a large overestimation of free energy. At the same time, the conformational analysis indicated that the faster pulling velocity may lead to the conformational change of NHSP90, which was proved to be false by the slower pulling velocity and the umbrella sampling method.

Million-atom molecular dynamics simulations reveal the interfacial interactions and assembly of plant PSII-LHCII supercomplex.

Protein-protein interface interactions dictate efficient excitation energy transfer from light-harvesting antennas to the photosystem II (PSII) core. In this work, we construct a 1.2 million atom-scale model of plant C2S2-type PSII-LHCII supercomplex and perform microsecond-scale molecular dynamics (MD) simulations to explore the interactions and assembly mechanisms of the sizeable PSII-LHCII supercomplex. We optimize the nonbonding interactions of the PSII-LHCII cryo-EM structure using microsecond-scale MD simulations. Binding free energy calculations with component decompositions reveal that hydrophobic interactions predominantly drive antenna-core association and the antenna-antenna interactions are relatively weak. Despite the positive electrostatic interaction energies, hydrogen bonds and salt bridges mainly provide directional or anchoring forces for interface binding. Analysis of the roles of small intrinsic subunits of PSII suggests that LHCII and CP26 first interact with small intrinsic subunits and then bind to the core proteins, whereas CP29 adopts a one-step binding process to the PSII core without the assistance of other factors. Our study provides insights into the molecular underpinnings of the self-organization and regulation of plant PSII-LHCII. It lays the framework for deciphering the general assembly principles of photosynthetic supercomplexes and possibly other macromolecular structures. The finding also has implications for repurposing photosynthetic systems to enhance photosynthesis.

Insights into the binding mode of AS1411 aptamer to nucleolin.

AS1411 aptamer can function as a recognition probe to detect the cell surface nucleolin overexpressed in cancer cells, however, little is known about their binding process. This study proposed a feasible binding mode for the first time and provided atomic-level descriptions for the high affinity and specific binding of AS1411. The binding pose predicted by docking was screened using knowledge-based criteria, and a microsecond molecular dynamics (MD) simulation showed the stable existence of the predicted structure in the solution. Structural analysis shows that the unique capping of the 5' end of AS1411 provides the specific binding with RBD1, and the interactions of hydrogen bond, salt bridge, and water-mediated network between AS1411 and RBD1,2 stabilize the binding. The calculation of per-residue decomposition emphasizes the dominant contribution of van der Waals energy and critical residues are screened. Our study provides the molecular basis of this specific binding and can guide rational AS1411-based aptamers design. Further insights require tight collaborations between the experiments and in silico studies.

Rapid and sensitive multiplex detection of COVID-19 antigens and antibody using electrochemical immunosensor-/aptasensor-enabled biochips.

We report protein- and aptamer-based electrochemical biochips for low-cost, one-step, sensitive and accurate multiplex detection of SARS-CoV-2 spike (S) and nucleocapsid (N) proteins, and IgG antibody in unprocessed clinical samples, allowing citizens to achieve rapid diagnosis at home or in community settings.

Structural Insights into the Cofactor Role of Heparin/Heparan Sulfate in Binding between the SARS-CoV-2 Spike Protein and Host Angiotensin-Converting Enzyme II.

The viral entry process of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) requires heparin and heparan sulfates from the cell surface, functioning as a cofactor for human angiotensin-converting enzyme 2 (ACE2) for recognizing the receptor-binding domain (RBD) of the spike (S) protein on the surface of the virion. In the present study, the binding poses of an oligosaccharide with four repeating units of GlcNS6S-IdoA2S (octa) predicted by Vina-Carb in the RBD binding site were employed in molecular dynamics (MD) simulations to provide atomic details for studying the cofactor mechanism. The molecular model in the MD simulations reproduced the length- and sequence-dependent behavior observed from the microarray experiments and revealed an important planar U-turn shape for HP/HS binding to RBD. The model for octa with this shape in the ACE2-RBD complex enhanced the interactions in the binding interface. The comparisons with the ACE2-RBD complex suggested that the presence of octa in the RBD binding site blocked the movements in a loop region at the distal end of the RBD binding interface and promoted the contacts of this loop region with the ACE2 N-terminus helix. This study shed light on the atomic and dynamic details for HP/HS interacting with RBD and provided insights into their cofactor role in the ACE2-RBD interactions.

Insights into Selectin Inhibitor Design from Endogenous Isomeric Ligands of SLea and SLex.

Selectins interact with cell-surface glycans to promote the initial tethering and rolling of leukocytes, and these interactions are targets for designs of inhibitors to neutralize diseases related to excessive inflammatory responses in many cardiovascular and immune dysfunctions, as well as tumor markers in different cancers. The isomeric endogenous tetrasaccharides, sialyl Lewis X (sLex) and sialyl Lewis A (sLea), are minimal sugar structures required for selectin binding. Understanding their subtle structural variances and significant advanced binding strengths of sLea over sLex could benefit the rational designs for selectin inhibitors. Modeling based on the E-selectin-sLex crystal structure in the present study demonstrated that the N-acetyl group of GlcNAc in sLex could form steric hindrances in the E-selectin-sLex complex, but the hydroxy methylene group of GlcNAc in sLea at the same position allows for stronger binding interactions. The subsequent designed inhibitor with a synthetic accessible linker molecule that has no exo-cyclic moieties replacing GlcNAc displayed comparable dynamic and energetic binding features to sLea. The present study deciphered the clues from endogenous isomeric sLea and sLex and provided insights into designing selectin inhibitors with simplified synthesis.

Molecular mechanism of methyl-dependent and spatial-specific DNA recognition of c-Jun homodimer.

DNA methylation is important in regulation of gene expression and normal development because it alters the interplay between protein and DNA. Experiments have shown that a single 5-methylcytosine at different CpG sites (mCpG) might have different effects on specific recognition, but the atomistic origin and dynamic details are largely unclear. In this work, we investigated the mechanism of monomethylation at different CpG sites in the cognate motif and the cooperativity of full methylation. By constructing four models of c-Jun/Jun protein binding to the 5[Formula: see text]-XGAGTCA-3[Formula: see text] (X represents C or methylated C) motif, we characterized the dynamics of the contact interface using the all-atom molecular dynamics method. Free energy analysis of MM/GBSA suggests that regardless of whether the C12pG13 site of the bottom strand is methylated, the effects from mC25 of the top strand are dominant and can moderately enhance the binding by [Formula: see text] 31 kcal/mol, whereas mC12 showed a relatively small contribution, in agreement with the experimental data. Remarkably, we found that this spatial-specific influence was induced by different regulatory rules. The influence of the mC25 site is mainly mediated by steric hindrance. The additional methyl group leads to the conformational changes in nearby residues and triggers an obvious structural bending in the protein, which results in the formation of a new T-Asn-C triad that enhances the specific recognition of TCA half-sites. The substitution of the methyl group at the mC12 site of the bottom strand breaks the original H-bonds directly. Such changes in electrostatic interactions also lead to the remote allosteric effects of protein by multifaceted interactions but have negligible contributions to binding. Although these two influence modes are different, they can both fine-tune the local environment, which might produce remote allosteric effects through protein-protein interactions. Further analysis reveals that the discrepancies in these two modes are primarily due to their location. Moreover, when both sites are methylated, the major determinant of binding specificity depends on the context and the location of the methylation site, which is the result of crosstalk and cooperativity.

Antiviral drug design based on the opening mechanism of spike glycoprotein in SARS-CoV-2.

The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters the host cell after the receptor binding domain (RBD) of the virus spike (S) glycoprotein binds to the human angiotensin-converting enzyme 2 (hACE2). This binding requires the RBD to undergo a conformational change from a closed to an open state. In the present study, a key pair of salt bridges formed by the side chains of K537 and E619, residues at the interfaces of SD1 and SD2, respectively, was identified to promote the opening of the RBD. Mutations of K537Q and E619D reduced their side chain lengths and eliminated this pair of salt bridges; as a result, the opening of the RBD was not observed in the MD simulations. Thus, blocking the formation of this pair of salt bridges is a promising approach for treating novel coronavirus disease 2019 (COVID-19). FDA approved drug molecules were screened by their capabilities of blocking the formation of the key pair of salt bridges, achieved by their positional stabilities in the cavity containing the side chains of K537 and E619 formed in the interface between SD1 and SD2. Simeprevir, imatinib, and naldemedine were identified to possess the desired capability with the most favorable interaction energies.

Correction: Antiviral drug design based on the opening mechanism of spike glycoprotein in SARS-CoV-2.

Correction for 'Antiviral drug design based on the opening mechanism of spike glycoprotein in SARS-CoV-2' by Ruichao Mao et al., Phys. Chem. Chem. Phys., 2021, DOI: 10.1039/d1cp01045j.

Insights into the Loop at the E-Selectin Binding Site: From Open to Close Conformation.

The loop at the E-selectin binding site displayed open and close conformations observed in crystal structures before and after complexing with sialyl lewis x (sLex), respectively, and these different conformations were less studied and could affect the binding and dissociation of selectin/sLex that are essential for the recruitment of leukocytes and early inflammatory response. Hereby, we studied the roles of different loop conformations by performing molecular dynamics simulations, including adaptive steered MD simulations and energy calculations. Results suggested that the loop in the E-selectin binding site could switch from open to close conformation after the binding of sLex spontaneously, and the close conformation enhanced the binding by making sLex immersed slightly deeper in the binding site. Potential mean force calculations showed that more work was required for sLex to dissociate when the loop was in the close conformation, benefiting the formation of the catch bonds and prolonging the bonding lifetime by having more durable interactions between sLex and the loop residues in the rebinding step. This study provided atomic and dynamic details of the influence of the loop conformations on E-selectin/sLex interactions and further elucidated their mechanisms.

Toward the Prediction of Multi-Spin State Charges of a Heme Model by Random Forest Regression.

The random forest regression (RFR) model was introduced to predict the multiple spin state charges of a heme model, which is important for the molecular dynamic simulation of the spin crossover phenomenon. In this work, a multiple spin state structure data set with 39,368 structures of the simplified heme-oxygen binding model was built from the non-adiabatic dynamic simulation trajectories. The ESP charges of each atom were calculated and used as the real-valued response. The conformational adapted charge model (CAC) of three spin states was constructed by an RFR model using symmetry functions. The results show that our RFR model can effectively predict the on the fly atomic charges with the varying conformations as well as the atomic charge of different spin states in the same conformation, thus achieving the balance of accuracy and efficiency. The average mean absolute error of the predicted charges of each spin state is <0.02 e. The comparison studies on descriptors showed a maximum 0.06 e improvement in prediction of the charge of Fe 2+ by using 11 manually selected structural parameters. We hope that this model can not only provide variable parameters for developing the force field of the multi-spin state but also facilitate automation, thus enabling large-scale simulations of atomistic systems.

How a single 5-methylation of cytosine regulates the recognition of C/EBPß transcription factor: a molecular dynamic simulation study.

CpG methylation can regulate gene expression by altering the specific binding of protein and DNA. In order to understand how a single 5mC regulates protein-DNA interactions, we have compared the structures and dynamics of CEBP/ßprotein-DNA complexes before and after methylation, and the results indicate that even a single 5mC can regulate protein-DNA recognition by steric-hindrance effect of methyl group and changing the hydrogen bond interactions. The interactions between the methyl group, mCpG motif, and the conserved residue arginine make the protein read out the variation of local environment, which further enhances the specific recognition and affects the base pair stacking. The stacking interactions can propagate along the backbone of DNA and lead to long-range allosteric effects, including obvious conformational variations for DNA base pairs.

Spin crossover dynamics studies on the thermally activated molecular oxygen binding mechanism on a model copper complex.

The theoretical description of the primary dioxygen (O2) binding and activation step in many copper or iron enzymes, suffers from the intrinsically electronic non-adiabaticity of the spin flip events of the triplet dioxygen molecule (3O2), mediated by spin-orbit couplings. In this work, we presented the early-stage ultrafast spin flip dynamics of O2 binding for a simplified monocopper complex, involving the coupled singlet and triplet electronic states. The on-the-fly trajectory surface hopping (TSH) simulations have identified the dynamical effects that may influence the mode of O2 coordination (end-on vs. side-on), and the electronic structures can be viewed as complexes of molecular O2 with Cu(i) or as Cu(ii)-superoxide compounds. In addition, significant spin flip events are obversed within the sub-picosecond regime. We hope this work may provide complimentary insights on the traditional interpretation of O2 binding on copper complexes and subsequent catalytic reaction mechanisms.

Trinucleotide Base Pair Stacking Free Energy for Understanding TF-DNA Recognition and the Functions of SNPs.

Single nucleotide polymorphisms (SNPs) affect base pair stacking, which is the primary factor for maintaining the stability of DNA. However, the mechanism of how SNPs lead to phenotype variations is still unclear. In this work, we connected SNPs and base pair stacking by a 3-mer base pair stacking free energy matrix. The SNPs with large base pair stacking free energy differences led to phenotype variations. A molecular dynamics (MD) simulation was then applied. Our results showed that base pair stacking played an important role in the transcription factor (TF)-DNA interaction. Changes in DNA structure mainly originate from TF-DNA interactions, and with the increased base pair stacking free energy, the structure of DNA approaches its free type, although its binding affinity was increased by the SNP. In addition, quantitative models using base pair stacking features revealed that base pair stacking can be used to predict TF binding specificity. As such, our work combined knowledge from bioinformatics and structural biology and provided a new understanding of the relationship between SNPs and phenotype variations. The 3-mer base pair stacking free energy matrix is useful in high-throughput screening of SNPs and predicting TF-DNA binding affinity.