Presented is an extension of the CHARMM General force field (CGenFF) to enable the modeling of sulfonyl-containing compounds. Model compounds containing chemical moieties such as sulfone, sulfonamide, sulfonate and sulfamate were used as the basis for the parameter optimization. Targeting high-level quantum mechanical and experimental crystal data, the new parameters were optimized in a hierarchical fashion designed to maintain compatibility with the remainder of the CHARMM additive force field. The optimized parameters satisfactorily reproduced equilibrium geometries, vibrational frequencies, interactions with water, gas phase dipole moments and dihedral potential energy scans. Validation involved both crystalline and liquid phase calculations showing the newly developed parameters to satisfactorily reproduce experimental unit cell geometries, crystal intramolecular geometries and pure solvent densities. The force field was subsequently applied to study conformational preference of a sulfonamide based peptide system. Good agreement with experimental IR/NMR data further validated the newly developed CGenFF parameters as a tool to investigate the dynamic behavior of sulfonyl groups in a biological environment. CGenFF now covers sulfonyl group containing moieties allowing for modeling and simulation of sulfonyl-containing compounds in the context of biomolecular systems including compounds of medicinal interest.

The General AMBER Force Field (GAFF) has been broadly used by researchers all over the world to perform in silico simulations and modelings on diverse scientific topics, especially in the field of computer-aided drug design whose primary task is to accurately predict the affinity and selectivity of receptor-ligand binding. The atomic partial charges in GAFF and the second generation of GAFF (GAFF2) were originally developed with the quantum mechanics derived restrained electrostatic potential charge, but in practice, users usually adopt an efficient charge method, Austin Model 1-bond charge corrections (AM1-BCC), based on which, without expensive ab initio calculations, the atomic charges could be efficiently and conveniently obtained with the ANTECHAMBER module implemented in the AMBER software package. In this work, we developed a new set of BCC parameters specifically for GAFF2 using 442 neutral organic solutes covering diverse functional groups in aqueous solution. Compared to the original BCC parameter set, the new parameter set significantly reduced the mean unsigned error (MUE) of hydration free energies from 1.03 kcal/mol to 0.37 kcal/mol. More excitingly, this new AM1-BCC model also showed excellent performance in the solvation free energy (SFE) calculation on diverse solutes in various organic solvents across a range of different dielectric constants. In this large-scale test with totally 895 neutral organic solventsolute systems, the new parameter set led to accurate SFE predictions with the MUE and the root-mean-square-error of 0.51 kcal/mol and 0.65 kcal/mol, respectively. This newly developed charge model, ABCG2, paved a promising path for the next generation GAFF development.

We investigated the effects of 17 widely used atomistic molecular mechanics force fields (MMFFs) on the structures and kinetics of amyloid peptide assembly. To this end, we performed large-scale all-atom molecular dynamics simulations in explicit water on the dimer of the sevenresidue fragment of the Alzheimer's amyloid-β peptide, Aβ16−22, for a total time of 0.34 ms. We compared the effects of these MMFFs by analyzing various global reaction coordinates, secondary structure contents, the fibril population, the in-register and out-of-register architectures, and the fibril formation time at 310 K. While the AMBER94, AMBER99, and AMBER12SB force fields do not predict any β-sheets, the seven force fields, AMBER96, GROMOS45a3, GROMOS53a5, GROMOS53a6, GROMOS43a1, GROMOS43a2, and GROMOS54a7, form β-sheets rapidly. In contrast, the following five force fields, AMBER99-ILDN, AMBER14SB, CHARMM22*, CHARMM36, and CHARMM36m, are the best candidates for studying amyloid peptide assembly, as they provide good balances in terms of structures and kinetics. We also investigated the assembly mechanisms of dimeric Aβ16−22 and found that the fibril formation rate is predominantly controlled by the total β-strand content.

