Insight into the dynamic interaction of different carbohydrates with human surfactant protein D: molecular dynamics simulations

Comparative Assessment of Water Models in Protein-Glycan Interaction: Insights from Alchemical Free Energy Calculations and Molecular Dynamics Simulations

Accurate computational simulations of protein-glycan dynamics are crucial for a comprehensive understanding of critical biological mechanisms, including host-pathogen interactions, immune system defenses, and intercellular communication. The accuracy of these simulations, including molecular dynamics (MD) simulation and alchemical free energy calculations, critically relies on the appropriate parameters, including the water model, because of the extensive hydrogen bonding with glycan hydroxyl groups. However, a systematic evaluation of water models' accuracy in simulating protein-glycan interaction at the molecular level is still lacking. In this study, we used full atomistic MD simulations and alchemical absolute binding free energy (ABFE) calculations to investigate the performance of five distinct water models in six protein-glycan complex systems. We evaluated water models' impact on structural dynamics and binding affinity through over 5.8 μs of simulation time per system. Our results reveal that most protein-glycan complexes are stable in the overall structural dynamics regardless of the water model used, while some show obvious fluctuations with specific water models. More importantly, we discover that the stability of the binding motif's conformation is dependent on the water model chosen when its residues form weak hydrogen bonds with the glycan. The water model also influences the conformational stability of the glycan in its bound state according to density functional theory (DFT) calculations. Using alchemical ABFE calculations, we find that the OPC water model exhibits exceptional consistency with experimental binding affinity data, whereas commonly used models such as TIP3P are less accurate. The findings demonstrate how different water models affect protein-glycan interactions and the accuracy of binding affinity calculations, which is crucial in developing therapeutic strategies targeting these interactions.

Elucidating the enhanced binding affinity of a double mutant SP-D with trimannose on the influenza A virus using molecular dynamics

•

The Asp325Ala mutation in SP-D promotes a trimannose conformational change to a more stable state.

•

The Arg343Val mutation in SP-D reduces its interaction with Glu333 to increase the binding affinity with trimannose.

•

The Arg343Val mutation contributes more to the increase of SP-D’s binding affinity with trimannose than Asp325Ala.

Surfactant protein D (SP-D) is an essential component of the human pulmonary surfactant system, which is crucial in the innate immune response against glycan-containing pathogens, including Influenza A viruses (IAV) and SARS-CoV-2. Previous studies have shown that wild-type (WT) SP-D can bind IAV but exhibits poor antiviral activities. However, a double mutant (DM) SP-D consisting of two point mutations (Asp325Ala and Arg343Val) inhibits IAV more potently. Presently, the structural mechanisms behind the point mutations’ effects on SP-D’s binding affinity with viral surface glycans are not fully understood. Here we use microsecond-scale, full-atomistic molecular dynamics (MD) simulations to understand the molecular mechanism of mutation-induced SP-D’s higher antiviral activity. We find that the Asp325Ala mutation promotes a trimannose conformational change to a more stable state. Arg343Val increases the binding with trimannose by increasing the hydrogen bonding interaction with Glu333. Free energy perturbation (FEP) binding free energy calculations indicate that the Arg343Val mutation contributes more to the increase of SP-D’s binding affinity with trimannose than Asp325Ala. This study provides a molecular-level exploration of how the two mutations increase SP-D binding affinity with trimannose, which is vital for further developing preventative strategies for related diseases.

MD Simulation Investigation on the Binding Process of Smoke-Derived Germination Stimulants to Its Receptor

Karrikins (KARs) are a class of smoke-derived seed germination stimulants with great significance in both agriculture and plant biology. By means of direct binding to the receptor protein KAI2, the compounds can initiate the KAR signal transduction pathway, hence triggering germination of the dormant seeds in the soil. In the research, several molecular dynamics (MD) simulation techniques were properly integrated to investigate the binding process of KAR₁ to KAI2 and reveal the details of the whole binding event. The calculated binding free energy, -7.00 kcal/mol, is in good agreement with the experimental measurement, -6.83 kcal/mol. The obtained PMF profile indicates the existence of three intermediate states in the binding process. The analysis of the simulation trajectories demonstrates that, in the intermediate structures, KAR₁ is stabilized by some hydrophobic residues (Phe26, Phe134, Leu142, Trp153, Phe157, Leu160, Phe194), along with several bridging water molecules, and meanwhile, the significant shifting occurs in the local conformation of the protein as the ligand's binding. A series of the residues (Gln141-Phe157) on the so-called "cap domain" are proposed to be responsible for capturing the ligand at the initial stage of the binding. Besides, the changes of the ligand's poses are also quantitatively characterized by the proper choice of the coordinate system. Our work will contribute to the more penetrating understanding of the ligand binding process and the receptor affinity difference between several members in the KAR family and help design new, more effective germination stimulants.

