Large-scale conformational changes in proteins are often associated with the binding of a substrate. Because conformational changes may be related to the function of an enzyme, understanding the ...kinetics and energetics of these motions is very important. We have delineated the atomically detailed conformational transition pathway of the phosphotransferase enzyme adenylate kinase (AdK) in the absence and presence of an inhibitor. The computed free energy profiles associated with conformational transitions offer detailed mechanistic insights into, as well as kinetic information on, the ligand binding mechanism. Specifically, potential of mean force calculations reveal that in the ligand-free state, there is no significant barrier separating the open and closed conformations of AdK. The enzyme samples near closed conformations, even in the absence of its substrate. The ligand binding event occurs late, toward the closed state, and transforms the free energy landscape. In the ligand-bound state, the closed conformation is energetically most favored with a large barrier to opening. These results emphasize the underlying dynamic nature of the enzyme and indicate that the conformational transitions in AdK are more intricate than a mere two-state jump between the crystal-bound and -unbound states. Based on the existence of the multiple conformations of the enzyme in the open and closed states, a different viewpoint of ligand binding is presented. Our estimated activation energy barrier for the conformational transition is also in reasonable accord with the experimental findings.
Prolactin and the prolactin receptors are members of a family of hormone/receptor pairs which include GH, erythropoietin, and other ligand/receptor pairs. The mechanisms of these ligand/receptor ...pairs have broad similarities, including general structures, ligand/receptor stoichiometries, and activation of several common signaling pathways. But significant variations in the structural and mechanistic details are present among these hormones and their type 1 receptors. The prolactin receptor is particularly interesting because it can be activated by three sequence-diverse human hormones: prolactin, GH, and placental lactogen. This system offers a unique opportunity to compare the detailed molecular mechanisms of these related hormone/receptor pairs. This review critically evaluates selected literature that informs these mechanisms, compares the mechanisms of the three lactogenic hormones, compares the mechanism with those of other class 1 ligand/receptor pairs, and identifies information that will be required to resolve mechanistic ambiguities. The literature describes distinct mechanistic differences between the three lactogenic hormones and their interaction with the prolactin receptor and describes more significant differences between the mechanisms by which other related ligands interact with and activate their receptors.
Buried ionizable groups in proteins often play important structural and functional roles. However, it is generally challenging to study the detailed molecular mechanisms solely based on experimental ...measurements. Free energy calculations using atomistic simulations, on the other hand, complement experimental studies and can provide high temporal and spatial resolution information that can lead to mechanistic insights. Nevertheless, it is also well recognized that sufficient sampling of such atomistic simulations can be challenging, considering that structural changes related to the buried charges may be very slow. In the present study, we describe a simple but effective enhanced sampling technique called replica exchange with charge tempering (REChgT) with a novel free energy method, multisite λ dynamics (MSλD), to study two systems containing buried charges, pK a prediction of a small molecule, orotate, in complex with the dihydroorotate dehydrogenase, and relative stability of a Glu–Lys pair buried in the hydrophobic core of two variants of Staphylococcal nuclease. Compared to the original MSλD simulations, the usage of REChgT dramatically increases sampling in both conformational and alchemical spaces, which directly translates into a significant reduction of wall time to converge the free energy calculations. This study highlights the importance of sufficient sampling toward developing improved free energy methods.
Allostery is involved in innumerable biological processes and plays a fundamental role in human disease. Thus, the exploration of allosteric modulation is crucial for research on biological ...mechanisms and in the development of novel therapeutics. The development of small-molecule allosteric effectors can be used as tools to probe biological mechanisms of interest. One of the main limitations in targeting allosteric sites is the difficulty in uncovering them for specific receptors. Furthermore, upon discovery of novel allosteric modulation, early lead generation is made more difficult as compared to that at orthosteric sites because there is likely no information about the types of molecules that can bind at the site. In the work described here, we present a novel drug discovery pipeline, FASTDock, which allows one to uncover ligandable sites as well as small molecules that target the given site without requiring pre-existing knowledge of ligands that can bind in the targeted site. By using a hierarchical screening strategy, this method has the potential to enable high-throughput screens of an exceptionally large database of targeted ligand space.
Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) ...in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.
Growing evidence suggests that the β-amyloid (Aβ) peptides of Alzheimer's disease are generated in early endosomes and that small oligomers are the principal toxic species. We sought to understand ...whether and how the solution pH, which is more acidic in endosomes than the extracellular environment, affects the conformational processes of Aβ. Using constant pH molecular dynamics simulations of two model peptides, Aβ(1-28) and Aβ(10-42), we found that the folding landscape of Aβ is strongly modulated by pH and is most favorable for hydrophobically driven aggregation at pH 6. Thus, our theoretical findings substantiate the possibility that Aβ oligomers develop intracellularly before secretion into the extracellular milieu, where they may disrupt synaptic activity or act as seeds for plaque formation.
As most relevant motions in biomolecular systems are inaccessible to conventional molecular dynamics simulations, algorithms that enhance sampling of rare events are indispensable. Increasing ...interest in intrinsically disordered systems and the desire to target ensembles of protein conformations (rather than single structures) in drug development motivate the need for enhanced sampling algorithms that are not limited to “two-basin” problems, and can efficiently determine structural ensembles. For systems that are not well-studied, this must often be done with little or no information about the dynamics of interest. Here we present a novel strategy to determine structural ensembles that uses dynamically defined sampling regions that are organized in a hierarchical framework. It is based on the weighted ensemble algorithm, where an ensemble of copies of the system (“replicas”) is directed to new regions of configuration space through merging and cloning operations. The sampling hierarchy allows for a large number of regions to be defined, while using only a small number of replicas that can be balanced over multiple length scales. We demonstrate this algorithm on two model systems that are analytically solvable and examine the 10-residue peptide chignolin in explicit solvent. The latter system is analyzed using a configuration space network, and novel hydrogen bonds are found that facilitate folding.
Display omitted
A halogen bond is a highly directional, non-covalent interaction between a halogen atom and another electronegative atom. It arises due to the formation of a small region of positive ...electrostatic potential opposite the covalent bond to the halogen, called the ‘sigma hole.’ Empirical force fields in which the electrostatic interactions are represented by atom-centered point charges cannot capture this effect because halogen atoms usually carry a negative charge and therefore interact unfavorably with other electronegative atoms. A strategy to overcome this problem is to attach a positively charged virtual particle to the halogen. In this work, we extend the additive CHARMM General Force Field (CGenFF) to include such interactions in model systems of phenyl-X, with X being Cl, Br or I including di- and trihalogenated species. The charges, Lennard-Jones parameters, and halogen-virtual particle distances were optimized to reproduce the orientation dependence of quantum mechanical interaction energies with water, acetone, and N-methylacetamide as well as experimental pure liquid properties and relative hydration free energies with respect to benzene. The resulting parameters were validated in molecular dynamics simulations on small-molecule crystals and on solvated protein–ligand complexes containing halogenated compounds. The inclusion of positive virtual sites leads to better agreement across experimental observables, including preservation of ligand binding poses as a direct result of the improved representation of halogen bonding.