Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Analysis of Conformational Variation in Macromolecular Structural Models

Experimental conditions or the presence of interacting components can lead to variations in the structural models of macromolecules. However, the role of these factors in conformational selection is often omitted by in silico methods to extract dynamic information from protein structural models. Structures of small peptides, considered building blocks for larger macromolecular structural models, can substantially differ in the context of a larger protein. This limitation is more evident in the case of modeling large multi-subunit macromolecular complexes using structures of the individual protein components. Here we report an analysis of variations in structural models of proteins with high sequence similarity. These models were analyzed for sequence features of the protein, the role of scaffolding segments including interacting proteins or affinity tags and the chemical components in the experimental conditions. Conformational features in these structural models could be rationalized by conformational selection events, perhaps induced by experimental conditions. This analysis was performed on a non-redundant dataset of protein structures from different SCOP classes. The sequence-conformation correlations that we note here suggest additional features that could be incorporated by in silico methods to extract dynamic information from protein structural models.

``Conformational Simulation'' of Sulfamethizole by Molecular Complexation and Insights from Charge Density Analysis: Role of Intramolecular S center dot center dot center dot O Chalcogen Bonding

Intramolecular S center dot center dot center dot O chalcogen bonding and its potential to lock molecular conformation have been examined in the crystal forms of sulfamethizole, a sulfonamide antibiotic. Molecular complexes of sulfamethizole, including salts and cocrystal, have been synthesized, and their crystal structures were analyzed in order to examine the possible conformational preferences of the molecule in various ionic states and supramolecular environments (neutral/cocrystal, anionic salt, and cationic salt forms). The electrostatic potential mapped on Hirshfeld surfaces generated for these crystal forms provides insights into the possible binding modes of the drug in different environments. Further, the observed conformation locking feature has been rationalized in terms of the experimental charge density features of the intramolecular S center dot center dot O chalcogen bonding in sulfamethizole. The study quantitatively illustrates and rationalizes an intriguing case of a local minimum of molecular conformation being exclusively preferred over the global minimum, as it facilitates more efficient intermolecular interactions in a supramolecular environment.

Conformational flexibility of DNA: polymorphism and handedness

It is shown that left-handed duplexes are possible for A, B, and D forms of DNA. These duplexes are stereochemically satisfactory and are consistent with the observed x-ray intensity data. On scrutiny the refined right-handed models of B and D DNA by Arnott and coworkers are found to be stereochemically unacceptable. It was possible to formulate a stereochemical guideline for molecular model building based on theory and analysis of single-crystal structure data of dinucleoside monophosphate and higher oligomers. This led to both right- and left-handed DNA duplexes. The right-handed B and D DNA duplexes so obtained are stereochemically superior to earlier models and agree well with the observed x-ray intensity data. The observation that DNA can exist in either handedness for all the polymorphous forms of DNA at once explained A in equilibrium B and B in equilibrium D transitions. Hence it is confirmed that polymorphism of DNA is a reflection on the conformational flexibility inherent in DNA, the same cause that ultimately allows DNA in either handedness. The possibility of various types of right- and left-handed duplexes generated by using dinucleoside monophosphate and trinucleoside diphosphate as repeating units resulted in a variety of models, called RL models. All these models have alternating right and left helical segments and inverted stacking at the bend region as suggested by us earlier. It turns out that the B-Z DNA model of Wang et al. is only an example of RL models.

Conformational Dynamics in Alkyl Chains of an Anchored Bilayer: A Molecular Dynamics Study

Molecular dynamics (MD) simulations are reported for an anchored bilayer formed by the intercalation of cetyl trimethyl ammonium (CTA) and CH3(CH2)15N+(CH3) ions in a layered solid, CdPS3. The intercalated CTA ions are organized with the cationic headgroups tethered to the inorganic sheet and the hydrocarbon tails arranged as bilayers. Simulations were performed at three temperatures, 65, 180, and 298 K, using an isothermal−isobaric ensemble that was subsequently switched once macroscopic parameters had converged to a canonical isothermal−isochoric ensemble. The simulations are able to reproduce the experimental features of this system, including the formation of the bilayer and layer-to-layer separation distance. An analysis of the conformation of the chains showed that at all three temperatures a fraction of the alkyl chains retained a planar all-trans conformation, and that gauche bonds occurred as part of a “kink” (gauche+−trans−gauche−) sequence and not as isolated gauche bonds. Trans−gauche isomerization rates for the alkyl chains in the anchored bilayer are slower than those in lipid bilayers at the same temperature and show a progressive increase as the torsion numbers approach the tail. A two-dimensional periodic Voronoi tessellation analysis was performed to obtain the single-molecular area of an alkyl chain in the bilayer. The single-molecular area relaxation times are an order of magnitude longer than the trans−gauche isomerization times. The results indicate that the trans−gauche isomerization is associated with the creation and annihilation of a kink defect sequence. The results of the present MD simulation explain the apparent conflicting estimates of the gauche disorder in this system as obtained from infrared and 13C nuclear magnetic resonance measurements.

Conformational Characterization of S-Methyl Dithiocarbazate by Infrared Spectra and Vibrational Assignments