Effect of External Electric Field on Substrate Transport of a Secondary Active Transporter

Substrate transport across a membrane accomplished by a secondary active transporter (SAT) is essential to the normal physiological function of living cells. In the present research, a series of all-atom molecular dynamics (MD) simulations under different electric field (EF) strengths was performed to investigate the effect of an external EF on the substrate transport of an SAT. The results show that EF both affects the interaction between substrate and related protein's residues by changing their conformations and tunes the timeline of the transport event, which collectively reduces the height of energy barrier for substrate transport and results in the appearance of two intermediate conformations under the existence of an external EF. Our work spotlights the crucial influence of external EFs on the substrate transport of SATs and could provide a more penetrating understanding of the substrate transport mechanism of SATs.

Insight into the interactive residues between two domains of human somatic Angiotensin-converting enzyme and Angiotensin II by MM-PBSA calculation and steered molecular dynamics simulation

Angiotensin-converting enzyme (ACE), a membrane-bound zinc metallopeptidase, catalyzes the formation of Angiotensin-II (AngII) and the deactivation of bradykinin in the renin-angiotensin-aldosterone and kallikrein-kinin systems. As a hydrolysis product of ACE, AngII is regarded as an inhibitor and displays stronger competitive inhibition in the C-domain than the N-domain of ACE. However, the AngII binding differences between the two domains and the mechanisms behind AngII dissociation from the C-domain are rarely explored. In this work, molecular docking, Molecular Mechanics/Poisson-Boltzmann Surface Area calculation, and steered molecular dynamics (SMD) are applied to explore the structures and interactions in the binding or unbinding of AngII with the two domains of human somatic ACE. Calculated free energy values suggest that the C-domain-AngII complex is more stable than the N-domain-AngII complex, consistent with available experimental data. SMD simulation results imply that electrostatic interaction is dominant in the dissociation of AngII from the C-domain. Moreover, Gln106, Asp121, Glu123, and Tyr213 may be the key residues in the unbinding pathway of AngII. The simulation results in our work provide insights into the interactions between the two domains of ACE and its natural peptide inhibitor AngII at a molecular level. Moreover, the results provide theoretical clues for the design of new inhibitors.

Molecular mechanisms of inhibition of influenza by surfactant protein D revealed by large-scale molecular dynamics simulation

Surfactant protein D (SP-D), a mammalian C-type lectin, is the primary innate inhibitor of influenza A virus (IAV) in the lung. Interactions of SP-D with highly branched viral N-linked glycans on hemagglutinin (HA), an abundant IAV envelope protein and critical virulence factor, promote viral aggregation and neutralization through as yet unknown molecular mechanisms. Two truncated human SP-D forms, wild-type (WT) and double mutant D325A+R343V, representing neck and carbohydrate recognition domains are compared in this study. Whereas both WT and D325A+R343V bind to isolated glycosylated HA, WT does not inhibit IAV in neutralization assays; in contrast, D325A+R343V neutralization compares well with that of full-length native SP-D. To elucidate the mechanism for these biochemical observations, we have determined crystal structures of D325A+R343V in the presence and absence of a viral nonamannoside (Man9). On the basis of the D325A+R343V-Man9 structure and other crystallographic data, models of complexes between HA and WT or D325A+R343V were produced and subjected to molecular dynamics. Simulations reveal that whereas WT and D325A+R343V both block the sialic acid receptor site of HA, the D325A+R343V complex is more stable, with stronger binding caused by additional hydrogen bonds and hydrophobic interactions with HA residues. Furthermore, the blocking mechanism of HA differs for WT and D325A+R343V because of alternate glycan binding modes. The combined results suggest a mechanism through which the mode of SP-D-HA interaction could significantly influence viral aggregation and neutralization. These studies provide the first atomic-level molecular view of an innate host defense lectin inhibiting its viral glycoprotein target.

