We describe the development, current features, and some directions for future development of the Amber package of computer programs. This package evolved from a program that was constructed in the late 1970s to do Assisted Model Building with Energy Refinement, and now contains a group of programs embodying a number of powerful tools of modern computational chemistry, focused on molecular dynamics and free energy calculations of proteins, nucleic acids, and carbohydrates.

A new derivation of the GLYCAM06 force field, which removes its previous specificity for carbohydrates, and its dependency on the AMBER force field and parameters, is presented. All pertinent force field terms have been explicitly specified and so no default or generic parameters are employed. The new GLYCAM is no longer limited to any particular class of biomolecules, but is extendible to all molecular classes in the spirit of a small-molecule force field. The torsion terms in the present work were all derived from quantum mechanical data from a collection of minimal molecular fragments and related small molecules. For carbohydrates, there is now a single parameter set applicable to both α-and β-anomers and to all monosaccharide ring sizes and conformations. We demonstrate that deriving dihedral parameters by fitting to QM data for internal rotational energy curves for representative small molecules generally leads to correct rotamer populations in molecular dynamics simulations, and that this approach removes the need for phase corrections in the dihedral terms. However, we note that there are cases where this approach is inadequate. Reported here are the basic components of the new force field as well as an illustration of its extension to carbohydrates. In addition to reproducing the gas-phase properties of an array of small test molecules, condensed-phase simulations employing GLYCAM06 are shown to reproduce rotamer populations for key small molecules and representative biopolymer building blocks in explicit water, as well as crystalline lattice properties, such as unit cell dimensions, and vibrational frequencies.

Summary The SARS-CoV-2 betacoronavirus uses its highly glycosylated trimeric Spike protein to bind to the cell surface receptor angiotensin converting enzyme 2 (ACE2) glycoprotein and facilitate host cell entry. We utilized glycomics-informed glycoproteomics to characterize site-specific microheterogeneity of glycosylation for a recombinant trimer Spike mimetic immunogen and for a soluble version of human ACE2. We combined this information with bioinformatics analyses of natural variants and with existing 3D structures of both glycoproteins to generate molecular dynamics simulations of each glycoprotein both alone and interacting with one another. Our results highlight roles for glycans in sterically masking polypeptide epitopes and directly modulating Spike-ACE2 interactions. Furthermore, our results illustrate the impact of viral evolution and divergence on Spike glycosylation, as well as the influence of natural variants on ACE2 receptor glycosylation. Taken together, these data can facilitate immunogen design to achieve antibody neutralization and inform therapeutic strategies to inhibit viral infection.

The relationship between the three-dimensional structures of oligosaccharides and polysaccharides and their biological properties has been the focus of many recent studies. The overall conformation of an oligosaccharide depends primarily on the orientation of the torsion angles (, , and ) between glycosyl residues. Numerous experimental studies have shown that in glucopyranosides the -torsion angle (O6-C6-C5-O5) displays a preference for gauche orientations, in disagreement with predictions based on gas-phase quantum mechanics calculations. In contrast, the -angle in galactopyranosides displays a high proportion of the anti-orientation. For oligosaccharides containing glycosidic linkages at the 6-position (136 linked), variations in rotamer population have a direct effect on the oligosaccharides' structure and function, and yet the physical origin of these conformational preferences remains unclear. Although it is generally recognized that the gauche effect in carbohydrates is a solvent-dependent phenomenon, the mechanism through which solvent induces the gauche preference is not understood. In the present work, quantum mechanics and solvated molecular dynamics calculations were performed on two representative carbohydrates, methyl ␣-D-glucopyranoside and methyl ␣-Dgalactopyranoside. We show that correct reproduction of the experimental rotamer distributions about the -angles is obtained only after explicit water is included in the molecular dynamics simulations. The primary role of the water appears to be to disrupt the hydrogen bonding within the carbohydrate, thereby allowing the rotamer populations to be determined by internal electronic and steric repulsions between the oxygen atoms. The results reported here provide a quantitative explanation of the conformational behavior of (136)-linked carbohydrates.T he phenomenon of carbohydrate recognition is critical to many biological functions, including the response of the immune system to bacterial pathogens, the attachment of viral inf luenza particles to host cells, and the hyperacute rejection of tissue transplants from nonhuman sources. Understanding the conformational properties of carbohydrates is essential for the elucidation of their mechanisms of action, which may aid in the design of carbohydrate-based vaccines, antiviral drugs, and other therapeutic agents. A significant amount of research has gone into the study of the three-dimensional structures and dynamics of oligosaccharides and polysaccharides for these reasons (1-18).Unlike polypeptides and proteins, oligosaccharides do not contain secondary-structural motifs and do not form wellorganized tertiary structures in solution. Rather, oligosaccharides often populate multiple conformational families, thus requiring both temporal and spatial descriptors to quantify their conformational properties. Experimental methods such as NMR spectroscopy and x-ray crystallography have a long history of application to carbohydrates. However, these methods generally result in a single three-dimensional model for th...

Severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) has resulted in a pandemic and continues to spread around the globe at an unprecedented rate. To date, no effective therapeutic is available to fight its associated disease, COVID-19. Our discovery of a novel insertion of glycosaminoglycan (GAG)-binding motif at S1/S2 proteolytic cleavage site (681-686 (PRRARS)) and two other GAG-binding-like motifs within SARS-CoV-2 spike glycoprotein (SGP) led us to hypothesize that host cell surface GAGs may interact SARS-CoV-2 SGPs to facilitate host cell entry. Using a surface plasmon resonance direct binding assay, we found that both monomeric and trimeric SARS-CoV-2 SGP bind more tightly to immobilized heparin (K D ¼ 40 pM and 73 pM, respectively) than the SARS-CoV and MERS-CoV SGPs (500 nM and 1 nM, respectively). In competitive binding studies, the IC 50 of heparin, tri-sulfated non-anticoagulant heparan sulfate, and non-anticoagulant low molecular weight heparin against SARS-CoV-2 SGP binding to immobilized heparin were 0.056 μM, 0.12 μM, and 26.4 μM, respectively. Finally, unbiased computational ligand docking indicates that heparan sulfate interacts with the GAG-binding motif at the S1/S2 site on each monomer interface in the trimeric SARS-CoV-2 SGP, and at another site (453-459 (YRLFRKS)) when the receptor-binding domain is in an open conformation. The current study serves a foundation to further investigate biological roles of GAGs in SARS-CoV-2 pathogenesis. Furthermore, our findings may provide additional basis for further heparin-based interventions for COVID-19 patients exhibiting thrombotic complications.