Infrared spectra are recorded for S-methyl dithiocarbazate and its N-deuterated compound in two molecular conformations in the solid state and in solution between 4000 and 30 cm−1. The assignments have been supported from a complete normal coordinate analysis; the conformation sensitive bands of the –CSNHNH2 grouping are discussed. The assignments are compared with those of related molecules to check the internal consistency and to obtain the pattern of the characteristic bands of thiocarbazoyl (–CSNHNH2) group. The magnitudes of the C–N and S–CH3 torsional barriers are estimated from the force constants.

Conformationally Restricted Nucleocyclitols: a Study into their Conformational Preferences and Supramolecular Architecture in the Solid State

With the intent of probing the feasibility of employing annulation as a tactic to engender axial rich conformations in nucleoside analogues, two adenine-derived, ``conformationally restricted'' nucleocylitols, 9 and 10, have been conceptualized as representatives of a hitherto unexplored class of nucleic acid base-cyclitol hybrids. A general synthetic strategy, with an inherent scope for diversification, allowed rapid functionalization of indane and tetralin to furnish 9 and 10 respectively in fair yield. Single-crystal X-ray diffraction analysis revealed that the two nucleocyclitols under study, though homologous, present completely dissimilar modes of molecular packing, marked, in particular, by the nature of involvement of the adenynyl NH2 group in the supramolecular assembly. In addition, the crystal structures of 9 and 10 also exhibit two different conformations of the functionalized cyclohexane ring. Thus, while the six-membered carbocycle in cyclopenta-annulated 9 exists in the expected chair (C) conformation that in cyclohexaannulated 10, which crystallizes as a dihydrate, shows an unusual twist-boat (TB) conformation. From a close analysis of the (HNMR)-H-1 spectroscopic data recorded for 9 and 10 in CD3OD, it was possible to put forth a putative explanation for the uncanny conformational preferences of crystalline 9 and 10.

Conformational Choices for the Stereochemically Constrained gamma-Amino Acid Residue Gabapentin: Theoretical Studies and Correlation With Experimental Results

Gabapentin (1-aminomethylcyclohexaneacetic acid, Gpn) is an achiral, conformationally constrained gamma amino acid residue. A survey of available crystal structures of Gpn peptides reveals that the torsion angles about the C-gamma-C-beta (theta(1)) and C-beta-C-alpha(theta(2)) bonds are overwhelmingly limited to gauche, gauche (g(+)g(+)/g(-)g(-)) conformations. The Gpn residue forms C-7 and C-9 hydrogen bonds in which the donor and acceptor atoms come from the flanking peptide units. In combination with alpha amino acid residues alpha gamma and gamma alpha segments can adopt C-12 hydrogen bonded structures. The conformational choices available to the Gpn residue have been probed using energy calculations, adopting a grid search strategy. Ramachandran phi-psi maps have been constructed for fixed values of theta(1) and theta(2), corresponding to the gauche and trans conformations. The sterically allowed and energetically favorable regions of conformational space have been defined and experimental observations compared. C-7 and C-9 hydrogen bonded conformational families have been identified using a grid search approach in which theta(1) and theta(2) values are varied over a range of +/- 10 degrees about ideal values at 1 degrees intervals. The theoretical analysis together with experimental observations for 59 Gpn residues from 35 crystal structures permits definition of the limited range of conformational possibilities at this gamma amino acid residue. .

Conformational studies on cyclic dipeptides

An analysis of 11 crystal structures of cyclic dipeptides so far reported in the literature is made, with main reference to the internal parameters of these molecules. Preferred conformations of the side chains of cyclic dipeptides with different α-amino acid residues have been studied by classical energy calculations. The possible conformations of the DKP ring are also studied. The significance of the non-bonded interaction in deciding the pathway for conformational change has also been investigated. The agreement between theoretical results and experimental observations is quite good, both with respect to the conformation of these molecules as well as the enthalpy difference as estimated from n.m.r. studies between different conformers.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.

Conformational characteristics of asparaginyl residues in proteins

Backbone conformations at 1064 asparaginyl residues in 123 non-homologous, high-resolution X-ray structures of proteins were analysed. Asn adopts conformations in left-handed x-helical region and other partially allowed regions in the Ramachandran map more readily than any other non-glycyl residue. Asn conformational clusters in the (phi,psi) regions of left-handed alpha-helix, right-handed alpha-helix and extended (beta) strands were investigated in detail for their occurrence in various secondary structures, especially in beta-turn regions. Preferences were observed for Asn conformations in different positions in various beta-turn types, including the first and fourth positions of the turn. Asparaginyl residues with extended conformations are found to occur frequently in irregular regions, although they are expected to occur predominantly in extended strands or in the third position of type II beta-turns. Asn conformations at the N-cap positions of helices strongly prefer extended conformation than alpha(L), which seems to be characteristic of non-glycyl residues at that position. In the linkers connecting two extended strands and those connecting an alpha-helix and an extended strand, Asn with alpha(L) or alpha(R) conformation is more favoured than Asn with the beta-conformation. Analysis of Asn-Asn doublets and Asn-X-Asn triplets permitted identification of conformational families in such sequences. Results of this investigation provide useful hints in modelling Asn-rich regions in proteins such as malaria parasite coat protein. (C) Munksgaard 1994.