Release of cytochrome c from the mitochondrial intermembrane space is critical to apoptosis induced by a variety of death stimuli. Bid is a BH3-only prodeath Bcl-2 family protein that can potently activate this efflux. In the current study, we investigated the mitochondrial localization of Bid and its interactions with mitochondrial phospholipids, focusing on their relationships with Bid-induced cytochrome c release. We found that Bid binding to the mitochondria required only three of its eight helical structures (alpha4-alpha6), but not the BH3 domain, and the binding could not be inhibited by the antideath molecule Bcl-x(L). Membrane fractionations indicated that tBid bound to mitochondrial outer membranes at both contact and noncontact sites. Bid could interact with specific cardiolipin species on intact mitochondria as identified by mass spectrometry. Like the binding to the mitochondria, this interaction could not be blocked by the mutation in the BH3 domain or by Bcl-x(L.) However, a cardiolipin-specific dye, 10-N-nonyl acridine orange, could preferentially suppress Bid binding to the mitochondrial contact site and inhibit Bid-induced mitochondrial cristae reorganization and cytochrome c release. These findings thus suggest that interactions of Bid with mitochondrial cardiolipin at the contact site can contribute significantly to its functions.

Additive force fields are designed to account for induced electronic polarization in a mean-field average way using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parameterized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-Dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parameterization strategy based on a restrained ensemble - maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field that includes water, ions, proteins, DNA and selected carbohydrates.

The mitochondrial pathway is critical for the efficient execution of death receptor-initiated apoptosis in certain cell types. Questions remain as to why the mitochondria are required in that scenario. We investigated the molecular events that determined the need for the mitochondria by using an in vivo model of anti-Fas-induced hepatocyte apoptosis. In wild-type mice, Fas stimulation resulted in normal activation of caspase-3, with the generation of the active p19-p12 complex. In biddeficient mice, caspase-3 activation was arrested after the initial cleavage at Asp 175 . This allowed the generation of the p12 small subunit, but the p20 large subunit could not be further processed to the p19 subunit. The p20-p12 complex generated by Fas stimulation in biddeficient hepatocytes was inactive, arresting the death program. Failure of p20/p12 caspase-3 to mature and to exhibit activity was because of the inhibition by the inhibitor-of-apoptosis proteins (IAPs), such as XIAP, and also to a low caspase-8 activity. This block could be overcome in wild-type mice by two mechanisms. Smac was released from mitochondria early following Fas activation and was competitively bound to the IAPs to reverse their effects. XIAP could also be cleaved, and this occurred later and was likely mediated by enhanced caspase activities. Both mechanisms were dependent on Bid and thus were not operative in bid-deficient hepatocytes. In conclusion, mitochondrial activation by Bid is required for reversing the IAP inhibition through Smac release. It is also required for the alternative activation of caspases through cytochrome c release, as demonstrated previously. Together, these events ensure a successful progression of the death program initiated by the death receptor activation in the hepatocyte.

A polarizable empirical force field based on the classical Drude oscillator is presented for the hexopyranose form of selected monosaccharides. Parameter optimization targeted quantum mechanical (QM) dipole moments, solute - water interaction energies, vibrational frequencies, and conformational energies. Validation of the model was based on experimental data on crystals, densities of aqueous-sugar solutions, diffusion constants of glucose, and rotational preferences of the exocylic hydroxymethyl of D-glucose and D-galactose in aqueous solution as well as additional QM data. Notably, the final model involves a single electrostatic model for all sixteen diastereomers of the monosaccharides, indicating the transferability of the polarizable model. The presented parameters are anticipated to lay the foundation for a comprehensive polarizable force field for saccharides that will be compatible with the polarizable Drude parameters for lipids and proteins, allowing for simulations of glycolipids and glycoproteins.

Accurate prediction of the absolute or relative protein−ligand binding affinity is one of the major tasks in computer-aided drug design projects, especially in the stage of lead optimization. In principle, the alchemical free energy (AFE) methods such as thermodynamic integration (TI) or free-energy perturbation (FEP) can fulfill this task, but in practice, a lot of hurdles prevent them from being routinely applied in daily drug design projects, such as the demanding computing resources, slow computing processes, unavailable or inaccurate force field parameters, and difficult and unfriendly setting up and post-analysis procedures. In this study, we have exploited practical protocols of applying the CPU (central processing unit)-TI and newly developed GPU (graphic processing unit)-TI modules and other tools in the AMBER software package, combined with ff14SB/GAFF1.8 force fields, to conduct efficient and accurate AFE calculations on protein−ligand binding free energies. We have tested 134 protein−ligand complexes in total for four target proteins (BACE, CDK2, MCL1, and PTP1B) and obtained overall comparable performance with the commercial Schrodinger FEP+ program (Wang et al.