How does (E)-2-(acetamidomethylene)succinate bind to its hydrolase? From the binding process to the final result

The binding of (E)-2-(acetamidomethylene)succinate (E-2AMS) to E-2AMS hydrolase is crucial for biological function of the enzyme and the last step reaction of vitamin B(6) biological degradation. In the present study, several molecular simulation methods, including molecular docking, conventional molecular dynamics (MD), steered MD (SMD), and free energy calculation methods, were properly integrated to investigate the detailed binding process of E-2AMS to its hydrolase and to assign the optimal enzyme-substrate complex conformation. It was demonstrated that the substrate binding conformation with trans-form amide bond is energetically preferred conformation, in which E-2AMS's pose not only ensures hydrogen bond formation of its amide oxygen atom with the vicinal oxyanion hole but also provides probability of the hydrophobic interaction between its methyl moiety and the related enzyme's hydrophobic cavity. Several key residues, Arg146, Arg167, Tyr168, Arg179, and Tyr259, orientate the E-2AMS's pose and stabilize its conformation in the active site via the hydrogen bond interaction with E-2AMS. Sequentially, the binding process of E-2AMS to E-2AMS hydrolase was studied by SMD simulation, which shows the surprising conformational reversal of E-2AMS. Several important intermediate structures and some significant residues were identified in the simulation. It is stressed that Arg146 and Arg167 are two pivotal residues responsible for the conformational reversal of E-2AMS in the binding or unbinding. Our research has shed light onto the full binding process of the substrate to E-2AMS hydrolase, which could provide more penetrating insight into the interaction of E-2AMS with the enzyme and would help in the further exploration on the catalysis mechanism.

Molecular dynamics simulations suggest ligand's binding to nicotinamidase/pyrazinamidase

The research on the binding process of ligand to pyrazinamidase (PncA) is crucial for elucidating the inherent relationship between resistance of Mycobacterium tuberculosis and PncA’s activity. In the present study, molecular dynamics (MD) simulation methods were performed to investigate the unbinding process of nicotinamide (NAM) from two PncA enzymes, which is the reverse of the corresponding binding process. The calculated potential of mean force (PMF) based on the steered molecular dynamics (SMD) simulations sheds light on an optimal binding/unbinding pathway of the ligand. The comparative analyses between two PncAs clearly exhibit the consistency of the binding/unbinding pathway in the two enzymes, implying the universality of the pathway in all kinds of PncAs. Several important residues dominating the pathway were also determined by the calculation of interaction energies. The structural change of the proteins induced by NAM’s unbinding or binding shows the great extent interior motion in some homologous region adjacent to the active sites of the two PncAs. The structure comparison substantiates that this region should be very important for the ligand’s binding in all PncAs. Additionally, MD simulations also show that the coordination position of the ligand is displaced by one water molecule in the unliganded enzymes. These results could provide the more penetrating understanding of drug resistance of M. tuberculosis and be helpful for the development of new antituberculosis drugs.

Theoretical improvement of the specific inhibitor of human carbonic anhydrase VII

The selectivity of a known arylsulfonamides inhibitor for two isozymes II and VII of human carbonic anhydrases (hCAs) was studied by homology modeling, molecular docking and molecular dynamics methods. The results show that the selectivity of the inhibitor for two isozymes is due to the different side chain lengths between N67 of hCA II and Q64 of hCA VII. One more methene group in the side chain of Q64 of hCA VII makes it possible to form the hydrogen bond with the bromide atom of the known inhibitor. From the point of view, the modification to the known inhibitor was performed to obtain an inhibitor with higher selectivity. The complex conformations of the new designed inhibitor and two isozymes designate the formation of the hydrogen bond between the newly added group (hydroxypropyl group) and Q64 of hCA VII but N67 of hCA II. The results of the binding free energy from the MM/PBSA approach also prove the selectivity improvement of the new inhibitor in comparison with the known inhibitor. The work will help the design of the isozyme-specific inhibitors of hCA VII.