Stereochemical analysis of the antigenic tip of the V3 loop peptide of HIV-1 gp120

A novel multiple turn conformation has been observed for a segment GPGRAFY in the crystal structure of a complex of HIV-1 gp120 V3 loop peptide with the Fab fragment of a neutralizing antibody [Ghiara ct al. (1994) Science 264, 82-85]. A structural motif has been defined for the peptide segment, employing idealized backbone conformations characterized by ranges of virtual C-alpha torsion angles and bond angles. A search of 122 high-resolution protein crystal structures has permitted identification of 24 examples of similar structural motifs. Two major conformational families have been identified, which differ primarily in the conformation at residue 3. The observed conformation at residue 3 in family 1 is left-handed helical (alpha(L)) and that in family 2 is right-handed helical (alpha(R)). Of the 10 examples in family 1, 9 examples have Gly residues at position 3. Of the 12 examples in family 2, 7 examples have Asn/Asp at position 3. Computer modeling of the V3 loop tip sequence using the two backbone conformational families as starting points leads to minimum-energy conformations in which antigenically important side-chains occupy similar spatial arrangements. This stereochemical analysis of the V3 loop tip sequence suggests a rational basis for the design of synthetic analog peptides for use as viral antagonists or synthetic antigens. (C) Munksgaard 1995.

Stereochemical Analysis of Higher α,α-Dialkylglycine Containing Peptides. Characterization of Local Helical Conformations at Dipropylglycine Residues and Observation of a Novel Hydrated Multiple β-Turn Structure in Crystals of a Glycine Rich Peptide

The peptide Boc-Gly-Dpg-Gly-Gly-Dpg-Gly-NHMe (1) has been synthesized to examine the conformational preferences of Dpg residues in the context of a poor helix promoting sequence. Single crystals of 1 were obtained in the space group P21/c with a = 13.716(2) Å, b = 12.960(2) Å, c = 22.266(4) Å, and β = 98.05(1)°; R = 6.3% for 3660 data with |Fo| > 4σ. The molecular conformation in crystals revealed that the Gly(1)-Dpg(2) segment adopts φ, ψ values distorted from those expected for an ideal type II‘ β-turn (φGly(1) = +72.0°, ψGly(1) = −166.0°; φDpg(2) = −54.0°, ψDpg(2) = −46.0°) with an inserted water molecule between Boc-CO and Gly(3)NH. The Gly(3)-Gly(4) segment adopts φ, ψ values which lie broadly in the right handed helical region (φGly(3) = −78.0°, ψGly(3) = −9.0°; φGly(4) = −80.0°, ψGly(4) = −18.0°). There is a chiral reversal at Dpg(5) which takes up φ, ψ values in the left handed helical region. The Dpg(5)-Gly(6) segment closely resembles an ideal type I‘ β-turn (φDpg(5) = +56.0°, ψDpg(5) = +32.0°; φGly(6) = +85.0°, ψGly(6) = −3.0°). Molecules of both chiral senses are found in the centrosymmetric crystal. The C-terminus forms a hydrated Schellman motif, with water insertion into the potential 6 → 1 hydrogen bond between Gly(1)CO and Gly(6)NH. NMR studies in CDCl3 suggest substantial retention of the multiple turn conformation observed in crystals. In solution the observed NOEs support local helical conformation at the two Dpg residues.

Mechanism of lectin-cell interactions: thermodynamic and kinetic analysis of the binding of winged bean acidic lectin (WBA II) to human erythrocytes

In order to identify the forces involved in the binding and to understand the mechanism involved, equilibrium and kinetic studies were performed on the binding of the winged bean acidic lectin to human erythrocytes. The magnitudes of delta S and delta H were positive and negative respectively, an observation differing markedly from the lectin-simple sugar interactions where delta S and delta H are generally negative. Analysis of the sign and magnitudes of these values indicate that ionic and hydrogen bonded interactions prevail over hydrophobic interactions resulting in net -ve delta H (-37.12 kJ.mol-1) and +ve delta S (14.4 J.mole-1 K-1 at 20 degrees C), thereby suggesting that this entropy driven reaction also reflects conformational changes in the lectin and/or the receptor. Presence of two kinds of receptors for WBA II on erythrocytes, as observed by equilibrium studies, is consistent with the biexponential dissociation rate constants (at 20 degrees C K1 = 1.67 x 10(-3) M-1 sec-1 and K2 = 11.1 x 10(-3) M-1 sec-1). These two rate constants differed by an order of magnitude accounting for the difference in the association constants of the two receptors of WBA II. However, the association process remains monoexponential suggesting no observable difference in the association rates of the lectin molecule with both the receptors, under the experimental conditions studied. The thermodynamic parameters calculated from kinetic data correlate well with those observed by equilibrium. A two-step binding mechanism is proposed based on the kinetic parameters for WBA II-receptor interaction